Displaying 1 result from an estimated 1 matches for "npod".
Did you mean:
ipod
2008 Sep 26
0
Confidence interval for binomial variance
...program to measure coverage:
#09.26.08 03.10 tbinomVarCI.r
#copyright 2008 by Robert A LaBudde, all 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)...