Displaying 2 results from an estimated 2 matches for "bootout".
2011 Feb 16
2
boot.ci error with large data sets
...n(E$A, E$B))}
#Some fake data:
test <- data.frame(rnorm(1000, 5), rnorm(1000, 3))
test1 <- data.frame(rnorm(10000, 5), rnorm(10000, 3))
names(test) <- c("A", "B")
names(test1) <- c("A", "B")
# Getting the boot object and the CI, seem to works fine
bootout <- boot(test, w.mean, R=1000, stype="i")
(bootci <- boot.ci(bootout, conf = 0.95, type = "bca"))
# Now with a bigger data set, boot.ci returns an error.
bootout1 <- boot(test1, w.mean, R=1000, stype="i")
(bootci1 <- boot.ci(bootout1, conf = 0.95, type = &qu...
2004 Sep 21
2
Bootstrap ICC estimate with nested data
...theta<-function(x,DATA){tmod<-lme(CINISMO~1,random=~1|IDGRUP,data=DATA[x,])
OUT<-as.numeric(VarCorr(tmod)[[1]])/(as.numeric(VarCorr(tmod)[[1]])+as.numer
ic(VarCorr(tmod)[[2]]))
return(OUT)}
Finally, I can run the bootstrap-t confidence limit function (or other
functions) as follows:
> bootout<-boott(1:nrow(TDAT),100,theta,data=TDAT)
This seems to work, but the estimates also seem strange. For intance, the
observed ICC value is larger than the 95% confidence intervals provided by
the bootstrap. It occurs to me that the results might be strange because
the sampling with replacement...