我有一个X的整数数(非常大)和一个概率p,我想用它从服从Poisson分布的s中抽取一个样本X(一个数)。例如,如果X = 10^8p=0.05,我希望s是我们得到的头数。
我可以很容易地用随机二项式做这个:

s=np.random.binomial(n=X, p=p)

如何使用random.poisson应用相同的想法?

最佳答案

只需乘以pX

np.random.poisson(10**8 * 0.05)

得到大于10*8的概率在数值上为零。
pjs教授强调,我们把概率和数结合成一个速率,这个速率是泊松过程的参数。
更值得一提的是,对于如此大的数量,你会发现二项式和泊松的pmf彼此非常相似,而且(使用概率函数或工程师称之为“cdf”)是高斯函数。

08-24 13:19