if (FALSE) {
# Only run example if R package 'Rgraphviz' from Bioconductor is available
#
if (requireNamespace("Rgraphviz", quietly = TRUE)) {
# Define order of treatments in depression dataset dat.linde2015
#
trts <- c("TCA", "SSRI", "SNRI", "NRI",
"Low-dose SARI", "NaSSa", "rMAO-A", "Hypericum", "Placebo")
# Outcome labels
#
outcomes <- c("Early response", "Early remission")
# (1) Early response
#
pw1 <- pairwise(treat = list(treatment1, treatment2, treatment3),
event = list(resp1, resp2, resp3),
n = list(n1, n2, n3),
studlab = id, data = dat.linde2015, sm = "OR")
#
net1 <- netmeta(pw1, common = FALSE,
seq = trts, ref = "Placebo", small.values = "undesirable")
# (2) Early remission
#
pw2 <- pairwise(treat = list(treatment1, treatment2, treatment3),
event = list(remi1, remi2, remi3),
n = list(n1, n2, n3),
studlab = id, data = dat.linde2015, sm = "OR")
#
net2 <- netmeta(pw2, common = FALSE,
seq = trts, ref = "Placebo", small.values = "undesirable")
# Partial order of treatment rankings
#
po <- netposet(netrank(net1), netrank(net2), outcomes = outcomes)
# Hasse diagram
#
hasse(po)
}
}
Run the code above in your browser using DataLab