I understand that when applying GLMMs (e.g. in logistic mixed effects regression), the interpretation of the coefficients for the fixed effects is that they are also conditional on the random effects (e.g. a unit change in a predictor marks the log odds change in an outcome variable only for subjects with the same random effects). These are very likely quite different from the marginal change for the population. I am looking for a way to estimate and test this marginal change.
I read this answer and the related links. It is great, but the package used here is restricted to using a single grouping variable. My model has random effects both for subject and for item (subjects are measured over different items), so I believe it does not apply.
I would also be happy using GEEs if it can help me, but my specification includes also random slopes (and 2 grouping variables), so I think they do not apply either. Can someone recommend a way to estimate the marginal effects for such data? I am less interested in the conditional effects.