跳到正文
格致开物MATHWIKI

局部渐近正态性

局部渐近正态性(local asymptotic normality,LAN)描述统计实验在真参数附近的似然形状:把参数差缩小到样本量平方根的倒数尺度后,对数似然比可近似为一个随机线性项减去确定的二次项。它是模型或实验序列的性质,并不是“所有估计量都近似正态”的同义词。

先看独立伯努利观测 Xi∼Bernoulli⁡(p0),其中 0<p0<1。若 n=100、候选基础概率 p0=0.4,一个局部变化 h=0.5 对应概率 p0+h/n=0.45。样本量增大时,比较的两个参数越来越近,而非始终相差0.05。

伯努利模型的直接展开

记成功总数为 Sn、d=h/n。对数似然比为 Λn(h)=Snlog⁡(1+d/p0)+(n−Sn)log⁡{1−d/(1−p0)}. 固定 h,充分大的 n 使两个对数的自变量远离零。用 log⁡(1+z)=z−z2/2+O(z3),整理线性项: Sndp0−(n−Sn)d1−p0=hΔn,Δn=Sn−np0np0(1−p0). 二次项则为 −h22{Sn/np02+1−Sn/n(1−p0)2}. 由大数律,括号趋于 I(p0)=1/[p0(1−p0)]。三阶余项的绝对值可界为常数乘 n|d|3=|h|3/n,故趋零。伯努利中心极限定理又给出 Δn→dN(0,I(p0))。因此完整得到 Λn(h)=hΔn−12h2I(p0)+oPp0(1). 这里的随机线性系数来自样本波动,负二次项控制局部似然的弯曲程度。

一次样本的精确值与二次值

设百次试验观察到四十六次成功。则 Δn=46−4010×0.24=2.5,I(p0)=25/6. 取 h=0.5,二次展开值为 0.5×2.5−(0.5)2(25/6)/2=35/48≈0.7292;精确值可直接计算为 46log⁡(0.45/0.4)+54log⁡(0.55/0.6)≈0.7194. 两者接近但不同,差来自有限样本展开余项与样本二阶系数。此处的基础概率是比较的参考值,并非从样本得到了“已知真值”。

百次伯努利四十六次成功的局部精确对数似然比与二次近似在参考值附近接近
横轴h对应p等于0.4加h除以十。远离局部中心时,有限样本高阶项更加明显。

实验序列中的定义

设 Pθ,n 是 n 规模实验的分布,θ0 是参数空间内点。在通常的非退化定义中,若存在正定信息矩阵 I 和随机向量 Δn,对每个固定局部向量 h 有 log⁡dPθ0+h/n,ndPθ0,n=h𝖳Δn−12h𝖳Ih+oP(1),Δn→dN(0,I)在 Pθ0,n 下, 就称该实验在 θ0 局部渐近正态。似然比在参考分布下解释;需要讨论更一般支配关系时,要相应处理测度的绝对连续部分。Duchi讲义给出该定义及二次均值可微性推出LAN的定理。[1]

对于独立同分布的二次均值可微密度族,得分二阶可积时可建立相应LAN展开;信息正定再保证通常的非退化局部比较。一般定理需控制密度平方根的余项,本篇已完整证明的是内部概率的伯努利特例。

变式。 若 Xi∼N(θ,4),令 Δn=n(X‾−θ0)/4,则直接展开平方差就精确得到 Λn(h)=hΔn−h2/8. 参考分布下 Δn∼N(0,1/4),没有余项,信息为 1/4。这展示了“局部像高斯实验”的一个本身就是高斯的基准。

如果伯努利基础概率为零,分布退化、两侧局部扰动不都合法,前述含 1/[p0(1−p0)] 的证明失效。参数边界、不可识别性或不同收敛尺度都应单独分析,不能仅凭似然能被软件求导便套用LAN。

参考文献

  1. ↑ John Duchi,Local Asymptotics and Optimality,Stanford Stats300b,2021,第6—12张,QMD、LAN及高斯平移实验。