使用scipy进行的Poisson分布拟合优度检验

时间:2019-05-01 16:57:22

标签: scipy poisson hypothesis-test

在测试泊松分布拟合优度时,我得到的结果与理论相反

从理论上我们知道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

0 个答案:

没有答案