Модель линейного смешанного эффекта в python с более чем несколькими независимыми переменными

Question

Я должен получить случайный наклон и перехват для моих данных. Ниже приведено уравнение, которое я в итоге написал, это правильный способ сделать это?

md = smf.mixedlm("""UNITS_sold ~ Price_1 + Price_2 + Price_3 + Price_3 + 
                    Coupon_1 +  Coupon_2 + Coupon_3 + Coupon_4""", 
                 df, 
                 groups = df["id"], 
                 re_formula = """~Price_1 + Price_2 + Price_3 + Price_3 + 
                                  Coupon_1 +  Coupon_2 + Coupon_3 + Coupon_4""") 

mdf = md.fit() 
print(mdf.summary())

Здесь купон является категориальной переменной, id - это моя иерархия уровня 1, и я должен предсказать единицы измерения. Товар продается.

Спасибо!