课程笔记

8.3两个正态总体的假设检验

阅读入口

设定与记号

设两组样本分别独立取自正态总体 N(μ1,σ12)N(\mu_1,\sigma_1^2)、N(μ2,σ22)N(\mu_2,\sigma_2^2),样本量为 n1,n2n_1,n_2,样本均值和方差分别记为 Xˉ,Yˉ,S12,S22\bar X,\bar Y,S_1^2,S_2^2。检验以第一组减第二组或第一组除第二组的顺序书写。

本节统一使用 Xˉ−Yˉ\bar X-\bar Y 的顺序,下面三种情形的拒绝域表都按这个顺序书写;换组时要同步改写整段公式。

均值差的检验

标准误怎么算,取决于两个总体的方差是否已知、是否假定相等。三种情形依次削弱条件。

两总体方差已知

两个 σ\sigma 都是已知常数,用

Z=Xˉ−Yˉ−Δ0σ12/n1+σ22/n2∼N(0,1).Z=\frac{\bar X-\bar Y-\Delta_0} {\sqrt{\sigma_1^2/n_1+\sigma_2^2/n_2}}\sim N(0,1).

分母是两个样本均值各自标准误的平方和再开方:Xˉ\bar X 的标准误为 σ1/n1\sigma_1/\sqrt{n_1},Yˉ\bar Y 的标准误为 σ2/n2\sigma_2/\sqrt{n_2},独立样本的方差可以直接相加。

方差未知但假定相等

此时两个总体方差相等,记公共方差为 σ2\sigma^2,可以用两组信息合成合并方差

Sp2=(n1−1)S12+(n2−1)S22n1+n2−2,S_p^2=\frac{(n_1-1)S_1^2+(n_2-1)S_2^2}{n_1+n_2-2},

再用它代替未知的 σ2\sigma^2:

T=Xˉ−Yˉ−Δ0Sp1n1+1n2∼t(n1+n2−2).T=\frac{\bar X-\bar Y-\Delta_0} {S_p\sqrt{\dfrac1{n_1}+\dfrac1{n_2}}}\sim t(n_1+n_2-2).

方差未知且不假定相等:Welch检验

不假定方差相等时,分母只能直接拼两个样本方差:

TW=Xˉ−Yˉ−Δ0S12/n1+S22/n2≈t(ν),T_W=\frac{\bar X-\bar Y-\Delta_0} {\sqrt{S_1^2/n_1+S_2^2/n_2}}\approx t(\nu),ν=(S12/n1+S22/n2)2(S12/n1)2n1−1+(S22/n2)2n2−1.\nu=\frac{\left(S_1^2/n_1+S_2^2/n_2\right)^2} {\dfrac{\left(S_1^2/n_1\right)^2}{n_1-1} +\dfrac{\left(S_2^2/n_2\right)^2}{n_2-1}}.

例如 n1=n2=10n_1=n_2=10、S12=8S_1^2=8、S22=2S_2^2=2 时,S12/n1=0.8S_1^2/n_1=0.8、S22/n2=0.2S_2^2/n_2=0.2,于是

ν=(0.8+0.2)20.82/9+0.22/9=10.68/9≈13.2,\nu=\frac{(0.8+0.2)^2}{0.8^2/9+0.2^2/9} =\frac{1}{0.68/9}\approx13.2,

比方差相等时假定的 n1+n2−2=18n_1+n_2-2=18 小得多。两组样本方差和样本量共同影响这个近似自由度。

拒绝域

三种情形只是统计量和自由度不同,拒绝域的形状一致。下表以 t(ν)t(\nu) 为例,ν\nu 取 n1+n2−2n_1+n_2-2(合并 tt)或Welch自由度;方差已知的 ZZ 情形把 tα(ν)t_\alpha(\nu) 换成 zαz_\alpha,ν\nu 不再出现。

备择假设拒绝域
μ1−μ2≠Δ0\mu_1-\mu_2\ne\Delta_0∣T∣>tα/2(ν)\lvert T\rvert>t_{\alpha/2}(\nu)
μ1−μ2>Δ0\mu_1-\mu_2>\Delta_0T>tα(ν)T>t_{\alpha}(\nu)
μ1−μ2<Δ0\mu_1-\mu_2<\Delta_0T<t1−α(ν)=−tα(ν)T<t_{1-\alpha}(\nu)=-t_{\alpha}(\nu)

这一段的三种情形合成一张拒绝域表,选型时按“题目给了多少信息”从上往下取第一行适用的,不必三种都算一遍。

方差比的F检验

检验 H0:σ12/σ22=r0H_0:\sigma_1^2/\sigma_2^2=r_0,使用

F=S12/S22r0∼F(n1−1,n2−1).F=\frac{S_1^2/S_2^2}{r_0}\sim F(n_1-1,n_2-1).

