### load BCG vaccine data
data(dat.bcg)
### meta-analysis of the log risk rates using a random-effects model
res <- rma(ai=tpos, bi=tneg, ci=cpos, di=cneg,
measure="RR", data=dat.bcg, method="DL")
res
### average risk rate with 95\% CI
predict(res, transf=exp)
### mixed-effects model with absolute latitude as a moderator
res <- rma(ai=tpos, bi=tneg, ci=cpos, di=cneg,
mods=ablat, measure="RR", data=dat.bcg, method="DL")
res
### predicted average risk rates for 10-60 degrees absolute latitude
predict(res, newmods=c(10, 20, 30, 40, 50, 60), transf=exp)
Run the code above in your browser using DataLab