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")

但是，出现以下错误：

unhashable type: 'numpy.ndarray'

您能帮我解决这个问题吗？该模型可以在R中轻松完成，但我无法在python中找到兼容的示例。