library(boot) aircondit$hours t.test(aircondit$hours) samplemeans=replicate(1000,mean(sample(aircondit$hours,12,replace=T))) quantile(samplemeans,c(.025,.975))