如何计算 statsmodels 线性混合效应模型的 AIC 值？

Question

线性混合效应模型 (link) 的 Statsmodels 文档页面声称“statsmodels LME 框架目前支持通过 Wald 检验和系数的置信区间、轮廓似然分析、似然比检验和AIC。”

我已尝试按如下方式检索 AIC：

model = smf.mixedlm(formula=formula,data=data_subset,groups=data_subset["m#"]).fit()
model.aic

但结果值是 'nan'。

我还查看了有关混合lm() 和混合lm.fit() 函数的文档，但无法找到任何可以表明他们希望计算 AIC 的参数。

感谢任何帮助！

Answer 1

调用 model.fit() 时应使用 reml=False。