f <- c(15, 24, 28, 23)
n <- c(50, 50, 50, 50)
ci.2x2.prop.bs(.05, f, n)
# Should return:
# Estimate SE z p LL UL
# AB: -0.27450980 0.13692496 -2.0048 0.04498 -0.54287780 -0.00614181
# A: -0.11764706 0.06846248 -1.7184 0.08572 -0.25183106 0.01653694
# B: -0.03921569 0.06846248 -0.5728 0.56678 -0.17339968 0.09496831
# A at b1: -0.25000000 0.09402223 -2.6589 0.00784 -0.43428019 -0.06571981
# A at b2: 0.01923077 0.09787658 0.1965 0.84423 -0.17260380 0.21106534
# B at a1: -0.17307692 0.09432431 -1.8349 0.06652 -0.35794917 0.01179533
# B at a2: 0.09615385 0.09758550 0.9853 0.32446 -0.09511021 0.28741790
Run the code above in your browser using DataLab