y <- c(30, 20, 15, 10, 10, 60, 20, 25, 20, 30, 10, 5, 50, 40,
20, 10, 0, 20, 50)
ci.cod(.05, y)
# Should return:
# Estimate SE LL UL
# 0.5921053 0.1814708 0.3813259 1.092679
Run the code above in your browser using DataLab