Как перевести модель lme4 в statsmodels в python?

Question

Как перевести следующую трехуровневую модель в statsmodel в Python?

fit = lmer(Y~ 1 + I(x1/100) + x2 + x3 + I(x1/x4) + x5 +
                   (1|level1:level2:level3) + 
                   (1|level1:level2) +
                   (1+ I(x1/100) + x2+ x3|level1),
                 data=df,
                 na.action = na.exclude,
                 REML=FALSE,
                 control = lmerControl(optimizer = "optimx", calc.derivs = FALSE,
     optCtrl = list(method = "nlminb", starttests = FALSE, kkt = FALSE)))