本节对正态总体参数 μ 和 σ² 的各种检验分别进行讨论.
7.2.1 单个正态总体均值的检验
Section titled “7.2.1 单个正态总体均值的检验”设 x₁, x₂, …, xₙ 是来自 N(μ, σ²) 的样本,考虑如下三种关于 μ 的检验问题:
I H _ 0: μ ≤slant μ_ 0 vs H _ 1: μ > μ_ 0, (7.2.1) II H _ 0: μ ≥slant μ_ 0 text vs H _ 1: μ < μ_ 0, (7.2.2) III H _ 0: μ = μ_ 0 vs H _ 1: μ ≠ μ_ 0, (7.2.3)其中 μ₀ 是已知常数.由于正态总体含两个参数,总体方差 σ² 已知与否对检验有影响.下面我们分 σ 已知和未知两种情况叙述.
一、 σ = σ₀ 已知时的 u 检验
Section titled “一、 σ = σ₀σ=σ0 已知时的 uu 检验”对于(7.2.1)式所示的单侧检验问题Ⅰ,由于 μ 的点估计是 barx ,且 barx∼ N(μ,σ₀²/n) ,故选用检验统计量
u = frac overline x - μ_ 0σ_ 0 / sqrt n (7.2.4)是恰当的.直觉告诉我们: 当样本均值 barx 不超过设定均值 μ₀ 时,应倾向于接受原假设; 当样本均值 barx 超过 μ₀ 时,应倾向于拒绝原假设.可是,在有随机性存在的场合,如果 barx 比 μ₀ 大一点就拒绝原假设似乎不当,只有当 barx 比 μ₀ 大到一定程度时拒绝原假设才是恰当的.这就存在一个临界值 c,拒绝域为
W _ 1 = {(x _ 1, x _ 2, … , x _ n): u ≥slant c }, (7.2.5)常简记为 {u≥slant c}.若要求检验的显著性水平为 α ,则 c 满足
P _ μ_ 0 (u ≥slant c) = α .由于在 μ = μ₀ 时, u ∼ N(0,1) ,故 c = u₁ - α (见图 7.2.1(a)),最后的拒绝域为

(a) H₁ : μ >μ₀
(7.2.6)
(b) H₁ : μ < μ₀
(c) H₁ : μ ≠ μ₀
该检验用的检验统计量是 u 统计量,故一般称为 u 检验.该检验的势函数是 μ 的函数,它可用正态分布写出,具体如下: 对 μ ∈ (-∞, ∞) ,
r l g (μ) & = P _ μ (X ∈ W _ 1) = P _ μ (u ≥slant u _ 1 - α) & = P _ μ (frac bar x - μ_ 0σ_ 0 / sqrt n ≥slant u _ 1 - α) & = P _ μ (frac bar x - μ + μ - μ_ 0σ_ 0 / sqrt n ≥slant u _ 1 - α) & = P _ μ (frac bar x - μσ_ 0 / sqrt n ≥slant frac μ_ 0 - μσ_ 0 / sqrt n + u _ 1 - α) & = 1 - Φ (sqrt n (μ_ 0 - μ) / σ_ 0 + u _ 1 - α).由此可见,势函数是 μ 的增函数,其图形见图 7.2.2(a).由增函数性质知,只要 g(μ₀) = α ,就可保证在 μ ≤slant μ₀ 时有 g(μ)≤slant α . 所以上述求出的检验是显著性水平为 α 的检验.
下面我们讲述用 p 值进行检验的方法.
类似于7.1.3节的讲述,对给定的样本观测值,可以计算出相应的检验统计量u的值,记为 u₀=frac√(n)(barx-μ₀)σ₀,这里的 barx 是样本观测值.因在 μ=μ₀ 时,u是服从标准正态分布

(a) H₁ : μ >μ₀

