python中的伽马分布

Question

我们在 R 中有 qgamma，在 excel 中有 gamm.inv，我无法在 python 中使用 invgamma 函数获得相同的结果。例如在excel中 GAMMA.INV(0.99,35,0.08)=4.01，我可以使用 qgamma 和 probability=0.99, alpha(scale)=35.0, shape=0.08 从 R 获得相同的值。

但是当我使用invgamma.pdf(0.99,35,0.8) 时，它会返回1.61787636512e-15。

正确的做法是什么？我用pdf、ppf、cdf 尝试了gamma 和invgamma，但没有什么与我从R 中的qgamma 或excel 中的GAMMA.INV 获得的值相匹配

Answer 1

qgamma 和 GAMMA.INV 给出 inverse cumulative distribution function 的 gamma distribution。

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

伽马分布的SciPy分布对象是scipy.stats.gamma，逆累积分布函数的方法是ppf，是“百分点函数”的简称（逆CDF的别称）。因此，您应该使用scipy.stats.gamma.ppf。确保传递正确的参数。