r-glmmadaptive

Generalized linear mixed models can be fitted in R with adaptive Gaussian quadrature. It supports longitudinal and clustered data analyses where random effects and non-Gaussian outcomes appear together.