f <- c(111, 161, 132)
n <- c(200, 200, 200)
ci.pairs.prop.bs(.05, f, n)
# Should return:
# Estimate SE z p LL UL
# 1 2 -0.2475248 0.04482323 -5.522243 3.346989e-08 -0.35483065 -0.14021885
# 1 3 -0.1039604 0.04833562 -2.150803 3.149174e-02 -0.21967489 0.01175409
# 2 3 0.1435644 0.04358401 3.293968 9.878366e-04 0.03922511 0.24790360
Run the code above in your browser using DataLab