第 17 章 · 点估计与最优性

17.4 Fisher信息与Cramér–Rao界

Fisher Information and Information Bounds
Fisher信息把模型随参数的局部变化与估计量的精度联系起来。本节明确共同支持和求导积分换序条件,推导得分均值、信息恒等式及独立样本的可加性,再证明无偏、有偏和多参数Cramér–Rao界与达界条件。伯努利、正态和指数算例展示精确比较,Uniform端点反例说明忽略支持变化为何会得到错误下界。
先修与阅读位置

主要先修:9.3 勒贝格积分与积分收敛定理、16.1 统计模型、似然与可识别性、17.1 偏差、风险与一致性。

本节以求导换序和Cauchy–Schwarz建立信息界,接19.3的正则似然理论。17.3的无偏最优估计提供另一种比较,不是本节证明的先修。 查看全书依赖图与分支路线。

17.4.1 得分函数之前,先检查正则性

在17.3节,无偏最优估计通过条件期望和充分性构造。本节从模型本身出发,为所有满足条件的估计量给出方差下界。Cramér–Rao界是有限样本不等式,不需要先把抽样分布近似为正态;它的关键在于参数变化能否通过密度导数表达。

定义17.4-1 得分函数与Fisher信息

先设 \(\Theta\) 为实数上的开区间,\(P_\theta\) 相对于不依赖参数的 \(\sigma\) 有限测度 \(\nu\) 有密度 \(p_\theta\),其正密度区域 \(A\) 不随参数改变。对完整数据 \(X\) 定义

\[\ell_\theta(x)=\log p_\theta(x),\qquad U_\theta(X)=\partial_\theta\ell_\theta(X),\qquad I(\theta)=\mathbb E_\theta\{U_\theta(X)^2\}.\]
(17.4-1)

\(U_\theta\) 称为得分函数(score)。本节 \(\log\) 使用自然对数。这里的 \(I(\theta)\) 是完整数据的信息,讨论单个观测时另记 \(I_1(\theta)\),不能遗漏样本量因子。

以下采用一组容易核查的充分条件。对每个内点 \(\theta_0\),在其某邻域内,\(p_\eta(x)\) 对 \(\eta\) 在 \(\nu\) 几乎处处二次连续可微;一阶、二阶导数分别被 \(\nu\) 可积的函数 \(H_1(x),H_2(x)\) 一致控制,即 \(|\partial_\eta p_\eta(x)|\le H_1(x)\)、\(|\partial_\eta^2 p_\eta(x)|\le H_2(x)\)。同时要求 \(I(\theta_0)<\infty\)。控制收敛定理允许相应的求导与积分换序;只需要一阶结论时,不必额外要求二阶条件。

对具体估计量 \(T(X)\),还需有限方差,以及估计量期望可微所需的换序条件。例如,邻域内 \(|T(x)\partial_\eta p_\eta(x)|\le H_T(x)\),其中 \(H_T\) 对 \(\nu\) 可积,便是一个充分条件。密度本身可换序,不自动保证乘上任意 \(T\) 后也可换序。更一般的正则性可用密度比在平方可积空间中的可微性表述,见van de Geer讲义§3.3。

得分在给定参考参数处可以从数据算出,但代入未知真参数后的 \(U_\theta(X)\) 不是不依赖未知参数的统计量。下文所有期望都使用与得分求导位置相同的 \(\theta\);在其他参考位置取值,均值一般不为零。

17.4.2 零均值得分、曲率与信息可加性

定理17.4-1 信息的两种表达

满足上一节相应的一阶、二阶条件时,有

\[\mathbb E_\theta U_\theta =\int_A\partial_\theta p_\theta\,d\nu =\partial_\theta\int_Ap_\theta\,d\nu=0,\qquad I(\theta)=\operatorname{Var}_\theta(U_\theta) =-\mathbb E_\theta(\partial_\theta^2\ell_\theta).\]
(17.4-2)

证明。第一式来自密度归一化的一阶求导。对归一化等式再求一次导数,并用正密度区域内的恒等式,得到

