代码之家 › 专栏 › 技术社区 › Aaron

scipy gamma分布与维基百科上的公式不匹配

gamma-function scipy python

0

Aaron · 技术社区 · 7 年前

根据维基百科,伽马分布的预期值为 k * theta . 然而,当我将其插入伽马分布的CDF时,我并没有得到预期的0.5。

k = 1.5
theta = 2.1
expected_value = k * theta
scipy.stats.gamma.cdf(expected_value, k, scale=theta)

最后一行的结果是0.6083,而不是预期的0.5。

1 回复 | 直到 7 年前

1

2

Warren Weckesser 7 年前

CDF在 median 为0.5。 k*theta 给出平均值(即预期值),而不是中值。

In [8]: from scipy.stats import gamma

In [9]: k = 1.5

In [10]: theta = 2.1

In [11]: med = gamma.median(k, scale=theta)

In [12]: med
Out[12]: 2.4842725785941049

In [13]: gamma.cdf(med, k, scale=theta)
Out[13]: 0.49999999999999994