M06 · 第 1 章 · 第一编 抽样与估计

样本、统计量与抽样分布:从重复抽样到正态枢轴

从统计模型、随机样本与经验分布出发,推导次序统计量、样本均值和均匀最大值的分布,建立正态样本下卡方、t、F 枢轴及充分性因子分解,并说明非独立抽样改变方差的边界。

报告页面错误
预备知识概率模型综合复习大数定律

本章目标

  1. 区分总体模型、随机样本、观测值、统计量与统计量的抽样分布。
  2. 由经验分布和次序统计量计算样本比例、极值概率与均匀最大值分布。
  3. 推导样本均值的期望、方差和标准误,并区分有限样本精确分布与渐近近似。
  4. 在正态 iid 条件下正确使用卡方、Student t 与 F 分布,并保留自由度。
  5. 用因子分解定理判断充分性,并识别无放回、分层和整群抽样造成的依赖。
页面阅读位置0% · 仅保存在此浏览器
章节未开始
本册完成进度0/6 章 · 0%
本页目录

统计模型把重复抽样写成概率对象

统计推断从观测数字追溯可能产生它们的概率模型族。记统计模型为 {Pθ:θΘ}\{P_\theta:\theta\in\Theta\},其中参数 θ\theta 在频率学解释下固定但未知。一次实验得到 x1,,xnx_1,\ldots,x_n;为了讨论同一实验重复进行时结果怎样波动,先把这些数还原为随机变量 X1,,XnX_1,\ldots,X_n。小写字母表示已经观察到的实现值,大写字母保留尚未实现的随机性。

随机样本、统计量与抽样分布

X1,,XnX_1,\ldots,X_n 相对于同一 PθP_\theta 独立同分布,则称它们为来自 PθP_\theta 的容量为 nn 的随机样本。统计量是形如

T=T(X1,,Xn)T=T(X_1,\ldots,X_n)

的可测函数,它只能含样本和已知常数,不能含未知参数。把观测值代入后得到统计量的实现值 t=T(x1,,xn)t=T(x_1,\ldots,x_n)。当同一抽样机制反复执行时,TTPθP_\theta 下形成的概率分布称为 TT 的抽样分布;该分布通常随 θ\theta 与样本量 nn 改变。

样本均值

X=1ni=1nXi\overline X=\frac1n\sum_{i=1}^nX_i

是统计量,因为公式中没有未知量。若总体均值 μ\mu 未知, i(Xiμ)2\sum_i(X_i-\mu)^2 不是统计量。无偏样本方差采用约定

S2=1n1i=1n(XiX)2.S^2=\frac1{n-1}\sum_{i=1}^n(X_i-\overline X)^2.

一次样本的直方图描述样本内部的数值分布;X\overline X 的抽样分布描述重新抽取整组样本后均值怎样变化。二者的横轴、随机单位和重复机制都不同。

经验分布把样本变成随机阶梯

经验分布函数

给定随机样本 X1,,XnX_1,\ldots,X_n,经验分布函数为

Fn(t)=1ni=1n1{Xit}.F_n(t)=\frac1n\sum_{i=1}^n\mathbf 1_{\{X_i\le t\}}.

对固定样本,FnF_n 是右连续的非降阶梯函数;每个观测贡献 1/n1/n 的跳跃,重复值使同一点的跳幅相加。

Fn(t)F_n(t) 同时是一项统计量和一条随机函数。固定 tt 时,指标变量 1{Xit}\mathbf 1_{\{X_i\le t\}} 是成功概率 F(t)F(t) 的 Bernoulli 变量,因此

Eθ[Fn(t)]=F(t),Varθ(Fn(t))=F(t)[1F(t)]n.E_\theta[F_n(t)]=F(t), \qquad \operatorname{Var}_\theta(F_n(t))=\frac{F(t)[1-F(t)]}{n}.

这个逐点结论说明经验分布在每个阈值处无偏,并给出 n1/2n^{-1/2} 量级的波动。它没有自动给出整条曲线的一致误差界;研究 suptFn(t)F(t)\sup_t|F_n(t)-F(t)| 需要更强的经验过程结论。

五次读数形成的经验阶梯

