search for: bcaresult

...i Misses cnt = 1000; # reps x = rnorm(10000) # create population m = mean(x) for (i in 1:cnt) { s = sample(x,10,replace=F) # sample population tresults = t.test(s) attach(tresults) C[i] = CIr(conf.int[1],conf.int[2]) if (CImiss(m,conf.int[1],conf.int[2])) Ccov = Ccov + 1 detach(tresults) bcaresults <- bcanon(s,5000,mean,alpha=c(.025,.975)) attach(bcaresults) B[i] = CIr(confpoints[1,2],confpoints[2,2]) if (CImiss(m,confpoints[1,2],confpoints[2,2])) Bcov = Bcov + 1 detach(bcaresults) } print(summary (C)) print(Ccov/cnt) print(summary (B)) print(Bcov/cnt) [[alternative HTML versi...