library(boot) aircondit$hours t.test(aircondit$hours) b=sample(aircondit$hours,12000,replace=T) y=rep(1:1000,12) samplemeans=tapply(b,y,mean) quantile(samplemeans,c(.025,.975))