y1 <- c(21, 4, 9, 12, 35, 18, 10, 22, 24, 1, 6, 8, 13, 16, 19)
y2 <- c(67, 28, 30, 28, 52, 40, 25, 37, 44, 10, 14, 20, 28, 40, 51)
ci.ratio.median.ps(.05, y1, y2)
# Should return:
# Median1 Median2 Median1/Median2 LL UL
# 13 30 0.4333333 0.3094838 0.6067451
Run the code above in your browser using DataLab