Is there a Python module that can fit a generalized linear mixed model? [closed]

I’ve built a mixed model in R using glmer from the lme4 package. 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 statsmodels module, but it is not clear to me if statsmodels will be able to fit the same model that glmer was able to.


Source : Link , Question Author : OuijiBoard , Answer Author : Community

Leave a Comment