python - scipy：如何使用 weibull_min.pdf？

Question

我发现scipy's 的weibull_min pdf功能不直观且难以使用。我想生成比例= 30 和形状= 2.5 的weibull PDF。这应该是这样的：

这是我到目前为止所尝试的：

import matplotlib.pyplot as plt
from scipy.stats import weibull_max, uniform
import numpy as np
aoas = np.linspace(0, 8, 1000)
speeds = np.linspace(1, 80, 1000)
plt.fill_between(speeds, speeds * 0, weibull_max.pdf(speeds, 2.5, 30), facecolor='k')
plt.ylabel('Probability Density')
plt.savefig('speedDist.pdf')
plt.clf()

当然，我的代码是错误的。文档说，weibull_minexpects是分位数，pdf(x, c, loc=0, scale=1)是形状因子，是比例因子。但是，当我将代码更改为 时，输出 pdf 仅由零值组成。所以，我很困惑。如何生成所需的 Weibull PDF？xcscaleweibull_max.pdf(speeds, 2.5, 0, 30)

标签： pythonscipystatisticsprobabilitydistribution

Answer 1

您的代码中有两个问题。

• weibull_min并且weibull_max分布不一样。您weibull_min在问题的标题和描述中引用，但您weibull_max在代码中使用。要匹配您显示的图，请使用weibull_min.

• 该pdf方法的参数是(x, shape, loc, scale)。你写了pdf(speeds, 2.5, 30)30 给loc参数。相反，您应该使用pdf(speeds, 2.5, scale=30)

这是您的脚本的变体：

import numpy as np
import matplotlib.pyplot as plt
from scipy.stats import weibull_min


speeds = np.linspace(0, 80, 1000)    
p = weibull_min.pdf(speeds, 2.5, scale=30)

plt.plot(speeds, p, 'b', linewidth=1)
plt.fill_between(speeds, speeds * 0, p, facecolor='b', alpha=0.1)
plt.ylabel('Probability Density')
plt.show()

它生成这个图：