统计各个函数的耗时_数理统计第31讲(Kolmogorov- Smirnov检验，列联表独立性检验)...

6.4.2分布拟合的卡方检验（续）

例 测得200件混泥土制件的抗压强度，按区间分布如下:

试在

水平下检验接受抗压强度的分布服从正态分布.

解 这个例子中没有原始数据，我们就用组中值来代替原始数据.

即

.

我们已经计算过正态分布参数的极大似然估计为

参数的极大似然估计值分别为

进而得到分布函数的估计，而每个区间中概率的估计就是在两区间端点的分布函数值之差，即

列出计算检验统计量的表格

每一类的频数都大于等于5，不需要再进行合并。

现在

，临界值为

所以拒绝域为

，而

，因此保留原假设，认为服从正态.

---

对于分布拟合的卡方检验问题，我们想要检验的是总体的分布是否为已知的分布

但实际上我们只检验了总体分布落入给定区间的概率是否与已知分布相同，即

这样即使

是成立的，

也未必成立，特别是对于连续总体，我们仅仅只用了有限个值，就会导致第二类错误，因此

卡方拟合优度检验一般不用来检验连续总体！

• Kolmogorov- Smirnov检验

思想：用经验分布函数取近似总体分布

这样检验问题就转化成

那么当原假设成立的时候，两者之差应该比较小，所以我们取检验统计量为两者之差的最大值

在

比较大的时候我们拒绝原假设.

因为经验分布是阶梯函数，且已知分布

是单调增的，所以它们图象之间的关系分为四种情形

而无论是哪一种情形，在任何一个区间上，

的最大值都只可能出现在区间端点，因此在整个区间的上确界可以写成

个值得最大值

其拒绝域的形式为

，当

时，

可以通过书后附表13得到，当

时，计算出

，然后在附表14中找到

（或者与之最接近的），表中给出了对应的

，然后用

得到临界值.

例 检验下面10个数是否来自标准正态总体.

解 取检验统计量为

这里

，查附表13得临界值，进而拒绝域为

计算统计量的值如下表

得出此表的过程：

• 首先将数据进行排列，得到各个次序统计量的值
• 将各个次序统计量的值代入
  中（这里就是标准正态分布函数）
• 写出各个区间的起止端点
• 得出
  ，
  就是
  与相应区间端点之差的较大者
• 取
  中最大的作为检验统计量的值

这里

，没有落入拒绝域，所以保留原假设，认为这10个数据来自标准正态总体.

对这种方法进行推广就得到Smirnov检验：

Smirnov检验的问题是两个总体分布是否相同，那么自然就转化成两个样本的经验分布函数是否相同，即假设为

，取检验统计量为

也是当这个统计量的值比较大是拒绝原假设，拒绝域形式为

如果是单边问题

，那么检验统计量就不加绝对值，即

.

6.5列联表独立性检验（检验两个特性是否独立）

例 根据男女色盲人数

判断性别与色盲是否独立.

上面这个例子就是检验两个特性（性别与色盲）之间是否独立，其中性别分为男女两类，色盲分为正常与色盲两类.

一般而言，假设两个特性

分别可以分成

类，即

与

,其中属于

的个体数为

，根据样本可以给出

的列联表如下

要检验的是两个特性

是否独立.

下面将独立这个条件转化成具体的含义：

设

当两个特性独立时，根据定义，这些概率应该满足

如果

都是

已知的，类似分类数据检验，用实际频数近似理论频数.这样就得到了检验统计量为 ( 注意这里利用了

)

而当原假设成立的时候，这个统计量近似卡方分布

拒绝域为

如果

都是

未知的（如色盲的例子），那么我们有

个参数需要估计，但是这些参数并不完全独立，而被

限制，所以自由的参数有

个，而它们的极大似然估计分别为（

记住）

所以这时我们选取的检验统计量为

此时自由度为

，所以拒绝域为

继续看色盲的例子，选取检验统计量为

现

，设显著性水平为0.01，则

，拒绝域为

现利用

，构造理论频数的表格

然后得到检验统计量的计算表格

因为

，落入拒绝域，所以拒绝原假设，认为性别和色盲不独立.

解题步骤：

• 构造检验统计量
• 列出理论频数表格
• 列出检验统计量的计算表格
• 判断有无落入拒绝域

另外，在

都是

未知的情况下，下面这个等价计算式子在计算器上敲比较方便，考试时可以采用