f <- c(11, 15, 20, 27)
n <- c(60, 60, 60, 60)
x <- c(10, 20, 30, 40)
ci.slope.prop.bs(.05, f, n, x)
# Should return:
# Estimate SE z p LL UL
# 0.008688525 0.002566409 3.3855 0.00071 0.003658456 0.01371859
Run the code above in your browser using DataLab