# percentile CI for mean of Mouse data control group.

library(bootstrap)

B <- 2000

results <- bootstrap(mouse.c, B, mean)
low <- quantile(results$thetastar, 0.05)
high <- quantile(results$thetastar, 0.95)
cat("percentile CI: (", low, ",", high, ")\n")
hist(results$thetastar)