#Compute Clopper-Pearson 95% CI (n=25) for input y cpci=function(y){ fl=qf(.025,2*y,2*(25-y+1)) lci=1/(1+(25-y+1)/(y*fl)) fu=qf(.975,2*(y+1),2*(25-y)) uci=1/(1+(25-y)/((y+1)*fu)) ci=c(lci,uci) return(ci) } # Store confidence bounds for outcomes y=0 through y=25 cimat=matrix(cpci(0:25),ncol=2) cimat[1,1]=0 cimat[26,2]=1 # Single p0 p0=.4 cicov=(p0>cimat[,1])&(p0cimat[,1])&(p0[i]