python中多随机效应的混合效应建模

Question

我有以下数据集，我想运行一个混合效应线性模型，如下所示：

 Type    Season    Price     Sales
  A       Fall      34        500
  A       Winter    20        300
  B       Fall      15        100
  B       Winter    10        150

fm = model = smf.mixedlm("Sales ~ Price", df, groups = df[["Type","Season"]], 
                re_formula = "~Price")

但是，我得到了以下错误：

^{pr2}$

你能帮我解决这个问题吗。这个模型在R中很容易实现，但是我在python中找不到兼容的示例。在