8._总体分布的卡方检验法

总体分布的卡方检验法

前面讨论的参数检验问题是在总体分布类型已知的情况下，对其中的未知参数进行检验。在实际问题中，有时并不知道总体的具体分布情况，这时就需要从总体中抽取的样本对总体的分布进行推断，以判断总体服从何种分布，这类统计检验称为非参数检验。例如检验假设＂总体服从正态分布＂等。本节仅介绍 $\chi^2$ 拟合优度检验，又简称 $\chi^2$ 检验法，它是英国统计学家老皮尔逊（K．Pearson，1857－1936年）于1900年提出的，不少人把此项工作视为近代统计学的开端。 $\chi^2$ 检验法是在总体 $X$ 的分布未知时，根据来自总体的样本，检验总体分布的假设的一种检验方法．设 $X_1, X_2, \cdots, X_n$ 是来自总体 $F(x)$ 的样本，$F_0(x)$ 是理论分布，检验问题的原假设是

该分布的检验问题是根据样本的观察值 $x_1, x_2, \cdots, x_n$ 的数据判断是否与理论分布相合．样本量较大时，可以用 $\chi^2$ 拟合优度检验。这类问题可以分以下两种情况来讨论。

总体 $X$ 为离散型分布

设总体 $X$ 是取值为有限个或可列个 $b_1, b_2, \cdots$ 的离散随机变量，将相邻的某些 $b_i$ 合并为一类，且样本观察值 $x_1, x_2, \cdots, x_n$ 落入每一个类内的个数不小于 5 ，记 $B_1, B_2, \cdots, B_k$ 为 $b_1, b_2, \cdots$ 被分的 $k$ 个类，$n_i$ 为每一个 $B_i$ 内的个数．记

则假设

可以转化为如下假设

即

现对总体做了 $n$ 次观察，各类出现的观察频数分别为 $n_1, n_2, \cdots, n_k$ ，且

当 $H_0$ 为真时，则各概率 $p_i$ 与频率 $\frac{n_i}{n}$ 相差应该不大，或者各观察频数 $n_i$ 对期望频数 $n p_i$的偏差（ $n_i-n p_i$ ）不大．英国统计学家老皮尔逊提出了基于观察频数和期望频数之差的检验统计量

其中取偏差平方是为了把偏差积累起来，每项除以 $n p_i$ 是要求期望频数 $n p_i$ 较小时，偏差平方 $\left(n_i-n p_i\right)^2$ 更小才合理。并证明了下列结论。

定理 当 $n$ 充分大 $(n \geqslant 50)$ 时，则统计量 $\chi^2=\sum_{i=1}^k \frac{\left(n-n p_i\right)^2}{n p_i}$ 近似服从 $\chi^2(k-1)$ 分布． 证明 略。 根据该定理，对给定的显著性水平 $\alpha$ ，确定 $l$ 值，使

求得 $l=\chi_\alpha^2(k-1)$ ，所以拒绝域为

若由所给的样本值 $x_1, x_2, \cdots, x_n$ 算得统计量 $\chi^2$ 的实测值落入拒绝域，则拒绝原假设 $H_0$ ，否则就认为差异不显著而接受原假设 $H_0$ ．

例 将一颗骰子郑 120 次，所得数据为

问这颗骰子是否均匀，对称（取 $\alpha=0.05$ ）？ 解 若这颗骰子是均匀的，对称的，则 $1 \sim 6$ 点中每点出现的可能性相同，都为 $1 / 6$ ．如果用 $A_i$ 表示第 $i$ 点出现 $(i=1,2, \cdots, 6)$ ，则待检验假设 $H_0: P\left(A_i\right)=1 / 6, i=1,2 \cdots, 6$ 。

在 $H_0$ 成立的条件下，理论概率 $p_i=p\left(A_i\right)=1 / 6$ ，由 $n=120$ 得频率 $n p_i=20$ ． 计算结果如下．

因此分布不含未知参数，又 $k=6, \alpha=0.05$ ，查表得 $\chi_\alpha^2(k-1)=\chi_{0.05}^2(5)=11.071$ ． 由上表知 $\chi^2=\sum_{i=1}^6 \frac{\left(n_i-n p_i\right)^2}{n p_i}=4.8<11.071$ ，故接受 $H_0$ ，认为这颗骰子是均匀对称的．

总体 $X$ 为连续型分布

