y2 <- c(21, 4, 9, 12, 35, 18, 10, 22, 24, 1, 6, 8, 13, 16, 19)
y1 <- c(67, 28, 30, 28, 52, 40, 25, 37, 44, 10, 14, 20, 28, 40, 51)
ci.ratio.mad.ps(.05, y1, y2)
# Should return:
# MAD1 MAD2 MAD1/MAD2 LL UL
# 12.71429 7.5 1.695238 1.109176 2.590961
Run the code above in your browser using DataLab