> r1=c(3,0)
> r2=c(27,7)
> n1=389
> n2=385
> cdf1=beta.diff(r1[1],r2[1],n1,n2)$cdf
> cdf1(0)
[1] 0.9986711
> cdf2=beta.diff(r1[2],r2[2],n1,n2)$cdf
> cdf2(0)
[1] 0.9714325