\[\partial_\theta^2p_\theta =p_\theta\{(\partial_\theta\ell_\theta)^2+\partial_\theta^2\ell_\theta\}, \qquad 0=\int_A\partial_\theta^2p_\theta\,d\nu =\mathbb E_\theta(U_\theta^2)+\mathbb E_\theta(\partial_\theta^2\ell_\theta).\]
(17.4-3)

平方得分期望有限,加上二阶密度导数的可积控制,也保证这里对数二阶导数的期望存在。结合零均值即得两种表达。证毕。

因此Fisher信息衡量的是模型在该参数附近的变化,并依赖参数化和计量单位。它不是已经观察到的某条似然曲线的高度。观测信息 \(-\partial_\theta^2\ell_\theta(X)\) 是随机量;在上述条件下,其期望才等于 \(I(\theta)\),单次取值也不一定非负。

若 \(X_1,\ldots,X_n\) 相互独立,每个观测的模型均满足一阶条件并有有限信息,联合密度的对数为各对数密度之和,因此

\[U_\theta^{(n)}=\sum_{i=1}^n U_{\theta,i},\qquad I_n(\theta)=\sum_{i=1}^n I_i(\theta); \qquad \text{iid情形:}\ I_n(\theta)=nI_1(\theta).\]
(17.4-4)

最后一步使用独立性与每项零均值,使交叉项期望消失。非同分布的独立观测仍可逐项相加;有关联的数据不能直接套用 \(nI_1\)。这与16.5节不放回抽样必须保留协方差的原因一致。

17.4.3 无偏界、有偏界与等号条件

设 \(T\) 是不含未知参数的实值统计量,用来估计可微目标 \(g(\theta)\)。记 \(m(\theta)=\mathbb E_\theta T\),偏差为 \(b(\theta)=m(\theta)-g(\theta)\)。仅在某个参数点碰巧无偏,不足以令 \(m'=g'\);无偏性必须在邻域内成立,通常要求在整个参数空间成立。

定理17.4-2 标量信息不等式

在一阶正则性、估计量期望换序和有限方差条件下,若 \(0<I(\theta)<\infty\),则

\[\operatorname{Var}_\theta(T) \ge\frac{\{m'(\theta)\}^2}{I(\theta)} =\frac{\{g'(\theta)+b'(\theta)\}^2}{I(\theta)}. \quad b\equiv0\text{时:}\ \operatorname{Var}_\theta(T)\ge\frac{\{g'(\theta)\}^2}{I(\theta)}.\]
(17.4-5)

平方损失下,相应均方误差界为

\[\mathbb E_\theta\{T-g(\theta)\}^2 \ge b(\theta)^2+\frac{\{g'(\theta)+b'(\theta)\}^2}{I(\theta)}.\]
(17.4-6)

证明。按假设将 \(m(\theta)=\int_A T(x)p_\theta(x)\,d\nu(x)\) 求导,利用 \(\partial_\theta p_\theta=p_\theta U_\theta\) 和零均值得分,可得

\[m'(\theta)=\mathbb E_\theta(TU_\theta) =\operatorname{Cov}_\theta(T,U_\theta),\qquad \{m'(\theta)\}^2\le\operatorname{Var}_\theta(T)\operatorname{Var}_\theta(U_\theta).\]
(17.4-7)

末式是Cauchy–Schwarz不等式,得分方差等于信息,除以正的 \(I(\theta)\) 得式(17.4-5)。再用均方误差等于方差加偏差平方,得式(17.4-6)。证明同时覆盖无偏和有偏估计,无需对有偏情形另加正态近似。证毕。

固定参数点处达界,当且仅当Cauchy–Schwarz取等号,即

