y1 <- c(32.1, 39.8, 26.3, 35.0, 43.1, 27.0, 40.9, 37.4, 34.0, 29.2)
y2 <- c(36.8, 44.0, 47.1, 42.7, 49.0, 39.6, 46.2, 31.6, 33.1, 48.4)
ci.median2(.05, y1, y2)
# Should return:
# Median1 Median2 Median1-Median2 SE LL UL
# 34.5 43.35 -8.85 4.494993 -17.66002 -0.0399751
Run the code above in your browser using DataLab