# NOT RUN {
### meta-analysis of the log risk ratios using a random-effects model
res1 <- rma(measure="RR", ai=tpos, bi=tneg, ci=cpos, di=cneg,
data=dat.bcg, method="ML")
### mixed-effects model with two moderators (latitude and publication year)
res2 <- rma(measure="RR", ai=tpos, bi=tneg, ci=cpos, di=cneg, mods = ~ ablat + year,
data=dat.bcg, method="ML")
fitstats(res1, res2)
logLik(res1)
logLik(res2)
deviance(res1)
deviance(res2)
AIC(res1, res2)
AIC(res1, res2, correct=TRUE)
BIC(res1, res2)
# }
Run the code above in your browser using DataCamp Workspace