如何在大熊猫中建立具有两个固定效应和一个随机效应的混合线性模型？

Question

编码新手，所以我不确定该怎么做。
我有一个看起来像这样的df：

        ID        A          B
0       1         20         4.0
1       1         24         5.0
2       2         24         0.0
3       2         9          2.0
4       3         24         1.0
5       3         24         0.0

我想制作一个混合线性模型，使我可以将A和B相关以显示ID内的可变性。 A和B是我的固定效果，而ID是我的随机效果。我希望它能使A依赖B。
我尝试使用此代码来完成此操作：

import statsmodels.api as sm
model = sm.MixedLM.from_formula("B~ A", df1, groups=df1["ID"])
result = model.fit()
print(result.summary())

但是，我一直收到无效的语法错误。
我还想知道如何在散点图上绘制该图，因为我想根据每个ID的每条线的斜率找到A和B之间的相关性。