如何解决R 中的 lmerTest 与 statsmodels/Python 中的 lme
我确实有一个复杂的线性混合效应模型。
model_lme = MixedLM.from_formula("ratings ~ LCDev455*ST0LT1*targettype*Female0Male1",groups=data["ResponseID"],data=data)
model_lme = model_lme.fit()
model_lme.summary()
这在 R、lmertest 中非常有效。甚至没有收敛警告。但是在python statsmodels中,它给出了错误:奇异矩阵。
我认为这两个包(R 中的 lmertest,Python 中的 lme)给出了相同的输出。
有人知道为什么会发生这种情况吗?我更像是一个初学者,所以我无法在这里真正模拟我的数据。
谢谢!
版权声明:本文内容由互联网用户自发贡献,该文观点与技术仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 dio@foxmail.com 举报,一经查实,本站将立刻删除。