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

Question

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

The Statsmodels LME framework currently supports post-estimation inference via Wald tests and confidence intervals on the coefficients, profile likelihood analysis, likelihood ratio testing, and AIC. [emphasis added]

我注意到了MixedLM.loglike方法，但是我似乎找不到一个运行似然比测试的函数/方法。在

有人能帮我指一下正确的方向吗？在

Answer 1

我正在运行一个开发分支，因此情况可能发生了变化，但是混合lm.fit（）应具有名为“llf”的属性。这是对数似然函数在估计参数下的值。如果你有两个嵌套的模型并且取-2倍于它们的llf值的差，在一个更简单的模型为真的零假设下，这将是一个chi^2随机变量，其自由度等于两个模型的自由度差。在

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