Are there any functions in R for running bootstraps with clustered (as opposed to stratified) data? I can't seem to find anything obvious in boot or Bootstrap, though I imagine boot can be manipulated to resample from clusters. Is that what people use? I do see some cluster bootstrap resampling in glsD, but I need non-parameteric resampling & the capacity to run a sem model. Thanks in advance, Peter
Looks like I may have found a function that addresses my needs. Bootcov in Design handles bootstrapping from clustered data and will save the coefficients. I'm not entirely sure it handles clusters the way I'd like, but I'm going through the code. If it doesn't, it looks easily re-writeable. As far as I can tell, boot in package boot would do clusters only if the estimation function passed to it pastes together data based on the clusters boot selects.