title: Is there a way of getting "marginal effects" from a `glmer` object
body: <p>I am estimating random effects logit model using <code>glmer</code> and I would like to report Marginal Effects for the independent variables. For <code>glm</code> models, package <code>mfx</code> helps compute marginal effects. Is there any package or function for <code>glmer</code> objects?</p>
  <p>Thanks for your help.</p>
  <p>A reproducible example is given below</p>
  <pre><code>## mydata &lt;- read.csv(&quot;http://www.ats.ucla.edu/stat/data/binary.csv&quot;)
  ## as of 2020-08-24:
  mydata &lt;- read.csv(&quot;https://stats.idre.ucla.edu/stat/data/binary.csv&quot;)
  mydata$rank &lt;- factor(mydata$rank) #creating ranks
  id &lt;- rep(1:ceiling(nrow(mydata)/2), times=c(2)) #creating ID variable
  mydata &lt;- cbind(mydata,data.frame(id,stringsAsFactors=FALSE)) 
  set.seed(12345)
  mydata$ran &lt;- runif(nrow(mydata),0,1) #creating a random variable
  
  library(lme4)
  cfelr &lt;- glmer(admit ~ (1 | id) + rank + gpa + ran + gre, data=mydata ,family = binomial)
  summary(cfelr)
  </code></pre>