library(bootstrap) m <- mean(mouse.c) s <- sd(mouse.c) low <- m - qt(0.95, df=length(mouse.c)-1)*s/sqrt(length(mouse.c)) high <- m - qt(0.05, df=length(mouse.c)-1)*s/sqrt(length(mouse.c)) cat("normal theory CI: (", low, ",", high, ")\n")