第三章:完备充分性与最优无偏估计
本章研究如何在无偏估计量这一受限类别中寻找最优估计量。核心工具是完备性(completeness)、Rao-Blackwell 定理与 Lehmann-Scheffé 定理。最后介绍 Basu 定理,并通过 Bernoulli、Poisson、Normal 与 Uniform 模型说明这些结果的具体应用。
1. 为什么研究无偏估计量
设目标参数为 \(g(\theta)\),估计量为 \(\delta(X)\),损失函数为 \(L\{g(\theta),a\}\)。估计量的风险定义为
在平方损失下,风险可以分解为
其中
1.1 为什么通常不存在无约束的一致最优估计量
命题 1.1(无约束的一致最优估计量通常不存在)
假设概率测度族 \(\{P_\theta:\theta\in\Theta\}\) 两两相互绝对连续,并且 \(g(\theta)\) 不是常数。在平方损失下,不存在估计量 \(\delta_0\),使其风险在每个参数点都不大于所有其他估计量的风险。
命题 1.1 的证明(点击展开)
固定 \(\theta_0\),并将 \(\delta_0\) 与常数估计量
比较。后者在 \(\theta_0\) 处的风险为零。若 \(\delta_0\) 对所有估计量都一致最优,则
因此
由于各 \(P_\theta\) 相互绝对连续,上式在每个 \(P_\theta\) 下都几乎处处成立。再选择另一个参数值 \(\theta_1\),同理可得
这与 \(g\) 不是常数矛盾。\(\square\)
为了得到有意义的比较,可以限制估计量的类别,例如考虑不变估计、Bayes 估计、minimax 估计或无偏估计。本章研究最后一种选择。
1.2 无偏估计、MRU 与 UMVU
定义 1.2(无偏估计)
如果估计量 \(\delta(X)\) 满足
则称 \(\delta\) 是 \(g(\theta)\) 的无偏估计量。
定义 1.3(MRU 与 UMVU)
若无偏估计量 \(\delta_0\) 在每个参数点的风险都不大于任何其他无偏估计量的风险,则称 \(\delta_0\) 是最小风险无偏估计量(minimum-risk unbiased estimator, MRU)。
在平方损失下,MRU 称为一致最小方差无偏估计量(uniformly minimum-variance unbiased estimator, UMVU)。
对无偏估计量,偏差项为零,因此
这里的最优性是受限最优性:UMVU 估计量只在无偏估计量类中最优;扩大比较范围后,它仍可能被某个有偏估计量支配。
1.3 无偏估计量是否存在
无偏估计本质上是一个逆问题:寻找函数 \(\delta\),使得
并非每个目标函数 \(g\) 都属于这个期望算子的值域。
命题 1.4(Bernoulli 样本中的可估目标)
设
实值函数 \(g(p)\) 存在基于该样本的无偏估计量,当且仅当 \(g\) 是次数不超过 \(n\) 的关于 \(p\) 的多项式。
命题 1.4 的证明(点击展开)
记
对任意统计量 \(D(X)\),令
按照 \(T\) 的取值对样本点分组,可得
右端是次数不超过 \(n\) 的多项式,因此必要性成立。
反之,Bernstein 多项式
构成次数不超过 \(n\) 的多项式空间的一组基。将 \(g\) 在这组基下展开,并把相应系数定义为 \(d_t\),取 \(D=d_T\),即可得到 \(g(p)\) 的无偏估计量。\(\square\)
因此,当 \(r\leq n\) 时,\(p^r\) 可以被无偏估计;而 \(1/p\)、\(\log p\) 以及多数非线性标准化形状参数都不能由固定大小为 \(n\) 的 Bernoulli 样本无偏估计。
2. 完备性:由期望得到唯一性
2.1 完备与有界完备统计量
定义 2.1(完备统计量)
统计量 \(T\) 称为关于模型族 \(\mathcal{P}\) 完备,如果对任意满足可积性的函数 \(a(T)\),只要
就有
定义 2.2(有界完备统计量)
如果上述结论只要求对有界函数 \(a\) 成立,则称 \(T\) 有界完备(boundedly complete)。
完备性蕴含有界完备性。它与充分性描述的是不同性质:充分性表示 \(T\) 之外不再含有参数信息;完备性表示模型族下的期望能够唯一识别 \(T\) 的函数。
引理 2.3(无偏函数的唯一性)
若 \(T\) 完备,并且 \(a_1(T)\) 与 \(a_2(T)\) 都是同一目标函数 \(g(\theta)\) 的无偏估计量,则
引理 2.3 的证明(点击展开)
差值 \(a_1(T)-a_2(T)\) 可积,并且对所有 \(\theta\) 都满足
由完备性的定义,差值几乎处处为零。\(\square\)
2.2 Binomial 分布族的完备性
例 2.4(Binomial 完备性)
设
若对所有 \(p\) 都有 \(E_pa(T)=0\),则
两边除以 \((1-p)^n\),并令
得到
一个多项式若在区间上恒为零,则所有系数均为零。因此 \(a(t)=0\) 对 \(t=0,\ldots,n\) 成立,故 Binomial 分布族完备。
2.3 Poisson 分布族的完备性
例 2.5(Poisson 完备性)
设
若 \(E_\lambda a(T)=0\),则
可积性保证该幂级数对每个 \(\lambda>0\) 绝对收敛。因为它在一个区间上恒为零,所以所有系数都为零,即 \(a(t)=0\) 对所有 \(t\) 成立。因此 Poisson 分布族完备。
3. 满指数族的完备充分性
3.1 完备充分性定理
定理 3.1(满指数族中的完备充分统计量)
设模型具有指数族形式
并且自然参数空间 \(H\) 包含 \(\mathbb{R}^k\) 中的一个非空开集。则 \(T\) 对由该开集索引的模型族是充分且完备的。因此,对于任何包含该子族的更大模型族,\(T\) 仍然完备。
定理 3.1 的证明(点击展开)
由因子分解定理,\(T\) 是充分统计量。
令 \(\nu\) 为基准测度 \(h(x)\mu(dx)\) 在映射 \(x\mapsto T(x)\) 下的像测度,则 \(T\) 的诱导分布为
假设 \(a(T)\) 可积,并且其期望在一个开集上恒为零。选择该开集的内点 \(\eta_0\)。对原点附近的任意 \(u\),
定义有限符号测度
则 \(M\) 的双边 Laplace 变换满足
对原点某个邻域中的所有 \(u\) 成立。由下面的局部 Laplace 变换唯一性引理,\(M=0\)。因为 \(e^{\eta_0^Tt}>0\),可得
每个 \(P_\eta^T\) 关于 \(\nu\) 都有严格正的密度,因此 \(a(T)=0\) 在每个 \(P_\eta\) 下几乎处处成立。故 \(T\) 完备。\(\square\)
引理 3.2(局部 Laplace 变换的唯一性)
设 \(M\) 是 \(\mathbb{R}^k\) 上的有限符号测度,其指数矩在原点的某个邻域内有限。若
在该邻域内成立,则 \(M=0\)。
引理 3.2 的证明(点击展开)
对实部落在指数矩存在邻域内的复向量 \(z\),定义
用稍大的实指数矩控制被积函数,可以证明 \(L_M\) 是解析函数,并允许在积分号下求导。
从 \(L_M\) 在实开集上恒为零出发,逐个坐标应用单复变量的恒等定理,可以把零恒等式延拓到包含虚轴的连通复管状区域。因此,对每个 \(s\in\mathbb{R}^k\),
由有限符号测度的 Fourier 变换唯一性,得到 \(M=0\)。\(\square\)
3.2 为什么需要“开集”条件
自然参数空间具有满仿射张成足以保证最小充分性,因为有限个似然比就可以恢复自然统计量。但是,完备性要求 Laplace 变换在一个邻域上相等。弯曲参数空间可能具有满仿射张成,却不包含任何开集。因此,最小充分性本身并不能推出完备性。
3.3 完整正态分布族
例 3.3(完整正态分布族)
设
自然参数为
自然参数空间为
它是开集。因此
是 \((\mu,\sigma^2)\) 的完备充分统计量。
当 \(n\geq2\) 时,它与
一一对应。因此 \((\overline X,S^2)\) 也是完备充分统计量。
一个常见误解
当 \(\mu\) 与 \(\sigma^2\) 都未知时,\((\overline X,S^2)\) 对完整正态分布族仍然完备。可能失去完备性的是自然参数空间不含开集的弯曲正态子模型。“多参数”并不意味着“不完备”。
4. Rao-Blackwell 改进
4.1 凸损失下的 Rao-Blackwell 定理
定理 4.1(凸损失下的 Rao-Blackwell 定理)
设 \(T\) 是充分统计量,\(\delta(X)\) 是可积估计量,并定义
充分性保证右端可以选为一个不依赖于未知参数 \(\theta\) 的 \(T\) 的函数。若 \(a\mapsto L\{g(\theta),a\}\) 是凸函数,则
如果 \(\delta\) 是 \(g(\theta)\) 的无偏估计量,则 \(\delta_T\) 仍然无偏。
定理 4.1 的证明(点击展开)
充分性提供了一个与 \(\theta\) 无关的条件分布核,因此条件期望可以视为同一个 \(T\) 的可测函数。由条件 Jensen 不等式,
两边取期望便得到风险不等式。
另外,由塔式法则,
因此无偏性得以保持。\(\square\)
4.2 平方损失下的精确改进量
推论 4.2(平方损失下的精确风险改进)
在平方损失下,
除非 \(\delta\) 已经是 \(T\) 的函数(忽略 \(P_\theta\)-零测集),否则改进是严格的。
推论 4.2 的证明(点击展开)
写成
第一项关于 \(T\) 的条件期望为零,因此展开平方后交叉项的期望为零。于是
其中
整理即可。\(\square\)
这个等式正是条件方差分解,也解释了 Rao-Blackwell 化与条件 Monte Carlo、现代方差缩减方法之间的联系:它平均掉了不影响推断目标的额外随机性。
4.3 Bernoulli 幂函数的 Rao-Blackwell 化
例 4.3(Bernoulli 幂函数)
设
当 \(r\leq n\) 时,\(X_1\cdots X_r\) 是 \(p^r\) 的无偏估计量。给定 \(T=t\) 后,含有 \(t\) 个 \(1\) 的所有排列等可能,因此
其中下降阶乘定义为
因而 \((T)_r/(n)_r\) 仍然无偏,并且方差不大于原始乘积估计量。
5. Lehmann-Scheffé 定理
定理 5.1(凸损失下的 Lehmann-Scheffé 定理)
设 \(T\) 是完备充分统计量。若 \(\delta_0(T)\) 是 \(g(\theta)\) 的无偏估计量,并且风险有限,则 \(\delta_0\) 在任何关于行动 \(a\) 为凸函数的损失下都是 MRU 估计量。
在平方损失下,\(\delta_0\) 是唯一的 UMVU 估计量。更一般地,在严格凸损失下,\(\delta_0\) 是唯一的 MRU 估计量,其中唯一性忽略各模型共同的零测集。
定理 5.1 的证明(点击展开)
任取 \(g(\theta)\) 的无偏估计量 \(D(X)\)。对其进行 Rao-Blackwell 化:
则 \(D_T\) 仍然无偏,并且风险不大于 \(D\)。由于 \(D_T(T)\) 与 \(\delta_0(T)\) 都是完备统计量 \(T\) 的无偏函数,引理 2.3 给出
因此
从而证明 MRU 最优性。
在平方损失下,推论 4.2 表明等号成立只能发生在 \(D=D_T=\delta_0\) 几乎处处时。严格凸损失下,由 Jensen 不等式的等号条件可得同样结论。\(\square\)
UMVU 的标准构造步骤
(i) 找到一个完备充分统计量 \(T\);
(ii) 找到目标函数的任意无偏估计量 \(D(X)\);
(iii) 计算 \(E(D\mid T)\),或者直接找到一个关于 \(T\) 的无偏函数;
(iv) 使用 Lehmann-Scheffé 定理得到唯一性与 UMVU 最优性。
6. 四类代表性 UMVU 构造
6.1 Bernoulli 幂函数
统计量
由因子分解定理知是充分统计量,由例 2.4 知是完备统计量。因此
6.2 Poisson 点概率
设
对固定的 \(k\geq0\),示性函数 \(\mathbf{1}\{X_1=k\}\) 无偏估计
给定 \(T=t\) 后,\((X_1,\ldots,X_n)\) 服从总数为 \(t\)、各格概率为 \(1/n\) 的 multinomial 分布。因此
由 Poisson 完备性和 Lehmann-Scheffé 定理,\(\widehat g_k(T)\) 是 \(g_k(\lambda)\) 的 UMVU 估计量。
6.3 正态总体的均值与方差
设
由例 3.3,\((\overline X,S^2)\) 是完整正态分布族的完备充分统计量。又因为
所以 \(\overline X\) 与 \(S^2\) 分别是 \(\mu\) 与 \(\sigma^2\) 的 UMVU 估计量。
此外,
因此
是 \(\mu^2\) 的无偏估计量。它是完备充分统计量的函数,所以也是 \(\mu^2\) 的 UMVU 估计量。
6.4 Uniform 单端点模型
设
由因子分解定理,\(M\) 是充分统计量。若对每个 \(\theta>0\) 都有 \(E_\theta a(M)=0\),则
因此,对每个 \(\theta\),上式中的不定积分均为零。由 Lebesgue 意义下的微分,
几乎处处成立,所以 \(M\) 完备。
另一方面,当 \(r>-n\) 时,
因此,只要目标函数有定义,
这个非正则例子表明,完备充分性并不局限于具有固定支撑集的指数族。
7. Basu 定理:完备性与辅助统计量
7.1 辅助统计量
定义 7.1(辅助统计量)
如果统计量 \(A(X)\) 的分布不依赖于参数,则称 \(A\) 是辅助统计量(ancillary statistic)。
辅助统计量仍然可以反映估计精度、设计不平衡或抽样几何结构。“分布与参数无关”并不等于“没有用”。
7.2 Basu 定理
定理 7.2(Basu 定理)
若 \(T\) 是充分且有界完备的统计量,\(A\) 是辅助统计量,则在每个 \(P_\theta\) 下,\(T\) 与 \(A\) 相互独立。
定理 7.2 的证明(点击展开)
固定 \(A\) 的值域中的可测集 \(B\)。由充分性,可以选择一个与 \(\theta\) 无关的版本
由辅助性,
也与 \(\theta\) 无关。因此
函数 \(q_B-c_B\) 有界。由有界完备性,
对 \(T\) 的值域中的任意可测集 \(C\),
因此 \(T\) 与 \(A\) 独立。\(\square\)
7.3 为什么 \(\overline X\) 与 \(S^2\) 独立
例 7.3(正态样本均值与样本方差的独立性)
固定 \(\sigma^2\),考虑只由 \(\mu\) 索引的子模型
此时 \(\overline X\) 是关于 \(\mu\) 的完备充分统计量,而
的分布与 \(\mu\) 无关,所以 \(S^2\) 是该子模型中的辅助统计量。由 Basu 定理,
由于该结论对每个固定的 \(\sigma^2\) 都成立,所以它在完整正态分布族中成立。
8. 函数型目标:估计正态分布函数
设
其中 \(\sigma^2\) 已知、\(n\geq2\),并固定 \(y\in\mathbb{R}\)。示性函数
无偏估计
此时完备充分统计量为 \(\overline X\)。由联合正态分布的条件分布公式,
因此,Rao-Blackwell 化给出
类似的 Gaussian 卷积恒等式给出总体密度在 \(y\) 处的 UMVU 估计量:
对上式取期望,相当于把 \(N(\mu,\sigma^2/n)\) 的密度与 \(N(0,\sigma^2(1-1/n))\) 的密度卷积,结果正是 \(N(\mu,\sigma^2)\) 在 \(y\) 处的密度。完备性进一步保证该估计量是唯一的 UMVU 估计量。
9. Uniform 双端点模型
设
并记
9.1 极端次序统计量的完备充分性
定理 9.1(极端次序统计量的完备充分性)
统计量 \((U,V)\) 是参数 \((a,b)\) 的完备充分统计量。
定理 9.1 的证明(点击展开)
样本联合密度为
因此由因子分解定理,\((U,V)\) 充分。
\((U,V)\) 的联合密度为
假设对所有 \(a<b\) 都有
两边乘以 \((b-a)^n\),得到
令
绝对可积性允许在 Lebesgue 意义下求导。对几乎处处的 \(a<b\),
并且
因为 \(G\equiv0\),所以 \(q(a,b)=0\) 几乎处处。权重 \(n(n-1)(b-a)^{n-2}\) 在 \(a<b\) 时为正,因此 \(c(a,b)=0\) 几乎处处。再由 \((U,V)\) 的联合密度可知
对每个参数对 \((a,b)\) 都成立。因此 \((U,V)\) 完备。\(\square\)
9.2 两个端点的 UMVU 估计
对标准 Uniform\((0,1)\) 样本,
经线性变换可得
解这两个线性方程,得到
由定理 9.1 与 Lehmann-Scheffé 定理,它们分别是 \(a\) 与 \(b\) 的 UMVU 估计量。
10. 凸损失下不需要随机化
设随机化估计量在观察到 \(X=x\) 后,按照分布 \(Q_x(da)\) 选择行动。定义该随机行动的重心
若随机化估计量无偏,则
所以 \(\overline\delta\) 是普通的非随机化无偏估计量。
由凸性和 Jensen 不等式,对每个 \(x\),
再对 \(X\) 积分,得到非随机化重心估计量的风险不大于原随机化估计量。因此,在凸损失下,随机化不会给无偏估计带来额外收益。
这与离散假设检验不同:在离散检验中,行动空间与约束结构不同,随机化仍可能十分重要。
11. 本章证明核查清单
在断言某个估计量为 UMVU,或者使用独立性定理之前,应检查以下问题:
- 目标函数在整个参数空间上是否存在无偏估计量?
- 用来条件化的统计量是否充分,从而使 \(E_\theta(\delta\mid T)\) 能够选成同一个与参数无关的统计量?
- 该统计量是完备的,还是仅仅最小充分?
- 使用指数族完备性时,自然参数空间是否包含开集,而不仅仅具有满仿射张成?
- 期望、幂级数或 Laplace 变换所需的可积性条件是否成立?
- Rao-Blackwell 定理所用的损失函数是否凸?只有排除等号情形后,才能声称改进是严格的。
- 在引用 Lehmann-Scheffé 定理之前,候选估计量是否既无偏,又是完备充分统计量的函数?
- 使用 Basu 定理时,第二个统计量是否确实关于同一个参数族辅助,而 \(T\) 是否关于该参数族有界完备?
- 最终结论只是无偏估计量类中的最优性,还是错误地声称了在所有估计量中的可容许性?
12. 本章小结
本章的逻辑可以概括为:
完备充分统计量既提供了降维结构,又保证无偏函数的唯一性。Rao-Blackwell 定理把任意无偏估计量改进为充分统计量的函数,Lehmann-Scheffé 定理进一步把该函数识别为唯一的 UMVU 估计量。Basu 定理则揭示了完备充分统计量与辅助统计量之间的独立关系。