在测试泊松分布拟合优度时,我得到的结果与理论相反
从理论上我们知道X〜Poisson(1)和Y〜Poisson(7),然后是X + Y〜Poisson(1 + 7)。当我尝试通过下面的代码使用scipy进行检查时,我得到的p值<0.05,这意味着拒绝原假设[X + Y遵循Poisson(8)]
import numpy as np
from scipy import stats
x = np.random.poisson(1.0,1000)
y = np.random.poisson(7.0,1000)
print stats.chisquare(x+y,stats.poisson.rvs(8,size=1000),ddof=2)
实际结果:Power_divergenceResult(statistic = 2748.7129256364551,pvalue = 1.525087375864822e-163)。 p值<0.05
预期结果:p值>> 0.5