使用python的混合效果模型（statemodels）

Question

我有一个包含以下列的数据集：

“sampno_vehno”、“车辆年龄”、“车辆车身类型”、“车辆燃油类型”、“行程”、“高峰时间”、“周末”、“持续时间”和“波动性”

我想在satemodels（一个python库）中使用mixedlm的混合效果模型，但我不知道如何编写代码

我尝试了以下代码，我知道它有一些问题

import statsmodels.formula.api as smf
import pandas as pd

df = pd.DataFrame(myData, columns=['sampno_vehno', 'vehicle_age', 'vehicle_body_type',
                                   'vehicle_fuel_type', 'trip', 'rushHours', 'weekend',
                                   'duration_min', 'volatility'], dtype=float)

model = 'volatility ~ vehicle_age + C(vehicle_body_type) + C(vehicle_fuel_type) + C(rushHours) + C(weekend) + duration_min + trip'
md = smf.mixedlm(model, df, groups=df['sampno_vehno'],re_formula = "~trip")
mdf = md.fit()
print(mdf.summary())