n <- c(65, 30, 29, 45, 50)
cor <- c(.87, .92, .85, .90, .88)
m1 <- c(20.1, 20.5, 19.3, 21.5, 19.4)
m2 <- c(10.4, 10.2, 8.5, 10.3, 7.8)
sd1 <- c(9.3, 9.9, 10.1, 10.5, 9.8)
sd2 <- c(7.8, 8.0, 8.4, 8.1, 8.7)
x1 <- c(2, 3, 3, 4, 4)
X <- matrix(x1, 5, 1)
meta.lm.mean.ps(.05, m1, m2, sd1, sd2, cor, n, X)
# Should return:
# Estimate SE t p LL UL df
# b0 8.00 1.2491990 6.404104 0.000 5.5378833 10.462117 217
# b1 0.85 0.3796019 2.239188 0.026 0.1018213 1.598179 217
Run the code above in your browser using DataLab