如何在Python statsmodels线性混合效果模型中拥有多个组?

I am trying to use the Python statsmodels linear mixed effects model to fit a model that has two random intercepts, e.g. two groups. I cannot figure out how to initialize the model so that I can do this. I am trying to use the Python statsmodels linea