泊松分布

泊松分布

import numpy as np
from scipy import stats
import matplotlib.pyplot as plt

rate = 2   #均值
n = np.arange(0,11)
n
y = stats.poisson.pmf(n,rate)
y
sum = stats.poisson.cdf(10,n,y) 
sum
plt.plot(y,'o-')
plt.show()

data = stats.poisson.rvs(mu = 2,loc = 0,size = 1000)
print ("Mean:%g" % np.mean(data))
print ("SD:%g" % np.std(data,ddof = 1))
plt.hist(data,bins = 10,normed = True)
plt.xlabel("numbers of accidents") 
plt.ylabel("simulating poisson random variable")
plt.show()

运行结果（举例）：