f <- c(26, 24, 38)
n <- c(60, 60, 60)
v1 <- c(-.5, -.5, 1)
ci.lc.prop.bs(.05/2, f, n, v1)
# Should return:
# Estimate SE z p LL UL
# 0.2119565 0.07602892 2.7878 0.00531 0.0415451 0.3823679
v2 <- c(1, -1, 0)
ci.lc.prop.bs(.05/2, f, n, v2)
# Should return:
# Estimate SE z p LL UL
# 0.03225806 0.08857951 0.36417 0.71573 -0.1662843 0.2308004
Run the code above in your browser using DataLab