Skip to content

Cramér–Rao 下界

Cramér–Rao lower bound · Information inequality

正则模型中估计量方差由目标导数、偏差导数与 Fisher 信息共同限制的信息下界。

形式陈述

设标量参数 θ 的支配模型满足共同支持、密度可微、微分可移入积分等正则条件,单样本 Fisher 信息 I(θ) 有限且正。若估计量 T 二阶可积,其均值为

m(θ)=Eθ[T],

m 可微,则 n 个独立样本下

Varθ(T)m(θ)2nI(θ).

T 无偏估计 g(θ),则 m=g,得到

Varθ(T)g(θ)2nI(θ).

若偏差 b(θ)=m(θ)g(θ),则带偏版本为

Varθ(T)(g(θ)+b(θ))2nI(θ).

结合MSE 分解还需加 b(θ)2;下界不会自动惩罚偏差大小。

推导

令样本 score 为 Sn(θ)=θn(θ)。正则性给

Eθ[Sn]=0,Varθ(Sn)=nI(θ).

m(θ)=T(x)pθ(x)dμ(x) 求导并移入积分:

m(θ)=Eθ[TSn]=Covθ(T,Sn).

Cauchy–Schwarz 不等式

m(θ)2Varθ(T)Varθ(Sn),

代入 score 方差即得结论。等号当且仅当中心化估计误差与 score 几乎处处线性相关:

Tm(θ)=a(θ)Sn(θ).

这一条件必须由同一个不含未知参数的统计量在相关参数范围内满足,故信息下界未必可达。

例子与计算

XiN(μ,σ2)σ 已知,单样本信息 I(μ)=1/σ2。对无偏估计 μ

Varμ(T)σ2n.

样本均值方差恰为 σ2/n,并且

X¯μ=σ2nSn(μ),

所以达到等号。

Bernoulli(p)I(p)=1/[p(1p)],无偏估计 p 的下界为 p(1p)/n,样本比例同样达到。估计 g(p)=p2 时下界为 4p3(1p)/n,但达到它需要额外检查,不由代入公式保证。

边界与失败情形

Uniform(0,θ) 支持依参数变化,score 零均值推导失败。样本最大值的误差尺度为 1/n,看似“突破”常规 1/n 界,实际是正则假设不成立。

Cramér–Rao 是点态方差下界,不是任意估计量的 minimax 风险下界。允许偏差后,可取常数估计量在某个参数点方差为零;完整 MSE 比较必须计入偏差并考虑参数范围。

信息矩阵奇异、参数在边界或 nuisance 参数存在时,需使用矩阵 Schur complement、约束切空间或其他下界。直接取奇异矩阵逆没有意义。

界可达性不是定理结论。若等号条件要求估计量显含未知参数,就不存在全参数有效的 efficient unbiased estimator。

带 nuisance 参数时,目标坐标的有效信息不是信息矩阵相应对角元素,而是消去 nuisance 后的 Schur complement。忽略交叉信息会给出过度乐观的界。若参数受约束,还需把 score 投影到允许切空间,普通矩阵逆公式不再完整。

有限样本下界也可严格却不可达。达到等号要求中心化估计误差与 score 在几乎处处意义下成比例,而且同一估计规则不能依赖未知参数;仅算出一个漂亮的界值不表示存在估计量实现它。

推论与应用

向量参数版本为协方差矩阵的半正定序下界

Covθ(T)G(θ)In(θ)1G(θ)T,

其中 G 是目标 Jacobian。比较标量方差时需指定方向或损失权重。

正则 MLE 的渐近协方差常达到信息逆,形成渐近效率;这与有限样本无偏并达到 Cramér–Rao 界是不同命题。

参考资料
  • Harald Cramér, Mathematical Methods of Statistics, Princeton University Press, 1946,Ch. 32。
  • C. R. Rao, “Information and the Accuracy Attainable in the Estimation of Statistical Parameters,” 1945。
  • Erich L. Lehmann and George Casella, Theory of Point Estimation, 2nd ed., Springer, 1998,Ch. 2。