“Slutsky 定理是 plug in 标准误、Wald 统计量和 nuisance 参数替换的基本接口。它只保证极限分布不变,不提供替换造成的有限样本误差大小。”
形式陈述
检验
其中
在内部真参数、可识别、二次可微和非奇异信息等条件下,三者在
直觉
三种检验都读取零假设附近同一个二次 likelihood 几何,只是取样位置不同:Wald 看完整估计量离约束面多远,Score 在受限点看上坡方向有多陡,LR 直接比较受限与完整最优值的高度差。局部曲率近似良好时三幅图首阶相同,强曲率、边界或优化失败时则会分开。
局部等价
对 score 在
因此三者在
例子与边界
取
若
向量约束与可复核数值
对
其中
Poisson 例取
这个例子中三者分别使用估计点的尺度、零点的尺度与两个最优值之间的似然差,因此得到不同的有限样本数值。
边界与失败情形
Wald 使用参数坐标中的距离,非线性重参数化会改变其有限样本数值,所得区间还可能越出参数空间。Score 依赖零假设处的信息矩阵;LR 比较同一模型中的两个最优似然值,因此对一一重参数化不变。
参数在边界、信息奇异或 nuisance 只在备择出现时,普通卡方极限失败。此时要按具体模型重新确定极限分布。
推论与应用
模型失配时,sandwich 协方差可以校准某些 Wald 或 Score 统计量;似然比的极限则要分析目标函数及其曲率,普通 Wilks 定理与信息等式可能不再适用。
计算上,Score 只拟合零模型,适合高维备择难拟合的筛查;Wald 需要完整 MLE 与协方差;LR 需要两次优化但沿 likelihood 几何比较。若完整 MLE 不存在,Wald 和 LR 可能失效,而某些有限的 score 检验仍可计算,但其极限条件仍需检查。
对 nuisance 参数,Score 必须用 efficient score
信息为相应 Schur complement。这一修正扣除了干扰参数估计与目标 score 的相关部分;直接使用目标信息块会错估相应方差。
三种检验的p 值差异明显时,可以检查参数是否接近边界、profile 曲线是否偏离二次形状,以及样本量和优化收敛情况。这些诊断帮助解释预先选定的主要检验为何与其他方法不同。正则零假设下的共同极限由卡方分布校准,但边界和奇异模型需要别的极限定理。
置信集合也对应不同反演:Wald 给局部椭球,LR 给 profile likelihood level set,Score 给零点可接受域。它们首阶覆盖相同,却可能在强曲率模型中有不同不对称边界;参数空间受限时 LR 区域通常更自然。
若只获得 approximate MLE,完整模型与受限模型的优化目标都应距各自上确界
参考资料
- A. Wald, “Tests of Statistical Hypotheses Concerning Several Parameters,” Transactions of the AMS 54, 1943。
- C. R. Rao, “Large Sample Tests of Statistical Hypotheses,” 1948。
- A. W. van der Vaart, Asymptotic Statistics, Cambridge University Press, 1998,Ch. 16。