[R] bootstap and lme4

From: FabbroThomas <thomas.fabbro_at_unifr.ch>
Date: Mon 09 May 2005 - 22:54:02 EST


Hi,

I am trying to get bootstrap confidence intervals on variance components and related statistics. To calculate the variance components I use the package lme4.

> off.fun <- function(data, i){
   d <- data[i,]
   lme1<- lmer(y ~ trt + (trt-1|group), d)    VarCorr(lme1)@reSumry$group[2,1] #just as an example
}
> off.boot <- boot(data=data.sim, statistic=off.fun, R=100)

If I choose small values of R (<10) then I get very reasonable results but for large R I get the following error massage:

Can anyone tell me why I get this error massage and what I can do to avoid this problem?

Thank you very much for your help!

Thomas

I am working with:
R 2.1.0
boot 1.2-22

lme4 0.95-6
Matrix 0.95-7
lattice 0.11-6
latticeExtra 0.1-1

On Mac OS X 10.3.9



R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html Received on Mon May 09 23:02:11 2005

This archive was generated by hypermail 2.1.8 : Fri 03 Mar 2006 - 03:31:40 EST