m <- c(33.5, 37.9, 38.0, 44.1)
sd <- c(3.84, 3.84, 3.65, 4.98)
q <- c(.5, .5, -.5, -.5)
ci.lc.stdmean.ws(.05, m, sd, .672, 20, q)
# Should return:
# Estimate adj Estimate SE LL UL
# Unweighted standardizer: -1.3013 -1.2666 0.31479 -1.9182 -0.6843
# Level 1 standardizer: -1.3932 -1.3375 0.36618 -2.1109 -0.6755
Run the code above in your browser using DataLab