Dear all!

 

Probably you are all sick of this topic but i haven't found a solution for this issue clicking through older posts. I want to know your opinion about the following approach using two great functions of the lme4 package: simulate() and refit().

 

With these functions I can easily generate, say x, "bootstrap-models" from an underlying model (of class mer) fitted to my data

 

for(i in 1:x) refit(model, simulate(model))

 

For each of these models, i can perform prediction (e.g. with fitted()) leading to x predictions. For a large x, the 95% percentile confidence interval should be a 95% prediction interval. Did I overlook something?

 

I would be very happy to read your comments.

 

Thanks a lot!

 

Kind regards,

Alex
 		 	   		  
	[[alternative HTML version deleted]]

