Displaying 1 result from an estimated 1 matches for "binwilson".
Did you mean:
benwilson
2008 Sep 26
0
Confidence interval for binomial variance
...ert A LaBudde, all rights reserved
#CI for binomial sample variance
#created: 09.26.08 by r.a. labudde
#changes:
require('binGroup')
binomVarCI<- function (n, x, conf=0.95) {
p<- x/n #proportion
if (p<0.25 | p>0.75 | x==0 | x==n) { #use proportion-based CI
pCI<- binWilson(n, x, conf.level=conf) #CI for proportion
vCI<- sort(c(pCI[1]*(1-pCI[1])/(n-1), pCI[2]*(1-pCI[2])/(n-1)))
} else { #use chi-square-based CI
phiL<- qchisq(0.025, n-1)/(n-1)
phiU<- qchisq(0.975, n-1)/(n-1)
#vest<- p*(1-p)/(n-1)) #variance estimate
vCI<- c...