如何在Python中实现R Poisson分布？

Question

我已经用R编写了一些代码，但似乎无法在Python中做同样的事情。

下面是代码-它绝对可以在R中工作。

我在使用numpy实现Python语法方面遇到麻烦。

myMaxAC = qpois(p=as.numeric(0.95),
                lambda=(121412)*(0.005))

为清楚起见，0.95是置信区间，121412是我的人口规模，0.005是人口中的频率。

我只想知道如何在Python中获得相同的答案，偶然是648。

Answer 1

您可以使用poisson.ppf：

from scipy.stats import poisson
myMaxAC = poisson.ppf(0.95, (121412)*(0.005))
print(myMaxAC)
648.0