如何对线性混合效应模型进行似然比检验?

Question

Statsmodels的线性混合效应模型的文档声称

statmodels LME框架目前支持通过Wald检验和系数置信区间、剖面似然分析、似然比检验和AIC进行后估计推断。(重点补充道)

我已经注意到MixedLM.loglike方法，但我似乎找不到运行似然比检验的函数/方法。

谁能给我指个正确的方向?

Answer 1

我正在运行一个开发分支，所以事情可能已经改变，但MixedLM.fit()返回的结果类应该有一个名为'llf'的属性。这是对数似然函数在估计参数处的值。如果你有两个嵌套模型，取-2倍的llf值之差，在零假设下，更简单的模型成立，这将是一个chi^2随机变量，其自由度等于两个模型的自由度之差。

注意，当使用LR测试时，许多人认为您应该将估计器切换到ML(而不是默认的REML)。