跳到正文
格致开物MATHWIKI

似然比检验

似然比检验(likelihood-ratio test)比较零假设允许的最佳拟合,与整个模型允许的最佳拟合。若零假设受到的限制使最大似然下降很多,数据便倾向于拒绝该限制。统计量的大小要结合零假设下的参照分布判断,不能仅凭“拟合改善了”作结论。

例如独立进行二十次伯努利试验,观察到十四次成功,要检验成功概率是否为 p0=1/2。零假设只有一个候选值,整个模型 0<p<1 的最大似然估计则为 p^=14/20=0.7。两个模型都解释相同的二十次试验,似然可以直接比较。

受限最大值与统计量

设零假设参数空间 Θ0⊆Θ,定义 Λ=supθ∈Θ0L(θ)supθ∈ΘL(θ),WLR=−2log⁡Λ. 只要两个上确界正且有限,就有 0<Λ≤1、WLR≥0。需要的是各参数空间上的全局上确界;数值优化停在局部点,可能产生错误统计量。某些边界例子只存在上确界而没有达到它的参数,因此定义写作上确界更稳妥。

伯努利样本的对数似然略去与参数无关的项为 ℓ(p)=slog⁡p+(n−s)log⁡(1−p). 在 0<s<n 时,导数为 s/p−(n−s)/(1−p),二阶导数为负,所以唯一极大值是 s/n。本例由此得到 WLR=2[14log⁡0.70.5+6log⁡0.30.5]≈3.29131514. 即对数似然改善约 1.64566,再乘二得到检验统计量。

二十次中十四次成功的相对似然在0.7达到峰值,0.5的对数似然比峰值低约1.646
纵轴为对数似然减去最大值。似然比使用两个峰值之差;参数轴上的距离只是其中一个因素。

伯努利情形的卡方近似推导

写 q=p^,有 WLR=2nD(q‖p0), 其中 D(q‖p0)=qlog⁡qp0+(1−q)log⁡1−q1−p0. 固定内部点 0<p0<1。把它看成 q 的函数,一阶导数在 q=p0 为零,二阶导数为 1/[q(1−q)],所以在 p0 邻域 D(q‖p0)=(q−p0)22p0(1−p0)+O(|q−p0|3). 零假设下,大数定律给出 q→p0,中心极限定理给出 Zn=n(q−p0)p0(1−p0) →d N(0,1). 于是 q−p0=Op(n−1/2),三阶余项乘 n 后为 Op(n−1/2),从而 WLR=Zn2+op(1) →d χ12. 这给出了本模型的一项完整渐近证明。一般的威尔克斯定理把它推广到光滑的嵌套参数模型:在可识别、内部真参数、非奇异信息矩阵、适当光滑与积分条件以及一致极大似然估计等正则条件下,极限自由度为两模型维数之差;光滑约束还需具有相应满秩。一般定理的条件与局部二次证明见Geyer讲义。[1]

二十次试验的有限样本校准

用 χ12 近似,本例尾概率约为 0.06965,五百分比临界值约 3.84146,因而没有在该近似规则下拒绝零假设。

二十次试验不必依赖渐近近似。零假设下 S∼Binomial⁡(20,1/2),似然比统计量关于 S=10 对称,并随 |S−10| 增大。以不小于观察统计量为尾部,精确概率为 pexact=Pr0(S≤6)+Pr0(S≥14)=2∑k=1420(20k)2−20=0.1153182983…. 精确值比卡方近似大,说明小样本、离散分布的校准可能有明显区别。它们使用同一个似然比统计量,但不是同一个参照分布。

成功率二分之一的二项分布在不多于六次或不少于十四次成功处构成似然比尾部
红色端点包含与十四次成功同样极端的六次成功。尾部按离散概率求和,不是连续卡方曲线下面积。

已知正态方差的变式

若 Xi∼N(μ,4) 独立,n=9、X‾=0.8,检验 μ=0,配方直接给出 WLR=nX‾24=1.44. 零假设下 nX‾/2 恰好服从标准正态,所以统计量恰好服从 χ12,尾概率约 0.23014。这里不是仅在大样本时成立。

一一参数变换不会改变两个模型中的最大似然值,因此似然比统计量具有这种不变性。相反,方差为零的边界、混合模型的不可识别点或奇异信息矩阵可能破坏通常的卡方极限;不能只按“参数少了几个”指定自由度。复合似然也不是原始联合密度,其似然比常需额外校准。瓦尔德检验和得分检验在正则大样本中与似然比检验密切相关,但有限样本统计量可以不同。

参考文献

  1. ↑ Charles J. Geyer,The Wilks, Wald, and Rao Tests,§2—4,正则模型、光滑子模型及威尔克斯定理;§8讨论以期望信息代替观测信息。