oldset <- settings.meta("RevMan5")
data(caffeine)
head(caffeine)
ma1 <- metabin(h.caf, n.caf, h.decaf, n.decaf, sm = "OR",
data = caffeine, studlab = paste(study, year))
# Add risk of bias assessment to meta-analysis
ma2 <- rob(D1, D2, D3, D4, D5, overall = rob, data = ma1, tool = "rob2")
# Print risk of bias assessment
rob(ma2)
# Forest plot with risk of bias assessment
forest(ma2)
settings.meta(oldset)
Run the code above in your browser using DataLab