### Show results of first three studies (first study has three treatment
### arms)
head(dat.linde2016, 5)
if (FALSE) {
### Load netmeta package
suppressPackageStartupMessages(library("netmeta"))
### Print odds ratios and confidence limits with two digits
oldset <- settings.meta(digits = 2)
### Define order of treatments in printouts and forest plots
trts <- c("SSRI",
"Face-to-face CBT", "Face-to-face interpsy", "Face-to-face PST",
"Face-to-face CBT + SSRI", "Face-to-face interpsy + SSRI",
"Face-to-face PST + SSRI",
"Face-to-face psychodyn", "Other face-to-face",
"TCA", "SNRI", "NRI", "Low-dose SARI", "NaSSa", "rMAO-A", "Ind drug",
"Hypericum",
"Remote CBT", "Self-help CBT", "No contact CBT",
"UC", "Placebo")
### Conduct random effects network meta-analysis
net <- netmeta(lnOR, selnOR, treat1, treat2, id,
data = dat.linde2016, reference.group = "placebo",
seq = trts, sm = "OR", fixed = FALSE)
### Network graph
netgraph(net, seq = "o", number = TRUE)
### Show results
net
forest(net, xlim = c(0.2, 50))
### Additive component network meta-analysis with placebo as inactive
### treatment
nc <- netcomb(net, inactive = "placebo")
nc
forest(nc, xlim = c(0.2, 50))
### Use previous settings
settings.meta(oldset)
}
Run the code above in your browser using DataLab