Robust Standard Errors for Mixed-Effects Models in Lme4 Package of R

Huber-White robust standard errors for a GLMM - R

This looks like a bug in the package, as far as I can tell (the bread.glmerMod function was calling estfun.lmerMod rather than estfun.glmerMod; there's a broader question here about the design of the generic functions, but never mind ...)

You should be able to install a fixed version from my fork via remotes::install_github("bbolker/merDeriv"), then reload the package and try again.

Alternately, download the tarball, change vcov.lmerMod to vcov.glmerMod in the last line of R/bread.glmerMod.R, and re-install the package ...

R code for pulling out fixed effect standard errors in lme4 package

summary(model)$coef[, 2, drop = FALSE]

Or

sqrt(diag(vcov(model)))


Related Topics



Leave a reply



Submit