m <- c(12.86, 24.57, 36.29, 53.21)
se <- c(2.85, 2.99, 3.73, 3.88)
test.mono.median.bs(.05, m, se)
# Should return:
# Estimate SE LL UL
# 1 2 -11.71 4.130690 -21.59879 -1.8212115
# 2 3 -11.72 4.780481 -23.16438 -0.2756247
# 3 4 -16.92 5.382128 -29.80471 -4.0352947
Run the code above in your browser using DataLab