y1 <- c(21.1, 4.9, 9.2, 12.4, 35.8, 18.1, 10.7, 22.9, 24.0, 1.2, 6.1, 8.3, 13.1, 16.2)
y2 <- c(67.0, 28.1, 30.9, 28.6, 52.0, 40.8, 25.8, 37.4, 44.9, 10.3, 14.9, 20.2, 28.8, 40.6)
ci.median.ps(.05, y1, y2)
# Should return:
# Median1 Median2 Median1-Median2 SE LL UL
# 12.75 29.85 -17.1 3.704248 -24.36019 -9.839807
# SE1 SE2 COV
# 3.379695 4.968956 11.1957
Run the code above in your browser using DataLab