bootMer compatibility?
AndrewLawrence opened this issue · comments
Hey, thanks for your work on this really useful software!
I was wondering if blme models are compatible with lme4::bootMer
?
lme4::bootMer(my_bmerMod, my_icc_function, nsim = 1000, use.u = FALSE, type = "parametric")
gives no obvious errors.- I get sensible looking results for my purpose (bootstrapping ICC).
lme4::bootMer
uses a generic refit method (which blme defines for bmerMods) so it isn't just running lmers.
Thanks in advance!
I'm not entirely sure - I've only briefly reviewed the code for bootMer
- but from my understanding the function works for blme
objects but understates uncertainty by simulating responses conditioned on the estimates of the model parameters. If your perspective on the use of priors is to regularize the estimates, then this is fine. However, if you're looking for a fully Bayesian approach the correct thing to do would be to simulate from the posterior predictive distribution. This would would require a full Bayes solution, like MCMCglmm
or rstanarm
.
A little bit of extra uncertainty can be added for the parameters by conditioning on the covariances of the random effects. That would require rewriting lme4:::.simulateFun
. The truly hard bit is simulating for the random effects themselves.
Thanks, this is very helpful! I am approaching this from a regularisation perspective (to add confidence intervals to Chen et al's "RME" method).