统计模型把重复抽样写成概率对象
统计推断从观测数字追溯可能产生它们的概率模型族。记统计模型为
{Pθ:θ∈Θ},其中参数 θ 在频率学解释下固定但未知。一次实验得到
x1,…,xn;为了讨论同一实验重复进行时结果怎样波动,先把这些数还原为随机变量
X1,…,Xn。小写字母表示已经观察到的实现值,大写字母保留尚未实现的随机性。
随机样本、统计量与抽样分布
若 X1,…,Xn 相对于同一 Pθ 独立同分布,则称它们为来自
Pθ 的容量为 n 的随机样本。统计量是形如
T=T(X1,…,Xn) 的可测函数,它只能含样本和已知常数,不能含未知参数。把观测值代入后得到统计量的实现值
t=T(x1,…,xn)。当同一抽样机制反复执行时,T 在 Pθ 下形成的概率分布称为
T 的抽样分布;该分布通常随 θ 与样本量 n 改变。
样本均值
X=n1i=1∑nXi
是统计量,因为公式中没有未知量。若总体均值 μ 未知,
∑i(Xi−μ)2 不是统计量。无偏样本方差采用约定
S2=n−11i=1∑n(Xi−X)2.
一次样本的直方图描述样本内部的数值分布;X 的抽样分布描述重新抽取整组样本后均值怎样变化。二者的横轴、随机单位和重复机制都不同。
经验分布把样本变成随机阶梯
经验分布函数
给定随机样本 X1,…,Xn,经验分布函数为
Fn(t)=n1i=1∑n1{Xi≤t}. 对固定样本,Fn 是右连续的非降阶梯函数;每个观测贡献 1/n 的跳跃,重复值使同一点的跳幅相加。
Fn(t) 同时是一项统计量和一条随机函数。固定 t 时,指标变量
1{Xi≤t} 是成功概率 F(t) 的 Bernoulli 变量,因此
Eθ[Fn(t)]=F(t),Varθ(Fn(t))=nF(t)[1−F(t)].
这个逐点结论说明经验分布在每个阈值处无偏,并给出 n−1/2 量级的波动。它没有自动给出整条曲线的一致误差界;研究
supt∣Fn(t)−F(t)∣ 需要更强的经验过程结论。
五次读数形成的经验阶梯
观测值为 2.1,1.4,2.1,3.0,0.9。排序后为
0.9,1.4,2.1,2.1,3.0,所以
Fn(t)=⎩⎨⎧0,1/5,2/5,4/5,1,t<0.9,0.9≤t<1.4,1.4≤t<2.1,2.1≤t<3.0,t≥3.0. 在 2.1 处有两个观测,跳幅为 2/5。例如样本中落在区间
(1.4,3.0] 的比例是
Fn(3.0)−Fn(1.4)=1−52=53, 对应 2.1,2.1,3.0 三项。直接计数与阶梯函数增量一致,也核对了左端开、右端闭的约定。
次序统计量记录位置而非原编号
将样本从小到大排列,记为
X(1)≤⋯≤X(n)。X(1) 与 X(n) 分别是样本最小值和最大值;
X(k) 是第 k 个次序统计量。排序会丢掉原来的观测编号,却保留分位数、极差和端点信息。
若公共分布连续且具有密度 f,第 k 个次序统计量的密度为
fX(k)(x)=(k−1)!(n−k)!n!F(x)k−1[1−F(x)]n−kf(x).
组合因子选择哪 k−1 个观测落在 x 左侧、哪一个落入 x 附近以及其余观测落在右侧。连续性排除了同点并列造成的额外分配;离散样本有正概率出现并列,应改用事件计数或 CDF 差分,不能直接沿用这个密度。
样本最大值的分布
设 X1,…,Xn iid,公共分布函数为 F,并令
M=X(n)。则对每个实数 m,
Pθ(M≤m)=F(m)n. 若 F 在区间内具有密度 f,则 M 在相应区间内的密度为
fM(m)=nF(m)n−1f(m).
证明
事件 {M≤m} 当且仅当每个观测都不超过 m。iid 假设给出
Pθ(M≤m)=Pθ(X1≤m,…,Xn≤m)=i=1∏nPθ(Xi≤m)=F(m)n. 在 F 可微处对 m 求导,链式法则给出密度公式。若样本相互依赖,乘积步骤失效;边缘分布都等于
F 仍不足以推出 FM=Fn。
样本均值的中心、尺度和形状
若 Xi iid,且 E[Xi]=μ、Var(Xi)=σ2<∞,则独立性使协方差交叉项消失:
E[X]=μ,Var(X)=n21i=1∑nVar(Xi)=nσ2.
σ 描述单个观测的标准差,σ/n 是样本均值的标准误。未知
σ 时常以 S/n 估计标准误,但这个随机估计值不等于真实抽样标准差。
均值的精确形状取决于总体分布。若总体正态,任意有限 n 都有
X∼N(μ,nσ2).
若总体非正态但满足经典 iid 中心极限定理的有限正方差条件,则
σn(X−μ)dN(0,1).
这是 n→∞ 的渐近结论,不会把有限样本均值精确变成正态。偏斜、重尾或强依赖会改变近似速度,甚至破坏所用定理的条件。
正态样本把平方和分解为卡方变量
正态样本均值与样本方差的精确分布
若 X1,…,Xn∼iidN(μ,σ2),其中
n≥2、σ>0,则
Z=σn(X−μ)∼N(0,1),Q=σ2(n−1)S2∼χn−12, 且 Z⊥Q,等价地 X⊥S2。因此
T=Sn(X−μ)=Q/(n−1)Z∼tn−1.
证明
将标准化样本向量写成
Y=((X1−μ)/σ,…,(Xn−μ)/σ)T,则
Y∼Nn(0,In)。把 Y 正交投影到由
(1,…,1)T/n 张成的一维子空间,投影坐标正是 Z;投影到其正交补得到残差向量。球对称标准正态在正交变换后仍有相互独立的标准正态坐标,所以一维投影与其余
n−1 个坐标独立,而残差平方长度等于
i=1∑nσ2(Xi−X)2=σ2(n−1)S2. 它是 n−1 个独立标准正态平方之和,故服从 χn−12。独立比值的定义随即给出
tn−1 分布。
正态条件承担了关键作用。一般非正态总体下,X 与 S2 通常相关或至少不独立,
Q 也未必服从卡方分布。大样本 t 近似属于另一层渐近论证,不能改写成上述有限样本精确结论。
五个正态读数形成的 t 统计量
假设 Xi∼iidN(μ,σ2),观测为
8.2,9.1,10.0,10.7,12.0。计算得
x=10.0,i=1∑5(xi−x)2=3.24+0.81+0+0.49+4.00=8.54, 所以 s2=8.54/4=2.135、s≈1.461164。若考察中心值
μ0=9,则
t=1.4611645(10−9)≈1.5303. 在模型条件成立时,这个统计量于 μ=9 下服从 t4。分母使用
s/5,不能误用单个观测的标准差 s;自由度是五个残差受到和为零约束后留下的四维。
两个独立残差平方和产生 F 比值
若 Q1∼χν12、Q2∼χν22 且相互独立,则
F=Q2/ν2Q1/ν1∼Fν1,ν2.
设两组样本分别来自 N(μ1,σ12) 与
N(μ2,σ22),样本量为 n1,n2≥2,两组彼此独立。各组卡方结论给出
S22/σ22S12/σ12∼Fn1−1,n2−1.
在 σ12=σ22 的条件下,S12/S22 直接服从该 F 分布。分子、分母互换时,自由度也必须互换:若
F∼Fν1,ν2,则 1/F∼Fν2,ν1。两组相关、总体非正态或方差估计方式改变时,这个有限样本结论不再自动成立。
因子分解检验样本压缩是否充分
充分统计量保留关于参数的全部样本信息,这句话必须相对于指定模型族解释。它不表示统计量能够还原每个观测,也不保证它适合估计所有模型外的量。
Neyman–Fisher 因子分解判据
设模型族由同一 σ-有限测度支配,联合密度或联合 PMF 为
fθ(x)。统计量 T(X) 对 θ 充分,当且仅当存在非负函数
gθ 与不依赖 θ 的非负函数 h,使几乎处处有
fθ(x)=gθ(T(x))h(x).
证明
若存在分解,则在固定 T=t 的纤维上,样本点之间的相对权重只由
h(x) 决定;正规化后得到的条件分布不含 θ,故给定 T 后原样本不再提供参数信息。反向上,若 T 充分,则给定
T=t 的条件分布可选为与 θ 无关的版本。把联合分布写成
“T 的边缘分布乘以给定 T 的条件分布”,分别记为
gθ(t) 与 h(x),便得到所述分解。支配性和几乎处处约定保证离散、连续情形使用同一表述。
例如 Xi∼iidPoisson(λ),则
fλ(x1,…,xn)=e−nλλ∑ixii=1∏nxi!1.
参数只经 T=∑iXi 进入第一因子,所以总计数对 λ 充分。样本顺序对估计共同率不再提供额外信息,但若科学问题关心时间趋势,iid Poisson 模型本身就不再合适,充分性结论也随模型改变。
抽样设计决定独立性与标准误
“随机抽取”不等于 iid。有限总体中,从 N 个固定数值无放回简单随机抽取 n 个,样本位置之间负相关。若总体均值为
YN,有限总体方差定义为
SN2=N−11j=1∑N(Yj−YN)2,
则样本均值满足
E[X]=YN,Var(X)=(1−Nn)nSN2.
1−n/N 是有限总体修正。若仍套用 σ2/n,就忽略了不放回造成的协方差。
分层抽样先把总体划为若干层,再在各层内抽样;总体均值估计量要按层规模加权,方差由各层抽样率共同决定。整群抽样一次抽取班级、医院或家庭,群内单位往往共享环境,正相关会使有效信息量低于同样数量的 iid 个体。若进入样本的概率还依赖未观测结果,选择机制具有信息性,仅分析已入样本者甚至会产生偏差。抽样分布必须绑定真实设计,不能只凭每个观测的边缘直方图来认定 iid。
三个容易越过的条件边界
一份样本的直方图就是均值的抽样分布
直方图的单位是单个观测;均值的抽样分布以整份样本为重复单位。要近似后者,需按同一设计重复抽取许多份容量为
n 的样本并逐份计算均值,或从已声明模型推导。
样本量足够大时,所有 t 与卡方公式都精确
中心极限定理常使均值近似正态;对于非正态样本,残差平方和通常没有精确卡方分布,
X 与 S2 也通常不独立。近似方法还要说明矩条件、尾部形状和误差容许度。
相同边缘分布仍会给出不同标准误
令 Z∼N(0,1)。模型 A 取独立的 X1,X2∼N(0,1),则
Var((X1+X2)/2)=1/2。模型 B 取 X1=X2=Z,每个观测的边缘仍为
N(0,1),但均值就是 Z,方差为 1。共享来源把标准误放大,边缘正态性不能替代联合独立性。
练习:从经验阶梯到精确枢轴
练习 1:经验分布的跳幅与样本比例
- 所属知识
- 经验分布
- 难度
- 2/5
样本为 1,3,2,3,5,4。求 F6(3)、F6(3−) 以及样本落在
[2,4] 内的比例。
查看提示
先排序,再按不超过阈值的观测数除以六。
查看解答
不超过 3 的观测为 1,2,3,3,故 F6(3)=4/6=2/3;严格小于
3 的观测只有 1,2,所以 F6(3−)=2/6=1/3。区间比例为
F6(4)−F6(2−)=65−61=32, 与直接计数 2,3,3,4 四项一致。
查看解答
独立性给出
P(L>ℓ)=(1−ℓ)5(0≤ℓ≤1), 故 P(L>0.2)=0.85=0.32768,并且
FL(ℓ)=1−(1−ℓ)5。密度为
5(1−ℓ)4,于是
E[L]=∫01ℓ5(1−ℓ)4dℓ=61. 期望落在支持 [0,1] 内,且最小值随五次抽样明显靠近零。
练习 3:样本均值的标准误
- 所属知识
- 样本均值
- 难度
- 2/5
某总体均值为 12、标准差为 6。从中 iid 抽取 36 个观测。求
X 的期望、方差和标准误;说明是否已知它的精确分布形状。
查看提示
独立同分布下先把单个观测方差除以 n。
查看解答
E[X]=12,Var(X)=3662=1,se(X)=1. 仅凭均值和方差不能确定有限样本均值的精确形状。若总体正态,则
X∼N(12,1);若总体非正态,中心极限定理至多在附加条件下提供近似。
练习 4:同一正态样本中的卡方量与 t 量
- 所属知识
- 正态样本分解
- 难度
- 3/5
设 n=10 的样本来自 N(5,4),某次样本得到
x=5.6,s2=3.24。计算方差枢轴 q 和均值 t 统计量,并写出各自的零假设分布。
查看提示
分别代入
Q=(n−1)S2/σ2 与
T=n(均值
−μ)/S。
查看解答
总体方差 σ2=4,所以
q=49×3.24=7.29,t=1.810(5.6−5)≈1.0541. 在所给正态模型下,随机量 Q∼χ92;以未知方差版本检验
μ=5 时,T∼t9。两个自由度都来自 n−1=9。
练习 5:充分统计量与抽样设计的边界
- 所属知识
- 因子分解与依赖
- 难度
- 4/5
(1)若 Xi∼iidBernoulli(p),用因子分解说明
K=∑iXi 对 p 充分。(2)有限总体 N=100,其中成功对象 40 个;无放回抽取
n=20。求成功比例 p=K/20 的方差,并与独立有放回抽样比较。
查看提示
先分解 Bernoulli 联合 PMF,再比较有放回和无放回的方差公式。
查看解答
Bernoulli 联合 PMF 为
p∑ixi(1−p)n−∑ixi=gp(K)⋅1, 故 K 充分。无放回时超几何方差给出
Var(p)=200.4(0.6)100−1100−20≈0.00969697. 有放回 iid 抽样的方差为 0.4(0.6)/20=0.012。有限总体修正
80/99 使无放回方差更小;两种设计不能共用同一标准误。
概念之间的统计链条
抽样分布的核对资料
课程 · 2016MIT 18.650 Statistics for Applications
Philippe Rigollet
适合把随机样本、估计量评价、置信区间、检验和回归放在同一课程结构中学习,并比较有限样本结论与渐近方法。
打开官方来源
MIT 18.650《Statistics for Applications》的公开讲义把统计模型、充分统计量、点估计与渐近分布放在同一推断链条中。阅读时应逐项核对本章使用的 iid 条件、正态有限样本结论与一般渐近结论,不把不同层次的定理合并。
书籍 · 2002Statistical Inference, Second Edition
George Casella, Roger L. Berger
适合深入学习充分性、完备性、UMVU、似然、信息量、枢轴量和 Neyman–Pearson 框架的严格条件。
打开官方来源
Casella 与 Berger 的《Statistical Inference》第二版系统处理随机样本、次序统计量、充分性以及常见抽样分布。本章的最大值分布、因子分解判据和正态样本分解可沿教材相应章节复核定义、自由度与支配条件。
书籍 · 2023Introductory Statistics 2e
Barbara Illowsky, Susan Dean
提供大量分步例题和练习,适合核对分布、样本统计量、置信区间、假设检验和简单线性回归的基础计算与解释。
打开官方来源
OpenStax《Introductory Statistics》提供样本均值、标准误、t 与卡方分布的分步数值案例,适合复算本章练习中的尺度和自由度。其入门叙述不替代这里对正态性、独立性和抽样设计的条件说明。
三项资料分别承担理论课程、严格教材与数值入门的角色。对任何外部公式,仍应检查随机单位、样本设计、参数化方式与自由度是否和当前问题一致。
下一章:估计量怎样比较
抽样分布说明统计量在重复抽样中的波动,却没有决定哪个统计量最适合估计未知参数。下一章以
点估计、充分性与信息量
为中心,比较矩估计与最大似然估计,按偏差、方差和均方误差评价方案,并在明确正则条件后讨论 Fisher
信息与 Cramér–Rao 下界。