f <- c(26, 24, 38)
n <- c(60, 60, 60)
v <- c(-.5, -.5, 1)
ci.lc.prop.bs(.05, f, n, v)
# Should return:
# Estimate SE z p LL UL
# 0.2119565 0.07602892 2.787841 0.005306059 0.06294259 0.3609705
Run the code above in your browser using DataLab