m <- c(12.86, 24.57, 36.29, 53.21)
sd <- c(13.185, 12.995, 14.773, 15.145)
n <- c(20, 20, 20, 20)
test.mono.mean.bs(.05, m, sd, n)
# Should return:
# Estimate SE LL UL
# 1 2 -11.71 4.139530 -22.07803 -1.3419744
# 2 3 -11.72 4.399497 -22.74731 -0.6926939
# 3 4 -16.92 4.730817 -28.76921 -5.0707936
Run the code above in your browser using DataLab