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