### load BCG vaccine data
data(dat.bcg)
### meta-analysis of the log relative risks using a random-effects model
res <- rma(ai=tpos, bi=tneg, ci=cpos, di=cneg, data=dat.bcg,
measure="RR", method="REML")
plot(res, qqplot=TRUE)
### mixed-effects model with two moderators (absolute latitude and publication year)
res <- rma(ai=tpos, bi=tneg, ci=cpos, di=cneg, mods = ~ ablat + year,
data=dat.bcg, measure="RR", method="REML")
plot(res, qqplot=TRUE)
Run the code above in your browser using DataLab