第八章 假设检验
第八章 假设检验
§1 假设检验
统计推断的另一类重要问题是假设检验问题。在总体的分布函数完全未知或只知其形式、但不知其参数的情况,为了推断总体的某些未知特性,提出某些关于总体的假设。例如,提出总体服从泊松分布的假设,又如,对于正态总体提出数学期望等于 $ \mu_{0} $ 的假设等。我们要根据样本对所提出的假设作出是接受,还是拒绝的决策。假设检验是作出这一决策的过程。这里,先结合例子来说明假设检验的基本思想和做法。
例 1 某车间用一台包装机包装葡萄糖. 袋装糖的净重是一个随机变量, 它服从正态分布. 当机器正常时, 其均值为 0.5 kg, 标准差为 0.015 kg. 某日开工后为检验包装机是否正常, 随机地抽取它所包装的糖 9 袋, 称得净重为 (kg)
0.497 0.506 0.518 0.524 0.498 0.511 0.520 0.515 0.512
问机器是否正常?
以 $ \mu $, $ \sigma $ 分别表示这一天袋装糖的净重总体 X 的均值和标准差。由于长期实践表明标准差比较稳定,我们就设 $ \sigma=0.015 $。于是 $ X\sim N(\mu,0.015^{2}) $,这里 $ \mu $ 未知。问题是根据样本值来判断 $ \mu=0.5 $ 还是 $ \mu\neq0.5 $。为此,我们提出两个相互对立的假设
$$ H_{0}:\mu=\mu_{0}=0.5 $$
和
$$ H_{1}:\mu\neq\mu_{0}. $$
然后,我们给出一个合理的法则,根据这一法则,利用已知样本作出决策是接受假设 $ H_{0} $(即拒绝假设 $ H_{1} $),还是拒绝假设 $ H_{0} $(即接受假设 $ H_{1} $)。如果作出的决策是接受 $ H_{0} $,则认为 $ \mu=\mu_{0} $,即认为机器工作是正常的,否则,则认为是不正常的。
由于要检验的假设涉及总体均值 $ \mu $,故首先想到是否可借助样本均值 $ \overline{X} $ 这一统计量来进行判断。我们知道, $ \overline{X} $ 是 $ \mu $ 的无偏估计, $ \overline{X} $ 的观察值 $ \overline{x} $ 的大小在一定程度上反映 $ \mu $ 的大小。因此,如果假设 $ H_{0} $ 为真,则观察值 $ \overline{x} $ 与 $ \mu_{0} $ 的偏差 $ \left|\overline{x}-\mu_{0}\right| $ 一般不应太大。若 $ \left|\overline{x}-\mu_{0}\right| $ 过分大,我们就怀疑假设 $ H_{0} $ 的正确性而拒绝 $ H_{0} $,
并考虑到当 $ H_{0} $ 为真时 $ \frac{\overline{X}-\mu_{0}}{\sigma/\sqrt{n}}\sim N(0,1) $. 而衡量 $ \left|\overline{x}-\mu_{0}\right| $ 的大小可归结为衡量 $ \frac{\left|\overline{x}-\mu_{0}\right|}{\sigma/\sqrt{n}} $ 的大小. 基于上面的想法, 我们可适当选定一正数 k, 使当观察值 $ \overline{x} $ 满足 $ \frac{\left|\overline{x}-\mu_{0}\right|}{\sigma/\sqrt{n}}\geqslant k $ 时就拒绝假设 $ H_{0} $, 反之, 若 $ \frac{\left|\overline{x}-\mu_{0}\right|}{\sigma/\sqrt{n}} 然而,由于作出决策的依据是一个样本,当实际上 $ H_{0} $为真时仍可能作出拒绝 $ H_{0} $的决策(这种可能性是无法消除的),这是一种错误,犯这种错误的概率记为 P{当 $ H_0 $ 为真拒绝 $ H_0\} $ 或 $ P_{\mu_0} $ {拒绝 $ H_0\} $ 或 $ P_{\mu \in H_0} $ {拒绝 $ H_0 $}。 记号 $ P_{\mu_0}\{\cdot\} $ 表示参数 $ \mu $ 取 $ \mu_0 $ 时事件 $ \{\cdot\} $ 的概率, $ P_{\mu \in H_0}\{\cdot\} $ 表示 $ \mu $ 取 $ H_0 $ 规定的值时事件 $ \{\cdot\} $ 的概率。我们无法排除犯这类错误的可能性,因此自然希望将犯这类错误的概率控制在一定限度之内,即给出一个较小的数 $ \alpha $ ( $ 0 < \alpha < 1 $),使犯这类错误的概率不超过 $ \alpha $,即使得 $$ P\{ 当 H_{0} 为真拒绝 H_{0}\}\leq\alpha $$ 为了确定常数 k,我们考虑统计量 $ \frac{X - \mu_{0}}{\sigma/\sqrt{n}} $。由于只允许犯这类错误的概率最大为 $ \alpha $,令(1.1)式右端取等号,即令 $$ P\{当 H_{0} 为真拒绝 H_{0}\}=P_{\mu_{0}}\{|\frac{\overline{X}-\mu_{0}}{\sigma/\sqrt{n}}|\geqslant k\}=\alpha $$ 由于当 $ H_{0} $ 为真时, $ Z=\frac{X-\mu_{0}}{\sigma/\sqrt{n}}\sim N(0,1) $,由标准正态分布分位点的定义得(如图 8-1). $$ k=z_{a/2}. $$ 因而,若 Z 的观察值满足 $$ |z|=\left|\frac{\overline{x}-\mu_{0}}{\sigma/\sqrt{n}}\right|\geqslant k=z_{a/2}, $$ 则拒绝 $ H_{0} $,而若 $$ |z|=\left|\frac{\overline{x}-\mu_{0}}{\sigma/\sqrt{n}}\right| 则接受 $ H_{0} $ 例如,在本例中取 $ \alpha=0.05 $,则有 $ k=z_{0.05/2}=z_{0.025}=1.96 $,又已知 n=9, $ \sigma=0.015 $,再由样本算得 $ \overline{x}=0.511 $,即有 $$ \left|\frac{\bar{x}-\mu_{0}}{\sigma/\sqrt{n}}\right|=2,2>1.96, $$ 于是拒绝 $ H_{0} $,认为这天包装机工作不正常。 上例中所采用的检验法则是符合实际推断原理的。因通常 $ \alpha $ 总是取得较小,一般取 $ \alpha=0.01,0.05 $。因而若 $ H_{0} $ 为真,即当 $ \mu=\mu_{0} $ 时, $ \{|\frac{\overline{X}-\mu_{0}}{\sigma/\sqrt{n}}|\geqslant z_{\alpha/2}\} $ 是一个小概率事件,根据实际推断原理,就可以认为,如果 $ H_{0} $ 为真,则由一次试验得到的观察值 $ \overline{x} $,满足不等式 $ \left|\frac{\overline{x}-\mu_{0}}{\sigma/\sqrt{n}}\right|\geqslant z_{\alpha/2} $ 几乎是不会发生的。现在在一次观察中竟然出现了满足 $ \left|\frac{\overline{x}-\mu_{0}}{\sigma/\sqrt{n}}\right|\geqslant z_{\alpha/2} $ 的 $ \overline{x} $,则我们有理由怀疑原来的假设 $ H_{0} $ 的正确性,因而拒绝 $ H_{0} $。若出现的观察值 $ \overline{x} $ 满足 $ \left|\frac{\overline{x}-\mu_{0}}{\sigma/\sqrt{n}}\right| 在上例的做法中,我们看到当样本容量固定时,选定 $ \alpha $后,数k就可以确定,然后按照统计量 $ Z=\frac{\overline{X}-\mu_{0}}{\sigma/\sqrt{n}} $的观察值的绝对值 $ \left|z\right| $大于等于k还是小于k来作出决策。数k是检验上述假设的一个门槛值。如果 $ \left|z\right|=\left|\frac{\overline{x}-\mu_{0}}{\sigma/\sqrt{n}}\right|\geqslant k $,则称 $ \overline{x} $与 $ \mu_{0} $的差异是显著的,这时拒绝 $ H_{0} $;反之,如果 $ \left|z\right|=\left|\frac{\overline{x}-\mu_{0}}{\sigma/\sqrt{n}}\right| 统计量 $ Z=\frac{\overline{X}-\mu_{0}}{\sigma/\sqrt{n}} $ 称为检验统计量. 前面的检验问题通常叙述成:在显著性水平 $ \alpha $下,检验假设 $$ H_{0}:\mu=\mu_{0},\quad H_{1}:\mu\neq\mu_{0}. $$ 也常说成“在显著性水平 $ \alpha $下,针对 $ H_{1} $检验 $ H_{0} $”。 $ H_{0} $称为原假设或零假设, $ H_{1} $称为备择假设(意指在原假设被拒绝后可供选择的假设)。我们要进行的工作是,根据样本,按上述检验方法作出决策在 $ H_{0} $与 $ H_{1} $两者之间接受其一。 当检验统计量取某个区域 C 中的值时,我们拒绝原假设 $ H_{0} $,则称区域 C 为拒绝域,拒绝域的边界点称为临界点。如在上例中拒绝域为 $ |z| \geqslant z_{\alpha/2} $,而 $ z = -z_{\alpha/2} $, $ z = z_{\alpha/2} $ 为临界点。 由于检验法则是根据样本作出的,总有可能作出错误的决策。如上面所说的那样,在假设 $ H_{0} $ 实际上为真时,我们可能犯拒绝 $ H_{0} $ 的错误,称这类“弃真”的错误为第Ⅰ类错误。又当 $ H_{0} $ 实际上不真时,我们也有可能接受 $ H_{0} $。称这类“取伪”的错误为第Ⅱ类错误。犯第Ⅱ类错误的概率记为 P{当 $ H_0 $ 不真接受 $ H_0\} $ 或 $ P_{\mu \in H_1} $ {接受 $ H_0 $}。 为此,在确定检验法则时,我们应尽可能使犯两类错误的概率都较小。但是,进一步讨论可知,一般来说,当样本容量固定时,若减少犯一类错误的概率,则犯另一类错误的概率往往增大。若要使犯两类错误的概率都减小,除非增加样本容量。在给定样本容量的情况下,一般来说,我们总是控制犯第Ⅰ类错误的概率,使它不大于 $ \alpha $。 $ \alpha $的大小视具体情况而定,通常 $ \alpha $取0.1,0.05,0.01,0.005等值。这种只对犯第Ⅰ类错误的概率加以控制,而不考虑犯第Ⅱ类错误的概率的检验,称为显著性检验。 形如(1.2)式中的备择假设 $ H_{1} $,表示 $ \mu $可能大于 $ \mu_{0} $,也可能小于 $ \mu_{0} $,称为双边备择假设,而称形如(1.2)的假设检验为双边假设检验. 有时,我们只关心总体均值是否增大,例如,试验新工艺以提高材料的强度。这时,所考虑的总体的均值应该越大越好。如果我们能判断在新工艺下总体均值较以往正常生产的大,则可考虑采用新工艺。此时,我们需要检验假设 $$ H_{0}:\mu\leqslant\mu_{0},\quad H_{1}:\mu>\mu_{0}. $$ 形如(1.3)的假设检验,称为右边检验.类似地,有时我们需要检验假设 $$ H_{0}:\mu\geqslant\mu_{0},\quad H_{1}:\mu<\mu_{0}. $$ 形如(1.4)的假设检验,称为左边检验.右边检验和左边检验统称为单边检验. 下面来讨论单边检验的拒绝域. 设总体 $ X \sim N(\mu, \sigma^{2}) $, $ \mu $ 未知、 $ \sigma $ 为已知, $ X_{1}, X_{2}, \cdots, X_{n} $ 是来自 X 的样本。给定显著性水平 $ \alpha $。我们来求检验问题(1.3) $$ H_{0}:\mu\leqslant\mu_{0},\quad H_{1}:\mu>\mu_{0} $$ 的拒绝域. 因 $ H_{0} $ 中的全部 $ \mu $ 都比 $ H_{1} $ 中的 $ \mu $ 要小,当 $ H_{1} $ 为真时,观察值 $ \bar{x} $ 往往偏大,因此,拒绝域的形式为 $$ \overline{x}\geqslant k\quad(k 是某一正常数 ). $$ 下面来确定常数 k,其做法与例 1 中的做法类似。 $$ H_{0} $$ $$ \begin{aligned}H_{0}&=P_{\mu\in H_{0}}\{\overline{X}\geqslant k\}\\&=P_{\mu\leqslant\mu_{0}}\{\frac{\overline{X}-\mu_{0}}{\sigma/\sqrt{n}}\geqslant\frac{k-\mu_{0}}{\sigma/\sqrt{n}}\}\\&\leqslant P_{\mu\leqslant\mu_{0}}\{\frac{\overline{X}-\mu}{\sigma/\sqrt{n}}\geqslant\frac{k-\mu_{0}}{\sigma/\sqrt{n}}\}\\ \end{aligned} $$ $ \text{q} $ $$ P_{\mu\leq\mu_{0}}\{\frac{\overline{X}-\mu}{\sigma/\sqrt{n}}\geq\frac{k-\mu_{0}}{\sigma/\sqrt{n}}\}=\alpha. $$ 由于 $ \frac{\overline{X}-\mu}{\sigma/\sqrt{n}}\sim N(0,1) $,由(1.5)得到 $ \frac{k-\mu_{0}}{\sigma/\sqrt{n}}=z_{\alpha} $(如图8-2), $ k=\mu_{0}+\frac{\sigma}{\sqrt{n}}z_{\alpha} $,即得检验问题 (1.3)的拒绝域为 $$ \overline{x}\geqslant\mu_{0}+\frac{\sigma}{\sqrt{n}}z_{a}, $$ 即 $$ z=\frac{\overline{x}-\mu_{0}}{\sigma/\sqrt{n}}\geq z_{a}. $$ 类似地,可得左边检验问题(1.4) $$ H_{0}:\mu\geqslant\mu_{0},H_{1}:\mu<\mu_{0} $$ 的拒绝域为 $$ z=\frac{\overline{x}-\mu_{0}}{\sigma/\sqrt{n}}\leqslant-z_{a}. $$ 例2 公司从生产商购买牛奶。公司怀疑生产商在牛奶中掺水以谋利。通过测定牛奶的冰点,可以检验出牛奶是否掺水。天然牛奶的冰点温度近似服从正态分布,均值 $ \mu_{0} = -0.545\ ^{\circ}C $,标准差 $ \sigma = 0.008\ ^{\circ}C $。牛奶掺水可使冰点温度升高而接近于水的冰点温度( $ 0\ ^{\circ}C $)。测得生产商提交的5批牛奶的冰点温度,其均值为 $ \overline{x} = -0.535\ ^{\circ}C $,问是否可以认为生产商在牛奶中掺了水?取 $ \alpha = 0.05 $。 解 按题意需检验假设 $$ \begin{aligned}&H_{0}:\mu\leqslant\mu_{0}=-0.545( 即设牛奶未掺水 ),\\&H_{1}:\mu>\mu_{0}( 即设牛奶已掺水 )\\ \end{aligned} $$ 这是右边检验问题,其拒绝域如(1.6)式所示,即为 $$ z=\frac{\overline{x}-\mu_{0}}{\sigma/\sqrt{n}}\geqslant z_{0.05}=1.645. $$ 现在 $ z=\frac{-0.535-(-0.545)}{0.008/\sqrt{5}}=2.7951>1.645 $,z 的值落在拒绝域中,所以我们在显著性水平 $ \alpha=0.05 $ 下拒绝 $ H_{0} $,即认为牛奶商在牛奶中掺了水。 综上所述,可得处理参数的假设检验问题的步骤如下: $ 1^{\circ} $ 根据实际问题的要求,提出原假设 $ H_{0} $及备择假设 $ H_{1} $ $ 2^{\circ} $ 给定显著性水平 $ \alpha $ 以及样本容量 n; $ 3^{\circ} $ 确定检验统计量以及拒绝域的形式; $ 4^{\circ} $ 按 P{当 $ H_{0} $ 为真拒绝 $ H_{0} $}≤α求出拒绝域; $ 5^{\circ} $ 取样,根据样本观察值作出决策,是接受 $ H_{0} $ 还是拒绝 $ H_{0} $ 下面我们只讨论正态总体参数的假设检验问题. 在§1中已讨论过正态总体 $ N(\mu,\sigma^{2}) $当 $ \sigma^{2} $已知时关于 $ \mu $的检验问题(1.2),(1.3),(1.4).在这些检验问题中,我们都是利用统计量 $ Z=\frac{\overline{X}-\mu_{0}}{\sigma/\sqrt{n}} $来确定拒绝域的.这种检验法常称为Z检验法. 设总体 $ X \sim N(\mu, \sigma^{2}) $,其中 $ \mu, \sigma^{2} $ 未知,我们来求检验问题 $$ H_{0}:\mu=\mu_{0},\quad H_{1}:\mu\neq\mu_{0} $$ 的拒绝域(显著性水平为 $ \alpha $). 设 $ X_{1}, X_{2}, \cdots, X_{n} $ 是来自总体 X 的样本. 由于 $ \sigma^{2} $ 未知, 现在不能利用 $ \frac{\overline{X}-\mu_{0}}{\sigma/\sqrt{n}} $ 来确定拒绝域了. 注意到 $ S^{2} $ 是 $ \sigma^{2} $ 的无偏估计, 我们用 S 来代替 $ \sigma $, 采用 $$ t=\frac{\overline{X}-\mu_{0}}{S/\sqrt{n}} $$ 作为检验统计量. 当观察值 $ \left|t\right|=\left|\frac{\overline{x}-\mu_{0}}{s/\sqrt{n}}\right| $ 过分大时就拒绝 $ H_{0} $,拒绝域的形式为 $$ |t|=\left|\frac{\overline{x}-\mu_{0}}{s/\sqrt{n}}\right|\geqslant k. $$ 由第六章§3定理三知,当 $ H_{0} $为真时, $ \frac{X-\mu_{0}}{S/\sqrt{n}}\sim t(n-1) $,故由 $$ P\{ 当 H_{0} 为真拒绝 H_{0}\}=P_{\mu_{0}}\{|\frac{\overline{X}-\mu_{0}}{S/\sqrt{n}}|\geqslant k\}=\alpha $$ 得 $ k=t_{\alpha/2}(n-1) $,即得拒绝域为 $$ |t|=\left|\frac{\overline{x}-\mu_{0}}{s/\sqrt{n}}\right|\geqslant t_{a/2}(n-1). $$ 对于正态总体 $ N(\mu, \sigma^{2}) $,当 $ \sigma^{2} $ 未知时,关于 $ \mu $ 的单边检验的拒绝域在表 8.1 中给出. 上述利用 t 统计量得出的检验法称为 t 检验法. 在实际中,正态总体的方差常为未知,所以我们常用t检验法来检验关于正态总体均值的检验问题. 例1 某种元件的寿命 X(以 h 计)服从正态分布 $ N(\mu, \sigma^{2}) $, $ \mu, \sigma^{2} $ 均未知。现测得16只元件的寿命如下: $$ \begin{array}{l l l l l l}{159}&{280}&{101}&{212}&{224}&{379}&{179}&{264}\\ {222}&{362}&{168}&{250}&{149}&{260}&{485}&{170}\end{array} $$ 问是否有理由认为元件的平均寿命大于 $ 225 \, h $? 解 按题意需检验 $$ H_{0}:\mu\leqslant\mu_{0}=225,\quad H_{1}:\mu>225. $$ 取 $ \alpha=0.05 $。由表8.1知此检验问题的拒绝域为 $$ t=\frac{\overline{x}-\mu_{0}}{s/\sqrt{n}}\geqslant t_{a}(n-1). $$ 现在 n=16, $ t_{0.05}(15)=1.7531 $. 又算得 $ \overline{x}=241.5, s=98.7259 $, 即有 $$ t=\frac{\overline{x}-\mu_{0}}{s/\sqrt{n}}=0.6685<1.7531. $$ t 没有落在拒绝域中,故接受 $ H_{0} $,即认为元件的平均寿命不大于 225 h. ☐ 我们还可以用 t 检验法检验具有相同方差的两正态总体均值差的假设. 设 $ X_{1}, X_{2}, \cdots, X_{n_{1}} $ 是来自正态总体 $ N(\mu_{1}, \sigma^{2}) $ 的样本, $ Y_{1}, Y_{2}, \cdots, Y_{n_{2}} $ 是来自正态总体 $ N(\mu_{2}, \sigma^{2}) $ 的样本, 且设两样本独立. 又分别记它们的样本均值为 $ \overline{X}, \overline{Y} $, 记样本方差为 $ S_{1}^{2}, S_{2}^{2} $. 设 $ \mu_{1}, \mu_{2}, \sigma^{2} $ 均为未知, 要特别引起注意的是, 在这里假设两总体的方差是相等的. 现在来求检验问题: $$ H_{0}:\mu_{1}-\mu_{2}=\delta,\quad H_{1}:\mu_{1}-\mu_{2}\neq\delta $$ ( $ \delta $为已知常数)的拒绝域.取显著性水平为 $ \alpha $. 引用下述t统计量作为检验统计量: $$ t=\frac{(\overline{X}-\overline{Y})-\delta}{S_{w}\sqrt{\frac{1}{n_{1}}+\frac{1}{n_{2}}}}, $$ 其中 $$ S_{w}^{2}=\frac{(n_{1}-1)S_{1}^{2}+(n_{2}-1)S_{2}^{2}}{n_{1}+n_{2}-2},S_{w}=\sqrt{S_{w}^{2}}. $$ 当 $ H_{0} $ 为真时,由第六章 §3 定理四知 $ t \sim t(n_{1} + n_{2} - 2) $. 与单个总体的 t 检验法相仿,其拒绝域的形式为 $$ \left|\frac{(\overline{x}-\overline{y})-\delta}{s_{w}\sqrt{\frac{1}{n_{1}}+\frac{1}{n_{2}}}}\right|\geqslant k. $$ 由 $$ P\{ 当 H_{0} 为真拒绝 H_{0}\}=P_{\mu_{1}-\mu_{2}=\delta}\{|\frac{(\overline{X}-\overline{Y})-\delta}{S_{w}\sqrt{\frac{1}{n_{1}}+\frac{1}{n_{2}}}}|\geqslant k\}=\alpha $$ 可得 $ k=t_{a/2}(n_{1}+n_{2}-2) $. 于是得拒绝域为 $$ |t|=\frac{|(\overline{x}-\overline{y})-\delta|}{s_{w}\sqrt{\frac{1}{n_{1}}+\frac{1}{n_{2}}}}\geqslant t_{a/2}(n_{1}+n_{2}-2). $$ 关于均值差的两个单边检验问题的拒绝域在表8.1中给出. 常用的是 $ \delta=0 $的情况. 当两个正态总体的方差均为已知(不一定相等)时,我们可用Z检验法来检验两正态总体均值差的假设问题,见表8.1. 例2 用两种方法(A 和 B)测定冰自-0.72 ℃转变为0 ℃的水的融化热(以cal/g计).测得以下的数据: $$ \begin{aligned}& 方法 A:\ 79.98\quad80.04\quad80.02\quad80.04\quad80.03\quad80.03\\&\quad80.04\quad79.97\quad80.05\quad80.03\quad80.02\quad80.00\quad80.02\end{aligned} $$ 方法 B:80.02 79.94 79.98 79.97 79.97 80.03 79.95 78.97 设这两个样本相互独立,且分别来自正态总体 $ N(\mu_{1}, \sigma^{2}) $ 和 $ N(\mu_{2}, \sigma^{2}) $, $ \mu_{1}, \mu_{2}, \sigma^{2} $ 均未知。试检验假设(取显著性水平 $ \alpha=0.05 $) $$ H_{0}:\mu_{1}-\mu_{2}{\leqslant}0,H_{1}:\mu_{1}-\mu_{2}>0. $$ 解 分别画出对应于方法 A 和方法 B 的数据的箱线图,如图 8-3. 这两种方法所得的结果是有明显差异的,现在来检验上述假设. $$ n_{1}=13,\quad\overline{x}_{A}=80.02,\quad s_{A}^{2}=0.024^{2} $$ $$ n_{2}=8,\ \overline{x}_{B}=79.98,\ s_{B}^{2}=0.03^{2} $$ $$ s_{w}^{2}=\frac{12\times s_{A}^{2}+7\times s_{B}^{2}}{19}=0.000\ 717\ 8. $$ $$ t=\frac{\overline{x}_{A}-\overline{x}_{B}}{s_{w}\sqrt{1/13+1/8}}=3.33>t_{0.05}(13+8-2)=1.7291. $$ 故拒绝 $ H_{0} $,认为方法A比方法B测得的融化热要大. 有时为了比较两种产品、两种仪器、两种方法等的差异,我们常在相同的条件下做对比试验,得到一批成对的观察值。然后分析观察数据作出推断。这种方 法常称为逐对比较法. 例 3 有两台光谱仪 $ I_{x}, I_{y} $,用来测量材料中某种金属的含量,为鉴定它们的测量结果有无显著的差异,制备了 9 件试块(它们的成分、金属含量、均匀性等均各不相同),现在分别用这两台仪器对每一试块测量一次,得到 9 对观察值如下. 问能否认为这两台仪器的测量结果有显著的差异(取 $ \alpha=0.01 $)? 解 本题中的数据是成对的,即对同一试块测出一对数据。我们看到一对与另一对之间的差异是由各种因素,如材料成分、金属含量、均匀性等因素引起的。由于各试块的特性有广泛的差别,就不能将仪器 $ I_{x} $ 对9个试块的测量结果(即表中第一行)看成是同分布随机变量的观察值。因而表中第一行不能看成是一个样本的样本值。同样,表中第二行也不能看成是一个样本的样本值。再者,对于每一对数据而言,它们是同一试块用不同仪器 $ I_{x}, I_{y} $ 测得的结果,因此,它们不是两个独立的随机变量的观察值。综上所述,我们不能用表8.1中第4栏的检验法来作检验。而同一对中两个数据的差异则可看成是仅由这两台仪器性能的差异所引起的,这样,局限于各对中两个数据来比较就能排除种种其他因素,而只考虑单独由仪器的性能所产生的影响。从而能比较这两台仪器的测量结果是否有显著的差异。 一般,设有 n 对相互独立的观察结果: $ (X_{1},Y_{1}),(X_{2},Y_{2}),\cdots,(X_{n},Y_{n}) $,令 $ D_{1}=X_{1}-Y_{1},D_{2}=X_{2}-Y_{2},\cdots,D_{n}=X_{n}-Y_{n} $,则 $ D_{1},D_{2},\cdots,D_{n} $ 相互独立。又由于 $ D_{1},D_{2},\cdots,D_{n} $ 是由同一因素所引起的,可认为它们服从同一分布。今假设 $ D_{i}\sim N(\mu_{D},\sigma_{D}^{2}),i=1,2,\cdots,n $。这就是说 $ D_{1},D_{2},\cdots,D_{n} $ 构成正态总体 $ N(\mu_{D},\sigma_{D}^{2}) $ 的一个样本,其中 $ \mu_{D},\sigma_{D}^{2} $ 未知。我们需要基于这一样本检验假设: (1) $ H_{0}:\mu_{D}=0 $, $ H_{1}:\mu_{D}\neq0 $; (2) $ H_{0}: \mu_{D} \leqslant 0 $, $ H_{1}: \mu_{D} > 0 $; (3) $ H_{0}: \mu_{D} \geqslant 0 $, $ H_{1}: \mu_{D} < 0 $. 分别记 $ D_{1}, D_{2}, \cdots, D_{n} $ 的样本均值和样本方差的观察值为 $ \overline{d}, s_{D}^{2} $,按表 8.1 第 2 栏中关于单个正态总体均值的 t 检验.知检验问题(1),(2),(3)的拒绝域分别为(显著性水平为 $ \alpha $): $$ |t|=\left|\frac{\overline{d}}{s_{D}/\sqrt{n}}\right|\geqslant t_{a/2}(n-1), $$ $$ t=\frac{\overline{d}}{s_{D}/\sqrt{n}}\geqslant t_{a}(n-1), $$ $$ t=\frac{\overline{d}}{s_{D}/\sqrt{n}}\leqslant-t_{a}(n-1). $$ 现在回过来讨论本例的检验问题.先作出同一试块分别由仪器 $ I_{x} $, $ I_{y} $ 测得的结果之差,列于上表的第三行.按题意需检验假设 $$ H_{0}:\mu_{D}=0,\ H_{1}:\mu_{D}\neq0. $$ 现在 n=9, $ t_{\alpha/2}(8)=t_{0.005}(8)=3.3554 $ 即知拒绝域为 $$ |t|=\left|\frac{\overline{d}}{s_{D}/\sqrt{n}}\right|\geqslant3.3554. $$ 由观察值得 $ \overline{d}=0.06,s_{D}=0.1227,|t|=\frac{0.06}{0.1227/\sqrt{9}}=1.467<3.3554 $现 $ \left|t\right| $的值不落在拒绝域内,故接受 $ H_{0} $,认为两台仪器的测量结果并无显著差异. 例4 做以下的实验以比较人对红光或绿光的反应时间(以s计)。实验在点亮红光或绿光的同时,启动计时器,要求受试者见到红光或绿光点亮时,就按下按钮,切断计时器,这就能测得反应时间。测量的结果如下表: 设 $ D_{i}=X_{i}-Y_{i}\quad(i=1,2,\cdots,8) $ 是来自正态总体 $ N(\mu_{D},\sigma_{D}^{2}) $ 的样本, $ \mu_{D},\sigma_{D}^{2} $ 均未知。试检验假设(取显著性水平 $ \alpha=0.05 $) $$ H_{0}:\mu_{D}\geqslant0,\quad H_{1}:\mu_{D}<0. $$ 解 现在 n=8, $ \overline{x}_{d}=-0.0625 $, $ s_{d}=0.0765 $,而 $$ \frac{x_{d}}{s_{d}/\sqrt{8}}=-2.311<-t_{0.05}(7)=-1.8946 $$ 故拒绝 $ H_{0} $,认为 $ \mu_{D}<0 $,即认为人对红光的反应时间小于对绿光的反应时间,也就是人对红光的反应要比绿光快。 现在来讨论有关正态总体方差的假设检验问题.以下分单个总体和两个总体的情况来讨论. 设总体 $ X \sim N(\mu, \sigma^{2}) $, $ \mu, \sigma^{2} $ 均未知, $ X_{1}, X_{2}, \cdots, X_{n} $ 是来自 X 的样本。要求检 验假设(显著性水平为 $ \alpha $) $$ H_{0}:\sigma^{2}=\sigma_{0}^{2},\quad H_{1}:\sigma^{2}\neq\sigma_{0}^{2}, $$ $ \sigma_{0}^{2} $ 为已知常数. 由于 $ S^{2} $ 是 $ \sigma^{2} $ 的无偏估计,当 $ H_{0} $ 为真时,观察值 $ s^{2} $ 与 $ \sigma_{0}^{2} $ 的比值 $ \frac{s^{2}}{\sigma_{0}^{2}} $ 一般来说应在1附近摆动,而不应过分大于1或过分小于1。由第六章§3定理二知当 $ H_{0} $ 为真时 $$ \frac{(n-1)S^{2}}{\sigma_{0}^{2}}\sim\chi^{2}(n-1) $$ 我们取 $$ \chi^{2}=\frac{(n-1)S^{2}}{\sigma_{0}^{2}} $$ 作为检验统计量,如上所说知道上述检验问题的拒绝域具有以下的形式: $$ \frac{(n-1)s^{2}}{\sigma_{0}^{2}}\leqslant k_{1}\quad\text{ 或 }\quad\frac{(n-1)s^{2}}{\sigma_{0}^{2}}\geqslant k_{2}, $$ 此处 $ k_{1}, k_{2} $ 的值由下式确定: $$ \begin{aligned}&P\{ 当 H_{0} 为真拒绝 H_{0}\}\\ &=P_{\sigma_{0}^{2}}\{(\frac{(n-1)S^{2}}{\sigma_{0}^{2}}\leqslant k_{1})\cup(\frac{(n-1)S^{2}}{\sigma_{0}^{2}}\geqslant k_{2})\}=\alpha.\\ \end{aligned} $$ 为计算方便起见,习惯上取 $$ P_{\sigma_{0}^{2}}\{\frac{(n-1)S^{2}}{\sigma_{0}^{2}}\leqslant k_{1}\}=\frac{\alpha}{2},P_{\sigma_{0}^{2}}\{\frac{(n-1)S^{2}}{\sigma_{0}^{2}}\geqslant k_{2}\}=\frac{\alpha}{2}, $$ 故得 $ k_{1}=\chi_{1-\alpha/2}^{2}(n-1) $, $ k_{2}=\chi_{\alpha/2}^{2}(n-1) $。于是得拒绝域为 $$ \frac{(n-1)s^{2}}{\sigma_{0}^{2}}\leqslant\chi_{1-a/2}^{2}(n-1)\quad\text{ 或 }\quad\frac{(n-1)s^{2}}{\sigma_{0}^{2}}\geqslant\chi_{a/2}^{2}(n-1)\textcircled{1} $$ 下面来求单边检验问题(显著性水平为 $ \alpha $) $$ H_{0}:\sigma^{2}\leqslant\sigma_{0}^{2},\ H_{1}:\sigma^{2}>\sigma_{0}^{2} $$ 的拒绝域.因 $ H_{0} $中的全部 $ \sigma^{2} $都比 $ H_{1} $中的 $ \sigma^{2} $要小,当 $ H_{1} $为真时, $ S^{2} $的观察值 $ s^{2} $往往偏大,因此拒绝域的形式为 $$ s^{2}\geqslant k. $$ 下面来确定常数 k. $$ H_{0} $$ $$ \begin{aligned}H_{0}\}&=P_{\sigma^{2}\leqslant\sigma_{0}^{2}}\{S^{2}\geqslant k\}\\&=P_{\sigma^{2}\leqslant\sigma_{0}^{2}}\{\frac{(n-1)S^{2}}{\sigma_{0}^{2}}\geqslant\frac{(n-1)k}{\sigma_{0}^{2}}\}\end{aligned} $$ $$ \leqslant P_{\sigma^{2}}\leqslant\sigma_{0}^{2}\{\frac{(n-1)S^{2}}{\sigma^{2}}\geqslant\frac{(n-1)k}{\sigma_{0}^{2}}\}\text{。(因为}\sigma^{2}\leqslant\sigma_{0}^{2}) $$ 要控制 P{当 $ H_{0} $ 为真拒绝 $ H_{0} $} \leqslant \alpha $,只需令$ $$ P_{\sigma^{2}}\leq_{\sigma_{0}^{2}}\{\frac{(n-1)S^{2}}{\sigma^{2}}\geq\frac{(n-1)k}{\sigma_{0}^{2}}\}=\alpha. $$ 因 $ \frac{(n-1)S^{2}}{\sigma^{2}}\sim\chi^{2}(n-1) $,由(3.3)得 $ \frac{(n-1)k}{\sigma_{0}^{2}} $ $ =\chi_{a}^{2}(n-1) $(见图8-4)。于是 $ k=\frac{\sigma_{0}^{2}}{n-1}\chi_{a}^{2} $ (n-1),得检验问题(3.2)的拒绝域为 $ s^{2}\geq\sigma_{0}^{2} $ $ \frac{\sigma_{0}^{2}}{n-1}\chi_{\alpha}^{2}(n-1) $,即 $$ \chi^{2}=\frac{(n-1)s^{2}}{\sigma_{0}^{2}}\geqslant\chi_{a}^{2}(n-1). $$ 类似地,可得左边检验问题 $$ H_{0}:\sigma^{2}\geqslant\sigma_{0}^{2},\quad H_{1}:\sigma^{2}<\sigma_{0}^{2} $$ 的拒绝域为 $$ \chi^{2}=\frac{(n-1)s^{2}}{\sigma_{0}^{2}}\leqslant\chi_{1-a}^{2}(n-1). $$ 以上检验法称为 $ \chi^{2} $ 检验法. 例1 某厂生产的某种型号的电池,其寿命(以h计)长期以来服从方差 $ \sigma^{2}=5000 $的正态分布,现有一批这种电池,从它的生产情况来看,寿命的波动性有所改变.现随机取26只电池,测出其寿命的样本方差 $ s^{2}=9\ 200 $.问根据这一数据能否推断这批电池的寿命的波动性较以往的有显著的变化(取 $ \alpha=0.02 $)? 解 本题要求在水平 $ \alpha=0.02 $ 下检验假设 $$ H_{0}:\sigma^{2}=5\ 000,\quad H_{1}:\sigma^{2}\neq5\ 000. $$ 现在 n=26, $ \chi_{a/2}^{2}(n-1)=\chi_{0.01}^{2}(25)=44.314,\chi_{1-a/2}^{2}(25)=\chi_{0.99}^{2}(25)=11.524.\sigma_{0}^{2}=5000 $,由(3.1)式拒绝域为 $$ \frac{(n-1)s^{2}}{\sigma_{0}^{2}}\geqslant44.314\quad 或 \quad\frac{(n-1)s^{2}}{\sigma_{0}^{2}}\leqslant11.524. $$ 由观察值 $ s^{2}=9\ 200 $ 得 $ \frac{(n-1)s^{2}}{\sigma_{0}^{2}}=46>44.314 $,所以拒绝 $ H_{0} $,认为这批电池寿命的波动性较以往的有显著的变化. □ 设 $ X_{1}, X_{2}, \cdots, X_{n_{1}} $ 是来自总体 $ N(\mu_{1}, \sigma_{1}^{2}) $ 的样本, $ Y_{1}, Y_{2}, \cdots, Y_{n_{2}} $ 是来自总体 $ N(\mu_{2}, \sigma_{2}^{2}) $ 的样本,且两样本独立。其样本方差分别为 $ S_{1}^{2}, S_{2}^{2} $。且设 $ \mu_{1}, \mu_{2}, \sigma_{1}^{2}, \sigma_{2}^{2} $ 均为未知。现在需要检验假设(显著性水平为 $ \alpha $)。 $$ H_{0}:\sigma_{1}^{2}{\leqslant}\sigma_{2}^{2},\quad H_{1}:\sigma_{1}^{2}>\sigma_{2}^{2}. $$ 当 $ H_{0} $ 为真时, $ E(S_{1}^{2})=\sigma_{1}^{2}\leqslant\sigma_{2}^{2}=E(S_{2}^{2}) $,当 $ H_{1} $ 为真时 $ E(S_{1}^{2})=\sigma_{1}^{2}>\sigma_{2}^{2}=E(S_{2}^{2}) $。当 $ H_{1} $ 为真时,观察值 $ \frac{S_{1}^{2}}{S_{2}^{2}} $ 有偏大的趋势,故拒绝域具有形式 $$ \begin{array}{l}\frac{s_{1}^{2}}{s_{2}^{2}}\geqslant k,\end{array} $$ 常数 k 确定如下: $$ \begin{aligned}P\{ 当 H_{0} 为真拒绝 H_{0}\}=&P_{\sigma_{1}^{2}\leqslant\sigma_{2}^{2}}\{\frac{S_{1}^{2}}{S_{2}^{2}}\geqslant k\}\\ \leqslant&P_{\sigma_{1}^{2}\leqslant\sigma_{2}^{2}}\{\frac{S_{1}^{2}/S_{2}^{2}}{\sigma_{1}^{2}/\sigma_{2}^{2}}\geqslant k\}\text{(因为 }\sigma_{1}^{2}/\sigma_{2}^{2}\leqslant1)\text{}.\end{aligned} $$ 要控制 P{当 $ H_{0} $ 为真拒绝 $ H_{0} $} \leqslant \alpha $,只需令$ $$ P_{\sigma_{1}^{2}\leqslant\sigma_{2}^{2}}\{\frac{S_{1}^{2}/S_{2}^{2}}{\sigma_{1}^{2}/\sigma_{2}^{2}}\geqslant k\}=\alpha. $$ 由第六章§3定理四知 $ \frac{S_{1}^{2}/S_{2}^{2}}{\sigma_{1}^{2}/\sigma_{2}^{2}}\sim F(n_{1}-1,n_{2}-1) $,由(3.7)式得 $ k=F_{\alpha}(n_{1}-1,n_{2}-1) $。即得检验问题(3.6)的拒绝域为 $$ F=\frac{s_{1}^{2}}{s_{2}^{2}}\geqslant F_{a}(n_{1}-1,n_{2}-1). $$ 上述检验法称为 F 检验法. 关于 $ \sigma_{1}^{2} $, $ \sigma_{2}^{2} $ 的另外两个检验问题的拒绝域在表8.1中给出. 例2 设§2例2中的两个样本分别来自总体 $ N(\mu_{A},\sigma_{A}^{2}),N(\mu_{B},\sigma_{B}^{2}) $,且两样本独立。试检验 $ H_{0}:\sigma_{A}^{2}=\sigma_{B}^{2},H_{1}:\sigma_{A}^{2}\neq\sigma_{B}^{2} $,以说明我们假设 $ \sigma_{A}^{2}=\sigma_{B}^{2} $是合理的。(取显著性水平 $ \alpha=0.01 $) 解 此处 $ n_{1}=13, n_{2}=8, \alpha=0.01 $ ,拒绝域为 $$ \frac{s_{A}^{2}}{s_{B}^{2}}\geqslant F_{0.005}(12,7)=8.18, $$ 或 $$ \frac{s_{A}^{2}}{s_{B}^{2}}\leqslant F_{0.005}(12,7)=\frac{1}{F_{0.005}(7,12)}=\frac{1}{5.52}=0.18. $$ 现在 $ s_{A}^{2}=0.024^{2}, s_{B}^{2}=0.03^{2}, s_{A}^{2}/s_{B}^{2}=0.64 $ $$ 0.18<0.64<8.18, $$ 故接受 $ H_{0} $,认为两总体方差相等。两总体方差相等也称两总体具有方差齐性,这也表明§2例2中假设两总体方差相等是合理的。 置信区间与假设检验之间有明显的联系,先考察置信区间与双边检验之间的对应关系.设 $ X_{1},\cdots,X_{n} $是一个来自总体的样本, $ x_{1},\cdots,x_{n} $是相应的样本值, $ \Theta $是参数 $ \theta $的可能取值范围. 设 $ \left(\theta(X_{1},\cdots,X_{n}),\theta(X_{1},\cdots,X_{n})\right) $是参数 $ \theta $的一个置信水平为 $ 1-\alpha $的置信区间,则对于任意 $ \theta\in\Theta $,有 $$ P_{\theta}\{\theta\mid(X_{1},\cdots,X_{n})<\theta<\theta(X_{1},\cdots,X_{n})\}\geqslant1-\alpha $$ 考虑显著性水平为 $ \alpha $的双边检验 $$ H_{0}:\theta=\theta_{0},H_{1}:\theta\neq\theta_{0}. $$ 由(4.1)式 $$ P_{\theta_{0}}\{\underline{\theta}(X_{1},\cdots,X_{n})<\theta_{0}<\underline{\theta}(X_{1},\cdots,X_{n})\}\geqslant1-\alpha $$ 即有 $$ P_{\theta_{0}}\{(\theta_{0}\leqslant\theta(X_{1},\cdots,X_{n}))\cup(\theta_{0}\geqslant\theta(X_{1},\cdots,X_{n}))\}\leqslant\alpha. $$ 按显著性水平为 $ \alpha $ 的假设检验的拒绝域的定义,检验(4.2)的拒绝域为 $$ \theta_{0}\leqslant\underline{\theta}(x_{1},\cdots,x_{n})\quad\text{ 或 }\quad\theta_{0}\geqslant\overline{\theta}(x_{1},\cdots,x_{n}); $$ 接受域为 $$ \theta(x_{1},\cdots,x_{n})<\theta_{0}<\theta(x_{1},\cdots,x_{n}). $$ 这就是说,当我们要检验假设(4.2)时,先求出 $ \theta $的置信水平为 $ 1-\alpha $的置信区间 $ (\underline{\theta},\overline{\theta}) $,然后考察区间 $ (\underline{\theta},\overline{\theta}) $是否包含 $ \theta_{0} $,若 $ \theta_{0}\in(\underline{\theta},\overline{\theta}) $则接受 $ H_{0} $,若 $ \theta_{0}\notin(\underline{\theta},\overline{\theta}) $,则拒绝 $ H_{0} $。 反之,对于任意 $ \theta_{0}\in\Theta $,考虑显著性水平为 $ \alpha $的假设检验问题 $$ H_{0}:\theta=\theta_{0},H_{1}:\theta\neq\theta_{0}, $$ 假设它的接受域为 $$ \underline{\theta}(x_{1},\cdots,x_{n})<\underline{\theta}_{0}<\theta(x_{1},\cdots,x_{n}), $$ 即有 $$ P_{\theta_{0}}\{\theta(X_{1},\cdots,X_{n})<\theta_{0}<\theta(X_{1},\cdots,X_{n})\}\geqslant1-\alpha. $$ 由 $ \theta_{0} $ 的任意性,由上式知对于任意 $ \theta\in\Theta $,有 $$ P_{\theta}\{\theta(X_{1},\cdots,X_{n})<\theta<\overline{\theta}(X_{1},\cdots,X_{n})\}\geqslant1-\alpha. $$ 因此 $ \left(\theta\left(X_{1},\cdots,X_{n}\right),\theta\left(X_{1},\cdots,X_{n}\right)\right) $是参数 $ \theta $的一个置信水平为 $ 1-\alpha $的置信区间. 这就是说,为要求出参数 $ \theta $ 的置信水平为 $ 1-\alpha $ 的置信区间,我们先求出显著性水平为 $ \alpha $ 的假设检验问题: $ H_{0}:\theta=\theta_{0}, H_{1}:\theta\neq\theta_{0} $ 的接受域 $ \theta(x_{1},\cdots,x_{n})<\theta_{0}<\bar{\theta}(x_{1},\cdots,x_{n}) $,那么, $ (\theta(X_{1},\cdots,X_{n}),\bar{\theta}(X_{1},\cdots,X_{n})) $ 就是 $ \theta $ 的置信水平为 $ 1-\alpha $ 的置信区间。 还可验证,置信水平为 $ 1-\alpha $ 的单侧置信区间 $ (-\infty, \theta(X_1, \cdots, X_n)) $ 与显著性水平为 $ \alpha $ 的左边检验问题 $ H_0: \theta \geqslant \theta_0 $, $ H_1: \theta < \theta_0 $ 有类似的对应关系。即若已求得单侧置信区间 $ (-\infty, \bar{\theta}(X_1, \cdots, X_n)) $,则当 $ \theta_0 \in (-\infty, \bar{\theta}(x_1, \cdots, x_n)) $ 时接受 $ H_0 $,当 $ \theta_0 \notin (-\infty, \bar{\theta}(x_1, \cdots, x_n)) $ 时拒绝 $ H_0 $。反之,若已求得检验问题 $ H_0: \theta \geqslant \theta_0 $, $ H_1: \theta < \theta_0 $ 的接受域为 $ -\infty < \theta_0 \leqslant \bar{\theta}(x_1, \cdots, x_n) $,则可得 $ \theta $ 的一个单侧置信区间 $ (-\infty, \bar{\theta}(X_1, \cdots, X_n)) $。 置信水平为 $ 1-\alpha $ 的单侧置信区间 $ (\underline{\theta}(X_{1},\cdots,X_{n}),\infty) $ 与显著性水平为 $ \alpha $ 的右边检验问题 $ H_{0}:\theta\leq\theta_{0},H_{1}:\theta>\theta_{0} $ 也有类似的对应关系.即若已求得单侧置信区间 $ (\underline{\theta}(X_{1},\cdots,X_{n}),\infty) $ 则当 $ \theta_{0}\in(\underline{\theta}(x_{1},\cdots,x_{n}),\infty) $ 时接受 $ H_{0} $,当 $ \theta_{0}\notin(\underline{\theta}(x_{1},\cdots,x_{n}),\infty) $ 时拒绝 $ H_{0} $ 反之,若已求得检验问题 $ H_{0}:\theta\leq\theta_{0},H_{1}:\theta>\theta_{0} $ 的接受域为 $ \underline{\theta}(x_{1},\cdots,x_{n})\leq\theta_{0}<\infty $ ,则可得 $ \theta $ 的一个单侧置信区间 $ (\underline{\theta}(X_{1},\cdots,X_{n}),\infty) $ 例1 设 $ X \sim N(\mu, 1) $, $ \mu $ 未知, $ \alpha = 0.05 $,n = 16,且由一样本算得 $ \overline{x} = 5.20 $,于是得到参数 $ \mu $ 的一个置信水平为 0.95 的置信区间 $$ \begin{aligned}(\overline{x}-\frac{1}{\sqrt{16}}z_{0.025},\overline{x}+\frac{1}{\sqrt{16}}z_{0.025})&=(5.20-0.49,5.20+0.49)\\&=(4.71,5.69).\end{aligned} $$ 现在考虑检验问题 $ H_{0}:\mu=5.5,\quad H_{1}:\mu\neq5.5 $ 。由于 $ 5.5\in(4.71,5.69) $,故接受 $ H_{0} $ 。 例2 数据如上例.试求右边检验问题 $ H_{0}:\mu\leq\mu_{0},H_{1}:\mu>\mu_{0} $ 的接受域,并求 $ \mu $ 的单侧置信下限 $ (\alpha=0.05) $. 解 检验问题的拒绝域为 $ z=\frac{x-\mu_{0}}{1/\sqrt{16}}\geqslant z_{0.05} $,或即 $ \mu_{0}\leqslant4.79 $。于是检验问题的接受域为 $ \mu_{0}>4.79 $。这样就得到 $ \mu $ 的单侧置信区间 $ (4.79,\infty) $,单侧置信下限 $ \mu=4.79 $。 以上我们在进行假设检验时,总是根据问题的要求,预先给出显著性水平以控制犯第Ⅰ类错误的概率,而犯第Ⅱ类错误的概率则依赖于样本容量的选择.在 一些实际问题中,我们除了希望控制犯第Ⅰ类错误的概率外,往往还希望控制犯第Ⅱ类错误的概率。在这一节,我们将阐明如何选取样本的容量使得犯第Ⅱ类错误的概率控制在预先给定的限度之内。为此,我们引入施行特征函数。 定义 若 C 是参数 $ \theta $ 的某检验问题的一个检验法, $$ \beta(\theta)=P_{\theta}( 接受 H_{0}) $$ 称为检验法 C 的施行特征函数或 OC 函数,其图形称为 OC 曲线. 由定义知,若此检验法的显著性水平为 $ \alpha $,那么当真值 $ \theta \in H_{0} $ 时, $ \beta(\theta) $ 就是作出正确判断(即 $ H_{0} $ 为真时接受 $ H_{0} $)的概率,故此时 $ \beta(\theta) \geqslant 1 - \alpha $;而当 $ \theta \in H_{1} $ 时,则 $ \beta(\theta) $ 就是犯第Ⅱ类错误的概率,而 $ 1 - \beta(\theta) $ 是作出正确判断(即 $ H_{0} $ 为不真时拒绝 $ H_{0} $)的概率。函数 $ 1 - \beta(\theta) $ 称为检验法 C 的功效函数。当 $ \theta^{*} \in H_{1} $ 时,值 $ 1 - \beta(\theta^{*}) $ 称为检验法 C 在点 $ \theta^{*} $ 的功效。它表示当参数 $ \theta $ 的真值为 $ \theta^{*} $ 时,检验法 C 作出正确判断的概率。 本书只介绍正态总体均值的检验法的 OC 函数及其图形. 右边检验问题. $ H_{0}:\mu\leq\mu_{0},\quad H_{1}:\mu>\mu_{0} $ 的 OC 函数是 $$ \begin{aligned}\beta(\mu)=P_{\mu}( 接受 H_{0})&=P_{\mu}\{\frac{\overline{X}-\mu_{0}}{\sigma/\sqrt{n}} 其中 $ \lambda=\frac{\mu-\mu_{0}}{\sigma/\sqrt{n}} $ 。其图形如图 8-5 所示。此 OC 函数 $ \beta(\mu) $ 有如下性质: (1)它是 $ \lambda=\frac{\mu-\mu_{0}}{\sigma/\sqrt{n}} $ 的单调递减连续函数; (2) $ \lim_{\mu \to \mu_0^+} \beta(\mu) = 1 - \alpha $, $ \lim_{\mu \to \infty} \beta(\mu) = 0 $. 由 $ \beta(\mu) $ 的连续性可知,当参数的真值 $ \mu(\mu>\mu_{0}) $ 在 $ \mu_{0} $ 附近时,检验法的功效很低,即 $ \beta(\mu) $ 的值很大,亦即犯第Ⅱ类错误的概率很大。因为 $ \alpha $ 通常取得比较小,而不管 $ \sigma $ 多么小,n 多么大,只要 n 给定,总存在 $ \mu_{0} $ 附近的点 $ \mu(\mu>\mu_{0}) $ 使 $ \beta(\mu) $ 几乎等于 $ 1-\alpha $. 这表明,无论样本容量 n 多么大,要想对所有 $ \mu \in H_{1} $,即真值为 $ H_{1} $ 所规定的任一点时,控制犯第Ⅱ类错误的概率都很小是不可能的。但是我们可以使用 OC 函数 $ \beta(\mu) $ 以确定样本容量 n,使当真值 $ \mu \geqslant \mu_{0} + \delta $ ( $ \delta > 0 $ 为取定的值)时,犯第Ⅱ类错误的概率不超过给定的 $ \beta $。这是由于 $ \beta(\mu) $ 是 $ \mu $ 的递减函数,故当 $ \mu \geqslant \mu_{0} + \delta $ 时有 $$ \beta(\mu_{0}+\delta)\geqslant\beta(\mu). $$ 于是只要 $ \beta(\mu_{0}+\delta)=\Phi(z_{\alpha}-\sqrt{n}\delta/\sigma)\leqslant\beta $ ,亦即只要 n 满足 $$ z_{a}-\sqrt{n}\delta/\sigma\leqslant-z_{\beta} $$ 即可. 这就是说,只要 $$ \sqrt{n}\geqslant\frac{(z_{\alpha}+z_{\beta})\sigma}{\delta}, $$ 就能使当 $ \mu \in H_1 $ 且 $ \mu \geqslant \mu_0 + \delta $ 时(即真值 $ \mu \geqslant \mu_0 + \delta $ 时)犯第Ⅱ类错误的概率不超过 $ \beta $。 类似地,可得左边检验问题 $ H_{0}:\mu\geqslant\mu_{0},H_{1}:\mu<\mu_{0} $ 的 OC 函数为 $$ \beta(\mu)=\Phi(z_{a}+\lambda),\quad\lambda=\frac{\mu-\mu_{0}}{\sigma/\sqrt{n}}. $$ 当真值 $ \mu \geqslant \mu_{0} $ 时 $ \beta(\mu) $ 为作出正确判断的概率;当真值 $ \mu < \mu_{0} $ 时, $ \beta(\mu) $ 给出犯第Ⅱ类错误的概率。只要样本容量 n 满足 $$ \sqrt{n}\geqslant\frac{(z_{\alpha}+z_{\beta})\sigma}{\delta}, $$ 就能使当 $ \mu\in H_{1} $ 且 $ \mu\leqslant\mu_{0}-\delta(\delta>0 $ ,为取定的值)时,犯第Ⅱ类错误的概率不超过给定的值 $ \beta $ 双边检验问题 $ H_{0}:\mu=\mu_{0}, H_{1}:\mu\neq\mu_{0} $ 的 OC 函数是 $$ \begin{aligned}\beta(\mu)&=P_{\mu}(接受 H_{0})=P_{\mu}\{-z_{a/2}<\frac{\overline{X}-\mu_{0}}{\sigma/\sqrt{n}} OC 曲线如图 8-6 所示. 注意 $ \beta(\mu) $ 是 $ |\lambda| $ 的严格单调下降函数. 在双边检验问题中,若要求对 $ H_{1} $ 中满足 $ \left|\mu - \mu_{0}\right| \geqslant \delta > 0 $ 的 $ \mu $ 处的函数值 $ \beta(\mu) \leqslant \beta $,则需解超越方程 $$ \begin{aligned}\beta=&\Phi(z_{\alpha/2}-\sqrt{n}\delta/\sigma)+\\&\Phi(z_{\alpha/2}+\sqrt{n}\delta/\sigma)-1\end{aligned} $$ 才能确定 n. 通常因 n 较大,故总可以认为 $ z_{\alpha/2}+\sqrt{n}\delta/\sigma\geqslant4 $,于是 $ \Phi(z_{\alpha/2}+\sqrt{n}\delta/\sigma)\approx1 $,故近似地有 $$ \beta\approx\Phi(z_{a/2}-\sqrt{n}\delta/\sigma). $$ 由此知只要样本容量 n 满足 $$ z_{\alpha/2}-\sqrt{n}\delta/\sigma\leq-z_{\beta}, $$ 即只要 n 满足 $$ \sqrt{n}\geqslant(z_{a/2}+z_{\beta})\frac{\sigma}{\delta}, $$ 就能使当 $ \mu\in H_{1} $ 且 $ |\mu-\mu_{0}|\geqslant\delta $ ( $ \delta>0 $, 为取定的值) 时, 犯第Ⅱ类错误的概率不超过给定的值 $ \beta $. 例1(工业产品质量抽验方案)设有一大批产品,产品质量指标 $ X \sim N(\mu, \sigma^{2}) $。以 $ \mu $ 小者为佳,厂方要求所确定的验收方案对高质量的产品 $ (\mu \leqslant \mu_{0}) $ 能以高概率 $ 1 - \alpha $ 为买方所接受。买方则要求低质产品 $ (\mu \geqslant \mu_{0} + \delta, \delta > 0) $ 能以高概率 $ 1 - \beta $ 被拒绝。 $ \alpha, \beta $ 由厂方与买方协商给出。并采取一次抽样以确定该批产品是否为买方所接受。问应怎样安排抽样方案。已知 $ \mu_{0} = 120, \delta = 20 $,且由工厂长期经验知 $ \sigma^{2} = 900 $。又经商定 $ \alpha, \beta $ 均取为 0.05。 解 检验问题可表达为 $$ H_{0}:\mu\leqslant\mu_{0},\quad H_{1}:\mu>\mu_{0}, $$ 且要求当 $ \mu \geqslant \mu_{1} + \delta $ 时能以 $ 1 - \beta = 0.95 $ 的概率拒绝 $ H_{0} $. 由 Z 检验,拒绝域为 $$ \frac{\overline{x}-\mu_{0}}{\sigma/\sqrt{n}}\geq z_{a}, $$ 故 OC 函数为 $$ \begin{align*}\beta(\mu)=&P_{\mu}\{\frac{\overline{X}-\mu_{0}}{\sigma/\sqrt{n}} 现要求当 $ \mu \geqslant \mu_{0} + \delta $ 时 $ \beta(\mu) \leqslant \beta $. 因 $ \beta(\mu) $ 是 $ \mu $ 的递减函数,故只需 $ \beta(\mu_{0} + \delta) = \beta $ 即可. 此时,由 (5.9) 式可得 $$ \sqrt{n}\geqslant\frac{(z_{a}+z_{b})\sigma}{\delta}. $$ 按给定的数据算得 $ n \geqslant 24.35 $,故取 n=25。且当 $ \bar{x} $ 满足 $ \frac{x-\mu_0}{\sigma/\sqrt{n}} \geqslant z_\alpha = z_{0.05} = 1.645 $ 时,即当 $ \bar{x} \geqslant 129.87 $ 时,买方就拒绝这批产品,而当 $ \bar{x} < 129.87 $ 时,买方接受这批产品。 右边检验问题 $ H_{0}:\mu\leq\mu_{0},H_{1}:\mu>\mu_{0} $ 的 t 检验法的 OC 函数是 $$ \beta(\mu)=P_{\mu}(H_{0})=P_{\mu}\{\frac{X-\mu_{0}}{S/\sqrt{n}} 其中变量 $$ \frac{\overline{X}-\mu_{0}}{S/\sqrt{n}}=\left(\frac{\overline{X}-\mu}{\sigma/\sqrt{n}}+\lambda\right)/\left(\frac{S}{\sigma}\right),\quad\lambda=\frac{\mu-\mu_{0}}{\sigma/\sqrt{n}}. $$ 我们称变量 $ \frac{\overline{X}-\mu_{0}}{S/\sqrt{n}} $服从非中心参数为 $ \lambda $、自由度为n-1的非中心t分布.在 $ \lambda=0 $时,它是通常的 $ t(n-1) $变量. 若给定 $ \alpha, \beta $ 以及 $ \delta > 0 $,则可从书末附表7查得所需容量 n,使得当 $ \mu \in H_{1} $ 且 $ \frac{\mu - \mu_{0}}{\sigma} \geqslant \delta $ 时犯第Ⅱ类错误的概率不超过 $ \beta $. 若给定 $ \alpha, \beta $ 以及 $ \delta > 0 $,对于左边检验问题 $ H_{0}: \mu \geq \mu_{0}, H_{1}: \mu < \mu_{0} $ 的 t 检验法,也可从附表 7 查得所需容量 n,使得当 $ \mu \in H_{1} $ 且 $ \frac{\mu - \mu_{0}}{\sigma} \leq -\delta $ 时犯第Ⅱ类错误的概率不超过 $ \beta $。对于双边检验问题 $ H_{0}: \mu = \mu_{0}, H_{1}: \mu \neq \mu_{0} $ 的 t 检验法也可从附表 7 查得所需容量 n,使得当 $ \mu \in H_{1} $,且 $ \frac{|\mu - \mu_{0}|}{\sigma} \geq \delta $ 时犯第Ⅱ类错误的概率不超过 $ \beta $。 例2 考虑在显著性水平 $ \alpha=0.05 $ 下进行t检验 $$ H_{0}:\mu{\leqslant}68,\quad H_{1}:\mu>68. $$ (1)要求在 $ H_{1} $中 $ \mu\geqslant\mu_{1}=68+\sigma $时犯第Ⅱ类错误的概率不超过 $ \beta=0.05 $。求所需的样本容量. (2)若样本容量为 n=30,问在 $ H_{1} $ 中 $ \mu=\mu_{1}=68+0.75\sigma $ 时犯第Ⅱ类错误的概率是多少? 解 (1)此处 $ \alpha=\beta=0.05,\mu_{0}=68,\delta=\frac{\mu_{1}-\mu_{0}}{\sigma}=\frac{(68+\sigma)-68}{\sigma}=1 $ ,查附表7得 n=13. (2)现在 $ \alpha=0.05, n=30, \delta=\frac{\mu_1-\mu_0}{\sigma}=\frac{(68+0.75\sigma)-68}{\sigma}=0.75 $,查附表 7,得 $ \beta=0.01 $。 例3 考虑在显著性水平 $ \alpha=0.05 $ 下进行t检验 $$ H_{0}:\mu=14,\quad H_{1}:\mu\neq14. $$ 要求在 $ H_{1} $ 中 $ \frac{|\mu-14|}{\sigma} \geqslant 0.4 $ 时犯第Ⅱ类错误的概率不超过 $ \beta = 0.1 $,求所需样本容量. 解 此处 $ \alpha=0.05,\beta=0.1,\delta=0.4 $ ,查附表7得 n=68. 在实际问题中,有时只给出 $ \alpha,\beta $ 及 $ \left|\mu_{1}-\mu_{0}\right| $ 的值,而需要确定所需的样本容量 n. 这时由于 $ \sigma $ 未知,不能确定 $ \delta=\left|\mu_{1}-\mu_{0}\right|/\sigma $ 的值,因而不能直接查表以确定 样本容量. 此时可采用下述近似方法. 先适当取一值 $ n_{1} $,抽取容量为 $ n_{1} $ 的样本,根据这一样本计算 $ s^{2} $ 的值,以 $ s^{2} $ 作为 $ \sigma^{2} $ 的估计,算出 $ \delta $ 的近似值. 由 $ \alpha, \beta, \delta $ 的值查附表7定出样本的容量,记为 $ n_{2} $. 若 $ n_{1} \geqslant n_{2} $,则取 $ n_{1} $ 作为所求的容量,即取 $ n = n_{1} $. 否则,再抽 $ n_{2} - n_{1} $ 个独立观察值与原来抽得的观察值合并,重新计算 $ \delta $ 的近似值. 然后用 $ \delta $ 的新近似值和 $ \alpha, \beta $ 查附表7,再次定出样本容量. 记为 $ n_{3} $. 若 $ n_{2} \geqslant n_{3} $,则取 $ n = n_{2} $,否则再按上法重复进行. 一般,只需试少数几次就可得到所求的样本容量 n. 下面考虑两个正态总体均值差的 t 检验. 若两个正态总体 $ N(\mu_1, \sigma_1^2) $, $ N(\mu_2, \sigma_2^2) $ 中 $ \sigma_1^2 = \sigma_2^2 = \sigma^2 $ 而 $ \sigma^2 $ 未知。在均值差 $ \mu_1 - \mu_2 $ 的检验问题 $ H_0: \mu_1 - \mu_2 = 0 $, $ H_1: \mu_1 - \mu_2 \neq 0 $ (或 $ H_0: \mu_1 - \mu_2 \leq 0 $, $ H_1: \mu_1 - \mu_2 > 0 $ 或 $ H_0: \mu_1 - \mu_2 \geq 0 $, $ H_1: \mu_1 - \mu_2 < 0 $)的 t 检验法中,当分别自两个总体取得的相互独立的样本其容量 $ n_1 = n_2 = n $ 时,给定 $ \alpha, \beta $ 以及 $ \delta = |\mu_1 - \mu_2| / \sigma $ 的值后可以查附表 8 得到所需样本容量,使当 $ |\mu_1 - \mu_2| / \sigma \geq \delta $ 时犯第Ⅱ类错误的概率小于或等于 $ \beta $。当仅给出 $ \alpha, \beta $ 以及 $ |\mu_1 - \mu_2| $ 的值时,可按类似于上面所说的方法处理。 燃料的辛烷值越高,燃料质量越好。因燃料B较燃料A价格便宜,因此,若两者辛烷值相同时,则使用燃料B;但若含量的均值差 $ \mu_{A}-\mu_{B}\geqslant5 $,则使用燃料A。设两总体的分布均可认为是正态的,而两个样本相互独立。问应采用哪种燃料(取 $ \alpha=0.01,\beta=0.01 $)? 解 按题意需要在显著性水平 $ \alpha=0.01 $ 下检验假设 $$ H_{0}:\mu_{A}-\mu_{B}{\leqslant}0,\quad H_{1}:\mu_{A}-\mu_{B}>0, $$ 并要求在 $ \mu_{A}-\mu_{B}\geqslant5 $ 时,犯第Ⅱ类错误的概率不超过 $ \beta=0.01 $ 所取的样本容量为 $ n_{A}=n_{B}=12 $,且有 $ \overline{x}_{A}=80.83 $, $ \overline{x}_{B}=78.67 $, $ s_{A}^{2}=5.61 $, $ s_{B}^{2}=6.06 $。经显著性水平为 0.1 的 F 检验知,可认为两总体的方差相等,即有 $ \sigma_{A}^{2}=\sigma_{B}^{2} $,记为 $ \sigma^{2} $。因 $ n_{1}=n_{2} $,取 $ \hat{\sigma}^{2}=(s_{A}^{2}+s_{B}^{2})/2=5.835 $ 作为 $ \sigma^{2} $ 的点估计,取 $ \hat{\sigma}=\sqrt{\hat{\sigma}^{2}} $,于是 $ \delta=5/\hat{\sigma}=2.07 $,查表,当 $ \alpha=0.01 $, $ \beta=0.01 $, $ \delta=2.07 $ 时 $ n\geqslant8 $。现 n=12,故已近似地满足要求。而右边检验的拒绝域为 $$ t=\frac{\overline{x}_{A}-\overline{x}_{B}}{s_{w}\sqrt{1/n_{A}+1/n_{B}}}\geqslant t_{0.01}(n_{1}+n_{2}-2)=2.5083. $$ 由样本观察值算得 t=2.19<2.5083,故接受 $ H_{0} $,即采用燃料 B. 上面介绍的各种检验法都是在总体分布形式为已知的前提下进行讨论的. 但在实际问题中,有时不能知道总体服从什么类型的分布,这时就需要根据样本来检验关于分布的假设.本节介绍 $ \chi^{2} $ 拟合检验法.它可以用来检验总体是否具有某一个指定的分布或属于某一个分布族,还介绍专用于检验分布是否为正态的“偏度、峰度检验法”. 设总体 X 的分布未知, $ x_{1}, x_{2}, \cdots, x_{n} $ 是来自 X 的样本值。我们来检验假设 $$ H_{0}: 总体 X 的分布函数为 F(x). $$ $$ H_{1}: 总体 X 的分布函数不是 F(x)\textcircled{1} $$ 其中设 $ F(x) $ 不含未知参数.(也常以分布律或概率密度代替 $ F(x) $.) 下面来定义检验统计量。将在 $ H_{0} $ 下 X 可能取值的全体 $ \Omega $ 分成互不相交的子集 $ A_{1}, A_{2}, \cdots, A_{k} $,以 $ f_{i} $ (i=1,2,\cdots,k) 记样本观察值 $ x_{1}, x_{2}, \cdots, x_{n} $ 中落在 $ A_{i} $ 的个数,这表示事件 $ A_{i}=\{X $ 的值落在子集 $ A_{i} $ 内\} 在 n 次独立试验中发生 $ f_{i} $ 次,于是在这 n 次试验中事件 $ A_{i} $ 发生的频率为 $ f_{i}/n $。另一方面,当 $ H_{0} $ 为真时,我们可以根据 $ H_{0} $ 中所假设的 X 的分布函数来计算事件 $ A_{i} $ 的概率,得到 $ p_{i}=P(A_{i}), i=1,2,\cdots,k $。频率 $ f_{i}/n $ 与概率 $ p_{i} $ 会有差异,但一般来说,当 $ H_{0} $ 为真,且试验的次数又甚多时,这种差异不应太大,因此 $ \left(\frac{f_{i}}{n}-p_{i}\right)^{2} $ 不应太大。我们采用形如 $$ \sum_{i=1}^{k}C_{i}\left(\frac{f_{i}}{n}-p_{i}\right)^{2} $$ 的统计量来度量样本与 $ H_{0} $ 中所假设的分布的吻合程度,其中 $ C_{i}\quad(i=1,2,\cdots,k) $ ②为给定的常数.皮尔逊证明,如果选取 $ C_{i}=n/p_{i}\quad(i=1,2,\cdots,k) $,则由(6.2)定义的统计量具有下述定理中所述的简单性质.于是我们就采用 $$ \chi^{2}=\sum_{i=1}^{k}\frac{n}{p_{i}}\left(\frac{f_{i}}{n}-p_{i}\right)^{2}=\sum_{i=1}^{n}\frac{f_{i}^{2}}{n p_{i}}-n $$ 作为检验统计量. 定理 若 n 充分大 $ (n \geqslant 50) $,则当 $ H_{0} $ 为真时统计量 (6.3) 近似服从 $ \chi^{2}(k-1) $ 分布。(证略) 据以上的讨论,当 $ H_{0} $ 为真时,(6.3) 式中的 $ \chi^{2} $ 不应太大,如 $ \chi^{2} $ 过分大就拒绝 $ H_{0} $,拒绝域的形式为 $$ \chi^{2}\geqslant G\quad(G 为正常数 ). $$ 对于给定的显著性水平 $ \alpha $,确定G使 P{当 $ H_{0} $ 为真时拒绝 $ H_{0} $}= $ P_{H_{0}} $ $ \{\chi^{2} \geqslant G\}=\alpha $ 由上述定理得 $ G = \chi_{a}^{2}(k-1) $. 即当样本观察值使 (6.3) 式中的 $ \chi^{2} $ 的值有 $$ \chi^{2}\geqslant\chi_{a}^{2}(k-1) $$ 则在显著性水平 $ \alpha $ 下拒绝 $ H_{0} $; 否则就接受 $ H_{0} $. 这就是单个分布的 $ \chi^{2} $ 拟合检验法. $ \chi^{2} $ 拟合检验法是基于上述定理得到的,所以使用时必须注意 n 不能小于 50。另外 $ np_{i} $ 不能太小,应有 $ np_{i} \geqslant 5 $,否则应适当合并 $ A_{i} $,以满足这个要求(见下例)。 例1. 下表列出了某一地区在夏季的一个月中由100个气象站报告的雷暴雨的次数. 其中 $ f_{i} $ 是报告雷暴雨次数为 i 的气象站数。试用 $ \chi^{2} $ 拟合检验法检验雷暴雨的次数 X 是否服从均值 $ \lambda=1 $ 的泊松分布(取显著性水平 $ \alpha=0.05 $)。 解 按题意需检验假设 $$ H_{0}:\quad P\{X=i\}=\frac{\lambda^{i}\mathrm{e}^{-\lambda}}{i!}=\frac{\mathrm{e}^{-1}}{i!},\quad i=0,1,\cdots. $$ 在 $ H_{0} $ 下 X 所有可能取的值为 $ \Omega=\{0,1,2,\cdots\} $,将 $ \Omega $ 分成如表所示的两两不相交的子集 $ A_{0}, A_{1}, \cdots, A_{6} $,则有 $ P\{X=i\} $ 为 $$ p_{i}=P\{X=i\}=\frac{\mathrm{e}^{-1}}{i!},i=0,1,\cdots,5. $$ 例如 $$ p_{0}=P\{X=0\}=e^{-1}=0.36788, $$ $$ p_{3}=P\{X=3\}=\frac{\mathrm{e}^{-1}}{3!}=0.06131, $$ $$ p_{6}=P\{X\geqslant6\}=1-\sum_{i=0}^{5}p_{i}=0.059. $$ $$ n=100. $$ $$ \Sigma=127.04 $$ 计算结果如表8-2所示,其中有些 $ np_{i}<5 $的组予以适当合并,使得每组均有 $ np_{i}\geqslant5 $,如表中第4列花括号所示。并组后 $ k=4,\chi^{2} $的自由度为k-1=4-1=3. $ \chi_{0.05}^{2}(k-1)=\chi_{0.05}^{2}(3)=7.815 $。现在 $ \chi^{2}=127.04-100=27.04>7.815 $,故在显著性水平0.05下拒绝 $ H_{0} $,认为样本不是来自均值 $ \lambda=1 $的泊松分布。 例2 在研究牛的毛色与牛角的有无,这样两对性状分离现象时,用黑色无角牛与红色有角牛杂交,子二代出现黑色无角牛192头,黑色有角牛78头,红色无角牛72头,红色有角牛18头,共360头,问这两对性状是否符合孟德尔遗传规律中9:3:3:1的遗传比例? 以 X 记各种牛的序号,按题意需检验各类牛的头数符合比例 9:3:3:1,即 $ (9/16):(3/16):(3/16):(1/16) $. 需检验假设: $ H_{0}:X $ 的分布律为 取显著性水平为0.1. 所需计算列在表8-3中(n=360). $$ \Sigma=363.37 $$ 现在 $ \chi^{2}=363.37-360=3.37, k=4, \chi_{0.1}^{2}(4-1)=6.251>3.37 $ ,故接受 $ H_{0} $,认为两性状符合孟德尔遗传规律中 9:3:3:1 的遗传比例. □ 在(一)中要检验的原假设是 $ H_{0} $:总体X的分布函数是 $ F(x) $,其中 $ F(x) $是已知的,这种情况是不多的.我们经常遇到的所需检验的原假设是 $$ H_{0}: 总体 X 的分布函数是 F\left(x;\theta_{1},\cdots,\theta_{r}\right), $$ 其中 F 的形式已知,而 $ \boldsymbol{\theta} = (\theta_{1}, \theta_{2}, \cdots, \theta_{r}) $ 是未知参数,它们在某一个范围取值。在 $ F(x; \theta_{1}, \theta_{2}, \cdots, \theta_{r}) $ 中当参数 $ \theta_{1}, \theta_{2}, \cdots, \theta_{r} $ 取不同的值时,就得到不同的分布,因而 $ F(x; \theta_{1}, \theta_{2}, \cdots, \theta_{r}) $ 代表一族分布。(6.4) 中的 $ H_{0} $ 表示总体 X 的分布属于分 布族 $ F(x;\theta_{1},\theta_{2},\cdots,\theta_{r}) $. 采用类似(一)中的方法来定义检验统计量,将在 $ H_{0} $ 下 X 可能取值的全体 $ \Omega $ 分成 k (k>r+1) 个互不相交的子集 $ A_{1}, A_{2}, \cdots, A_{k} $,以 $ f_{i} $ (i=1,2,\cdots,k) 记样本观察值 $ x_{1}, x_{2}, \cdots, x_{n} $ 落在 $ A_{i} $ 的个数,则事件 $ A_{i}=\{X $ 的值落在 $ A_{i} $ 内\} 的频率为 $ f_{i}/n $. 另一方面,当 $ H_{0} $ 为真时,由 $ H_{0} $ 所假设的分布函数来计算 $ P(A_{i}) $,得到 $ P(A_{i})=p_{i}(\theta_{1},\theta_{2},\cdots,\theta_{r})=p_{i}(\pmb{\theta})=p_{i} $. 此时,需先利用样本求出未知参数的最大似然估计(在 $ H_{0} $ 下),以估计值作为参数值,求出 $ p_{i} $ 的估计值 $ \hat{p}_{i}=\hat{P}(A_{i}) $,在 (6.3) 式中以 $ \hat{p}_{i} $ 代替 $ p_{i} $,取 $$ \chi^{2}=\sum_{i=1}^{k}\frac{f_{i}^{2}}{n\hat{p}_{i}}\cdots n $$ 作为检验假设 $ H_{0} $ 的统计量. 可以证明,在某些条件下,在 $ H_{0} $ 为真时近似地有 $$ \chi^{2}=\sum_{i=1}^{k}\frac{f_{i}^{2}}{n\hat{p}_{i}}-n\sim\chi^{2}(k-r-1) $$ 与在(一)中一样可得假设检验问题(6.4)的拒绝域为 $$ \chi^{2}\geqslant\chi_{a}^{2}(k-r-1) $$ $ \alpha $为显著性水平.以上就是用来检验分布族的 $ \chi^{2} $拟合检验法. 例 3 在一实验中,每隔一定时间观察一次由某种铀所放射的到达计数器上的 $ \alpha $ 粒子数 X,共观察了 100 次,得结果如下表: 其中 $ f_{i} $ 是观察到有 i 个 $ \alpha $ 粒子的次数. 从理论上考虑知 X 应服从泊松分布 $$ P\{X=i\}=\frac{\lambda^{i}\mathrm{e}^{-\lambda}}{i!},\quad i=0,1,2,\cdots. $$ 问(6.7)式是否符合实际(取 $ \alpha=0.05 $)?即在显著性水平0.05下检验假设 $ H_{0} $: 总体 X 服从泊松分布 $$ P\{X=i\}=\frac{\lambda^{i}\mathrm{e}^{-\lambda}}{i!},\quad i=0,1,2,\cdots. $$ 解 因在 $ H_{0} $ 中参数 $ \lambda $ 未具体给出,所以先估计 $ \lambda $. 由最大似然估计法得 $ \hat{\lambda} = \bar{x} = 4.2 $. 在 $ H_{0} $ 假设下,即在 X 服从泊松分布的假设下,X 所有可能取的值为 $ \Omega = \{0,1,2,\cdots\} $,将 $ \Omega $ 分成如表 8-4 所示的两两不相交的子集 $ A_{0}, A_{1}, \cdots, A_{12} $. 则 $ P\{X=i\} $ 有估计 $$ \hat{p}_{i}=\hat{P}\{X=i\}=\frac{4.2^{i}\mathrm{e}^{-4.2}}{i!},\quad i=0,1,\cdots. $$ 例如 $$ \hat{p}_{0}=\hat{P}\{X=0\}=\mathrm{e}^{-4.2}=0.015, $$ $$ \hat{p}_{3}=\hat{P}\{X=3\}=\frac{4.2^{3}\mathrm{e}^{-4.2}}{3!}=0.185, $$ $$ \hat{p}_{12}=\hat{P}\{X\geqslant12\}=1-\sum_{i=0}^{11}\hat{p}_{i}=0.002. $$ $$ \Sigma=106,281 $$ 计算结果如表8-5所示,其中有些 $ n\hat{p}_{i}<5 $的组予以适当合并,使得每组均有 $ n\hat{p}_{i}\geqslant5 $,如表中第四列花括号所示。此处,并组后k=8,但因在计算概率时,估计了一个参数 $ \lambda $,故 $ r=1,\chi^{2} $的自由度为8-1-1=6. $ \chi_{0.05}^{2}(k-r-1)=\chi_{0.05}^{2}(6)=12.592 $,现在 $ \chi^{2}=106.281-100=6.281<12.592 $,故在显著性水平0.05下接受 $ H_{0} $.即认为样本来自泊松分布总体.也就是说认为理论上的结论是符合实际的.□ 注意 本题答案是“接受 $ H_{0} $,认为总体 X 的分布属于泊松分布族,即认为 $ X \sim \pi(\lambda) $”,亦即“认为必有某一个参数 $ \lambda_{0}, X \sim \pi(\lambda_{0}) $”,而不能将答案误写成“X 服从以 $ \lambda = 4.2 $ 为参数的泊松分布”. 例4 自1965年1月1日至1971年2月9日共2231天中,全世界记录到里氏震级4级和4级以上地震计162次,统计如下: 试检验相继两次地震间隔的天数 X 服从指数分布 ( $ \alpha=0.05 $). 解 按题意需检验假设 $ H_{0} $: X 的概率密度为 $$ f(x)=\{\begin{aligned}&\frac{1}{\theta}\mathrm{e}^{-x/\theta},&x&>0,\\ &0,&x&\leqslant0.\end{aligned}. $$ 在这里, $ H_{0} $ 中的参数 $ \theta $ 未给出,先由最大似然估计法求得 $ \theta $ 的估计为 $ \theta = \overline{x} = \frac{223}{162} = 13.77 $。在 $ H_{0} $ 下,X 可能取值的全体 $ \Omega $ 为区间 $ [0, \infty) $。将区间 $ [0, \infty) $ 分为 k = 9 个互不重叠的小区间: $ A_{1} = [0, 4.5] $, $ A_{2} = (4.5, 9.5] $, $ \cdots $, $ A_{9} = (39.5, \infty) $,如表 8-6 第二列所示。若 $ H_{0} $ 为真,X 的分布函数的估计为 $$ \hat{F}(x)=\{\begin{aligned}&1-\mathrm{e}^{-x/13.77},&x&>0,\\ &0,&x&\leqslant0.\end{aligned}. $$ 由上式可得概率 $ p_{i}=P(A_{i}) $ 的估计: $$ \hat{p}_{i}=\hat{P}(A_{i})=\hat{P}\{a_{i} 例如 $$ \hat{p}_{2}=\hat{P}(A_{2})=\hat{P}\{4,5 而 $$ \hat{p}_{9}=\hat{P}(A_{9})=1-\sum_{i=1}^{8}\hat{F}(A_{i})=0.056\ 8. $$ 将计算结果列表如下. 现在 $ \chi^{2}=163.5633-162=1.5633 $ ,因为 $ \chi_{0.05}^{2}(k-r-1)=\chi_{0.05}^{2}(8-1-1)=\chi_{0.05}^{2}(6)=12.592>1.5633 $ ,故在显著性水平 0.05 下接受 $ H_{0} $,认为 X 服从指 $$ \Sigma=163,5633 $$ 数分布. 例 5 对于第六章 §2 例 1 中的数据,试验验它们是否来自正态总体 X(取显著性水平 $ \alpha=0.1 $). 解 需检验假设 $ H_{0} $: X 的概率密度为 $$ f(x)=\frac{1}{\sqrt{2\pi}\sigma}\mathrm{e}^{-\frac{(x-\mu)^2}{2\sigma^2}},\quad-\infty 因在 $ H_{0} $ 中未给出 $ \mu $, $ \sigma^{2} $ 的数值。需先估计 $ \mu $, $ \sigma^{2} $。由最大似然估计法得 $ \mu $, $ \sigma^{2} $ 的估计值为 $ \hat{\mu}=143.8 $, $ \hat{\mu}^{2}=(6.0)^{2} $。我们将在 $ H_{0} $ 下 X 可能取值的区间 $ (-\infty, \infty) $ 分为 7 个小区间,并取事件 $ A_{i} $ 如表 8-7 中第一列所示。若 $ H_{0} $ 为真 X 的概率密度的估计为 $$ \hat{f}(x)=\frac{1}{\sqrt{2\pi}\times6}e^{-\frac{(x-143.8)^{2}}{2\times6^{2}}},-\infty 按上式并查标准正态分布的分布函数表即可得概率 $ P(A_{i}) $ 的估计.例如 $$ \begin{aligned}\hat{p}_{2}&=\hat{P}(A_{2})=\hat{P}\{129.5 将计算结果列表如下: $$ \Sigma=87.67 $$ 现在 $ \chi^{2}=87.67-84=3.67 $,因为 $ \chi_{0.1}^{2}(k-r-1)=\chi_{0.1}^{2}(5-2-1)=\chi_{0.1}^{2}(2)=4.605>3.67 $,故在水平 0.1 下接受 $ H_{0} $,即认为数据来自正态分布总体。 根据第五章关于中心极限定理的论述知道,正态分布随机变量是较广泛地 存在的,因此,当研究一连续型总体时,人们往往先考察它是否服从正态分布。上面介绍的 $ \chi^{2} $ 拟合检验法虽然是检验总体分布的较一般的方法,但用它来检验总体的正态性时,犯第Ⅱ类错误的概率往往较大。为此,统计学家们对检验正态总体的种种方法进行了比较,根据奥野忠一等人在20世纪70年代进行的大量模拟计算的结果,认为正态性检验方法中,总的来说,以“偏度、峰度检验法”及“夏皮罗-威尔克法”较为有效。在这里我们仅介绍偏度、峰度检验法。 随机变量 X 的偏度和峰度指的是 X 的标准化变量 $ [X - E(X)] / \sqrt{D(X)} $ 的三阶矩和四阶矩: $$ \nu_{1}=E\left[\left(\frac{X-E(X)}{\sqrt{D(X)}}\right)^{3}\right]=\frac{E[(X-E(X))^{3}]}{(D(X))^{3/2}}, $$ $$ \nu_{2}=E\left[\left(\frac{X-E(X)}{\sqrt{D(X)}}\right)^{4}\right]=\frac{E\left[(X-E(X))^{4}\right]}{(D(X))^{2}}. $$ 当随机变量 X 服从正态分布时, $ \nu_{1}=0 $ 且 $ \nu_{2}=3 $. 设 $ X_{1}, X_{2}, \cdots, X_{n} $ 是来自总体 X 的样本,则 $ \nu_{1}, \nu_{2} $ 的矩估计量分别是 $$ G_{1}=B_{3}/B_{2}^{3/2},\quad G_{2}=B_{4}/B_{2}^{2}, $$ 其中 $ B_{k}\quad(k=2,3,4) $ 是样本 k 阶中心矩,并分别称 $ G_{1}, G_{2} $ 为样本偏度和样本峰度. 若总体 X 为正态变量,则可证当 n 充分大时,近似地有 $$ G_{1}\sim N\left(0,\frac{6(n-2)}{(n+1)(n+3)}\right), $$ $$ G_{2}\sim N\left(3-\frac{6}{n+1},\frac{24n(n-2)(n-3)}{(n+1)^{2}(n+3)(n+5)}\right). $$ 设 $ X_{1}, X_{2}, \cdots, X_{n} $ 是来自总体 X 的样本,现在来检验假设 $ H_{0} $: X 为正态总体. 记 $$ \sigma_{1}=\sqrt{\frac{6(n-2)}{(n+1)(n+3)}},\quad\sigma_{2}=\sqrt{\frac{24n(n-2)(n-3)}{(n+1)^{2}(n+3)(n+5)}}, $$ $ \mu_{2}=3-\frac{6}{n+1},U_{1}=G_{1}/\sigma_{1},U_{2}=(G_{2}-\mu_{2})/\sigma_{2} $. 当 $ H_{0} $ 为真且 n 充分大时,近似地有 $$ U_{1}\sim N(0,1),\quad U_{2}\sim N(0,1). $$ 由第六章§3知样本偏度 $ G_{1} $、样本峰度 $ G_{2} $分别依概率收敛于总体偏度 $ \nu_{1} $和总体峰度 $ \nu_{2} $.因此当 $ H_{0} $为真且n充分大时,一般来说, $ G_{1} $与 $ \nu_{1}=0 $的偏离不应太大,而 $ G_{2} $与 $ \nu_{2}=3 $的偏离不应太大.故从直观来看当 $ \left|U_{1}\right| $的观察值 $ \left|u_{1}\right| $或 $ \left|U_{2}\right| $的观察值 $ \left|u_{2}\right| $过大时就拒绝 $ H_{0} $.取显著性水平为 $ \alpha $, $ H_{0} $的拒绝域为 $$ u_{1}\mid\geqslant k_{1}\quad 或 \quad|u_{2}\mid\geqslant k_{2}, $$ 其中 $ k_{1}, k_{2} $ 由以下两式确定: $$ P_{H_{0}}\{|U_{1}|\geqslant k_{1}\}=\frac{\alpha}{2},\quad P_{H_{0}}\{|U_{2}|\geqslant k_{2}\}=\frac{\alpha}{2}. $$ 这里记号 $ P_{H_{0}}\{\cdot\} $ 表示当 $ H_{0} $ 为真时事件 $ \{\cdot\} $ 的概率,即有 $ k_{1}=z_{\alpha/4}, k_{2}=z_{\alpha/4} $. 于是得拒绝域为 $$ u_{1}\mid\geqslant z_{a/4}\quad 或 \quad|u_{2}\mid\geqslant z_{a/4}. $$ 下面来验证当 n 充分大时上述检验法近似地满足显著性水平为 $ \alpha $ 的要求. 事实上当 n 充分大时有 $$ \begin{aligned}P\{& 当 H_{0} 为真拒绝 H_{0}\}\\&=P_{H_{0}}\{(|U_{1}|\geqslant z_{\alpha/4})\bigcup(|U_{2}|\geqslant z_{\alpha/4})\}\\&\leqslant P_{H_{0}}\{|U_{1}|\geqslant z_{\alpha/4}\}+P_{H_{0}}\{|U_{2}|\geqslant z_{\alpha/4}\}=\frac{\alpha}{2}+\frac{\alpha}{2}=\alpha.\end{aligned} $$ 例 6 试用偏度、峰度检验法检验例 5 中的数据是否来自正态总体(取 $ \alpha = 0.1 $). $ H_{0} $: 数据来自正态总体. 这里 $ \alpha=0.1, n=84, \sigma_{1}=\sqrt{\frac{6(n-2)}{(n+1)(n+3)}}=0.2579 $ $$ \sigma_{2}=\sqrt{\frac{24n(n-2)(n-3)}{(n+1)^{2}(n+3)(n+5)}}=0.4892 $$ $ \mu_{2}=3-\frac{6}{n+1}=2.9294 $. 下面来计算样本中心矩 $ B_{2}, B_{3}, B_{4} $,计算时可利用以下关系式: $$ B_{2}=A_{2}-A_{1}^{2},B_{3}=A_{3}-3A_{2}A_{1}+2A_{1}^{3}, $$ $$ B_{4}=A_{4}-4A_{3}A_{1}+6A_{2}A_{1}^{2}-3A_{1}^{4}, $$ 其中 $ A_{k}=\frac{1}{n}\sum_{i=1}^{n}X_{i}^{k} $ (k=1,2,3,4)为 k 阶样本矩。经计算得 $ A_{k} $ (k=1,2,3,4), $ B_{k} $ (k=2,3,4)的观察值分别为 $$ A_{1}=143.7738,\quad A_{2}=20706.13,\quad A_{3}=2987099, $$ $$ A_{4}=4.316426\times10^{8},\quad B_{2}=35.2246,\quad B_{3}=-28.5,\quad B_{4}=3840. $$ 样本偏度和样本峰度的观察值分别为 $$ g_{1}=-0.1363,\quad g_{2}=3.0948. $$ 而 $ z_{a/4}=z_{0.025}=1.96 $ 。由(6.11)式,拒绝域为 $$ |u_{1}|=|g_{1}/\sigma_{1}|\geqslant1.96\quad\text{ 或 }\quad|u_{2}|=|g_{2}-\mu_{2}|/\sigma_{2}\geqslant1.96. $$ 现算得 $ \left|u_{1}\right|=0.5285<1.96,\left|u_{2}\right|=0.3381<1.96 $ ,故接受 $ H_{0} $,认为数据来自 正态分布的总体. 上述检验法称为偏度、峰度检验法。使用这一检验法时样本容量以大于100为宜。 本节介绍一种有效的、且使用方便的检验方法——秩和检验法. 设有两个连续型总体,它们的概率密度函数分别为 $ f_{1}(x) $, $ f_{2}(x) $,均为未知,但已知 $$ f_{1}(x)=f_{2}(x-a),\quad a 为未知常数 , $$ 即 $ f_{1} $ 与 $ f_{2} $ 至多只差一平移。我们要检验下述各项假设 $$ H_{0}:a=0,H_{1}:a<0. $$ $$ H_{0}:a=0,H_{1}:a>0. $$ $$ H_{0}:a=0,H_{1}:a\neq0. $$ 特别,若两总体的均值存在,分别记作 $ \mu_{1}, \mu_{2} $,则由于 $ f_{1}, f_{2} $ 至多只差一平移,故有 $$ \mu_{2}=\mu_{1}-a. $$ 此时,上述各项假设分别等价于 $$ H_{0}:\mu_{1}=\mu_{2},H_{1}:\mu_{1}<\mu_{2}. $$ $$ H_{0}:\mu_{1}=\mu_{2},H_{1}:\mu_{1}>\mu_{2}. $$ $$ H_{0}:\mu_{1}=\mu_{2},H_{1}:\mu_{1}\neq\mu_{2}. $$ 现在来介绍威尔柯克斯(Frank Wilcoxon)提出的秩和检验法以检验上述假设.为此,先引入秩的概念. 秩设X为一总体,将一容量为n的样本观察值按自小到大的次序编号排列成 $$ x_{(1)} 称 $ x_{(i)} $ 的足标 i 为 $ x_{(i)} $ 的秩, i=1,2,\cdots,n. 现设自1,2两总体分别抽取容量为 $ n_{1},n_{2} $的样本,且设两样本独立.这里总假定 $ n_{1}\leqslant n_{2} $.我们将这 $ n_{1}+n_{2} $个观察值放在一起,按自小到大的次序排列,求出每个观察值的秩,然后将属于第1个总体的样本观察值的秩相加,其和记为 $ R_{1} $,称为第1样本的秩和.其余观察值的秩的总和记作 $ R_{2} $,称为第2样本的秩和.显然 $ R_{1},R_{2} $是离散型的随机变量,且有 $$ R_{1}+R_{2}=\frac{1}{2}(n_{1}+n_{2})(n_{1}+n_{2}+1). $$ 所以, $ R_{1}, R_{2} $ 中的一个确定后另一个随之而定.这样,我们只要考虑统计量 $ R_{1} $ 即可. 现在来解决双边检验问题(7.4). 对此,先作直观分析. 当 $ H_{0} $ 为真时,即有 $ f_{1}(x)=f_{2}(x) $,这时两个独立样本实际上是来自同一个总体. 因而第1个样本中诸元素的秩应该随机地、分散地在自然数 $ 1 \sim n_{1} + n_{2} $ 中取值,一般来说不应过分集中取较小的或较大的值. 考虑到 $$ \frac{1}{2}n_{1}(n_{1}+1)\leqslant R_{1}\leqslant\frac{1}{2}n_{1}(n_{1}+2n_{2}+1) $$ 即知当 $ H_{0} $为真时秩和 $ R_{1} $一般来说不应取太靠近上述不等式两端的值.因而,当 $ R_{1} $的观察值 $ r_{1} $过分大或过分小时,我们都拒绝 $ H_{0} $. 据以上分析,对于双边检验(7.4),在给定显著性水平 $ \alpha $下, $ H_{0} $的拒绝域为 $$ r_{1}{\leqslant}C_{U}\left(\frac{\alpha}{2}\right)\quad 或 \quad r_{1}{\geqslant}C_{L}\left(\frac{\alpha}{2}\right), $$ 其中临界点 $ C_{U}\left(\frac{\alpha}{2}\right) $ 是满足 $ P_{a=0}\{R_{1}\leqslant C_{U}(\frac{\alpha}{2})\}\leqslant\frac{\alpha}{2} $ 的最大整数,而 $ C_{L}\left(\frac{\alpha}{2}\right) $ 是满足 $ P_{a=0}\{R_{1}\geqslant C_{L}(\frac{\alpha}{2})\}\leqslant\frac{\alpha}{2} $ 的最小整数。而犯第Ⅰ类错误的概率为 $$ P_{a=0}\{R_{1}\leqslant C_{U}(\frac{\alpha}{2})\}+P_{a=0}\{R_{1}\geqslant C_{L}(\frac{\alpha}{2})\}\leqslant\frac{\alpha}{2}+\frac{\alpha}{2}=\alpha. $$ 如果知道 $ R_{1} $ 的分布,则临界点 $ C_{U}\left(\frac{\alpha}{2}\right) $, $ C_{L}\left(\frac{\alpha}{2}\right) $ 是不难求得的。下面以 $ n_{1}=3, n_{2}=4 $ 的情况为例来说明求临界点的方法。 当 $ n_{1}=3,n_{2}=4 $时,第1个样本中各观察值的秩的不同取法共有 $ \binom{3+4}{3}=35 $种.现将这35种情况列表如下. 由于这35种情况的出现是等可能的,由上表容易求得 $ R_{1} $的分布律和分布函数如下: 于是,对于不同的 $ \alpha $值,容易写出检验问题(7.4)的临界点和拒绝域.例如,给定 $ \alpha=0.2 $.由上表知 $$ P_{a=0}\{R_{1}\leqslant7\}=2/35<0,1=\frac{\alpha}{2}, $$ $$ P_{a=0}\{R_{1}\geqslant17\}=2/35<0.1=\frac{\alpha}{2}. $$ 即有 $ C_{U}(0,1)=7, C_{L}(0,1)=17 $. 故当 $ n_{1}=3, n_{2}=4 $,在显著性水平0.2下检验问题(7.4)的拒绝域为 $$ r_{1}\leqslant7\quad 或 \quad r_{1}\geqslant17. $$ 此时,犯第Ⅰ类错误的概率为 $$ P_{a=0}\{R_{1}\leqslant7\}+P_{a=0}\{R_{1}\geqslant17\}=2/35+2/35=0.114. $$ 类似地可得左边检验(7.2)的拒绝域为(显著性水平为 $ \alpha $) $$ r_{1}\leqslant C_{U}(\alpha), $$ 此处,临界点 $ C_{U}(\alpha) $ 是满足 $ P_{a=0}\{R_{1}\leqslant C_{U}(\alpha)\}\leqslant\alpha $ 的最大整数。而右边检验问题的拒绝域为(显著性水平为 $ \alpha $) $$ r_{1}\geqslant C_{L}\left(\alpha\right), $$ 此处,临界点 $ C_{L}(\alpha) $ 是满足 $ P_{\alpha=0}\{R_{1}\geqslant C_{L}(\alpha)\}\leqslant\alpha $ 的最小整数. 例如,若给定 $ \alpha=0.1 $,抽取的样本容量为 $ n_{1}=3, n_{2}=4 $,则由上表知检验问题(7.3)的拒绝域为 $$ r_{1}\geqslant17. $$ 此时犯第Ⅰ类错误的概率为2/35<0.1. 书末附表9中列出了 $ n_{1} $和 $ n_{2} $自2到10为止的 $ n_{1},n_{2} $的各种组合的临界点,以及相应的犯第Ⅰ类错误的概率. 例1 为查明某种血清是否会抑制白血病,选取患白血病已到晚期的老鼠9只,其中有5只接受这种治疗,另4只则不做这种治疗。设两样本相互独立。从试验开始时计算,其存活时间(以月计)如下: 设治疗与否的存活时间的概率密度至多只差一个平移. 取 $ \alpha=0.05 $,问这种血清对白血病是否有抑制作用? 解 本题需检验接受治疗的老鼠的存活期是否有增长。分别以 $ \mu_{1}, \mu_{2} $ 表示不作治疗和接受治疗的老鼠的存活时间总体的均值。需要检验的假设是 $$ H_{0}:\mu_{1}=\mu_{2},\quad H_{1}:\mu_{1}<\mu_{2}. $$ 这里, $ n_{1}=4,n_{2}=5,\alpha=0.05 $。先计算对应于 $ n_{1}=4 $ 的一组观察值的秩和。将两组数据放在一起按自小到大的次序排列。对来自第1个总体 $ (n_{1}=4) $ 的数据下面加一表示之。即有 所以 $ R_{1} $ 的观察值为 $ r_{1}=1+2+4+5=12 $ 。查附表9知 $ C_{U}(0.05)=12 $ ,即拒绝域为 $ r_{1}\leq12 $ 。而现在 $ r_{1}=12 $ ,故拒绝 $ H_{0} $,即认为这种血清对白血病有抑制作用。 可以证明,当 $ H_{0} $为真时(即a=0时) $$ \begin{align*}\mu_{R_{1}}&=E(R_{1})=\frac{n_{1}(n_{1}+n_{2}+1)}{2},\\\sigma_{R_{1}}^{2}&=D(R_{1})=\frac{n_{1}n_{2}(n_{1}+n_{2}+1)}{12}.\end{align*} $$ 而当 $ n_{1}, n_{2} \geqslant 10 $ ,当 $ H_{0} $ 为真时,近似地有 $$ R_{1}\sim N(\mu_{R_{1}},\sigma_{R_{1}}^{2}). $$ 因此,当 $ n_{1},n_{2}\geqslant10 $时我们可以采用 $$ Z=\frac{R_{1}-\mu_{R_{1}}}{\sigma_{R_{1}}} $$ 作为检验统计量. 在显著性水平 $ \alpha $ 下双边检验、右边检验、左边检验的近似拒绝域分别为 $$ |z|\geqslant z_{a/2},\quad z\geqslant z_{a},\quad z\leqslant-z_{a}. $$ 这里 z 是 Z 的观察值. 例2 某商店为了确定向公司A或公司B购买某种商品,将A,B公司以往各次进货的次品率进行比较,数据如下,设两样本独立。问两公司的商品的质量有无显著差异。设两公司的商品的次品率的密度至多只差一个平移,取显著性水平 $ \alpha=0.05 $。 解 分别以 $ \mu_{A}, \mu_{B} $ 记公司 A, B 的商品次品率总体的均值. 所需检验的假设是 $$ H_{0}:\mu_{A}=\mu_{B},\quad H_{1}:\mu_{A}\neq\mu_{B}. $$ 先将数据按自小到大的次序排列,得到对应于 $ n_{1}=10 $的样本的秩和为 $$ r_{1}=1+3+5+8+12+14+15+17+20+21=116. $$ 又当 $ H_{0} $为真时, $$ \begin{align*}E(R_{1})&=\frac{1}{2}n_{1}(n_{1}+n_{2}+1)=\frac{1}{2}\times10(10+13+1)=120,\\&D(R_{1})=\frac{1}{12}n_{1}n_{2}(n_{1}+n_{2}+1)=260.\end{align*} $$ 故知当 $ H_{0} $为真时近似地有 $$ R_{1}\sim N(120,260). $$ 拒绝域为 $$ \frac{\left|r_{1}-120\right|}{\sqrt{260}}\geqslant z_{0.025}=1.96. $$ 现在 $ R_{1} $ 的观察值为 $ r_{1}=116 $ ,得 $ \left|r_{1}-120\right|/\sqrt{260}=0.25<1.96 $ ,故接受 $ H_{0} $ ,认为两个公司商品的质量无显著差异. □ 在实际问题中(7.5)式中会出现某些观察值相等的情况,对于这种观察值的秩定义为足标的平均值。例如,若抽得的样本按次序排成0,1,1,1,2,3,3,则三个1的秩均为 $ (2+3+4)/3=3 $,两个3的秩均为 $ (6+7)/2=6.5 $。 将两个样本 $ n_{1}+n_{2}=n $ 个元素按自小到大的次序排列,若出现 k 个秩相同的组,设其中有 $ t_{i} $ 个数的秩为 $ a_{i}, i=1,2,\cdots,k, a_{1}<\cdots $$ \sigma_{R_{1}}^{2}=\frac{n_{1}n_{2}\left[n(n^{2}-1)-\sum_{i=1}^{k}t_{i}(t_{i}^{2}-1)\right]}{12n(n-1)}. $$ 当 k 不大时,附表 9 仍能使用,但表载值为近似值。又当 $ n_{1}, n_{2} \geqslant 10, H_{0} $ 为真,且 k 不大时,近似地有 $$ \begin{array}{r}{R_{1}\sim N(\mu_{R_{1}},\sigma_{R_{1}}^{2}),}\end{array} $$ 其中 $ \mu_{R_{1}} = n_{1}(n_{1} + n_{2} + 1)/2 $,而 $ \sigma_{R_{1}}^{2} $ 由(7.9)式确定。这时我们就采用 $$ Z=\frac{R_{1}-\mu_{R_{1}}}{\sigma_{R_{1}}} $$ 作为检验统计量来检验假设检验问题(7.2)′~(7.4)′. 例 3 两位化验员各自读得某种液体黏度如下: 设数据可以认为分别来自仅均值可能有差异的两个总体的样本. 试在 $ \alpha=0.05 $ 下,检验假设 $$ H_{0}:\mu_{1}=\mu_{2},\quad H_{1}:\mu_{1}>\mu_{2}, $$ 其中 $ \mu_{1}, \mu_{2} $ 分别为两总体的均值. 解 将两个样本的元素混合,按自小到大次序排列。并求出各个元素的秩如下: 现在 $ n_{1}=10, n_{2}=11, n=21, \mu_{R_{1}}=10 \times 22/2 = 110, k=2, \sum_{i=1}^{c} t_{i} (t_{i}^{2} - 1) = 3 \times (9 - 1) + 2 \times (4 - 1) = 30 $,按(7.9)式得 $ \sigma_{R_{1}}^{2} = 201 $。当 $ H_{0} $ 为真时近似地有 $$ R_{1}\sim N(110,201). $$ 拒绝域为 $$ \frac{r_{1}-110}{\sqrt{201}}\geqslant z_{0.05}=1.645. $$ 现在 $ R_{1} $ 的观察值为 $ r_{1}=121 $ ,得 $ (r_{1}-110)/\sqrt{201}=0.776<1.645 $ 。故接受 $ H_{0} $ 认为两位化验员所测得的数据无显著差异。 以上讨论的假设检验方法称为临界值法。本节介绍另一种被称为 p 值检验法的检验方法。先从一个例题讲起。 例1 设总体 $ X \sim N(\mu, \sigma^{2}) $, $ \mu $ 未知, $ \sigma^{2} = 100 $,现有样本 $ x_{1}, x_{2}, \cdots, x_{52} $,算得 $ \overline{x} = 62.75 $。现在来检验假设 $$ H_{0}:\mu=\mu_{0}=60,\qquad H_{1}:\mu>60. $$ 采用 Z 检验法,检验统计量为 $$ Z=\frac{\overline{X}-\mu_{0}}{\sigma/\sqrt{n}}. $$ 以数据代入,得 Z 的观察值为 $$ z_{0}=\frac{62.75-60}{10/\sqrt{52}}=1.983. $$ 概率 $$ P\{Z\geqslant z_{0}\}=P\{Z\geqslant1,983\}=1-\Phi(1,983)=0,023,8. $$ 此即为图8-7中标准正态曲线下位于 $ z_{0} $右边的尾部面积. 此概率称为 Z 检验法的右边检验的 p 值. 记为 $$ P\{Z\geqslant z_{0}\}=p 值 (=0,023\ 7). $$ 若显著性水平 $ \alpha \geqslant p = 0.0238 $,则对应的临界值 $ z_{\alpha} \leqslant 1.983 $,这表示观察值 $ z_{0} = 1.983 $ 落在拒绝域内(图 8-7(1)),因而拒绝 $ H_{0} $;又若显著性水平 $ \alpha < p = 0.0238 $,则对应的临界值 $ z_{\alpha} > 1.983 $,这表示观察值 $ z_{0} = 1.983 $ 不落在拒绝域内(图 8-7(2)),因而接受 $ H_{0} $。 据此,p值 $ =P\{Z\geqslant z_{0}\}=0.023 $ 8是原假设 $ H_{0} $ 可被拒绝的最小显著性水平. 一般,p值的定义是: 定义 假设检验问题的 p 值(probability value)是由检验统计量的样本观察值得出的原假设可被拒绝的最小显著性水平. 任一检验问题的 p 值可以根据检验统计量的样本观察值以及检验统计量在 $ H_{0} $ 下一个特定的参数值(一般是 $ H_{0} $ 与 $ H_{1} $ 所规定的参数的分界点)对应的分布求出。例如在正态总体 $ N(\mu, \sigma^{2}) $ 均值的检验中,当 $ \sigma $ 未知时,可采用检验统计量 $ t = \frac{\overline{X} - \mu_{0}}{S/\sqrt{n}} $,在以下三个检验问题中,当 $ \mu = \mu_{0} $ 时 $ t \sim t(n-1) $。如果由样本求得统计量 t 的观察值为 $ t_{0} $,那么在检验问题 $ H_{0}:\mu \leqslant \mu_{0}, H_{1}:\mu > \mu_{0} $ 中, p 值 $ =P_{\mu_{0}}\{t\geqslant t_{0}\}=t_{0} $ 右侧尾部面积,如图(8-8(1)); $$ H_{0}:\mu\geqslant\mu_{0},H_{1}:\mu<\mu_{0} 中 , $$ p 值 $ =P_{\mu_{0}}\{t\leqslant t_{0}\}=t_{0} $ 左侧尾部面积,如图(8-8(2)) $ H_{0}:\mu=\mu_{0}, H_{1}:\mu\neq\mu_{0} $ 中, (i) 当 $ t_{0} \geq 0 $ 时 $$ \begin{aligned}p\ 值 &=P_{\mu_{0}}\{|t|\geqslant t_{0}\}=P_{\mu_{0}}\{(t\leqslant-t_{0})\bigcup(t\geqslant t_{0})\}\\&=2\times(t_{0}\ 右侧尾部面积 )( 如图 8-9(1)).\end{aligned} $$ (ii)当 $ t_{0}<0 $时 $$ \begin{aligned}p\ 值 &=P_{\mu_{0}}\{\mid t\mid\geqslant-t_{0}\}=P_{\mu_{0}}\{(t\leqslant t_{0})\cup(t\geqslant-t_{0})\}\\&=2\times(t_{0}\ 左侧尾部面积 )( 如图 8-9(2)).\end{aligned} $$ 综合(i)(ii),p值=2×(由 $ t_{0} $界定的尾部面积). 上述各图中的曲线均为 $ t(n-1) $ 分布的概率密度曲线. 在现代计算机统计软件中,一般都给出检验问题的p值. 按 p 值的定义,对于任意指定的显著性水平 $ \alpha $,就有 (1)若 p 值 $ \leqslant\alpha $,则在显著性水平 $ \alpha $ 下拒绝 $ H_{0} $ (2) 若 p 值 $ >\alpha $,则在显著性水平 $ \alpha $ 下接受 $ H_{0} $ 有了这两条结论就能方便地确定 $ H_{0} $ 的拒绝域. 这种利用 p 值来确定检验拒绝域的方法, 称为 p 值检验法. 用临界值法来确定 $ H_{0} $ 的拒绝域时,例如当取 $ \alpha=0.05 $ 时知道要拒绝 $ H_{0} $,再取 $ \alpha=0.01 $ 也要拒绝 $ H_{0} $,但不能知道将 $ \alpha $ 再降低一些是否也要拒绝 $ H_{0} $。而 p 值法给出了拒绝 $ H_{0} $ 的最小显著性水平。因此 p 值法比临界值法给出了有关拒绝域的拒绝值法。 绝域的更多的信息. 例2 用p值法检验本章§1例2的检验问题 $$ H_{0}:\mu\leqslant\mu_{0}=-0.545,\quad H_{1}:\mu>\mu_{0}\quad\alpha=0.05. $$ 解 用 Z 检验法,现在检验统计量 $ Z=\frac{\overline{X}-\mu_{0}}{\sigma/\sqrt{n}} $ 的观察值为 $$ z_{0}=\frac{-0.535-(-0.545)}{0.008/\sqrt{5}}=2.7955. $$ $$ p\ 值 =P\{Z\geqslant2,795\ 5\}=1-\Phi(2,795\ 5)=0,002\ 6. $$ p值 $ \leq\alpha=0.05 $,故拒绝 $ H_{0} $ 例 3 用 p 检验法检验本章 §2 例 1 的检验问题 $$ H_{0}:\mu\leqslant\mu_{0}=225,\quad H_{1}:\mu>225,\quad\alpha=0.05. $$ 解 用 t 检验法,现在检验统计量 $ t=\frac{\overline{X}-\mu_{0}}{S/\sqrt{n}} $ 的观察值为 $$ t_{0}=\frac{241.5-225}{98.725\ 9/\sqrt{16}}=0.668\ 5 $$ 由计算机算得 $$ p\ 值 =P\{t\geqslant0,668\ 5\}=0,257\ 0, $$ p值 $ >\alpha=0.05 $,故接受 $ H_{0} $ 例4 用p值检验法检验本章§3例1中检验问题 $$ H_{0}:\sigma^{2}=5\ 000,H_{1}:\sigma^{2}\neq5\ 000,\alpha=0.02. $$ 解 用 $ \chi^{2} $ 检验法. 现在检验统计量 $ \chi^{2}=\frac{(n-1)S^{2}}{\sigma_{0}^{2}} $ 的观察值为 $$ \chi^{2}=\frac{25\times9\ 200}{5\ 000}=46 $$ 由计算机得 $$ p\ 值 =2\times P\{\chi^{2}\geqslant46\}=0.012\ 8, $$ p值 $ <\alpha=0.02 $,故拒绝 $ H_{0} $ p 值表示反对原假设 $ H_{0} $ 的依据的强度,p 值越小,反对 $ H_{0} $ 的依据越强、越充分(譬如对于某个检验问题的检验统计量的观察值的 p 值 = 0.0009,p 值如此的小,以至于几乎不可能在 $ H_{0} $ 为真时出现目前的观察值,这说明拒绝 $ H_{0} $ 的理由很强,我们就拒绝 $ H_{0} $)。 一般,若 p 值 $ \leq 0.01 $,称推断拒绝 $ H_{0} $ 的依据很强或称检验是高度显著的;若 0.01 < p 值 $ \leq 0.05 $ 称推断拒绝 $ H_{0} $ 的依据是强的或称检验是显著的;若 0.05 < p 值 $ \leq 0.1 $ 称推断拒绝 $ H_{0} $ 的理由是弱的,检验是不显著的;若 p 值 > 0.1 一般来说没有理由拒绝 $ H_{0} $。基于 p 值,研究者可以使用任意希望的显著性水平来 作计算。在杂志上或在一些技术报告中,许多研究者在讲述假设检验的结果时,常不明显地论及显著性水平以及临界值,代之以简单地引用假设检验的p值,利用或让读者利用它来评价反对原假设的依据的强度,作出推断。 统计推断就是由样本来推断总体,它包括两个基本问题:统计估计和假设检验。上一章讲述了参数估计,本章讨论假设检验问题。有关总体分布的未知参数或未知分布形式的种种论断叫统计假设,人们要根据样本所提供的信息对所考虑的假设作出接受或拒绝的决策。假设检验就是作出这一决策的过程。 一般,人们总是对原假设 $ H_{0} $ 作出接受或拒绝的决策。由于作出判断原假设 $ H_{0} $ 是否为真的依据是一个样本,由于样本的随机性,当 $ H_{0} $ 为真时,检验统计量的观察值也会落入拒绝域,致使我们作出拒绝 $ H_{0} $ 的错误决策;而当 $ H_{0} $ 为不真时,检验统计量的观察值也会未落入拒绝域,致使我们作出接受 $ H_{0} $ 的错误决策。 我们使用“接受假设”或“拒绝假设”这样的术语。接受一个假设并不意味着确信它是真的,它只意味着决定采取某种行动(例如 A);拒绝一个假设也不意味着它是假的,这也仅仅是作出采取另一种不同的行动(例如 B)。不论哪种情况,都存在作出错误选择的可能性。 当样本容量 n 固定时,减小犯第 I 类错误的概率,就会增大犯第Ⅱ类错误的概率,反之亦然。我们的做法是控制犯第 I 类错误的概率,使 $$ P\{ 当 H_{0} 为真拒绝 H_{0}\}\leqslant\alpha, $$ 其中 $ 0 < \alpha < 1 $ 是给定的小的数。 $ \alpha $ 称为检验的显著性水平。这种只对犯第Ⅰ类错误的概率加以控制而不考虑犯第Ⅱ类错误的概率的检验称为显著性检验。 在进行显著性检验时,犯第Ⅰ类错误的概率是由我们控制的。 $ \alpha $取得小,则概率P{当 $ H_{0} $为真拒绝 $ H_{0} $}就小,这保证了当 $ H_{0} $为真时错误地拒绝 $ H_{0} $的可能性很小。这意味着 $ H_{0} $是受到保护的,也表明 $ H_{0} $、 $ H_{1} $的地位不是对等的。于是,在一对对立假设中,选哪一个作为 $ H_{0} $需要小心。例如,考虑某种药品是否为真,这里可能犯两种错误:(1)将假药误作为真药,则冒着伤害病人的健康甚至生命的风险;(2)将真药误作为假药,则冒着造成经济损失的风险。显然,犯错误(1)比犯错误(2)的后果严重,因此,我们选取“ $ H_{0} $:药品为假, $ H_{1} $:药品为真”,即是使得犯第Ⅰ类错误“当药品为假时错判药品为真”的概率 $ \leqslant\alpha $。就是说,选择 $ H_{0} $、 $ H_{1} $使得两类错误中后果严重的错误成为第Ⅰ类错误。这是选择 $ H_{0} $、 $ H_{1} $的一个原则。 如果在两类错误中,没有一类错误的后果严重更需要避免时,常常取 $ H_{0} $为维持现状,即 取 $ H_{0} $ 为“无效益”、“无改进”、“无价值”等等.例如,取 $ H_{0} $: 新技术未提高效益, $ H_{1} $: 新技术提高效益. 实际上,我们感兴趣的是 $ H_{1} $ “提高效益”,但对采用新技术应持慎重态度。选取 $ H_{0} $ 为“新技术未提高效益”,一旦 $ H_{0} $ 被拒绝了,表示有较强的理由去采用新技术。 在实际问题中,情况比较复杂,如何选取 $ H_{0} $、 $ H_{1} $ 只能在实践中积累经验,根据实际情况去判断了. 注意,拒绝域的形式是由 $ H_{1} $确定的. 我们还介绍了置信区间与假设检验的关系.知道了置信区间就能容易判明是否接受原假设;反之,知道了检验的接受域就得到了相应的置信区间. 原假设 备择假设 检验统计量 单边检验 双边检验 显著性水平 拒绝域 显著性检验 一个正态总体的参数的检验 两个正态总体均值差、方差比的检验 成对数据的检验 $ \chi^{2} $ 分布拟合检验 偏度、峰度检验 秩和检验 $$ \begin{array}{l} 3,25 \quad 3,27 \quad 3,24 \quad 3,26 \quad 3,24 \end{array} $$ 设测定值总体服从正态分布,但参数均未知,问在 $ \alpha=0.01 $ 下能否接受假设:这批矿砂的镍含量的均值为 3.25. $$ \begin{aligned}&0.693\quad0.749\quad0.654\quad0.670\quad0.662\quad0.672\quad0.615\quad0.606\quad0.690\quad0.628\\&0.668\quad0.611\quad0.606\quad0.609\quad0.601\quad0.553\quad0.570\quad0.844\quad0.576\quad0.933\end{aligned} $$ 设这一工厂生产的矩形的宽度与长度的比值总体服从正态分布,其均值为 $ \mu $,方差为 $ \sigma^{2} $, $ \mu $, $ \sigma^{2} $ 均未知。试检验假设(取 $ \alpha=0.05 $) $$ H_{0}:\mu=0.618,\quad H_{1}:\mu\neq0.618. $$ 设装配时间的总体服从正态分布 $ N(\mu,\sigma^{2}), \mu, \sigma^{2} $ 均未知. 是否可以认为装配时间的均值显著大于10(取 $ \alpha=0.05 $)? 16 25 21 20 23 21 19 15 13 23 17 20 29 18 22 16 22 设维生素含量服从正态分布 $ N(\mu, \sigma^{2}) $, $ \mu, \sigma^{2} $ 均未知,问这批罐头是否符合要求(取显著性水平 $ \alpha = 0.05 $). 设两组数据分别来自正态总体,且两总体方差相等,但参数均未知。两样本相互独立。问两位作家所写的小品文中包含由3个字母组成的单字的比例是否有显著的差异(取 $ \alpha=0.05 $)? 设两样本分别来自正态总体,且两总体的方差相等,但参数均未知。两样本独立。分别以 $ \mu_{1} $, $ \mu_{2} $ 记对应于老、新过程的总体的均值,试检验假设 ( $ \alpha=0.05 $) $$ H_{0}:\mu_{1}-\mu_{2}\leqslant2,\quad H_{1}:\mu_{1}-\mu_{2}>2. $$ 设各对数据的差 $ D_{i}=X_{i}-Y_{i}\quad(i=1,2,\cdots,8) $ 是来自正态总体 $ N(\mu_{D},\sigma_{D}^{2}) $ 的样本, $ \mu_{D},\sigma_{D}^{2} $ 均未知.问是否可以认为早晨的身高比晚上的身高要高(取 $ \alpha=0.05 $)? 设 $ D_{i}=X_{i}-Y_{i}\quad(i=1,2,\cdots,15) $ 是来自正态总体 $ N(\mu_{D},\sigma_{D}^{2}) $ 的样本, $ \mu_{D},\sigma_{D}^{2} $ 均未知.问是否可以认为以材料A制成的后跟比材料B的耐穿(取 $ \alpha=0.05 $)? 设 $ D_{i}=X_{i}-Y_{i}\quad(i=1,2,\cdots,10) $ 是来自正态总体 $ N(\mu_{D},\sigma_{D}^{2}) $ 的样本, $ \mu_{D},\sigma_{D}^{2} $ 均未知。问以这两种种子种植的谷物的产量是否有显著的差异(取 $ \alpha=0.05 $)? $$ \begin{array}{l} 90 \quad 105 \quad 101 \quad 95 \quad 100 \quad 100 \quad 101 \quad 105 \quad 93 \quad 97 \end{array} $$ 考察提纯后群体是否比原群体整齐?取显著性水平 $ \alpha=0.01 $,并设小麦株高服从 $ N(\mu,\sigma^{2}) $。 $$ H_{0}:\sigma^{2}=0.11^{2},\quad H_{1}:\sigma^{2}\neq0.11^{2}. $$ $$ H_{0}:\sigma\geqslant0.04\%,\quad H_{1}:\sigma<0.04\% $$ $$ H_{0}:\sigma_{1}^{2}=\sigma_{2}^{2},\quad H_{1}:\sigma_{1}^{2}\neq\sigma_{2}^{2}, $$ 以说明在第6题中我们假设 $ \sigma_{1}^{2}=\sigma_{2}^{2} $是合理的. $$ H_{0}:\sigma_{1}^{2}=\sigma_{2}^{2},\quad H_{1}:\sigma_{1}^{2}\neq\sigma_{2}^{2}. $$ 以说明在第7题中我们假设 $ \sigma_{1}^{2}=\sigma_{2}^{2} $是合理的. 设两样本依次来自正态总体 $ N(\mu_{1}, \sigma_{1}^{2}) $, $ N(\mu_{2}, \sigma_{2}^{2}) $, $ \mu_{i}, \sigma_{i} $ (i=1,2)均未知,两样本相互独立。 (1)试检验假设 $ H_{0}:\sigma_{1}^{2}=\sigma_{2}^{2}, H_{1}:\sigma_{1}^{2}\neq\sigma_{2}^{2} $ (取 $ \alpha=0.05 $). (2)若能接受 $ H_{0} $,接着检验假设 $ H_{0}^{\prime}:\mu_{1}=\mu_{2},H_{1}^{\prime}:\mu_{1}\neq\mu_{2} $(取 $ \alpha=0.05 $). 设数据Ⅰ,Ⅱ分别来自正态总体 $ N(\mu_{1},\sigma_{1}^{2}) $, $ N(\mu_{2},\sigma_{2}^{2}) $, $ \mu_{1},\mu_{2},\sigma_{1}^{2},\sigma_{2}^{2} $均未知,两样本独立. (1)试检验假设 $ H_{0}:\sigma_{1}^{2}=\sigma_{2}^{2}, H_{1}:\sigma_{1}^{2}\neq\sigma_{2}^{2} $ (取 $ \alpha=0.1 $). (2)若能接受 $ H_{0} $,接着检验假设 $ H_{0}^{\prime}:\mu_{1}=\mu_{2},H_{1}^{\prime}:\mu_{1}\neq\mu_{2} $(取 $ \alpha=0.1 $). $$ H_{0}:\sigma_{1}^{2}\leqslant\sigma_{2}^{2},\quad H_{1}:\sigma_{1}^{2}>\sigma_{2}^{2}. $$ $$ H_{0}:\mu{\geqslant}15,\quad H_{1}:\mu<15. $$ 已知 $ \sigma^{2}=2.5 $. 取 $ \alpha=0.05 $. 若要求当 $ H_{1} $ 中的 $ \mu\leq13 $ 时犯第Ⅱ类错误的概率不超过 $ \beta=0.05 $,求所需的样本容量. $$ \begin{array}{l} 108 \quad 124 \quad 124 \quad 106 \quad 138 \quad 163 \quad 159 \quad 134 \end{array} $$ 设数据来自正态总体 $ N(\mu,\sigma^{2}),\mu,\sigma^{2} $ 未知. (1)试检验假设 $ H_{0}:\mu\leq125,H_{1}:\mu>125 $ ,取 $ \alpha=0.05 $ (2)若要求在上述 $ H_{1} $中 $ (\mu-125)/\sigma\geqslant1.4 $时,犯第Ⅱ类错误的概率不超过 $ \beta=0.1 $,求所需的样本容量. $$ H_{0}:\mu_{1}\leqslant2\mu_{2},\quad H_{1}:\mu_{1}>2\mu_{2}. $$ 此处 $ \mu_{1},\mu_{2} $ 分别是服用原有止痛片和服用新止痛片后至起作用的时间间隔的总体的均值.设两总体均为正态且方差分别为已知值 $ \sigma_{1}^{2},\sigma_{2}^{2} $.现分别在两总体中取一样本 $ X_{1},X_{2},\cdots,X_{n_{1}} $ 和 $ Y_{1},Y_{2},\cdots,Y_{n_{2}} $,设两个样本独立.试给出上述假设 $ H_{0} $ 的拒绝域,取显著性水平为 $ \alpha $. 问能否认为一页的印刷错误的个数服从泊松分布(取 $ \alpha=0.05 $). 取 $ \alpha=0.05 $,试检验假设 $ H_{0} $: 灯泡寿命服从指数分布 $$ f(t)=\{\begin{aligned}&0.005\mathrm{e}^{-0.005t},&t\geqslant0,\\ &0,&t<0.\end{aligned}. $$ (1)画出数据的直方图. (2)试取 $ \alpha=0.1 $检验数据来自正态总体 $ N(60,15^{2}) $ 试取 $ \alpha=0.05 $检验假设 $ H_{0} $: X 服从超几何分布 $$ \mathrm{P}\{X=k\}=\binom{5}{k}\binom{3}{3-k}\bigg/\binom{8}{3},\quad k=0,1,2,3. $$ 即检验假设 $ H_{0} $: 红球的只数为 5. 试取 $ \alpha=0.05 $,检验各类鱼数量的比例较10年前是否有显著的改变. $$ P\{X=x\}=p^{x-1}(1-p),\quad x=1,2,\cdots. $$ 今获得一样本如下: (1) 求 p 的最大似然估计值. (2)取 $ \alpha=0.05 $,检验假设: $ H_{0} $:数据来自总体 $ P\{X=x\}=p^{x-1}(1-p) $, $ x=1,2,\cdots $ 设两样本独立且1,2两队队员行李重量总体的概率密度至多差一个平移。记两总体的均值分别为 $ \mu_{1},\mu_{2} $,且 $ \mu_{1},\mu_{2} $均未知。试检验假设: $ H_{0}:\mu_{1}=\mu_{2},H_{1}:\mu_{1}<\mu_{2} $(取 $ \alpha=0.05 $)。 设两样本独立且数据所属的两总体的概率密度至多差一个平移,试问能否认为型号 A 的计算器平均使用时间比型号 B 来得长 ( $ \alpha=0.01 $)? 设两样本独立且数据所属的两总体的概率密度至多差一个平移. 问能否认为工人 A、工人 B 平均每天完成的件数没有显著差异 ( $ \alpha=0.1 $)? (2)考察生长在老鼠身上的肿块的大小。以 X 表示在老鼠身上生长了 15 天的肿块的直径(以 mm 计),设 $ X \sim N(\mu, \sigma^{2}) $, $ \mu, \sigma^{2} $ 均未知。今随机地取 9 只老鼠(在它们身上的肿块都长了 15 天),测得 $ \overline{x} = 4.3, s = 1.2 $,试取 $ \alpha = 0.05 $,用 p 值检验法检验假设 $ H_{0}: \mu = 4.0, H_{1}: \mu \neq 4.0 $,求出 p 值。 (3)用 p 值检验法检验 2 例 4 的检验问题. (4)用 p 值检验法检验第 27 题中的检验问题.

§2 正态总体均值的假设检验
(一) 单个总体 $ N(\mu,\sigma^{2}) $ 均值 $ \mu $ 的检验
(二) 两个正态总体均值差的检验(t检验)

(三) 基于成对数据的检验(t检验)
x(%) 0.20 0.30 0.40 0.50 0.60 0.70 0.80 0.90 1.00 y(%) 0.10 0.21 0.52 0.32 0.78 0.59 0.68 0.77 0.89 d=x-y(%) 0.10 0.09 -0.12 0.18 -0.18 0.11 0.12 0.13 0.11 红光(x) 0.30 0.23 0.41 0.53 0.24 0.36 0.38 0.51 绿光(y) 0.43 0.32 0.58 0.46 0.27 0.41 0.38 0.61 d=x-y -0.13 -0.09 -0.17 0.07 -0.03 -0.05 0.00 -0.10 §3 正态总体方差的假设检验
(一) 单个总体的情况

原假设 $ H_0 $ 检验统计量 备择假设 $ H_1 $ 拒绝域 1 $ \mu \le \mu_0 $\n $ \mu \ge \mu_0 $\n $ \mu = \mu_0 $\n( $ \sigma^2 $ 已知) $ Z = \frac{\overline{X} - \mu_0}{\sigma/\sqrt{n}} $ $ \mu > \mu_0 $\n $ \mu < \mu_0 $\n $ \mu \ne \mu_0 $ $ z \ge z_a $\n $ z \le -z_a $\n $ z \ge z_{a/2} $ 2 $ \mu \le \mu_0 $\n $ \mu \ge \mu_0 $\n $ \mu = \mu_0 $\n( $ \sigma^2 $ 未知) $ t = \frac{\overline{X} - \mu_0}{S/\sqrt{n}} $ $ \mu > \mu_0 $\n $ \mu < \mu_0 $\n $ \mu \ne \mu_0 $ $ t \ge t_a $(n-1)\n $ t \le -t_a $(n-1)\n $ t \ge t_{a-2} $(n-1) 3 $ \mu_1 - \mu_2 \le \delta $\n $ \mu_1 - \mu_2 \ge \delta $\n $ \mu_1 - \mu_2 = \delta $\n( $ \sigma_1^2 $, $ \sigma_2^2 $ 已知) $ Z = \frac{\overline{X} - \overline{Y} - \delta}{\sqrt{\frac{\sigma_1^2}{n_1} + \frac{\sigma_2^2}{n_2}}} $ $ \mu_1 - \mu_2 > \delta $\n $ \mu_1 - \mu_2 < \delta $\n $ \mu_1 - \mu_2 \ne \delta $ $ z \ge z_a $\n $ z \le -z_a $\n $ z \ge z_{a/2} $ 续表 原假设 $ H_0 $ 检验统计量 备择假设 $ H_1 $ 拒绝域 4 $ \mu_1 - \mu_2 \le \delta $\n $ \mu_1 - \mu_2 \ge \delta $\n $ \mu_1 - \mu_2 = \delta $\n $ (\sigma_1^2 = \sigma_2^2 = \sigma^2 $ 未知) $ t = \frac{\overline{X} - \overline{Y} - \delta}{S_w \sqrt{\frac{1}{n_1} + \frac{1}{n_2}}} $\n $ S_w^2 = \frac{(n_1 - 1)S_1^2 + (n_2 - 1)S_2^2}{n_1 + n_2 - 2} $ $ \mu_1 - \mu_2 > \delta $\n $ \mu_1 - \mu_2 < \delta $\n $ \mu_1 - \mu_2 \ne \delta $ $ t \ge t_a (n_1 + n_2 - 2) $\n $ t \le -t_a (n_1 + n_2 - 2) $\n $ t \ge t_{a/2} (n_1 + n_2 - 2) $ 5 $ \sigma^2 \le \sigma_0^2 $\n $ \sigma^2 \ge \sigma_0^2 $\n $ \sigma^2 = \sigma_0^2 $\n( $ \mu $ 未知) $ \chi^2 = \frac{(n-1)S^2}{\sigma_0^2} $ $ \sigma^2 > \sigma_0^2 $\n $ \sigma^2 < \sigma_0^2 $\n $ \sigma^2 \ne \sigma_0^2 $ $ \chi^2 \ge \chi_a^2 (n-1) $\n $ \chi^2 \le \chi_{a/2}^2 (n-1) $\n $ \chi^2 \ge \chi_{a/2}^2 (n-1) $ 6 $ \sigma_1^2 \le \sigma_2^2 $\n $ \sigma_1^2 \ge \sigma_2^2 $\n $ \sigma_1^2 = \sigma_2^2 $\n( $ \mu_1, \mu_2 $ 未知) $ F = \frac{S_1^2}{S_2^2} $ $ \sigma_1^2 > \sigma_2^2 $\n $ \sigma_1^2 < \sigma_2^2 $\n $ \sigma_1^2 \ne \sigma_2^2 $ $ F \ge F_a (n_1 - 1, n_2 - 1) $\n $ F \le F_{1-a} (n_1 - 1, n_2 - 1) $\n $ F \ge F_{a/2} (n_1 - 1, n_2 - 1) $或 $ F \le F_{1-a/2} (n_1 - 1, n_2 - 1) $ 7 $ \mu_D \le 0 $\n $ \mu_D \ge 0 $\n $ \mu_D = 0 $\n(成对数据) $ t = \frac{\overline{D} - 0}{S_D / \sqrt{n}} $ $ \mu_D > 0 $\n $ \mu_D < 0 $\n $ \mu_D \ne 0 $ $ t \ge t_a (n-1) $\n $ t \le -t_a (n-1) $\n $ t \ge t_{a/2} (n-1) $ (二) 两个总体的情况
§4 置信区间与假设检验之间的关系
§5 样本容量的选取
1. Z 检验法的 OC 函数


2. t 检验法的 OC 函数
燃料 A 81 84 79 76 82 83 84 80 79 82 81 79 燃料 B 76 74 78 79 80 79 82 76 81 79 82 78 §6 分布拟合检验
(一) 单个分布的 $ \chi^{2} $ 拟合检验法
i 0 1 2 3 4 5 ≥6 $ f_{i} $ 22 37 20 13 6 2 0 $ A_{i} $ $ A_{0} $ $ A_{1} $ $ A_{2} $ $ A_{3} $ $ A_{4} $ $ A_{5} $ $ A_{6} $ $A_i$ $f_i$ $p_i$ $np_i$ $f_i^{\frac{2}{2}}/(np_i)$ $A_0:\{X=0\}$ 22 $e^{-1}$ 36.788 13.16 $A_1:\{X=1\}$ 37 $e^{-1}$ 36.788 37.21 $A_2:\{X=2\}$ 20 $e^{-1}/2$ 18.394 21.75 $A_3:\{X=3\}$ 13 $e^{-1}/6$ 6.131 $A_4:\{X=4\}$ 6 $e^{-1}/24$ 1.533 $A_5:\{X=5\}$ 2 $e^{-1}/120$ 0.307 54.92 $A_6:\{X\geq6\}$ 0 $1-\sum_{i=0}^{5}p_i$ 0.059 序号 1 2 3 4 种类 黑色无角 黑色有角 红色无角 红色有角 数量 192 78 72 18 $ A_{i} $ $ A_{1} $ $ A_{2} $ $ A_{3} $ $ A_{4} $ X 1 2 3 4 $ p_{k} $ 9/16 3/16 3/16 1/16 $ A_i $ $ f_i $ $ p_i $ $ np_i $ $ f_i^2/(np_i) $ $ A_1 $ 192 9/16 $ 360 \times 9/16 = 202.5 $ $ 192^2/202.5 = 182.04 $ $ A_2 $ 78 3/16 $ 360 \times 3/16 = 67.5 $ $ 78^2/67.5 = 90.13 $ $ A_3 $ 72 3/16 $ 360 \times 3/16 = 67.5 $ $ 72^2/67.5 = 76.8 $ $ A_4 $ 18 1/16 $ 360 \times 1/16 = 22.5 $ $ 18^2/22.5 = 14.4 $ (二) 分布族的 $ \chi^{2} $拟合检验
i 0 1 2 3 4 5 6 7 8 9 10 11 $ \geq 12 $ $ f_i $ 1 5 16 17 26 11 9 9 2 1 2 1 0 $ A_i $ $ A_0 $ $ A_1 $ $ A_2 $ $ A_3 $ $ A_4 $ $ A_5 $ $ A_6 $ $ A_7 $ $ A_8 $ $ A_9 $ $ A_{10} $ $ A_{11} $ $ A_{12} $ $ A_i $ $ f_i $ $ \hat{p}_i $ $ n\hat{p}_i $ $ f_i^2/n\hat{p}_i $ $ A_0 $ 1 0.015 1.5 4.615 $ A_1 $ 5 0.063 6.3 19.394 $ A_2 $ 16 0.132 13.2 15.622 $ A_3 $ 17 0.185 18.5 34.845 $ A_4 $ 26 0.194 19.4 7.423 $ A_5 $ 11 0.163 16.3 7.105 $ A_6 $ 9 0.114 11.4 11.739 $ A_7 $ 9 0.069 6.9 $ A_8 $ 2 0.036 3.6 $ A_9 $ 1 0.017 1.7 $ A_{10} $ 2 0.007 0.7 5.538 $ A_{11} $ 1 0.003 0.3 $ A_{12} $ 0 0.002 0.2 相继两次地震间隔天数 x 0~4 5~9 10~14 15~19 20~24 25~29 30~34 35~39 ≥40 出现的频数 50 31 26 17 10 8 6 6 8① $ A_i $ $ f_i $ $ \hat{p}_i $ $ n\hat{p}_i $ $ f_i^2/(n\hat{p}_i) $ $ A_1:0 \le x \le 4.5 $ 50 0.2788 45.1656 55.3519 $ A_2:4.5 < x \le 9.5 $ 31 0.2196 35.5752 27.0132 $ A_3:9.5 < x \le 14.5 $ 26 0.1527 24.7374 27.3270 $ A_4:14.5 < x \le 19.5 $ 17 0.1062 17.2044 16.7980 $ A_5:19.5 < x \le 24.5 $ 10 0.0739 11.9718 8.3530 $ A_6:24.5 < x \le 29.5 $ 8 0.0514 8.3268 7.6860 $ A_7:29.5 < x \le 34.5 $ 6 0.0358 5.7996 6.2073 $ A_8:34.5 < x \le 39.5 $ 6 0.0248 4.0176 14.8269 $ A_9:39.5 < x < \infty $ 8 0.0568 9.2016 13.2192 $ A_i $ $ f_i $ $ \hat{p}_i $ $ n\hat{p}_i $ $ f_i^2/n\hat{p}_i $ $ A_1:x \le 129.5 $ 1 0.0087 0.73 4.91 $ A_2:129.5 < x \le 134.5 $ 4 0.0519 4.36 6.79 $ A_3:134.5 < x \le 139.5 $ 10 0.1752 14.72 41.55 $ A_4:139.5 < x \le 144.5 $ 33 0.3120 26.21 24.40 $ A_5:144.5 < x \le 149.5 $ 24 0.2811 23.61 10.02 $ A_6:149.5 < x \le 154.5 $ 9 0.1336 11.22 14.37 $ A_7:154.5 < x < \infty $ 3 0.0375 3.15 (三) 偏度、峰度检验
解 现在来检验假设
§7 秩和检验
秩 $ R_1 $ 秩 $ R_1 $ 秩 $ R_1 $ 秩 $ R_1 $ 秩 $ R_1 $ 123 6 136 10 167 14 247 13 356 14 124 7 137 11 234 9 256 13 357 15 125 8 145 10 235 10 257 14 367 16 126 9 146 11 236 11 267 15 456 15 127 10 147 12 237 12 345 12 457 16 134 8 156 12 245 11 346 13 467 17 135 9 157 13 246 12 347 14 567 18 $R_{1}$ 6 7 8 9 10 11 12 $P\{R_{1}=r_{1}\}$ 1/35 1/35 2/35 3/35 4/35 4/35 5/35 $P\{R_{1}\leq r_{1}\}$ 1/35 2/35 4/35 7/35 11/35 15/35 20/35 $R_{1}$ 13 14 15 16 17 18 $P\{R_{1}=r_{1}\}$ 4/35 4/35 3/35 2/35 1/35 1/35 $P\{R_{1}\leq r_{1}\}$ 24/35 28/35 31/35 33/35 34/35 1 不作治疗 1.9 0.5 0.9 2.1 接受治疗 3.1 5.3 1.4 4.6 2.8 数据 $ \underline{\text{0.5}} $ $ \underline{\text{0.9}} $ 1.4 $ \underline{\text{1.9}} $ $ \underline{\text{2.1}} $ 2.8 3.1 4.6 5.3 秩 1 2 3 4 5 6 7 8 9 A 7.0 3.5 9.6 8.1 6.2 5.1 10.4 4.0 2.0 10.5 B 5.7 3.2 4.2 11.0 9.7 6.9 3.6 4.8 5.6 8.4 10.1 化验员 A 82 73 91 84 77 98 81 79 87 85 化验员 B 80 76 92 86 74 96 83 79 80 75 79 数据 $ \underline{\text{73}} $ 74 75 76 $ \underline{\text{77}} $ $ \underline{\text{79}} $ 79 79 80 80 $ \underline{\text{81}} $ $ \underline{\text{82}} $ 83 $ \underline{\text{84}} $ $ \underline{\text{85}} $ 86 $ \underline{\text{87}} $ $ \underline{\text{91}} $ 92 96 $ \underline{\text{98}} $ 秩 1 2 3 4 5 7 7 7 9.5 9.5 11 12 13 14 15 16 17 18 19 20 21 §8 假设检验问题的 p 值检验法




小结
假设检验的两类错误 真实情况(未知) 所作决策 接受 $ H_0 $ 拒绝 $ H_0 $ $ H_0 $为真 正确 犯第Ⅰ类错误 $ H_0 $不真 犯第Ⅱ类错误 正确 重要术语及主题
习题
9.8 10.4 10.6 9.6 9.7 9.9 10.9 11.1 9.6 10.2 10.3 9.6 9.9 11.2 10.6 9.8 10.5 10.1 10.5 9.7
马克·吐温 0.225 0.262 0.217 0.240 0.230 0.229 0.235 0.217 斯诺特格拉斯 0.209 0.205 0.196 0.210 0.202 0.207 0.224 0.223 0.201
老过程 6 4 5 5 6 5 5 6 4 6 7 4 新过程 2 1 2 2 1 0 3 2 1 0 1 3
序号 1 2 3 4 5 6 7 8 早上 $ x_i $ 172 168 180 181 160 163 165 177 晚上 $ y_i $ 172 167 177 179 159 161 166 175
男子 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 材料A( $ x_{i} $) 6.6 7.0 8.3 8.2 5.2 9.3 7.9 8.5 7.8 7.5 6.1 8.9 6.1 9.4 9.1 材料B( $ y_{i} $) 7.4 5.4 8.8 8.0 6.8 9.1 6.3 7.5 7.0 6.5 4.4 7.7 4.2 9.4 9.1
土地编号i 1 2 3 4 5 6 7 8 9 10 种子A( $ x_i $) 23 35 29 42 39 29 37 34 35 28 种子B( $ y_i $) 26 39 35 40 38 24 36 27 41 27
x 101 100 99 99 98 100 98 99 99 99 y 100 98 100 99 98 99 98 98 99 100
工程杂志上的论文(数据Ⅰ) 未出版的学术报告(数据Ⅱ) 1.79 1.75 1.67 1.65 2.39 2.51 2.86 1.87 1.74 1.94 2.56 2.29 2.49 1.62 2.06 1.33 2.36 2.58 1.96 1.69 1.70 2.62 2.41
错误个数 $ f_i $ 0 1 2 3 4 5 6 ≥7 含 $ f_i $ 个错误的页数 36 40 19 2 0 2 1 0
寿命 t(h) 0≤t≤100 100<t≤200 200<t≤300 t>300 灯泡数 121 78 43 58
分数 x $ 20 \leq x \leq 30 $ $ 30 < x \leq 40 $ $ 40 < x \leq 50 $ $ 50 < x \leq 60 $ 学生数 5 15 30 51 分数 x $ 60 < x \leq 70 $ $ 70 < x \leq 80 $ $ 80 < x \leq 90 $ $ 90 < x \leq 100 $ 学生数 60 23 10 6
x 0 1 2 3 次数 1 31 55 25
序号 1 2 3 4 种类 鲑鱼 鲈鱼 竹夹鱼 鲇鱼 数量(条) 132 100 200 168 $ \sum=600 $
x 1 2 3 4 5 6 7 8 9 10 11 12 ≥13 观察到x的次数 48 31 20 9 6 5 4 2 1 1 2 1 0
1 队 34 39 41 28 33 2 队 36 40 35 31 39 36
型号A 5.5 5.6 6.3 4.6 5.3 5.0 6.2 5.8 5.1 5.2 5.9 型号B 3.8 4.3 4.2 4.0 4.9 4.5 5.2 4.8 4.5 3.9 3.7 4.6
工人A 49 52 53 47 50 工人B 56 48 58 46 55