EM <- 900*1.06/2.06 EF <- 900*1/2.06 X <- (480-EM)^2/EM+(420-EF)^2/EF 1- pchisq(X,1) print(res <- binom.test(480,900,480/900)) res$conf.int/(1-res$conf.int)