I’ve built a mixed model in R using
lme4package. I now wish to incorporate this model in Python. My R code is of the form
fit_model <- glmer(A * B * C + (1|x) + (1|y) + (1|z), data=data, family=binomial, nAGQ=0).
Can this model fitting be replicated in Python? I’ve come across the
statsmodelsmodule, but it is not clear to me if
statsmodelswill be able to fit the same model that
glmerwas able to.