# bootstrap t CI for mean of Mouse data control group.

library(bootstrap)

B1 <- 1000
B2 <- 250
mean.reps <- rep(0, B1)
t.reps <- rep(0, B1)

for (i in 1:B1) {
  bsample <- sample(mouse.c, length(mouse.c), replace=T)
  mean.reps[i] <- mean(bsample)
  results <- bootstrap(bsample, B2, mean)
  t.reps[i] <- (mean(bsample) - mean(mouse.c))/sd(results$thetastar)
}
low <- mean(mouse.c) - quantile(t.reps, 0.95)*sd(mean.reps)
high <- mean(mouse.c) - quantile(t.reps, 0.05)*sd(mean.reps)
cat("bootstrap t CI: (", low, ",", high, ")\n")
hist(t.reps)