观测值为 2.1,1.4,2.1,3.0,0.92.1,1.4,2.1,3.0,0.9。排序后为 0.9,1.4,2.1,2.1,3.00.9,1.4,2.1,2.1,3.0,所以

Fn(t)={0,t<0.9,1/5,0.9t<1.4,2/5,1.4t<2.1,4/5,2.1t<3.0,1,t3.0.F_n(t)= \begin{cases} 0,&t<0.9,\\ 1/5,&0.9\le t<1.4,\\ 2/5,&1.4\le t<2.1,\\ 4/5,&2.1\le t<3.0,\\ 1,&t\ge3.0. \end{cases}

2.12.1 处有两个观测,跳幅为 2/52/5。例如样本中落在区间 (1.4,3.0](1.4,3.0] 的比例是

Fn(3.0)Fn(1.4)=125=35,F_n(3.0)-F_n(1.4)=1-\frac25=\frac35,

对应 2.1,2.1,3.02.1,2.1,3.0 三项。直接计数与阶梯函数增量一致,也核对了左端开、右端闭的约定。

次序统计量记录位置而非原编号

将样本从小到大排列,记为 X(1)X(n)X_{(1)}\le\cdots\le X_{(n)}X(1)X_{(1)}X(n)X_{(n)} 分别是样本最小值和最大值; X(k)X_{(k)} 是第 kk 个次序统计量。排序会丢掉原来的观测编号,却保留分位数、极差和端点信息。

若公共分布连续且具有密度 ff,第 kk 个次序统计量的密度为

fX(k)(x)=n!(k1)!(nk)!F(x)k1[1F(x)]nkf(x).f_{X_{(k)}}(x) =\frac{n!}{(k-1)!(n-k)!} F(x)^{k-1}[1-F(x)]^{n-k}f(x).

组合因子选择哪 k1k-1 个观测落在 xx 左侧、哪一个落入 xx 附近以及其余观测落在右侧。连续性排除了同点并列造成的额外分配;离散样本有正概率出现并列,应改用事件计数或 CDF 差分,不能直接沿用这个密度。

样本最大值的分布

X1,,XnX_1,\ldots,X_n iid,公共分布函数为 FF,并令 M=X(n)M=X_{(n)}。则对每个实数 mm

Pθ(Mm)=F(m)n.P_\theta(M\le m)=F(m)^n.

FF 在区间内具有密度 ff,则 MM 在相应区间内的密度为

fM(m)=nF(m)n1f(m).f_M(m)=nF(m)^{n-1}f(m).
证明

事件 {Mm}\{M\le m\} 当且仅当每个观测都不超过 mm。iid 假设给出

Pθ(Mm)=Pθ(X1m,,Xnm)=i=1nPθ(Xim)=F(m)n.P_\theta(M\le m) =P_\theta(X_1\le m,\ldots,X_n\le m) =\prod_{i=1}^nP_\theta(X_i\le m) =F(m)^n.

FF 可微处对 mm 求导,链式法则给出密度公式。若样本相互依赖,乘积步骤失效;边缘分布都等于 FF 仍不足以推出 FM=FnF_M=F^n

四个均匀观测的最大值

X1,,X4iidUnif(0,θ)X_1,\ldots,X_4\stackrel{\mathrm{iid}}\sim\operatorname{Unif}(0,\theta)θ>0\theta>0,并令 M=X(4)M=X_{(4)}。对 0mθ0\le m\le\theta

FM(m)=(mθ)4,fM(m)=4m3θ4.F_M(m)=\left(\frac m\theta\right)^4, \qquad f_M(m)=\frac{4m^3}{\theta^4}.

因此

Pθ(Mθ/2)=116,P_\theta(M\le\theta/2)=\frac1{16},

而期望为

Eθ[M]=0θm4m3θ4dm=4θ5.E_\theta[M] =\int_0^\theta m\frac{4m^3}{\theta^4}\,\mathrm dm =\frac{4\theta}{5}.

最大值平均低于总体上界,差额为 θ/5\theta/5。密度积分 0θ4m3/θ4dm=1\int_0^\theta4m^3/\theta^4\,\mathrm dm=1,并且 FM(0)=0,FM(θ)=1F_M(0)=0,F_M(\theta)=1,端点与归一化均成立。

样本均值的中心、尺度和形状

XiX_i iid,且 E[Xi]=μE[X_i]=\muVar(Xi)=σ2<\operatorname{Var}(X_i)=\sigma^2<\infty,则独立性使协方差交叉项消失:

E[X]=μ,Var(X)=1n2i=1nVar(Xi)=σ2n.E[\overline X]=\mu, \qquad \operatorname{Var}(\overline X) =\frac1{n^2}\sum_{i=1}^n\operatorname{Var}(X_i) =\frac{\sigma^2}{n}.

σ\sigma 描述单个观测的标准差,σ/n\sigma/\sqrt n 是样本均值的标准误。未知 σ\sigma 时常以 S/nS/\sqrt n 估计标准误,但这个随机估计值不等于真实抽样标准差。

均值的精确形状取决于总体分布。若总体正态,任意有限 nn 都有

XN ⁣(μ,σ2n).\overline X\sim N\!\left(\mu,\frac{\sigma^2}{n}\right).

若总体非正态但满足经典 iid 中心极限定理的有限正方差条件,则

n(Xμ)σdN(0,1).\frac{\sqrt n(\overline X-\mu)}{\sigma} \xrightarrow{d}N(0,1).

这是 nn\to\infty 的渐近结论,不会把有限样本均值精确变成正态。偏斜、重尾或强依赖会改变近似速度,甚至破坏所用定理的条件。

正态样本把平方和分解为卡方变量

正态样本均值与样本方差的精确分布

X1,,XniidN(μ,σ2)X_1,\ldots,X_n\stackrel{\mathrm{iid}}\sim N(\mu,\sigma^2),其中 n2n\ge2σ>0\sigma>0,则

Z=n(Xμ)σN(0,1),Q=(n1)S2σ2χn12,Z=\frac{\sqrt n(\overline X-\mu)}{\sigma}\sim N(0,1), \qquad Q=\frac{(n-1)S^2}{\sigma^2}\sim\chi^2_{n-1},

ZQZ\perp Q,等价地 XS2\overline X\perp S^2。因此

T=n(Xμ)S=ZQ/(n1)tn1.T=\frac{\sqrt n(\overline X-\mu)}{S} =\frac{Z}{\sqrt{Q/(n-1)}} \sim t_{n-1}.
证明

将标准化样本向量写成 Y=((X1μ)/σ,,(Xnμ)/σ)TY=((X_1-\mu)/\sigma,\ldots,(X_n-\mu)/\sigma)^{\mathsf T},则 YNn(0,In)Y\sim N_n(0,I_n)。把 YY 正交投影到由 (1,,1)T/n(1,\ldots,1)^{\mathsf T}/\sqrt n 张成的一维子空间,投影坐标正是 ZZ;投影到其正交补得到残差向量。球对称标准正态在正交变换后仍有相互独立的标准正态坐标,所以一维投影与其余 n1n-1 个坐标独立,而残差平方长度等于

i=1n(XiX)2σ2=(n1)S2σ2.\sum_{i=1}^n\frac{(X_i-\overline X)^2}{\sigma^2} =\frac{(n-1)S^2}{\sigma^2}.

它是 n1n-1 个独立标准正态平方之和,故服从 χn12\chi^2_{n-1}。独立比值的定义随即给出 tn1t_{n-1} 分布。

正态条件承担了关键作用。一般非正态总体下,X\overline XS2S^2 通常相关或至少不独立, QQ 也未必服从卡方分布。大样本 t 近似属于另一层渐近论证,不能改写成上述有限样本精确结论。

五个正态读数形成的 t 统计量

假设 XiiidN(μ,σ2)X_i\stackrel{\mathrm{iid}}\sim N(\mu,\sigma^2),观测为 8.2,9.1,10.0,10.7,12.08.2,9.1,10.0,10.7,12.0。计算得

x=10.0,i=15(xix)2=3.24+0.81+0+0.49+4.00=8.54,\overline x=10.0, \qquad \sum_{i=1}^5(x_i-\overline x)^2 =3.24+0.81+0+0.49+4.00=8.54,

所以 s2=8.54/4=2.135s^2=8.54/4=2.135s1.461164s\approx1.461164。若考察中心值 μ0=9\mu_0=9,则

t=5(109)1.4611641.5303.t=\frac{\sqrt5(10-9)}{1.461164} \approx1.5303.

在模型条件成立时,这个统计量于 μ=9\mu=9 下服从 t4t_4。分母使用 s/5s/\sqrt5,不能误用单个观测的标准差 ss;自由度是五个残差受到和为零约束后留下的四维。

两个独立残差平方和产生 F 比值

Q1χν12Q_1\sim\chi^2_{\nu_1}Q2χν22Q_2\sim\chi^2_{\nu_2} 且相互独立,则

F=Q1/ν1Q2/ν2Fν1,ν2.F=\frac{Q_1/\nu_1}{Q_2/\nu_2}\sim F_{\nu_1,\nu_2}.

设两组样本分别来自 N(μ1,σ12)N(\mu_1,\sigma_1^2)N(μ2,σ22)N(\mu_2,\sigma_2^2),样本量为 n1,n22n_1,n_2\ge2,两组彼此独立。各组卡方结论给出

S12/σ12S22/σ22Fn11,n21.\frac{S_1^2/\sigma_1^2}{S_2^2/\sigma_2^2} \sim F_{n_1-1,n_2-1}.

σ12=σ22\sigma_1^2=\sigma_2^2 的条件下,S12/S22S_1^2/S_2^2 直接服从该 F 分布。分子、分母互换时,自由度也必须互换:若 FFν1,ν2F\sim F_{\nu_1,\nu_2},则 1/FFν2,ν11/F\sim F_{\nu_2,\nu_1}。两组相关、总体非正态或方差估计方式改变时,这个有限样本结论不再自动成立。

因子分解检验样本压缩是否充分

充分统计量保留关于参数的全部样本信息,这句话必须相对于指定模型族解释。它不表示统计量能够还原每个观测,也不保证它适合估计所有模型外的量。

Neyman–Fisher 因子分解判据

设模型族由同一 σ\sigma-有限测度支配,联合密度或联合 PMF 为 fθ(x)f_\theta(x)。统计量 T(X)T(X)θ\theta 充分,当且仅当存在非负函数 gθg_\theta 与不依赖 θ\theta 的非负函数 hh,使几乎处处有

fθ(x)=gθ(T(x))h(x).f_\theta(x)=g_\theta(T(x))h(x).
证明

若存在分解,则在固定 T=tT=t 的纤维上,样本点之间的相对权重只由 h(x)h(x) 决定;正规化后得到的条件分布不含 θ\theta,故给定 TT 后原样本不再提供参数信息。反向上,若 TT 充分,则给定 T=tT=t 的条件分布可选为与 θ\theta 无关的版本。把联合分布写成 “TT 的边缘分布乘以给定 TT 的条件分布”,分别记为 gθ(t)g_\theta(t)h(x)h(x),便得到所述分解。支配性和几乎处处约定保证离散、连续情形使用同一表述。

例如 XiiidPoisson(λ)X_i\stackrel{\mathrm{iid}}\sim\operatorname{Poisson}(\lambda),则

fλ(x1,,xn)=enλλixii=1n1xi!.f_\lambda(x_1,\ldots,x_n) =e^{-n\lambda}\lambda^{\sum_i x_i} \prod_{i=1}^n\frac1{x_i!}.

参数只经 T=iXiT=\sum_iX_i 进入第一因子,所以总计数对 λ\lambda 充分。样本顺序对估计共同率不再提供额外信息,但若科学问题关心时间趋势,iid Poisson 模型本身就不再合适,充分性结论也随模型改变。

抽样设计决定独立性与标准误

“随机抽取”不等于 iid。有限总体中,从 NN 个固定数值无放回简单随机抽取 nn 个,样本位置之间负相关。若总体均值为 YN\overline Y_N,有限总体方差定义为

SN2=1N1j=1N(YjYN)2,S_N^2=\frac1{N-1}\sum_{j=1}^N(Y_j-\overline Y_N)^2,

则样本均值满足

E[X]=YN,Var(X)=(1nN)SN2n.E[\overline X]=\overline Y_N, \qquad \operatorname{Var}(\overline X) =\left(1-\frac nN\right)\frac{S_N^2}{n}.

1n/N1-n/N 是有限总体修正。若仍套用 σ2/n\sigma^2/n,就忽略了不放回造成的协方差。

分层抽样先把总体划为若干层,再在各层内抽样;总体均值估计量要按层规模加权,方差由各层抽样率共同决定。整群抽样一次抽取班级、医院或家庭,群内单位往往共享环境,正相关会使有效信息量低于同样数量的 iid 个体。若进入样本的概率还依赖未观测结果,选择机制具有信息性,仅分析已入样本者甚至会产生偏差。抽样分布必须绑定真实设计,不能只凭每个观测的边缘直方图来认定 iid。

三个容易越过的条件边界

一份样本的直方图就是均值的抽样分布

直方图的单位是单个观测;均值的抽样分布以整份样本为重复单位。要近似后者,需按同一设计重复抽取许多份容量为 nn 的样本并逐份计算均值,或从已声明模型推导。

样本量足够大时,所有 t 与卡方公式都精确

中心极限定理常使均值近似正态;对于非正态样本,残差平方和通常没有精确卡方分布, X\overline XS2S^2 也通常不独立。近似方法还要说明矩条件、尾部形状和误差容许度。

相同边缘分布仍会给出不同标准误

ZN(0,1)Z\sim N(0,1)。模型 A 取独立的 X1,X2N(0,1)X_1,X_2\sim N(0,1),则 Var((X1+X2)/2)=1/2\operatorname{Var}((X_1+X_2)/2)=1/2。模型 B 取 X1=X2=ZX_1=X_2=Z,每个观测的边缘仍为 N(0,1)N(0,1),但均值就是 ZZ,方差为 11。共享来源把标准误放大,边缘正态性不能替代联合独立性。

练习:从经验阶梯到精确枢轴

练习 1:经验分布的跳幅与样本比例

样本为 1,3,2,3,5,41,3,2,3,5,4。求 F6(3)F_6(3)F6(3)F_6(3^-) 以及样本落在 [2,4][2,4] 内的比例。

查看提示
先排序,再按不超过阈值的观测数除以六。
查看解答

不超过 33 的观测为 1,2,3,31,2,3,3,故 F6(3)=4/6=2/3F_6(3)=4/6=2/3;严格小于 33 的观测只有 1,21,2,所以 F6(3)=2/6=1/3F_6(3^-)=2/6=1/3。区间比例为

F6(4)F6(2)=5616=23,F_6(4)-F_6(2^-)=\frac56-\frac16=\frac23,

与直接计数 2,3,3,42,3,3,4 四项一致。

练习 2:均匀样本最小值的尾概率

X1,,X5iidUnif(0,1)X_1,\ldots,X_5\stackrel{\mathrm{iid}}\sim\operatorname{Unif}(0,1)L=X(1)L=X_{(1)}。求 P(L>0.2)P(L>0.2)FL()F_L(\ell)E[L]E[L]

查看提示
最小值大于 m 等价于每个观测都大于 m。
查看解答

独立性给出

P(L>)=(1)5(01),P(L>\ell)=(1-\ell)^5\quad(0\le\ell\le1),

P(L>0.2)=0.85=0.32768P(L>0.2)=0.8^5=0.32768,并且 FL()=1(1)5F_L(\ell)=1-(1-\ell)^5。密度为 5(1)45(1-\ell)^4,于是

E[L]=015(1)4d=16.E[L]=\int_0^1\ell\,5(1-\ell)^4\,\mathrm d\ell=\frac16.

期望落在支持 [0,1][0,1] 内,且最小值随五次抽样明显靠近零。

练习 3:样本均值的标准误

某总体均值为 1212、标准差为 66。从中 iid 抽取 3636 个观测。求 X\overline X 的期望、方差和标准误;说明是否已知它的精确分布形状。

查看提示
独立同分布下先把单个观测方差除以 n。
查看解答
E[X]=12,Var(X)=6236=1,se(X)=1.E[\overline X]=12, \qquad \operatorname{Var}(\overline X)=\frac{6^2}{36}=1, \qquad \operatorname{se}(\overline X)=1.

仅凭均值和方差不能确定有限样本均值的精确形状。若总体正态,则 XN(12,1)\overline X\sim N(12,1);若总体非正态,中心极限定理至多在附加条件下提供近似。

练习 4:同一正态样本中的卡方量与 t 量

n=10n=10 的样本来自 N(5,4)N(5,4),某次样本得到 x=5.6,s2=3.24\overline x=5.6,s^2=3.24。计算方差枢轴 qq 和均值 t 统计量,并写出各自的零假设分布。

查看提示
分别代入 Q=(n1)S2/σ2Q=(n-1)S^{2}/\sigma^{2}T=n(T=\sqrt{n}(均值μ)/S-\mu)/S
查看解答

总体方差 σ2=4\sigma^2=4,所以

q=9×3.244=7.29,t=10(5.65)1.81.0541.q=\frac{9\times3.24}{4}=7.29, \qquad t=\frac{\sqrt{10}(5.6-5)}{1.8}\approx1.0541.

在所给正态模型下,随机量 Qχ92Q\sim\chi^2_9;以未知方差版本检验 μ=5\mu=5 时,Tt9T\sim t_9。两个自由度都来自 n1=9n-1=9

练习 5:充分统计量与抽样设计的边界

(1)若 XiiidBernoulli(p)X_i\stackrel{\mathrm{iid}}\sim\operatorname{Bernoulli}(p),用因子分解说明 K=iXiK=\sum_iX_ipp 充分。(2)有限总体 N=100N=100,其中成功对象 4040 个;无放回抽取 n=20n=20。求成功比例 p^=K/20\widehat p=K/20 的方差,并与独立有放回抽样比较。

查看提示
先分解 Bernoulli 联合 PMF,再比较有放回和无放回的方差公式。
查看解答

Bernoulli 联合 PMF 为

pixi(1p)nixi=gp(K)1,p^{\sum_i x_i}(1-p)^{n-\sum_i x_i} =g_p(K)\cdot1,

KK 充分。无放回时超几何方差给出

Var(p^)=0.4(0.6)201002010010.00969697.\operatorname{Var}(\widehat p) =\frac{0.4(0.6)}{20}\frac{100-20}{100-1} \approx0.00969697.

有放回 iid 抽样的方差为 0.4(0.6)/20=0.0120.4(0.6)/20=0.012。有限总体修正 80/9980/99 使无放回方差更小;两种设计不能共用同一标准误。

概念之间的统计链条

抽样分布的核对资料

课程 · 2016

MIT 18.650 Statistics for Applications

Philippe Rigollet

适合把随机样本、估计量评价、置信区间、检验和回归放在同一课程结构中学习,并比较有限样本结论与渐近方法。

打开官方来源

MIT 18.650《Statistics for Applications》的公开讲义把统计模型、充分统计量、点估计与渐近分布放在同一推断链条中。阅读时应逐项核对本章使用的 iid 条件、正态有限样本结论与一般渐近结论,不把不同层次的定理合并。

书籍 · 2002

Statistical Inference, Second Edition

George Casella, Roger L. Berger

适合深入学习充分性、完备性、UMVU、似然、信息量、枢轴量和 Neyman–Pearson 框架的严格条件。

打开官方来源

Casella 与 Berger 的《Statistical Inference》第二版系统处理随机样本、次序统计量、充分性以及常见抽样分布。本章的最大值分布、因子分解判据和正态样本分解可沿教材相应章节复核定义、自由度与支配条件。

书籍 · 2023

Introductory Statistics 2e

Barbara Illowsky, Susan Dean

提供大量分步例题和练习,适合核对分布、样本统计量、置信区间、假设检验和简单线性回归的基础计算与解释。

打开官方来源

OpenStax《Introductory Statistics》提供样本均值、标准误、t 与卡方分布的分步数值案例,适合复算本章练习中的尺度和自由度。其入门叙述不替代这里对正态性、独立性和抽样设计的条件说明。

三项资料分别承担理论课程、严格教材与数值入门的角色。对任何外部公式,仍应检查随机单位、样本设计、参数化方式与自由度是否和当前问题一致。

下一章:估计量怎样比较

抽样分布说明统计量在重复抽样中的波动,却没有决定哪个统计量最适合估计未知参数。下一章以

点估计、充分性与信息量 为中心,比较矩估计与最大似然估计,按偏差、方差和均方误差评价方案,并在明确正则条件后讨论 Fisher 信息与 Cramér–Rao 下界。