binom.ci<-function(r,n,alpha=.05){
#Exact binom conf int from Armitage and Berry '94 (pg 121)
#if (r==0) lower<-0
lower<-ifelse(r==0,0,r/(r+(n-r+1)*qf((1-alpha/2),(2*(n-r+1)),(2*r))))
#if (r==n) upper<-1
upper<-ifelse(r==n,1,(r+1)/(r+1+(n-r)/qf((1-alpha/2),(2*(r+1)),(2*(n-r)))))
return(list(lower=lower,upper=upper))}