library(boot) data(aircondit) aircondit$hours t.test(aircondit$hours) samplemean=function (x,i) {mean(x[i])} b=boot(aircondit$hours,samplemean,1000) boot.ci(b)