Hi All, I would like to do double boostrapping to estimate 95% CI coverage. So, I can only need to estimate 95% confidence interval from each bootstrapped sample. Since we don't have a closed form of 95% CI, in order to get 95% CI for each sample, we need to use bootstrapping. For outer bootstrapping, I used boot function in library(boot). How can I add double bootstrap inside the function. Is there any good way? or reference to see? Thanks, Becky [[alternative HTML version deleted]]