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

Question

Statsmodels的线性混合效果模型的documentation声称

  

Statsmodels LME框架目前支持通过Wald检验的估计后推断和系数，轮廓似然分析，似然比检验和AIC的置信区间。 [强调补充]

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

有人可以指出我正确的方向吗？

Answer 1

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

请注意，许多人认为在使用LR测试时应将估算器切换为ML（而不是默认的REML）。