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