在 Python 中使用泊松-伽马分布？

问：

有了泊松分布，我想修改它并使用泊松-伽马分布（不幸的是，维基百科上似乎不存在英文链接），这是一种与泊松随机变量 Poiss（λ） 相关的概率分布，其中参数 λ 不是恒定的，而是作为 Gamma 随机变量变化的;泊松混合物，其中参数 λ 具有 Gamma 分布。

我怎样才能申请我的？（在 scipy、numpy 或其他中）Poisson-Gamma Distributionfrom_1_to_3

from scipy.stats import poisson

mu_X = 1.7

one = poisson.pmf(1, mu_X) 
two = poisson.pmf(2, mu_X) 
three = poisson.pmf(3, mu_X) 

from_1_to_3 = one + two + three

更新在评论中，他们向我指出，也许泊松-伽马分布实际上是伽马-泊松混合物（维基百科链接）。如果是这种情况（我不知道是否是这种情况，我不确定），这是否意味着失败和成功的数量也会发挥作用，因此答案将是这个 Stackoverflow 答案 Scipy.nbinom vs 负二项分布公式？Negative binomial distribution

我希望有人能帮助我。提前感谢大家！