n <- c(120, 170, 150, 135)
rel <- c(.89, .87, .73, .71)
group <- c(1, 1, 2, 2)
r <- 10
meta.sub.cronbach(.05, n, rel, r, group)
# Should return:
# Estimate SE LL UL
# Set A: 0.88 0.01068845 0.8581268 0.8999386
# Set B: 0.72 0.02515130 0.6684484 0.7668524
# Set A - Set B: 0.16 0.02732821 0.1082933 0.2152731
Run the code above in your browser using DataLab