设 $X_1, X_2, \cdots, X_n$ 是来自总体 $X$ 的一个样本，其观察值为 $x_1, x_2, \cdots, x_n$ ，总体分布未知，现想用一个已知连续分布函数 $F_0(x)$ 去拟合这批数据，故需要对如下假设做出检验：

这类问题称为连续分布的拟合优度检验问题。具体步骤如下． （1）把 $X$ 的取值范围划分为 $k$ 个互不相交的小区间，记为 $A_1, A_2, \cdots, A_k$ ，如可取为

其中 $a_0<a_1<\cdots<a_{k-1}<a_k, a_0$ 可取 $-\infty, a_k$ 可取 $+\infty$ ；区间的划分视具体情况而定，使每个小区间所含样本值个数不小于 5 ，而区间个数 $k$ 不要太大，也不要太小； （2）当落入第 $i$ 个区间内时，就把它看成属于 $A_i$ 类，因此这 $k$ 个区间相当于 $k$ 类．在 $H_0$ 为真时，可算出总体 $X$ 的值落入第 $i$ 个小区间 $A_i$ 的概率

样本观测值落入这 $k$ 个区间的频数分别为 $n_1, n_2, \cdots, n_k$ ．接下来的步骤与离散型分布类似． 在对总体分布的假设检验中，有时只知道总体 $X$ 的分布函数的形式，但其中还含有末知参数，即分布函数为

其中 $\theta_1, \theta_2, \cdots, \theta_r$ 为未知参数．设 $X_1, X_2, \cdots, X_n$ 是取自总体 $X$ 的样本，现要用此样本来检验假设：

此类情况可按如下步骤进行检验： （1）利用样本 $X_1, X_2, \cdots, X_n$ 求出 $\theta_1, \theta_2, \cdots, \theta_r$ 的最大似然估计 $\hat{\theta}_1, \hat{\theta}_2, \cdots, \hat{\theta}_r$ ； （2）用 $\hat{\theta}_i$ 代替 $F\left(x, \theta_1, \theta_2, \cdots, \theta_r\right)$ 中的 $\theta_i(i=1,2, \cdots, r)$ ； （3）利用 $F\left(x, \hat{\theta}_1, \hat{\theta}_2, \cdots, \hat{\theta}_r\right)$ 计算 $p_i$ 的估计值 $\hat{p}_i(i=1,2, \cdots, k)$ ； （4）计算要检验的统计量

当 $n$ 充分大时，统计量 $\chi^2$ 近似服从 $\chi_\alpha^2(k-r-1)$ 分布，其中 $k$ 为互不相交的小区间的个数， $r$ 是被估计的参数的个数。 （5）对给定的显著性水平 $\alpha$ ，得拒绝域

注意：在使用皮尔逊 $\chi^2$ 检验法时，要求 $n \geqslant 50$ ，以及每个理论频数 $n p_i \geqslant 5(i=1, \cdots, k)$ ，否则应适当地合并相邻的小区间，使 $n p_i$ 满足要求。

例研究混凝土抗压强度的分布． 200 件混凝土制件的抗压强度以分组形式列出（如表 8．4．1 所示）．$n=\sum_{i=1}^6 n_i=200$ ．要求在给定的检验水平 $\alpha=0.05$ 下检验假设

表 8．4．1 200 件混凝土制件的抗压强度 解 原假设所定的正态分布的参数是未知的，需先求 $\mu$ 和 $\sigma^2$ 的极大似然估计值．由第 7章知，$\mu$ 和 $\sigma^2$ 的极大似然估计值为

设 $x_i^*$ 为第 $i$ 组的组中值，有

原假设 $H_0$ 改写成 $X$ 是正态 $N\left(221,12.33^2\right)$ 分布，计算每个区间的理论概率值

其中

为了计算统计量 $\chi^2$ 的值，把需要进行的计算列表，如表 8．4．2 所示．

从上面计算得出 $\chi^2$ 的观察值为 1.33 ．在检验水平 $a=0.05$ 下，查自由度 $m=6-2-1=3$ 的 $\chi^2$分布表，得到临界值 $\chi_{0.05}^2(3)=7.815$ ．由于 $\chi^2=1.33<7.815=\chi_{0.05}^2(3)$ ，不能拒绝原假设，所以认为混凝土制件的抗压强度的分布是正态分布 $N(221,152)$ 。