python混合效果模型（statemodels）

Question

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

“ sampno_vehno”，“ vehicle_age”，“ vehicle_body_type”，“ vehicle_fuel_type”，“ trip”，“ rushHours”，“ weekend”，“ duration_min”和“ volatility”。

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

我尝试了以下代码，但我知道它存在一些问题。

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