Is there a way of getting “marginal effects” from a `glmer` object

Question

I am estimating random effects logit model using glmer and I would like to report Marginal Effects for the independent variables. For glm models, p

Answer 1

Here's an approach using the margins() package:

library(margins)
library(lme4)

gm1 <- glmer(cbind(incidence, size - incidence) ~ period +
                 (1 | herd),
             data = cbpp,
             family = binomial)

m <- margins(gm1, data = cbpp)
m