m1 <- c(53, 60, 53, 57)
m2 <- c(55, 62, 58, 61)
sd1 <- c(4.1, 4.2, 4.5, 4.0)
sd2 <- c(4.2, 4.7, 4.9, 4.8)
cor <- c(.72, .78, .81, .85)
n <- c(30, 50, 30, 70)
meta.ave.mean.ps(.05, m1, m2, sd1, sd2, cor, n, bystudy = TRUE)
# Should return:
# Estimate SE LL UL df
# Average -3.25 0.2340603 -3.713965 -2.7860352 107.657
# Study 1 -2.00 0.5672507 -3.160158 -0.8398421 29.000
# Study 2 -2.00 0.4227434 -2.849535 -1.1504653 49.000
# Study 3 -5.00 0.5335104 -6.091151 -3.9088487 29.000
# Study 4 -4.00 0.3023716 -4.603215 -3.3967852 69.000
Run the code above in your browser using DataLab