Displaying 1 result from an estimated 1 matches for "vtrue".
Did you mean:
true
2008 Sep 26
0
Confidence interval for binomial variance
...rights reserved
#test of CI for binomial sample variance
#created: 09.26.08 by r.a. labudde
#changes:
nReal <- 1000
for (POD in c(0.05, 0.10, 0.15, 0.20, 0.25, 0.30, 0.35, 0.50)) {
cat('\nPOD: ', sprintf('%8.4f', POD), '\n')
for (nRepl in c(6, 12, 20, 50)) {
vtrue<- POD*(1-POD)/nRepl
pcover<- 0
for (iReal in 1:nReal) {
x<- rbinom(1, nRepl, POD)
vCI<- binomVarCI(nRepl, x)
#vest<- (x/nRepl)*(1 - x/nRepl)/(nRepl-1)
if (vtrue >= vCI[1] & vtrue<= vCI[2]) pcover<- pcover + 1
} #iReal
pcove...