n <- c(583, 470, 546, 680)
rel <- c(.91, .89, .90, .89)
x1 <- c(1, 0, 0, 0)
X <- matrix(x1, 4, 1)
meta.lm.cronbach(.05, n, rel, 10, X)
# Should return:
# Estimate SE z p LL UL
# b0 -2.2408328 0.03675883 -60.960391 0.000 -2.3128788 -2.16878684
# b1 -0.1689006 0.07204625 -2.344336 0.019 -0.3101087 -0.02769259
Run the code above in your browser using DataLab