python - R 中的 lmerTest 与 statsmodels/Python 中的 lme

Question

我确实有一个复杂的线性混合效应模型。

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）给出了相同的输出。

有谁知道为什么会发生这种情况？我更像是一个初学者，所以我不能在这里真正模拟我的数据。

谢谢！