y <- c(21, 4, 9, 12, 35, 18, 10, 22, 24, 1, 6, 8, 13, 16, 19)
x <- c(67, 28, 30, 28, 52, 40, 25, 37, 44, 10, 14, 20, 28, 40, 51)
ci.spear(.05, y, x)
# Should return:
# Estimate SE LL UL
# 0.8699639 0.08241326 0.5840951 0.9638297
Run the code above in your browser using DataLab