python - Scipy: ワイブル分布に合わせるには?

Question

ワイブル分布をフィッティングすることで完全に混乱しました。

weibull_params = sp.stats.exponweib.fit(df.speed, floc=0, f0=1)
# Returns (1, 1.7358162061451901, 0, 9.4955614228786978)

これらのパラメーターはhttps://en.wikipedia.org/wiki/Weibull_distributionのワイブル分布にどのように対応していますか? 具体的にはa、c wiki のlambdaとはk何ですか?

http://docs.scipy.org/doc/scipy/reference/generated/scipy.stats.exponweib.html#scipy.stats.exponweibでは、pdfは次のように定義されています。

exponweib.pdf(x, a, c) =
    a * c * (1-exp(-x**c))**(a-1) * exp(-x**c)*x**(c-1)

しかしウィキペディアでは、PDF は

---

exponweib.pdfまた、返されるパラメーターを使用すると、異なる結果が得られます

df['speed'].hist(bins=arange(0, df.speed.max()), alpha=0.5, normed=True)

def weib(x,lamb,k):
    return (k / lamb) * (x / lamb)**(k-1) * np.exp(-(x/lamb)**k)
k_shape, lamb_scale = weibull_params[1], weibull_params[3]
plt.plot(x, weib(x, lamb_scale, k_shape), label='self-defined weibull')

plt.plot(x, sp.stats.exponweib.pdf(x, k_shape, lamb_scale, loc=0, scale=1),'--', label ='custom_order')
plt.legend()