Displaying 1 result from an estimated 1 matches for "bcaresults".
Did you mean:
pcaresults
2003 Oct 08
0
Bootstrap Question
...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 versio...