如何在混合效应线性回归（Python）中运行变量交互？

Question

我正在用 python 分析我的数据框（附上表格的图片）。我能够将所需的变量拟合到混合效应线性回归模型中：

model = smf.mixedlm("stage ~ overallscore + spatialreasoning + numericalmem",
                    data=mf, groups="group", re_formula="1")
result = model.fit()

print(result)

result.summary()

摘要打印。但是当我试图获得变量“stage”和“overallscore”之间的相互作用时，我遇到了一个问题。

model2 = smf.mixedlm("stage ~ overallscore + spatialreasoning + numericalmem:stage,overallscore(5)",
                    data=mf, groups="group", re_formula="1")
result2 = model2.fit()

print(result2)

result2.summary()

它说系列对象不可调用。

我什至在没有指定阶段 5 的情况下尝试过这个，但它仍然不起作用。然后我尝试输入其他变量而不是这两个变量，但仍然出现相同的错误。我不确定如何解决这个问题，因此将不胜感激 :) 。到目前为止，这是我一直在使用的指导（按照他们的方法）：https://www.pythonfordatascience.org/mixed-effects-regression-python/