Python中的二元泊松分布

Question

我想从二元Possion分布中画N次。R中是否有类似于bivpois包的Python模块?

在Python中，我只知道库scipy.stats.poisson和numpy.random.possion，它们允许我从单变量泊松分布中进行绘制。依赖于单个参数lambda，但不依赖于二元或多元参数

Answer 1

你可以很容易地自己做，因为我没有看到任何内置方法:

https://en.wikipedia.org/wiki/Poisson_distribution Bivariate_Poisson_distribution

步骤:

1. 生成3个独立泊松变量Z_i，参数为lambda_i
2. 在Poi(lambda_i + lambda_3)
之后为i = 1,2生成两个P_i = Z_i + Z_3代码:

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)