\[T-m(\theta)=\frac{m'(\theta)}{I(\theta)}\,U_\theta(X) \quad P_\theta\text{几乎必然}.\]
(17.4-8)

这也包括 \(m'=0\) 的情形:达界要求 \(T\) 在该分布下为常数。若要构造对所有参数都达界的估计量,等式必须由同一个不含未知参数的 \(T\) 满足。右边在化简前可以含 \(\theta\),不能把未知参数直接放进估计公式。达到无偏界便在相应无偏类中方差最小;反过来,无偏类中最优的估计量未必达到这个界。

17.4.4 精确达界与偏差的作用

例17.4-1 伯努利样本均值

设 \(X_i\) iid服从 \(\operatorname{Bernoulli}(p)\),且 \(0<p<1\)。两点支持集固定,每个内点的导数条件都可直接核查。计算得

\[U_p^{(n)}=\frac{\sum_iX_i-np}{p(1-p)},\qquad I_n(p)=\frac{n}{p(1-p)},\qquad \operatorname{Var}_p(\bar X)=\frac{p(1-p)}n=\frac1{I_n(p)}.\]
(17.4-9)

样本均值无偏,并且 \(\bar X-p=U_p^{(n)}/I_n(p)\),故精确达界。若 \(n=20,p=0.3\),信息为 \(2000/21\),方差为 \(0.0105\);这些是模型在假定参数点的数值,不是把一次样本比例当成真参数后的保证。端点 \(p=0,1\) 不在本定理的开参数空间内。

例17.4-2 有偏估计可有更小方差

设 \(X_i\) iid服从 \(N(\mu,4)\),\(n=16\),方差已知,目标是 \(\mu\)。联合得分为 \(4(\bar X-\mu)\),信息为四。样本均值无偏,方差等于 \(1/4\),达到无偏界。另取不含未知参数的收缩估计 \(T=0.8\bar X\),则

\[m(\mu)=0.8\mu,\quad b(\mu)=-0.2\mu,\quad \operatorname{Var}_\mu(T)=0.16=\frac{0.8^2}{4},\quad \operatorname{MSE}_\mu(T)=0.16+0.04\mu^2.\]
(17.4-10)

它达到与自身偏差函数相对应的有偏界,方差小于无偏界并不构成矛盾。均方误差在 \(|\mu|<1.5\) 时小于 \(0.25\),在 \(|\mu|>1.5\) 时更大。比如 \(\mu=0\) 时为 \(0.16\),\(\mu=2\) 时为 \(0.32\);降低方差的代价要与偏差一起比较。

正态样本均值与收缩估计的精确均方误差 样本量16,已知正态总体方差4。样本均值的均方误差恒为0.25;乘以0.8的收缩估计均方误差为0.16加0.04乘真实均值平方,两曲线在真实均值正负1.5处相交。真实均值零时收缩估计风险0.16,真实均值2时风险0.32。 正态样本:16 次独立观测,已知方差 4 均方误差 0.0 0.2 0.4 0.6 -3 -1.5 0 1.5 3 均值 0.16 样本均值:无偏,风险恒为 0.25 收缩估计:存在偏差
图17.4-1 例17.4-2的精确风险公式:\(\bar X\) 的均方误差为 \(0.25\),\(0.8\bar X\) 为 \(0.16+0.04\mu^2\)。阴影标出后者风险较低的参数区间;曲线由公式取点连线,没有模拟数据或抽样近似。
例17.4-3 指数均值达界,率参数不必达界

设 \(X_i\) iid服从率参数为 \(\lambda>0\) 的指数分布,密度为 \(\lambda e^{-\lambda x}\mathbf1_{\{x>0\}}\)。单个得分是 \(1/\lambda-X_i\),故 \(I_n(\lambda)=n/\lambda^2\)。估计均值 \(g(\lambda)=1/\lambda\) 时,信息界为 \(1/(n\lambda^2)\),样本均值恰有此方差,并满足式(17.4-8)。

若改为估计率参数,令 \(S=\sum_iX_i\)。它服从形状 \(n\)、率 \(\lambda\) 的Gamma分布。把密度乘以 \(s^{-k}\) 积分可得 \(\mathbb E_\lambda S^{-k}=\lambda^k\Gamma(n-k)/\Gamma(n)\),要求 \(n>k\)。于是当 \(n>2\) 时,

\[\widehat\lambda_U=\frac{n-1}{S},\qquad \mathbb E_\lambda\widehat\lambda_U=\lambda,\qquad \operatorname{Var}_\lambda(\widehat\lambda_U)=\frac{\lambda^2}{n-2} >\frac{\lambda^2}{n}.\]
(17.4-11)

例如 \(n=5,\lambda=2\),其方差为 \(4/3\),率参数无偏界为 \(4/5\);若观察到 \(S=2.5\),无偏率估计为 \(1.6\),最大似然率估计为二。后者期望为 \(n\lambda/(n-1)\),在有限样本中有偏。非线性变换目标以后,原来对均值达界的结论不能直接移用于率参数;也不能因一个估计是MLE便断言它精确达界。

17.4.5 多参数信息矩阵与目标梯度

设参数 \(\theta\in\Theta\subset\mathbb R^d\),\(\Theta\) 开放,\(d\) 固定;对各分量满足前述共同支持、局部可微和换序条件。得分变为列向量,信息变为矩阵:

\[U_\theta=\nabla_\theta\ell_\theta(X),\qquad \mathcal I(\theta)=\mathbb E_\theta(U_\theta U_\theta^\top) =\operatorname{Cov}_\theta(U_\theta) =-\mathbb E_\theta(\nabla_\theta^2\ell_\theta).\]
(17.4-12)

最后一种表达仍需二阶条件。矩阵总是半正定;以下要求它有限且正定,才使用通常的逆矩阵。信息奇异时,不能直接写出该逆,须另行分析参数化和可以估计的方向。

定理17.4-3 多参数信息界

实值统计量 \(T\) 的期望为 \(m(\theta)\),满足各分量的换序条件及有限方差,则

\[\operatorname{Var}_\theta(T)\ge \nabla m(\theta)^\top\mathcal I(\theta)^{-1}\nabla m(\theta). \quad \text{无偏估计 }g\text{ 时,取 }m=g.\]
(17.4-13)

证明。逐个分量求导给出 \(\operatorname{Cov}_\theta(U_\theta,T)=\nabla m\)。任取列向量 \(a\),对 \(T\) 与 \(a^\top U_\theta\) 应用Cauchy–Schwarz,得到

\[(a^\top\nabla m)^2\le \operatorname{Var}_\theta(T)\,a^\top\mathcal I a. \qquad a=\mathcal I^{-1}\nabla m.\]
(17.4-14)

当梯度不为零时,把这个 \(a\) 代入,两边分别出现 \((\nabla m^\top\mathcal I^{-1}\nabla m)^2\) 与同一正量的一次方,约去即得结论;梯度为零时下界为零。等号要求 \(T-m=(\nabla m)^\top\mathcal I^{-1}U_\theta\) 几乎必然成立。证毕。

对 \(k\) 维估计量 \(\mathbf T\),设 \(m=\mathbb E_\theta\mathbf T\),其Jacobian \(Dm\) 是 \(k\times d\) 矩阵。对每个线性组合 \(c^\top\mathbf T\) 使用上一定理,便得

\[\operatorname{Cov}_\theta(\mathbf T) \succeq Dm(\theta)\mathcal I(\theta)^{-1}Dm(\theta)^\top. \quad \mathbb E_\theta\mathbf T=\theta\text{ 时:}\ \operatorname{Cov}_\theta(\mathbf T)\succeq\mathcal I(\theta)^{-1}.\]
(17.4-15)

\(\succeq\) 表示两矩阵之差半正定,即每个线性组合的方差都满足相应下界,并非逐项比较非对角元素。若目标只是一维函数,其他未知分量仍是干扰参数,必须使用完整信息矩阵的逆。

以 iid的 \(N(\mu,v)\) 为例,参数取均值和方差 \((\mu,v)\),其中 \(\mu\in\mathbb R,v>0\)。单个得分的两分量是 \((X-\mu)/v\) 与 \(\{(X-\mu)^2-v\}/(2v^2)\)。正态三阶中心矩为零、四阶为 \(3v^2\),故

\[\mathcal I_n(\mu,v)= \begin{pmatrix}n/v&0\\0&n/(2v^2)\end{pmatrix}, \qquad \mathcal I_n^{-1}= \begin{pmatrix}v/n&0\\0&2v^2/n\end{pmatrix}.\]
(17.4-16)

样本均值仍达到均值的下界。对于 \(n\ge2\),无偏样本方差 \(s^2\) 的精确方差却是 \(2v^2/(n-1)\),严格大于方差参数的界 \(2v^2/n\),这里可调用16.2节的卡方分布。无偏最优与精确达界并不是同一个条件。

17.4.6 Uniform端点为何不满足这些条件

设 \(X_i\) iid服从 \(U(0,\theta)\),\(\theta>0\),最大值为 \(M\)。Lebesgue支配测度并不依赖参数,但正密度区域依赖端点。对一组严格位于内部的观测,忽略支持限制而求导,会得到

\[\widetilde U_\theta=-\frac n\theta,\qquad \mathbb E_\theta\widetilde U_\theta=-\frac n\theta\ne0,\qquad \mathbb E_\theta\widetilde U_\theta^2=\frac{n^2}{\theta^2}.\]
(17.4-17)

它是常数,方差为零,平方期望却不为零;甚至形式上的负二阶对数导数期望为 \(-n/\theta^2\)。本节的几个信息表达已经不一致。这里并非得分均值定理出了错,而是其假设没有成立。

单个观测的归一化积分,按Leibniz法则求导,具体显示遗漏的项:

\[\frac d{d\theta}\int_0^\theta\frac1\theta\,dx =\int_0^\theta-\frac1{\theta^2}\,dx+\frac1\theta =-\frac1\theta+\frac1\theta=0.\]
(17.4-18)

只积分内部导数会漏掉右端点的贡献;不能通过“固定观测几乎不会落在边界”把它消除。若在整个实轴上写带指示函数的密度,则相应的求导与积分换序同样不合法。

由独立性,\(M\) 的精确分布函数在 \(0\le m\le\theta\) 上为 \((m/\theta)^n\)。对其密度积分,有 \(\mathbb E_\theta M^k=n\theta^k/(n+k)\),\(k>0\)。因此

\[T_U=\frac{n+1}{n}M,\qquad \mathbb E_\theta T_U=\theta,\qquad \operatorname{Var}_\theta(T_U)=\frac{\theta^2}{n(n+2)} <\frac{\theta^2}{n^2}.\]
(17.4-19)

最右边恰是把式(17.4-17)的平方期望误用为正规信息后得到的“界”。例如 \(n=4,\theta=2\),真实方差为 \(1/6\),这个错误界为 \(1/4\)。若再误用可加性,还会得到另一错误界 \(\theta^2/n=1\)。低于错误界只说明公式使用失当。支持随参数变化的模型必须单独检查;这一个反例也不能代替对所有其他非正则模型的分析。

Uniform密度的支持集随端点改变 横轴从0到3.5,纵轴为密度。端点为2时,密度在0到2等于0.5;端点为3时,在0到3等于三分之一。蓝色阴影从2到3显示参数变化新增的正密度区域。只求原支持内部高度的导数,会遗漏变化边界贡献的概率质量。 端点改变时,正密度区域也改变 端点为 2 端点为 3 密度 0.0 0.2 0.4 0.6 0.5 1/3 0 1 2 3 3.5 观测值 阴影:新增的正密度区域;内部导数不能替代完整的积分求导
图17.4-2 \(U(0,2)\) 与 \(U(0,3)\) 的精确密度,阴影为 \((2,3)\) 内新增的概率质量。两者都可用Lebesgue测度作支配测度,但共同正密度支持的条件不成立;密度积分的导数必须计入边界变化。

本节的界约束满足所述条件的估计量,未承诺下界总能达到,更没有给出所有MLE的有限样本最优性。正则MLE的渐近有效性还需要一致性、内部真参数、可逆信息和余项控制等条件,相关论证见19.3节。统计推断中精确结论与渐近结论应分别使用。

17.4.7 习题与解答

习题17.4-1 参考点得分的均值

令 \(X\sim\operatorname{Bernoulli}(p)\),在固定参考点 \(p_0\in(0,1)\) 求得分 \(U_{p_0}(X)\)。计算它在真实参数 \(p\) 下的期望,并解释为什么“得分均值为零”不是对任意参考点的陈述。

参考点得分为 \((X-p_0)/\{p_0(1-p_0)\}\),真实分布给出期望 \((p-p_0)/\{p_0(1-p_0)\}\)。分母为正,故只有两参数相同时期望才为零。归一化积分求导的证明中,密度与导数都取在同一个参数位置;把期望密度换成另一个参数,就不能沿用那一步。固定参考点得分可以从数据计算,而代入未知真参数的表达主要用于理论分析。

习题17.4-2 指数率MLE的有偏信息界

令 \(n=6\),\(X_i\) iid服从率 \(\lambda\) 的指数分布,\(S=\sum_iX_i\),\(T=6/S\)。计算 \(T\) 的偏差、方差与有偏信息界,并在 \(\lambda=3\) 时给出均方误差和其下界。

由Gamma逆矩公式,\(\mathbb E_\lambda T=6\lambda/5\),偏差为 \(\lambda/5\),均值导数为 \(6/5\)。又有 \(\mathbb E S^{-2}=\lambda^2/(5\cdot4)\),故方差为 \(36\lambda^2/20-36\lambda^2/25=9\lambda^2/25\)。完整样本信息为 \(6/\lambda^2\),有偏方差界为 \(6\lambda^2/25\),严格小于实际方差。加偏差平方后,均方误差为 \(10\lambda^2/25\),下界为 \(7\lambda^2/25\)。当率为三时,二者分别为 \(3.6\) 与 \(2.52\);MLE的名称没有使不等式自动取等号。

习题17.4-3 不能用对角元倒数代替完整逆矩阵

某正则二参数模型在考察点的信息矩阵为 \(\begin{pmatrix}4&1\\1&1\end{pmatrix}\)。估计第一个参数时计算无偏界;若第二个参数已知,界有何变化?说明这一比较的条件。

矩阵的行列式为三且首个主子式为四,故正定。逆矩阵为 \(\frac13\begin{pmatrix}1&-1\\-1&4\end{pmatrix}\)。目标梯度为 \((1,0)^\top\),所以两参数都未知时下界为三分之一。第二个参数已知并固定时,限制为一参数子模型,在同一点的信息为四,下界为四分之一。知道干扰参数能提高可用精度;两种比较针对不同的参数模型与无偏约束,不能对未知干扰参数时直接取首个对角元的倒数。

习题17.4-4 找出Uniform证明断裂的位置

对本节的 \(U(0,\theta)\) 模型,取无偏估计 \(T_U=(n+1)M/n\)。比较 \(d\mathbb E_\theta T_U/d\theta\) 与 \(\operatorname{Cov}_\theta(T_U,\widetilde U_\theta)\)。这指出Cramér–Rao证明的哪一步不能使用?

由于最大值的精确一阶矩,估计量期望就是端点参数,导数为一。形式得分 \(\widetilde U_\theta=-n/\theta\) 在固定参数下是常数,与任意有限方差统计量的协方差都为零。这两个量不相等,说明期望导数等于与得分协方差的关键恒等式在这里失效。根源包括支持边界改变以及由此造成的换序和零均值条件失效;Cauchy–Schwarz不等式本身仍然正确。不能把这个反例解释为一个合法信息界被更好估计量突破。

关键术语

  • 正则支配族:本节要求共同正密度支持及适当的局部可微、可积控制条件。
  • 得分函数:对数密度关于参数的导数,在真参数处描述局部模型变化。
  • Fisher信息:正则模型下得分的协方差,标量情形为其方差。
  • 观测信息:给定数据时的负对数密度Hessian,其期望在二阶条件下等于Fisher信息。
  • Cramér–Rao界:由期望导数与得分协方差得到的方差下界。
  • 有偏信息界:使用估计量均值导数,必须把目标导数与偏差导数一起计入。
  • 达界条件:中心化估计量几乎必然等于得分的指定线性组合。
  • 干扰参数:模型中影响推断、但不是当前目标的其他未知参数。
  • 半正定比较:矩阵差对所有向量的二次型非负,表示所有线性组合的方差比较。
  • 非正则端点模型:本节Uniform例的支持边界变化,破坏通常的信息恒等式。

参考文献

  1. Peter Kempthorne,MIT 18.655,Lecture 13:Unbiased Estimation and Risk Inequalities,2016。讲义陈述期望为一般函数的标量信息不等式、得分恒等式、可加性和多参数扩展。
  2. UC Berkeley,Stat 210A Reader:Score Function and Fisher Information,2025。核对得分参考点、信息的两种表达、标量及多参数下界,以及有限样本达界的限制。
  3. Sara van de Geer,Mathematical Statistics,2010,ETH Zürich讲义。§3.3–3.4给出平方可积可微性的条件和信息不等式证明,§3.2.3讨论Uniform最大值。