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.5222 0.00000 -0.35483065 -0.14021885
# 1 3 -0.1039604 0.04833562 -2.1508 0.03149 -0.21967489 0.01175409
# 2 3 0.1435644 0.04358401 3.2940 0.00099 0.03922511 0.24790360
Run the code above in your browser using DataLab