Fred,

A simple thing to do to bootstrap the fixed effects is

z <- lmer(...)
simz <- data.frame(t(apply(simulate(z, nsim = 2000),2, function(y)  
fixef(refit(z, y)))),
                    check.names = FALSE)

This gives you nsim (=2000) bootstrap estimates that you can use.   
Presumably, you can do the same for other components of the model.

Cheers, Tony

On Jan 23, 2012, at 6:57 PM, Jianyun Wu wrote:

> Hi All,
>
> Is there any package in R for bootstrapping linear mixed effects  
> model?
>
> I came cross a  presentation describing a package merBoot with various
> bootstrapping procedures for mixed effect models. However the  
> package was
> never released.
>
> Any alternative out there?
>
> Thanks
>
> Fred
>
> 	[[alternative HTML version deleted]]
>
> _______________________________________________
> R-sig-mixed-models@r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-sig-mixed-models

Anthony Ragnar Ives
Department of Zoology
UW-Madison
(608) 262-1519


	[[alternative HTML version deleted]]


