data("antidepressants")
#
pw1 <- pairwise(studlab = studyid, treat = drug_name,
n = ntotal, event = responders,
data = antidepressants, sm = "OR")
# Use subset to reduce runtime
pw0 <- subset(pw1, studyid < 60)
#
net0 <- netmeta(pw0, reference.group = "tra")
ranks0 <- tcc(net0, swd = 1.20, small.values = "undesirable")
# Comparison other drugs vs trazodone
forest(ranks0,
label.left = "Favours trazodone",
label.right = "Favours other drug")
# Comparison escitalopram vs other drugs
forest(ranks0, reference.group = "esc", baseline = FALSE,
label.left = "Favours other drug",
label.right = "Favours escitalopram")
if (FALSE) {
# Store a PDF file in the current working directory showing all results
# (this is the default, i.e., if argument 'reference.group' is missing)
forest(ranks0, baseline = FALSE, reference.group = trts,
file = "forest_tcc_antidepressants.pdf")
# Run analysis with full data set
net1 <- netmeta(pw1, reference.group = "tra")
ranks1 <- tcc(net1, swd = 1.20, small.values = "undesirable")
# Comparison other drugs vs trazodone
forest(ranks1,
label.left = "Favours trazodone",
label.right = "Favours other drug")
# Comparison escitalopram vs other drugs
forest(ranks1, reference.group = "esc", baseline = FALSE,
label.left = "Favours other drug",
label.right = "Favours escitalopram")
}
Run the code above in your browser using DataLab