Python中的二元泊松分布

我想从二元泊松分布中抽取N次。 在R中是否有类似于包

bivpois

的Python模块

---

在Python中，我只知道库

scipy.stats.poisson

和

numpy.random.possion

，它们允许我根据单个参数lambda从单变量poisson分布进行绘制，但不是从双变量或多变量中。你可以很容易地自己做，因为我没有看到任何内置方法：

步骤：

用参数lambda_i生成3个独立的泊松变量Z_i

为i=1，2生成两个P_i=Z_i+Z_3，其跟随Poi（lambda_i+lambda_3）

代码：

---

这是一个只链接的评论。请在这里总结相关信息，以获得一个完整的答案。刚刚修复了我的答案，提供了一个关于doVery nice的总结。谢谢你的更新

import numpy
lam1 = 1
lam2 = 2
lam3 = 3
#wrap next part in a loop to generate more than 1 sample
a = np.random.poisson(lam1)
b = np.random.poisson(lam2)
c = np.random.poisson(lam3)
bivariate1 = a + c #follows Poi(lam1+lam3)
bivariate2 = b + c #follows Poi(lam2+lam3)