m <- c(33.5, 37.9, 38.0, 44.1)
sd <- c(3.49, 3.84, 3.65, 4.98)
n <- c(10, 10, 10, 10)
v <- c(.5, .5, -.5, -.5)
ci.lc.mean.scheffe(.05, m, sd, n, v)
# Should return:
# Estimate SE t p LL UL
# -5.35 1.275231 -4.1953 0.00228 -9.089451 -1.610549
Run the code above in your browser using DataLab