Прежде всего, вам нужно только найти вероятность больше семи, верно?
Если да, я думаю, что это один из подходов, поэтому просим следовать:
from scipy import stats
occur_past = 5
ask_current = 7
mu = occur_past
x = ask_current
upto7 = scipy.stats.poisson.pmf(x, mu)
above7 = 1 - upto7
Вывод будет 0,8955551370429461
Вы можете просто преобразовать в функцию следующим образом:
def poisson(mu, x):
uptocurrent = scipy.stats.poisson.pmf(x, mu)
abovecurrent = 1 - uptocurrent
return abovecurrent
Надеюсь, это поможет вам.