data(antidepressants)
#
ranks <- tcc(treat = drug_name, studlab = studyid,
event = responders, n = ntotal, data = antidepressants,
mcid = 1.25, sm = "OR", small.values = "undesirable")
# Visualize treatment choice criteron for treatment "escitalopram"
forest(ranks, treat = "escitalopram")
# Fit the model
fit <- mtrank(ranks)
# Visualize the results
forest(fit)
# Calculate pairwise probabilities
paired_pref(fit, type = "better",
treat1 = "bupropion", treat2 = "escitalopram")
# \donttest{
# Same results using pairwise object as input to tcc()
# (and running a network meta-analysis)
#
pw <- pairwise(treat = drug_name, studlab = studyid,
event = responders, n = ntotal, data = antidepressants,
sm = "OR")
#
netmeta(pw)
#
ranks2 <- tcc(pw, mcid = 1.25, small.values = "undesirable")
#
fit2 <- mtrank(ranks2)
#
paired_pref(fit2, type = "better",
treat1 = "bupropion", treat2 = "escitalopram")
# }
Run the code above in your browser using DataLab