我有一个X
的整数数(非常大)和一个概率p,我想用它从服从Poisson分布的s
中抽取一个样本X
(一个数)。例如,如果X = 10^8
和p=0.05
,我希望s是我们得到的头数。
我可以很容易地用随机二项式做这个:
s=np.random.binomial(n=X, p=p)
如何使用
random.poisson
应用相同的想法? 最佳答案
只需乘以p
和X
:
np.random.poisson(10**8 * 0.05)
得到大于10*8的概率在数值上为零。
pjs教授强调,我们把概率和数结合成一个速率,这个速率是泊松过程的参数。
更值得一提的是,对于如此大的数量,你会发现二项式和泊松的pmf彼此非常相似,而且(使用概率函数或工程师称之为“cdf”)是高斯函数。