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))}