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.301263 -1.266557 0.3147937 -1.918248 -0.6842788
# Level 1 standardizer: -1.393229 -1.337500 0.3661824 -2.110934 -0.6755248
Run the code above in your browser using DataLab