检验方差相等时 r0=1r_0=1,统计量就是两个样本方差之比。记 ν1=n1−1\nu_1=n_1-1、ν2=n2−1\nu_2=n_2-1,拒绝域为

  • 双侧 σ12/σ22≠r0\sigma_1^2/\sigma_2^2\ne r_0:F<F1−α/2(ν1,ν2)F<F_{1-\alpha/2}(\nu_1,\nu_2) 或 F>Fα/2(ν1,ν2)F>F_{\alpha/2}(\nu_1,\nu_2);
  • 右侧 σ12/σ22>r0\sigma_1^2/\sigma_2^2>r_0:F>Fα(ν1,ν2)F>F_{\alpha}(\nu_1,\nu_2);
  • 左侧 σ12/σ22<r0\sigma_1^2/\sigma_2^2<r_0:F<F1−α(ν1,ν2)F<F_{1-\alpha}(\nu_1,\nu_2)。

FF 分布不对称,双侧检验必须分别取左右两个临界值,不能写成“绝对值超过某个临界值”。若手头的表只给上侧分位点,用

F1−α(ν1,ν2)=1Fα(ν2,ν1)F_{1-\alpha}(\nu_1,\nu_2)=\frac{1}{F_{\alpha}(\nu_2,\nu_1)}

换算,两个自由度要同时交换位置。

方差比检验与8.2的方差检验是一对:单总体用 χ2\chi^2,两总体用 FF,自由度分别是 n−1n-1 与 (n1−1,n2−1)(n_1-1,n_2-1)。

学习安排

本节先讨论均值,再讨论方差。均值部分按条件强弱分三层——两个方差都已知、方差未知但假定相等、方差未知且不假定相等(Welch);随后把三种情形的拒绝域合并成一张表,再讲配对设计如何化为单个总体的差值检验;最后是方差比的 FF 检验。

选择哪一种统计量,取决于题目给了多少信息:给了两个 σ\sigma 就用 ZZ;只给样本方差并且明说方差相等,就用合并 tt;只说“方差未知”而不提相等,就用Welch。

配对样本检验

同一对象的前后测量、按特征匹配的个体对,都属于配对设计。两张表看起来是两组数据,但配对内的两个观测不独立,不能套用上面的独立两样本公式。

配对设计先逐对作差:

Di=Xi−Yi,i=1,…,n,D_i=X_i-Y_i,\qquad i=1,\ldots,n,

其中 nn 是对数。若 D1,…,DnD_1,\ldots,D_n 相互独立且都服从正态分布,则检验均值差 H0:μD=Δ0H_0:\mu_D=\Delta_0 就化为单个正态总体的均值检验:

T=Dˉ−Δ0SD/n∼t(n−1),T=\frac{\bar D-\Delta_0}{S_D/\sqrt n}\sim t(n-1),

其中 Dˉ\bar D、SDS_D 分别是差值样本的均值和标准差。

课程约定

  • 全章固定“第一组减第二组”:分子写 Xˉ−Yˉ−Δ0\bar X-\bar Y-\Delta_0,备择方向也按 μ1−μ2\mu_1-\mu_2 描述。
  • 合并方差的自由度记为 n1+n2−2n_1+n_2-2;Welch的自由度记为 ν\nu,手工查表时取相邻整数并注明是近似。
  • 配对检验先逐对作差,自由度按对数 n−1n-1 计。
  • 方差比检验默认 r0=1r_0=1;本课程的表只给上侧分位点,左侧临界值用 1/Fα(ν2,ν1)1/F_\alpha(\nu_2,\nu_1) 换算,两个自由度同时交换。

易错点

  • 把合并方差算成两个样本方差的算术平均。
  • 把配对检验的自由度写成 n1+n2−2n_1+n_2-2,或把两列数据的方差直接相加当作差值的方差。
  • 交换 FF 统计量的分子分母时只改了统计量,没有同步交换自由度与备择方向。
  • 双侧 FF 检验为了查表把较大的样本方差放到分子后,仍按 α\alpha 而不是 α/2\alpha/2 取临界值。
  • 把“先做 FF 检验、再按结果选 tt 检验”当成固定流程;FF 检验对非正态敏感,未拒绝也不等于方差相等。

复习要点

三种均值差情形的公式只差分母,复习时按“标准误怎么来”记:方差已知直接拼 σi2/ni\sigma_i^2/n_i;方差相等先按自由度加权合成 Sp2S_p^2,再乘 1/n1+1/n2\sqrt{1/n_1+1/n_2};方差不等直接拼 Si2/niS_i^2/n_i,自由度另外算出。Welch的自由度不超过 n1+n2−2n_1+n_2-2,两组容量相等且样本方差相等时取等号;自由度变小会让临界值变大,判别更保守。