内容正文:
《概率论与数理统计》
第四节 总体分布的假设检验
----拟合优度检验
在前面的课程中,我们已经了解了假设检验的基本思想,并讨论了当总体分布为正态时,关于其中未知参数的假设检验问题 .
然而可能遇到这样的情形,总体服从何种理论分布并不知道,要求我们直接对总体分布提出一个假设 .
例如,从1500到1931年的432年间,每年爆发战争的次数可以看作一个随机变量,椐统计,这432年间共爆发了299次战争,具体数据如下:
战争次数X
0
1
2
3
4
223
142
48
15
4
发生 X次战争的年数
在概率论中,大家对泊松分布产生的一般条件已有所了解,容易想到,每年爆发战争的次数,可以用一个泊松随机变量来近似描述.也就是说,我们可以假设每年爆发战争次数分布X近似泊松分布.
上面的数据能否证实X 具有
泊松分布的假设是正确的?
现在的问题是:
又如,某钟表厂对生产的钟进行精确性检查,抽取100个钟作试验,拨准后隔24小时以后进行检查,将每个钟的误差(快或慢)按秒记录下来.
问该厂生产的钟的误差是否服从正态
分布?
再如,某工厂制造一批骰子,声称它是均匀的.
为检验骰子是否均匀,要把骰子实地投掷若干次,统计各点出现的频率与1/6的差距.
也就是说,在投掷中,出现1点,2点,…,6点的概率都应是1/6.
得到的数据能否说明“骰子均匀”的假设是可信的?
问题是:
K.皮尔逊
这是一项很重要的工作,不少人把它视为近代统计学的开端.
2-检验法.
解决这类问题的工具是英国统计学家K.皮尔逊
在1900年发表的一篇文章中引进的所谓
2检验法是在总体X 的分布未知时,根据来自总体的样本,检验关于总体分布的假设的一种检验方法.
H0:总体X的分布函数为F(x)
然后根据样本的经验分布和所假设的理论分布之间的吻合程度来决定是否接受原假设.
这种检验通常称作拟合优度检验,它是一种非参数检验.
使用2-检验对总体分布进行检验时,我们先
提出原假设:
似然估计法估计参数,然后作检验.
类型已知,但其参数未知,这时需要先用极大
分布拟合的2-检验的基本原理和步骤如下:
在用2-检验法检验假设H0时,若在H0下分布
3.根据所假设的理论分布,可以算出总体X的值落入每个Ai的概率pi,于是npi就是落入Ai的样本值的理论频数.
1. 将总体X的取值范围分成k个互不重迭的小区间,记作A1, A2, …, Ak .
2.把落入第i个小区间Ai的样本值的个数记作ni , 称为实测频数. 所有实测频数之和
n1+ n2+ …+ nk等于样本容量n.
标志着经验分布与理论分布之间的差异的大小.
皮尔逊引进如下统计量表示经验分布
与理论分布之间的差异:
统计量2的分布是什么?
在理论分布
已知的条件下,
npi是常量
实测频数
理论频数
皮尔逊证明了如下定理:
若原假设中的理论分布F(x)已经完全给定,那么当n 时,统计量
如果理论分布F(x)中有r个未知参数需用相应的估计量来代替,那么当n 时,统计量2 的分布渐近服从 (k-r-1)个自由度的2分布.
的分布渐近k-1个自由度的2 分布.
为了便于理解,我们对定理作一点直观的说明.
在理论分布F(x)完全给定的情况下,每个pi 都是确定的常数. 由棣莫佛-拉普拉斯中心极限定理,当n充分大时,实测频数 fi 渐近正态.
因此
是k个近似正态的变量的平方和.
这些变量之间存在着一个制约关系:
故统计量2渐近服从(k-1)个自由度的2分布。
在F(x)尚未完全给定的情况下,每个未知参数用相应的估计量代替,就相当于增加一个制约条件,因此,自由度也随之减少一个.
若有r个未知参数需用相应的估计量来代替,自由度就减少r个.
故统计量2渐近服从(k-r-1)个自由度的2分布。
如果根据所给的样本值 X1,X2, …,Xn算得统计量2的实测值落入拒绝域,则拒绝原假设,否则就认为差异不显著而接受原假设.
得拒绝域:
(不需估计参数)
(估计r 个参数)
根据这个定理,对给定的显著性水平,
查2 分布表可得临界值2 ,使得
皮尔逊定理是在n无限增大时推导出来的,因而在使用时要注意n要足够大,以及npi 不太小这两个条件.
根据计算实践,要求n不小于50,以及npi 都不小于 5. 否则应适当合并区间,使npi满足这个要求 .
让我们回到开始的一个例子,检验每年爆发战争次数分布是否服从泊松分布.
将有关计算结果列表如下:
提出假设H0: X服从参数为的泊松分布
根据观察结果,得参数的极大似然估计为
按参数为0