(c) H₁ : μ ≠ μ₀
的随机变量,令
此即说明 u₀ = u₁ - p ,于是由正态分布函数的反函数的单调性有如下结论: - 当 p₁ ≤slant α 时,u₁ - α ≤slant u₀,于是观测值落在拒绝域里,应拒绝原假设. - 当 p₁ > α 时, u₁ - α > u₀ ,于是观测值不在拒绝域里,应接受原假设。 由此可以看出,(7.2.7)计算出的值就是该检验的 p 值. 对检验问题(7.2.2)所示的单侧检验问题Ⅱ的讨论是完全类似的.仍选用 u 作为检验统计量,考虑到(7.2.2)的备择假设 H₁ 在左侧,其拒绝域(见图 7.2.1(b))为
W _ II = {u ≤slant u _ α }. (7.2.8)而检验的 p 值为
p _ II = P (u ≤slant u _ 0) = Φ (u _ 0), (7.2.9)u₀, u 的含义同上,后面还会用到就不再一一指出了. 对检验问题(7.2.3)所示的双侧检验问题Ⅲ,也可类似进行讨论,只不过检验的p值稍有不同.仍选用u作为检验统计量,考虑到(7.2.3)的备择假设 H₁ 分散在两侧,故其拒绝域亦应在两侧,即拒绝域应有如下形式
W _ text III = {| u | ≥slant c }.对给定的显著性水平 α ( 0 < α < 1 ),由 P_μ₀(|u| ≥slant c) = α 可定出 c = u₁ - α / 2 (见图 7.2.1(c)),最后的拒绝域为
W _ text III = {| u | ≥slant u _ 1 - α / 2 }. (7.2.10)下面介绍双侧检验的 p 值的计算.在检验统计量分布对称场合,双侧检验的 p 值的计算与单侧检验是类似的,不对称场合我们在后面介绍. 仿上,令
p _ III = P (| u | ≥slant | u _ 0 |) = 2 (1 - Φ (| u _ 0 |)) , (7.2.11)此即说明 |u₀| = u₁ - p_(III) / 2 ,这里要用到 u₀ 的绝对值是因为对双侧假设检验,观测值可能为正,也可能为负,二者机会相同,于是有类似的结论: - 当 p_(III) ≤slant α 时,u₁ - α / 2 ≤slant |u₀|,于是观测值落在拒绝域里,应拒绝原假设. - 当 p_(III) > α 时,u₁ - α / 2 > |u₀|,于是观测值不在拒绝域里,应接受原假设. 由此可以看出,(7.2.11)计算出的值就是该检验的 p 值.
从甲地发送一个信号到乙地.设乙地接收到的信号值是一个服从正态分布 N(μ, 0.2²) 的随机变量,其中 μ 为甲地发送的真实信号值.现甲地重复发送同一信号 5 次,乙地接收到的信号值为
l l l l l 8. 0 5 & 8. 1 5 & 8. 2 & 8. 1 & 8. 2 5,设接收方有理由猜测甲地发送的信号值为8,问能否接受这猜测?
解
这是一个假设检验的问题,总体 X ∼ N(μ, 0.2²),待检验的原假设 H₀ 与备择假设 H₁ 分别为
H _ 0: μ = 8 vs H _ 1: μ ≠ 8.这是一个双侧检验问题,检验的拒绝域为 {|u|≥slant u₁-α/2}.取显著性水平 α=0.05,则查表知 u₀.975=1.96.由该例中观测值可计算得出
overline x = 8. 1 5, u _ 0 = sqrt 5 (8. 1 5 - 8) / 0. 2 = 1. 6 8,u₀ 值未落入拒绝域 {|u|≥slant 1.96} 内,故不能拒绝原假设,即接受原假设,可认为猜测成立. 我们也可以采用 p 值完成此次检验.此处 u₀ = 1.68,根据(7.2.11)式,
p = 2 (1 - Φ (1. 6 8)) = 0. 0 9 3.由于 p 值大于事先给定的水平 0.05, 故不能拒绝原假设,结论是相同的. 进一步,我们从 p 值还可以看到,只要事先给定的显著性水平不高于0.093,则都不能拒绝原假设;而若事先给定的显著性水平高于0.093,如事先给定的显著性水平为0.10,则检验就会作出拒绝原假设的结论. 说明:在实际中也经常会遇到如下两个检验问题:
l l l l IV & H _ 0: μ = μ_ 0 & text vs & H _ 1: μ > μ_ 0, V & H _ 0: μ = μ_ 0 & text vs & H _ 1: μ < μ_ 0.它仍可用检验统计量 u 施行检验.检验问题Ⅳ的拒绝域与检验问题Ⅰ的拒绝域相同,即 W_(N) = {u ≥slant u₁-α}.这是因为检验问题Ⅳ与Ⅰ的备择假设相同,而Ⅳ的原假设是Ⅰ的原假设的子集,由于此时 u 检验的势函数是 μ 的单调增函数,因此,检验问题Ⅳ的显著性水平为 α 的检验与检验问题Ⅰ的显著性水平为 α 的检验是相同的,从而拒绝域也相同,它们的检验的 p 值也相同.类似地,检验问题 V 与检验问题Ⅱ的拒绝域以及 p 值也是相同的.这个现象在以后其他检验中也会出现,结论是相似的.由此,本文中不再考虑诸如 IV, V 的检验问题 (若出现此类检验问题,可归结为检验问题 I, II 处理).
二、 σ 未知时的 t 检验
Section titled “二、 σσ 未知时的 tt 检验”对检验问题 I,由于 σ 未知,无法使用 (7.2.4) 式作检验.一个自然的想法是将 (7.2.4) 式中未知的 σ 替换成样本标准差 s,这就形成 t 检验统计量
t = frac sqrt n (bar x - μ_ 0)s. (7.2.12)由推论5.4.2 知,在 μ = μ₀ 时,t ∼ t(n-1),从而检验问题 I 的拒绝域为
W _ 1 = {t ≥slant t _ 1 - α (n - 1) }, (7.2.13)检验的 p 值是类似的,对给定的样本观测值,可以计算出相应的检验统计量 t 的值,记为 t₀ = frac√(n)(barx - μ₀)s,这里的 barx, s 可由样本观测值算得.因 t 是服从自由度是 n - 1 的 t 分布的随机变量,则
p _ I = P (t ≥slant t _ 0). (7.2.14)对另两组检验问题的讨论是完全类似于上一小节的,罗列结果如下:检验问题Ⅱ的拒绝域为
W _ II = {t ≤slant t _ α (n - 1) }, (7.2.15)p 值为
p _ text I I = P (t ≤slant t _ 0). (7.2.16)检验问题Ⅲ的拒绝域为
W _ III = {| t | ≥slant t _ 1 - α / 2 (n - 1) }, (7.2.17)p 值为
p _ text III = P (| t | ≥slant | t _ 0 |), (7.2.18)同样可证明这三个检验都是显著性水平为 α 的检验.
某厂生产的某种铝材的长度服从正态分布,其均值设定为240 cm. 现从该厂抽取5件产品,测得其长度(单位: cm)为
l l l l l 2 3 9. 7 & 2 3 9. 6 & 2 3 9 & 2 4 0 & 2 3 9. 2,试判断该厂此类铝材的长度是否满足设定要求? 这是一个关于正态均值的双侧假设检验问题.原假设是 H₀: μ = 240,备择假设是 H₁: μ ≠ 240.由于 σ 未知,故采用 t 检验,其拒绝域为 {|t| ≥slant t₁ - α / 2(n - 1)},若取 α = 0.05,则查表得 t₀.975(4) = 2.7764.现由样本计算得到 barx = 239.5, s = 0.4,故
t _ 0 = sqrt 5 (2 3 9. 5 - 2 4 0) / 0. 4 = - 2. 7 9 5,由于 |t₀| = 2.795 > 2.7764 ,故拒绝原假设,认为该厂生产的铝材的长度不满足设定要求. 下面用 p 值再作一次检验.此处 t₀ = -2.795,记 t 是服从自由度是 4 的 t 分布的随机变量,则根据 (7.2.18),
p = P (| t | ≥slant 2. 7 9 5) = 2 P (t ≥slant 2. 7 9 5),利用统计软件(如 R 或 Excel)可计算出具体 p 值为 0.0491, 由于 p 值小于事先给定的显著性水平 0.05, 故拒绝原假设,结论是相同的. 综上,关于单个正态总体的均值的检验问题可汇总成表7.2.1. 表 7.2.1 单个正态总体均值的假设检验
| 检验法 | H₀ | H₁ | 检验统计量 | 拒绝域 | p值 |
| u检验(σ=σ0已知) | μ≤μ0 | μ>μ0 | {u≥u1-α} | 1-Φ(u0) | |
| μ≥μ0 | μ<μ0 | u=x-μ0/σ0/√n | {u≤uα} | Φ(u0) | |
| μ=μ0 | μ≠μ0 | {|u|≥u1-α/2} | 2(1-Φ(|u0|)) | ||
| t检验(σ未知) | μ≤μ0 | μ>μ0 | {t≥t1-α(n-1)} | P(t≥t0) | |
| μ≥μ0 | μ<μ0 | t=x-μ0/s/√n | {t≤tα(n-1)} | P(t≤t0) | |
| μ=μ0 | μ≠μ0 | {|t|≥t1-α/2(n-1)} | P(|t|≥|t0|) |
u₀ = √(n)(overlinex -μ₀) / σ₀,t₀ = √(n)(overlinex -μ₀) / s,u 是服从 N(0,1) 的随机变量, t 是服从 t(n - 1) 的随机变量.
7.2.2 假设检验与置信区间的关系
Section titled “7.2.2 假设检验与置信区间的关系”细心的读者可能会发现,这里用的检验统计量与6.6.3节中置信区间所用的枢轴量很相似,这不是偶然的,两者之间存在非常密切的关系,现以标准差未知场合为例叙述如下. 设 x₁, x₂, …, xₙ 是来自正态总体 N(μ, σ²) 的样本,现讨论在 σ 未知场合关于均值 μ 的检验问题。分三种情况: 首先考虑双侧检验问题Ⅲ,显著性水平为 α 的检验的接受域为
overline W _ III = {| overline x - μ_ 0 | ≤slant frac ssqrt n t _ 1 - α / 2 (n - 1) },它可以改写为
overline W _ III = {overline x - frac ssqrt n t _ 1 - α / 2 (n - 1) ≤slant μ_ 0 ≤slant overline x + frac ssqrt n t _ 1 - α / 2 (n - 1) },这里 μ₀ 并无限制,若让 μ₀ 在 (- ∞, ∞) 内取值,就可得到 μ 的 1 - α 置信区间 [overlinex ± (s)/(√(n)) t₁ - α / 2(n - 1)].反之,若有一个如上的 1 - α 置信区间,也可获得关于 H₀: μ = μ₀ 的显著性水平为 α 的显著性检验.所以,“正态均值 μ 的 1 - α 置信区间”与“关于 H₀: μ = μ₀ vs H₁: μ ≠ μ₀ 的双侧检验问题的显著性水平为 α 的检验”是一一对应的. 类似地考虑单侧检验问题Ⅰ,显著性水平为 α 的检验的接受域为
overline W _ I = {bar x - μ_ 0 ≤slant frac ssqrt n t _ 1 - α (n - 1) } = {μ_ 0 ≥slant bar x - frac ssqrt n t _ 1 - α (n - 1) },这就给出了参数 μ 的 1 - α 置信下限.反之,对上述给定的 μ 的 1 - α 置信下限,我们也可以得到关于 H₀: μ ≤slant μ₀ 的单侧检验问题的显著性水平为 α 的检验,它们之间也是一一对应的.同样,对单侧检验问题 II, 其显著性水平为 α 的检验与参数 μ 的 1 - α 置信上限也是一一对应的.
7.2.3 两个正态总体均值差的检验
Section titled “7.2.3 两个正态总体均值差的检验”设 x₁, x₂, …, xₘ 是来自正态总体 N(μ₁, σ₁²) 的样本, y₁, y₂, …, yₙ 是来自另一个正态总体 N(μ₂, σ₂²) 的样本,两个样本相互独立。考虑如下三类检验问题:
I H _ 0: μ_ 1 - μ_ 2 ≤slant 0 text vs H _ 1: μ_ 1 - μ_ 2 > 0. (7.2.19) II H _ 0: μ_ 1 - μ_ 2 ≥slant 0 vs H _ 1: μ_ 1 - μ_ 2 < 0. (7.2.20) III H _ 0: μ_ 1 - μ_ 2 = 0 vs H _ 1: μ_ 1 - μ_ 2 ≠ 0. (7.2.21)这里对常用的两种情形进行讨论.
一、 σ₁, σ₂ 已知时的两样本 u 检验
Section titled “一、 σ₁, σ₂σ1,σ2 已知时的两样本 uu 检验”此时 μ₁-μ₂ 的点估计 barx-bary 的分布完全已知,
overline x - overline y ∼ N (μ_ 1 - μ_ 2, frac σ_ 1 ^ 2m + frac σ_ 2 ^ 2n).由此可采用 u 检验方法,检验统计量为
u = frac bar x - bar ysqrt frac σ_ 1 ^ 2m + frac σ_ 2 ^ 2n.在 μ₁ = μ₂ 时, u ∼ N(0,1).检验的拒绝域取决于备择假设的具体内容.对 (7.2.19) 所示的检验问题 I,检验的拒绝域与 p 值分别为
W _ I = {u ≥slant u _ 1 - α }, p _ I = 1 - Φ (u _ 0),其中 u₀ = fracoverlinex - overlineysqrtfracσ₁²m + fracσ₂²n 是由样本计算得到的检验统计量的值.对(7.2.20)所示的检验问题 Ⅱ,检验的拒绝域与 p 值分别为
W _ II = {u ≤slant u _ α }, p _ II = Φ (u _ 0).对(7.2.21)所示的检验问题Ⅲ,检验的拒绝域与p值分别为
W _ III = {| u | ≥slant u _ 1 - α / 2 }, p _ III = 2 (1 - Φ (| u _ 0 |)).二、 σ₁=σ₂=σ 但未知时的两样本 t 检验 在 σ₁² = σ₂² = σ² 但未知时,首先
overline x - overline y ∼ N (μ_ 1 - μ_ 2, (frac 1m + frac 1n) σ^ 2),其次,由于
frac 1σ^ 2 ∑_ i = 1 ^ m (x _ i - bar x) ^ 2 ∼ χ^ 2 (m - 1), frac 1σ^ 2 ∑_ i = 1 ^ n (y _ i - bar y) ^ 2 ∼ χ^ 2 (n - 1),故 (1)/(σ²) (∑ (xᵢ - barx)² +∑ (yᵢ - bary)²)∼ χ² (m + n - 2) ,记
s _ w ^ 2 = frac 1m + n - 2 [ ∑_ i = 1 ^ m (x _ i - bar x) ^ 2 + ∑_ i = 1 ^ n (y _ i - bar y) ^ 2 ],于是有
t = frac (bar x - bar y) - (μ_ 1 - μ_ 2)s _ w sqrt frac 1m + frac 1n ∼ t (m + n - 2).当 μ₁=μ₂ 时,检验统计量为
t = frac overline x - overline ys _ w sqrt frac 1m + frac 1n.对检验问题Ⅰ,检验的拒绝域与p值分别为
W _ 1 = {t ≥slant t _ 1 - α (m + n - 2) }, p _ 1 = P (t ≥slant t _ 0),其中 t₀ = fracoverlinex - overlineys_w√(1/m + 1/n) 是由样本计算得到的检验统计量的值, t 是服从自由度是 n + m - 2 的 t 分布的随机变量.对检验问题Ⅱ, 检验的拒绝域与 p 值分别为
W _ II = {t ≤slant t _ α (m + n - 2) }, p _ II = P (t ≤slant t _ 0).对检验问题Ⅲ,检验的拒绝域与 p 值分别为
W _ text III = {| t | ≥slant t _ 1 - α / 2 (m + n - 2) }, p _ text III = P (| t | ≥slant | t _ 0 |).某厂铸造车间为提高铸件的耐磨性而试制了一种镍合金铸件以取代铜合金铸件,为此,从两种铸件中各抽取一个容量分别为 8 和 9 的样本,测得其硬度(一种耐磨性指标)为 镍合金:76.43 76.21 73.58 69.69 65.29 70.83 82.75 72.34, 铜合金:73.66 64.27 69.34 71.37 69.77 68.12 67.27 68.07 62.61. 根据专业经验,硬度服从正态分布,且方差保持不变,试在显著性水平 α=0.05 下判断镍合金的硬度是否有明显提高.
解
用 X 表示镍合金的硬度,Y 表示铜合金的硬度,则由假定,X ∼ N(μ₁, σ²) , Y ∼ N(μ₂, σ²),要检验的假设是: H₀ : μ₁ = μ₂ vs H₁ : μ₁ > μ₂.由于两者方差未知但相等,故采用两样本 t 检验,经计算,
overline x = 7 3. 3 9, overline y = 6 8. 2 7 5 6,∑_ i = 1 ^ 8 (x _ i - overline x) ^ 2 = 1 9 1. 7 9 5 8, ∑_ i = 1 ^ 9 (y _ i - overline y) ^ 2 = 9 1. 1 5 4 8,从而 s_w = √((1)/(8 + 9 - 2) (191.7958 + 91.1548)) = 4.3432,
t _ 0 = frac 7 3 . 3 9 - 6 8 . 2 7 5 64 . 3 4 3 2 × sqrt frac 18 + frac 19 = 2. 4 2 3 4,查表知 t₀.95(15) = 1.7531 ,由于 t > t₀.95(15) ,故拒绝原假设,可判断镍合金硬度有显著提高. 下面用 p 值再作一次检验.此处 t₀ = 2.4234,因 t 是服从自由度是 15 的 t 分布的随机变量,则
p = P (t ≥slant 2. 4 2 3 4),利用任一款统计软件(如 Excel)可计算出具体 p 值为 0.0142, 由于 p 值小于事先给定的显著性水平 0.05, 故拒绝原假设,结论是相同的. 利用假设检验与置信区间的关系对其他情况下的检验问题可仿 6.6.6 节中两正态总体均值差的置信区间类似进行,我们下面以表格形式列出(见表 7.2.2),而不作推导. 表 7.2.2 两个正态总体均值的假设检验
| 检验法 | H₀ | H₁ | 检验统计量 | 拒绝域 | p值 |
| u检验 | μ₁ ≤slant μ₂ | μ₁ >μ₂ | u = fracbarx - bary√((σ₁²)/(m) + (σ₂²)/(n)) | {u ≥slant u₁-α} | 1-Φ(u₁) |
| (σ₁, σ₂) | μ₁ ≥slant μ₂ | μ₁ < μ₂ | {u ≤slant u_α} | Φ(u₁) | |
| 已知) | μ₁ = μ₂ | μ₁ ≠ μ₂ | √((σ₁²)/(m) + (σ₂²)/(n)) | {|u| ≥slant u₁-α/2} | 2(1-Φ(|u₁|)) |
| t检验 | μ₁ ≤slant μ₂ | μ₁ >μ₂ | t = fracbarx - barys_w √(1/m + 1/n) | {t ≥slant t₁-α(m+n-2)} | P(T₁ ≥slant t₁) |
| (σ₁ = σ₂) | μ₁ ≥slant μ₂ | μ₁ < μ₂ | {t ≤slant t_α(m+n-2)} | P(T₁ ≤slant t₁) | |
| 未知) | μ₁ = μ₂ | μ₁ ≠ μ₂ | {|t| ≥slant t₁-α/2(m+n-2)} | P(|T₁| ≥slant |t₁|) | |
| 大样本u检验(m,n充分大) | μ₁ ≤slant μ₂ | μ₁ > μ₂ | u = fracbarx - bary√((sₓ²)/(m) + (s_y²)/(n)) | {u ≥slant u₁-α} | 1-Φ(u₂) |
| μ₁ ≥slant μ₂ | μ₁ < μ₂ | {u ≤slant u_α} | Φ(u₂) | ||
| μ₁ = μ₂ | μ₁ ≠ μ₂ | {|u| ≥slant u₁-α/2} | 2(1-Φ(|u₂|)) | ||
| 近似t检验(m,n不很大) | μ₁ ≤slant μ₂ | μ₁ > μ₂ | t = fracbarx - bary√((sₓ²)/(m) + (s_y²)/(n)) | {t ≥slant t₁-α(l)} | P(T₂ ≥slant t₂) |
| μ₁ ≥slant μ₂ | μ₁ < μ₂ | {t ≤slant t_α(l)} | P(T₂ ≤slant t₂) | ||
| μ₁ = μ₂ | μ₁ ≠ μ₂ | {|t| ≥slant t₁-α/2(l)} | P(|T₂| ≥slant |t₂|) |
u₁ = fracoverlinex - overlineysqrtfracσ₁²m + fracσ₂²n, u₂ = fracoverlinex - overlineysqrtfracsₓ²m + fracs_(y)²n, t₁ = fracoverlinex - overlineys_(w)√(1/m + 1/n), t₂ = fracoverlinex - overlineysqrtfracsₓ²m + fracs_(y)²n, T₁ 是服从自由度为 n + m - 2 的 t 分布 的随机变量,T₂ 是服从自由度为l的t分布的随机变量,l与 s_(w) 的表达式见6.6.6节.
7.2.4 成对数据检验
Section titled “7.2.4 成对数据检验”在对两个总体均值进行比较时,有时数据是成对出现的,此时若采用二样本 t 检验所得出的结论有可能是不对的,下面看一个例子.
为了比较两种谷物种子的优劣,特选取 10 块土质不全相同的土地,并将每块土地分为面积相同的两部分,分别种植这两种种子,施肥与田间管理在 20 小块土地上都是一样,下面是各小块上的单位产量:
| 土地 | 1 | 2 | 3 | 4 | 5 | 6 | 7 | 8 | 9 | 10 |
| 种子一的单位产量x | 23 | 35 | 29 | 42 | 39 | 29 | 37 | 34 | 35 | 28 |
| 种子二的单位产量y | 30 | 39 | 35 | 40 | 38 | 34 | 36 | 33 | 41 | 31 |
| 差d=x-y | -7 | -4 | -6 | 2 | 1 | -5 | 1 | 1 | -6 | -3 |
解
假定 x ∼ N(μ₁, σ₁²), y ∼ N(μ₂, σ₂²) ,且 x 与 y 独立,这里假定两个总体的方差相等是合理的。我们先用二样本 t 检验讨论此问题。为此,记两种种子的单位产量的样本均值分别为 overlinex, overliney ,样本方差分别为 sₓ², s_y² 。如今要对如下检验问题:
H _ 0: μ_ 1 = μ_ 2 vs H _ 1: μ_ 1 ≠ μ_ 2作出判断.在假定 σ₁² = σ₂² = σ² 下,采用二样本 t 检验,检验统计量 t₁ 与拒绝域 W₁ 分别是
r l t _ 1 & = frac overline x - overline ys _ w / sqrt n / 2, W _ 1 & = {| t _ 1 | > t _ 1 - α / 2 (2 n - 2) },其中 s_w² = (sₓ² + s_y²) / 2, α 是给定的显著性水平.由给出的数据可算得
bar x = 3 3. 1, bar y = 3 5. 7, s _ x ^ 2 = 3 3. 2 1 1 1, s _ y ^ 2 = 1 4. 2 3 3 3, s _ w ^ 2 = 2 3. 7 2 2 2,从而可算得两样本的 t 检验统计量的值 (s_w = √(23.7222) = 4.8705)
t _ 1 0 = frac 3 3 . 1 - 3 5 . 74 . 8 7 0 5 / sqrt 1 0 / 2 = - 1. 1 9 3 7. (7.2.22)若给定 α = 0.05 ,查表得 t₀.975(18) = 2.1009 ,由于 |t₁0| < 2.1009 ,故不应拒绝原假设,即认为两种种子的单位产量平均值没有显著差别.此处检验的 p 值为0.2480.
下面我们换一个角度来讨论此问题.在这个问题中出现了成对数据,同一块土地上用两种种子得两个产量,其差 dᵢ=xᵢ-yᵢ(i=1,2,…,10) 排除了土质差异这个不可控因素的影响,主要反映两种种子的优劣.对这种信息我们应加以利用. 在正态性假定下,d=x-y∼ N(μ,σ_(d)²),其中 μ=μ₁-μ₂,σ_(d)²=σ₁²+σ₂².原先要比较 μ₁ 与 μ₂ 的大小,如今则转化为考察 μ 是否为零,即考察如下检验问题:
H _ 0: μ = 0 vs H _ 1: μ ≠ 0,即把双样本的检验问题转化为单样本 t 检验问题.这时检验的 t 统计量为
t _ 2 = overline d / (s _ d / sqrt n),其中
overline d = frac 1n ∑_ i = 1 ^ n d _ i, s _ d = (frac 1n - 1 ∑_ i = 1 ^ n (d _ i - overline d) ^ 2) ^ 1 / 2.在给定显著性水平 α 下,该检验问题的拒绝域是
W _ 2 = {| t _ 2 | ≥slant t _ 1 - α / 2 (n - 1) },这就是成对数据的 t 检验. 在本例中可算得
n = 1 0, overline d = - 2. 6, s _ d = 3. 5 0 2 4,于是
t _ 2 0 = frac - 2 . 63 . 5 0 2 4 / sqrt 1 0 = frac - 2 . 61 . 1 0 7 6 = - 2. 3 4 7 5. (7.2.23)对给定的显著性水平 α = 0.05 ,可查表得 t₀.975(9) = 2.2622. 由于 |t₂0| > 2.2622 ,故应拒绝原假设 H₀:μ = 0 ,即可认为两种种子的平均单位产量有显著差异,此处检验的 p 值为0.0435.进一步,平均单位产量差的估计量为 hatμ = overlinex -overliney = -2.6 ,可见种子 y 要比种子 x 的平均单位产量高. 本问题中两种处理方法得到完全不同的结论,我们指出成对数据 t 检验方法更加合理.这是因为成对数据的差 dᵢ 已消除了试验单元(如土质)之间的差别,从而用于检验的标准差 s_(d)=3.5024 (见(7.2.23))已排除土质差异的影响,只保留种子间的差异.而二样本 t 检验中用于检验的标准差 s_(w)=4.8705 还含有土质差异,从而使得标准差增大,导致因子不显著.所以成对数据场合化为单样本 t 检验所作的结论更可信些.假如上述问题中 10 块土地的土质完全一样,即参加比较的试验单元完全一样,则用二样本 t 检验会更好一些,因为它可提供更多的自由度去估计误差. 应注意,成对数据的获得事先要作周密的安排(即试验设计).在获得成对数据时不能发生“错位”,从而准确获得“成对数据”的信息.
7.2.5 正态总体方差的检验
Section titled “7.2.5 正态总体方差的检验”一、单个正态总体方差的 χ² 检验
Section titled “一、单个正态总体方差的 χ²χ2 检验”设 x₁, x₂, …, xₙ 是来自 N(μ, σ²) 的样本,对方差亦可考虑如下三个检验问题:
I H _ 0: σ^ 2 ≤slant σ_ 0 ^ 2 vs H _ 1: σ^ 2 > σ_ 0 ^ 2, (7.2.24) II H _ 0: σ^ 2 ≥slant σ_ 0 ^ 2 text vs H _ 1: σ^ 2 < σ_ 0 ^ 2, (7.2.25) III H _ 0: σ^ 2 = σ_ 0 ^ 2 vs H _ 1: σ^ 2 ≠ σ_ 0 ^ 2, (7.2.26)其中 σ₀² 是已知常数.此处通常假定 μ 未知,它们采用的检验统计量是相同的,均为
χ^ 2 = (n - 1) s ^ 2 / σ_ 0 ^ 2. (7.2.27)在 σ² = σ₀² 时, χ² ∼ χ²(n - 1) ,于是,若取显著性水平为 α ,则对应三个检验问题的显著性水平为 α 的检验的拒绝域依次为
r l & W _ I = {χ^ 2 ≥slant χ_ 1 - α ^ 2 (n - 1) }, & W _ II = {χ^ 2 ≤slant χ_ α ^ 2 (n - 1) }, & W _ III = {χ^ 2 ≤slant χ_ α / 2 ^ 2 (n - 1) text 或 χ^ 2 ≥slant χ_ 1 - α / 2 ^ 2 (n - 1) }.χ² 分布是偏态分布,三种拒绝域形式见图 7.2.3.
(b) H₁ : σ² < σ₀²

