Как реализовать это распределение R Пуассона в Python?

Question

Я что-то кодировал в R, но я не могу сделать то же самое в Python.

Ниже приведен код - он определенно работает в R.

У меня проблемыс синтаксисом Python для достижения того же с numpy.

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