Python 中的二元泊松分布

Question

我想從雙變量 Possion 分布中抽取 N 次。 R中是否有類似於package bivpois的Python模塊？

In Python, I only know the libraries scipy.stats.poisson and numpy.random.possion which allow me to make draws from a univariate Poisson distribution depending on a single parameter lambda, but not from a bivariate or multivariate.

Answer 1

您可以自己輕松完成，因為我沒有看到任何內置方法：

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

腳步：

1. 生成 3 個帶參數 lambda_i 的獨立泊松變量 Z_i
2. 為 i = 1, 2 生成兩個 P_i = Z_i + Z_3，它遵循 Poi(lambda_i + lambda_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)