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.theil(.05, y, x)
# Should return:
# Estimate SE LL UL
# 0.5 0.1085927 0.3243243 0.75
Run the code above in your browser using DataLab