python - python中的 Gamma 分布

Question

我们有qgamma在 R 和 gamm.inv在 excel 中，我无法使用 invgamma 获得相同的结果python中的函数。例如在excel中GAMMA.INV(0.99,35,0.08)=4.01 ，我可以使用 qgamma 从 R 获得相同的值与 probability=0.99, alpha(scale)=35.0, shape=0.08 .
但是当我使用 invgamma.pdf(0.99,35,0.8)它返回 1.61787636512e-15 .
正确的做法是什么？我都试过了gamma和 invgamma与 pdf , ppf , cdf但与我从 qgamma 获得的值不匹配在 R 或 GAMMA.INV在excel中

Answer 1

qgamma和 GAMMA.INV给 inverse cumulative distribution function的 gamma distribution .
invgamma代表 inverse gamma distribution ，这是完全不同的事情。解释差异超出了本网站的范围，因为这是一个统计问题，而不是编程问题，但在继续之前您应该真正理解它。

Gamma 分布的 SciPy 分布对象是 scipy.stats.gamma ，逆累积分布函数的方法是 ppf ，“百分点函数”的缩写(逆 CDF 的另一个名称)。因此，您应该使用 scipy.stats.gamma.ppf .确保传递正确的参数。