(a) H₁ : σ² > σ₀²

(c) H₁ :
图 7.2.3 三种
图 7.2.3 三种 χ² 检验的拒绝域 (n-1=6, α=0.05)
我们亦可给出检验的 p 值.对单侧检验,想法是类似的,记 χ₀² = (n - 1)s² / σ₀² 是由样本计算得到的检验统计量的值,χ² 表示服从自由度为 n - 1 的 χ² 分布的随机变量,则检验问题 I, II 的 p 值分别为 p_(I) = P(χ² ≥slant χ₀²), p_(II) = P(χ² ≤slant χ₀²).对双侧检验则稍稍复杂一点,事实上,双侧检验的拒绝域在两侧,用 χ₀² 可算得两个尾部概率 P(χ² ≤slant χ₀²) 和 P(χ² ≥slant χ₀²),其和为 1, 其中必有一个 ≤slant 0.5.检验的注意力总放在拒绝域上,故应从中选一个小的与 α / 2 比较,从而检验问题 III 的 p 值为
对这样定义的 p 值,不难发现下述结论仍然是成立的: - 当 p_(mathbbI) ≤slant α 时,应拒绝原假设. - 当 p_(III) > α 时,应接受原假设. 最后说明一点,由于方差与标准差是一一对应关系,不等式 σ² ≤slant σ₀² 等价于 σ ≤slant σ₀,故上述讨论一样适用于对标准差的检验问题.
某类钢板每块的重量 X 服从正态分布,其一项质量指标是钢板重量(单位: N)的方差不得超过 0.016. 现从某天生产的钢板中随机抽取 25 块,得其样本方差s²=0.025,问该天生产的钢板重量的方差是否满足要求?
解
这是一个关于正态总体方差的单侧检验问题.原假设为 H₀: σ² ≤slant 0.016,备择假设为 H₁: σ² > 0.016,此处 n = 25,若取 α = 0.05,则查表知 χ₀.95²(24) = 36.415,现计算可得
χ_ 0 ^ 2 = frac (n - 1) s ^ 2σ_ 0 ^ 2 = frac 2 4 × 0 . 0 2 50 . 0 1 6 = 3 7. 5 > 3 6. 4 1 5.由此,在显著性水平0.05下,我们拒绝原假设,认为该天生产的钢板重量不符合要求. 下面再用 p 值作此次检验.此处 χ₀² = 37.5,记 χ² 是自由度为 24 的 χ² 分布的随机变量,则
p = P (χ^ 2 ≥slant 3 7. 5).利用任一款统计软件(如 Excel)可计算出具体 p 值为 0.0390, 由于 p 值小于事先给定的显著性水平 0.05, 故拒绝原假设,结论自然也是相同的.
二、两个正态总体方差比的 F 检验
Section titled “二、两个正态总体方差比的 FF 检验”设 x₁, x₂, …, xₘ 是来自 N(μ₁, σ₁²) 的样本, y₁, y₂, …, yₙ 是来自 N(μ₂, σ₂²) 的样本。考虑如下三个假设检验问题:
I H _ 0: σ_ 1 ^ 2 ≤slant σ_ 2 ^ 2 vs H _ 1: σ_ 1 ^ 2 > σ_ 2 ^ 2, (7.2.29) II H _ 0: σ_ 1 ^ 2 ≥slant σ_ 2 ^ 2 text vs H _ 1: σ_ 1 ^ 2 < σ_ 2 ^ 2, (7.2.30) III H _ 0: σ_ 1 ^ 2 = σ_ 2 ^ 2 vs H _ 1: σ_ 1 ^ 2 ≠ σ_ 2 ^ 2. (7.2.31)此处 μ₁,μ₂ 均未知,记 sₓ²,s_(y)² 分别是由 x₁,x₂,…,xₘ 算得的 σ₁² 的无偏估计和由 y₁,y₂,…,yₙ 算得的 σ₂² 的无偏估计(两个都是样本方差),则可建立如下的检验统计量:
F = frac s _ x ^ 2s _ y ^ 2. (7.2.32)当 σ₁² = σ₂² 时, F ∼ F(m - 1, n - 1) ,由此给出三个检验问题对应的拒绝域依次为
r l & W _ I = {F ≥slant F _ 1 - α (m - 1, n - 1) }, & W _ II = {F ≤slant F _ α (m - 1, n - 1) }, & W _ III = {F ≤slant F _ α / 2 (m - 1, n - 1) text 或 F ≥slant F _ 1 - α / 2 (m - 1, n - 1) }.此时检验的 p 值的讨论与前述 χ² 是相似的,记 F₀ = sₓ² / s_y² 是由样本计算得到的检验统计量的值, F 表示服从 F(m - 1, n - 1) 分布的随机变量,则检验问题 I,II,III 的 p 值分别为
l p _ I = P ( F ≥slant F _ 0 ) , p _ II = P ( F ≤slant F _ 0 ) , p _ III = 2 min { P ( F ≥slant F _ 0 ) , P ( F ≤slant F _ 0 ) }.甲、乙两台机床加工某种零件,零件的直径服从正态分布,总体方差反映了加工精度,为比较两台机床的加工精度有无差别,现从各自加工的零件中分别抽取 7 件产品和 8 件产品,测得其直径为
X (text 机 床 甲 ): 1 6. 2 1 6. 8 1 5. 8 1 5. 5 1 6. 7 1 5. 6 1 5. 8,Y (text 机 床 乙 : 1 5. 9 1 6. 0 1 6. 4 1 6. 1 1 6. 5 1 5. 8 1 5. 7 1 5. 0.这就形成了一个双侧假设检验问题,原假设是 H₀:σ₁²=σ₂²,备择假设为 H₁:σ₁²≠σ₂².此处 m=7, n=8, 经计算 sₓ²=0.2729, s_(y)²=0.2164,于是 F₀=(0.2729)/(0.2164)=1.261,若取 α=0.05,查表知 F₀.975(6,7)=5.12, F₀.025(6,7)=frac1F₀.975(7,6)=(1)/(5.70)=0.175.其拒绝域为
W = {F ≤slant 0. 1 7 5 text 或 F ≥slant 5. 1 2 }.由此可见,样本未落入拒绝域,即在显著性水平0.05下可以认为两台机床的加工精度无显著差异. 下面再用 p 值做此次检验.此处 F₀ = 1.261,记 F 为服从 F(6,7) 分布的随机变量,由于是双侧检验问题,故
p = 2 min {P (F ≥slant 1. 2 6 1), P (F ≤slant 1. 2 6 1) } = 2 min {0. 3 8 0 3, 0. 6 1 9 7 } = 0. 7 6 0 6,由于 p 值大于事先给定的显著性水平 0.05, 故不能拒绝原假设. 关于正态总体方差的假设检验汇总列于表 7.2.3 中. 表 7.2.3 正态总体方差的假设检验
| 检验法 | H₀ | H₁ | 检验统计量 | 拒绝域 | p值 |
| χ² 检验 | σ² ≤slant σ₀² | σ² >σ₀² | χ² ≥slant χ₁-α² (n-1) | P(χ² ≥slant χ₀²) | |
| σ² ≥slant σ₀² | σ² < σ₀² | χ² = ((n-1)s²)/(σ₀²) | χ² ≤slant χ_α² (n-1) | P(χ² ≤slant χ₀²) | |
| σ² = σ₀² | σ² ≠ σ₀² | χ² ≤slant χ_(α/2)² (n-1) 或 χ² ≥slant χ₁-α/2² (n-1) | 2min {P(χ² ≤slant χ₀²), P(χ² ≥slant χ₀²)} | ||
| F检验 | σ₁² ≤slant σ₂² | σ₁² >σ₂² | F ≥slant F₁-α (m-1,n-1) | P(F ≥slant F₀) | |
| σ₁² ≥slant σ₂² | σ₁² < σ₂² | F = (sₓ²)/(s_y²) | F ≤slant F_α (m-1,n-1) | P(F ≤slant F₀) | |
| σ₁² = σ₂² | σ₁² ≠ σ₂² | F ≤slant F_(α/2) (m-1,n-1) 或 F ≥slant F₁-α/2 (m-1,n-1) | 2min {P(F ≤slant F₀), P(F ≥slant F₀)} |
说明:本节习题均采用拒绝域的形式完成,在可以计算检验的 p 值时要求计算出 p 值.
- 有一批枪弹,出厂时,其初速率 v ∼ N(950, 100) (单位: m/s ). 经过较长时间储存,取 9 发进行测试,得样本值 (单位: m/s ) 如下:
据经验,枪弹经储存后其初速率仍服从正态分布,且标准差保持不变,问是否可认为这批枪弹的初速率有显著降低 (α=0.05) ? 2. 已知某炼铁厂铁水含碳量服从正态分布 N(4.55, 0.108²).现在测定了 9 炉铁水,其平均含碳量为 4.484, 如果铁水含碳量的方差没有变化,可否认为现在生产的铁水平均含碳量仍为 4.55 ( α = 0.05)? 3. 由经验知某零件质量 X ∼ N(15, 0.05²) (单位: g), 技术革新后,抽出 6 个零件,测得质量为 14.7 15.1 14.8 15.0 15.2 14.6. 已知方差不变,问平均质量是否仍为 15 g (取 α = 0.05 )? 4. 化肥厂用自动包装机包装化肥,每包的质量服从正态分布,其平均质量为 100 kg, 标准差为 1.2 kg. 某日开工后,为了确定这天包装机工作是否正常,随机抽取 9 袋化肥,称得质量 (单元: kg) 如下:
l l l l l l l l l 9 9. 3 & 9 8. 7 & 1 0 0. 5 & 1 0 1. 2 & 9 8. 3 & 9 9. 7 & 9 9. 5 & 1 0 2. 1 & 1 0 0. 5.设方差稳定不变,问这一天包装机的工作是否正常(取 α=0.05 )? 5. 设需要对某正态总体的均值进行假设检验
H _ 0: μ = 1 5 vs H _ 1: μ < 1 5.已知 σ²=2.5 ,取 α=0.05 ,若要求当 H₁ 中的 μ≤slant13 时犯第二类错误的概率不超过 0.05,求所需的样本容量. 6. 从一批钢管中抽取 10 根,测得其内径 (单位: mm) 为
l l l l l 1 0 0. 3 6 & 1 0 0. 3 1 & 9 9. 9 9 & 1 0 0. 1 1 & 1 0 0. 6 4 1 0 0. 8 5 & 9 9. 4 2 & 9 9. 9 1 & 9 9. 3 5 & 1 0 0. 1 0.设这批钢管内径服从正态分布 N(μ, σ²) ,试分别在下列条件下: (1) 已知 σ = 0.5 ; (2) σ 未知, 检验假设 (α = 0.05)
H _ 0: μ = 1 0 0 vs H _ 1: μ > 1 0 0.- 假定考生成绩服从正态分布,在某地一次数学统考中,随机抽取了 36 位考生的成绩,算得平均成绩为 66.5 分,标准差为 15 分,问在显著性水平 0.05 下,是否可以认为这次考试全体考生的平均成绩为 70 分?
- 一位中学校长在报纸上看到这样的报道:“这一城市的初中学生平均每周看 8 ~h 电视。”她认为她所在学校的学生看电视的时间明显小于该数字。为此她在该校随机调查了100个学生,得知平均每周看电视的时间 barx = 6.5 ~h ,样本标准差为 s = 2 ~h 。问是否可以认为这位校长的看法是对的(取 α = 0.05 )?
- 设在木材中抽出 100 根,测其小头直径,得到样本平均数为 barx = 11.2 ~cm,样本标准差 s = 2.6 ~cm,问能否认为该批木材小头的平均直径不低于 12 ~cm (取 α = 0.05 )?
- 考察一鱼塘中鱼的含汞量,随机地取 10 条鱼测得各条鱼的含汞量(单位:mg)为 0.8 1.6 0.9 0.8 1.2 0.4 0.7 1.0 1.2 1.1. 设鱼的含汞量服从正态分布 N(μ, σ²) ,试检验假设 H₀: μ ≤slant 1.2 vs H₁: μ > 1.2 (取 α = 0.10 ).
- 如果一个矩形的宽度 w 与长度 l 的比 w/l = 1/2 (√(5) - 1) ≈ 0.618 ,这样的矩形称为黄金矩形。下面列出从某工艺品工厂随机抽取的20个矩形宽度与长度的比值: 0.693 0.749 0.654 0.670 0.662 0.672 0.615 0.606 0.690 0.628 0.668 0.611 0.606 0.609 0.553 0.570 0.844 0.576 0.933 0.630. 设这一工厂生产的矩形的宽度与长度的比值总体服从正态分布,其均值为 μ,试检验假设(取 α=0.05 )
- 下面给出两种型号的计算器充电以后所能使用的时间(单位:h)的观测值: 型号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的计算器平均使用时间比型号B来得长(取 α=0.01 )?
- 从某锌矿的东、西两支矿脉中,各抽取样本容量分别为9与8的样本进行测试,得样本含锌平均数及样本方差如下:
若东、西两支矿脉的含锌量都服从正态分布且方差相同,问东、西两支矿脉含锌量的平均值是否可以看作一样(取 α=0.05 )? 14. 在针织品漂白工艺过程中,要考察温度对针织品断裂强力(主要质量指标)的影响.为了比较 70^(circ) C 与 80^(circ) C 的影响有无差别,在这两个温度下,分别重复做了 8 次试验,得数据 (单位: N) 如下: 70^(circ) C 时的强力:20.5 18.8 19.8 20.9 21.5 19.5 21.0 21.2, 80^(circ)C 时的强力:17.7 20.3 20.0 18.8 19.0 20.1 20.0 19.1. 根据经验,温度对针织品断裂强度的波动没有影响.问在70℃时的平均断裂强力与80℃时的平均断裂强力间是否有显著差别(假定断裂强力服从正态分布,取 α=0.05 )? 15. 一药厂生产一种新的止痛片,厂方希望验证服用新药片后至开始起作用的时间间隔较原有止痛片至少缩短一半,因此厂方提出需检验假设
H _ 0: μ_ 1 = 2 μ_ 2 vs H _ 1: μ_ 1 > 2 μ_ 2.此处 μ₁, μ₂ 分别是服用原有止痛片和服用新止痛片后至开始起作用的时间间隔的总体的均值.设两总体均为正态分布且方差分别为已知值 σ₁², σ₂² ,现分别在两总体中取一样本 x₁, x₂, …, xₙ 和 y₁, y₂, …, yₘ ,设两个样本独立.试给出上述假设检验问题的检验统计量及拒绝域. 16. 对冷却到 -0.72^(circ) C 的样品用 A, B 两种测量方法测量其融化到 0^(circ) C 时的潜热,数据如下: 方法 A: 79.98 80.04 80.02 80.04 80.03 80.03 80.04 79.97 80.05 80.03 80.02 80.00 80.02, 方法 B: 80.02 79.94 79.98 79.97 80.03 79.95 79.97 79.97. 假设它们服从正态分布,方差相等,试检验:两种测量方法的平均性能是否相等(取 α=0.05 )? 17. 为了比较测定活水中氯气含量的两种方法,特在各种场合收集到 8 个污水水样,每个水样均用这两种方法测定氯气含量(单位: mg/l), 具体数据如下:
| 水样号 | 方法一(x) | 方法二(y) | 差(d=x-y) |
| 1 | 0.36 | 0.39 | -0.03 |
| 2 | 1.35 | 0.84 | 0.51 |
| 3 | 2.56 | 1.76 | 0.80 |
| 4 | 3.92 | 3.35 | 0.57 |
| 5 | 5.35 | 4.69 | 0.66 |
| 6 | 8.33 | 7.70 | 0.63 |
| 7 | 10.70 | 10.52 | 0.18 |
| 8 | 10.91 | 10.92 | -0.01 |
- 测得两批电子器件的样品的电阻(单位: Ω )为 A 批 (x) : 0.140 0.138 0.143 0.142 0.144 0.137; B 批(y): 0.135 0.140 0.142 0.136 0.138 0.140. 设这两批器材的电阻值分别服从分布 N(μ₁, σ₁²) , N(μ₂, σ₂²) ,且两样本独立. (1) 试检验两个总体的方差是否相等 (取 α=0.05 ). (2) 试检验两个总体的均值是否相等 (取 α=0.05 ).
- 某厂使用两种不同的原料生产同一类型产品,随机选取使用原料 A 生产的样品 22 件,测得平均质量为 2.36 kg, 样本标准差为 0.57 kg. 取使用原料 B 生产的样品 24 件,测得平均质量为 2.55 kg, 样本标准差为 0.48 kg. 设产品质量服从正态分布,两个样本独立.问能否认为使用原料 B 生产的产品平均质量较使用原料 A 显著大 (取 α = 0.05 )?