# NOT RUN {
# Network meta-analysis of count mortality statistics
#
data(Woods2010)
p0 <- pairwise(treatment, event = r, n = N,
studlab = author, data = Woods2010, sm = "OR")
net0 <- netmeta(p0)
cilayout(bracket = "(", separator = " - ")
oldopts <- options(width = 100)
# League table for fixed effect model
#
netleague(net0, digits = 2)
# League table for fixed effect and random effects model
#
netleague(net0, comb.random = TRUE, digits = 2)
# Change order of treatments according to treatment ranking
#
netleague(net0, comb.random = TRUE, digits = 2,
seq = netrank(net0))
#
print(netrank(net0), comb.random = TRUE)
# Use depression dataset
#
data(Linde2015)
cilayout()
#
# Define order of treatments
#
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
#
p1 <- pairwise(treat = list(treatment1, treatment2, treatment3),
event = list(resp1, resp2, resp3),
n = list(n1, n2, n3),
studlab = id, data = Linde2015, sm = "OR")
#
net1 <- netmeta(p1,
comb.fixed = FALSE, comb.random = TRUE,
seq = trts, ref = "Placebo")
#
# (2) Early remission
#
p2 <- pairwise(treat = list(treatment1, treatment2, treatment3),
event = list(remi1, remi2, remi3),
n = list(n1, n2, n3),
studlab = id, data = Linde2015, sm = "OR")
#
net2 <- netmeta(p2,
comb.fixed = FALSE, comb.random = TRUE,
seq = trts, ref = "Placebo")
options(width = 200)
netleague(net1, digits = 2)
netleague(net1, digits = 2, ci = FALSE)
netleague(net2, digits = 2, ci = FALSE)
netleague(net1, net2, digits = 2, ci = FALSE)
netleague(net1, net2, seq = netrank(net1, small = "bad"), ci = FALSE)
netleague(net1, net2, seq = netrank(net2, small = "bad"), ci = FALSE)
print(netrank(net1, small = "bad"), comb.random = TRUE)
print(netrank(net2, small = "bad"), comb.random = TRUE)
options(oldopts)
# }
# NOT RUN {
# Generate a partial order of treatment rankings
#
np <- netposet(net1, net2, outcomes = outcomes, small.values = rep("bad",2))
hasse(np)
plot(np)
# }
Run the code above in your browser using DataLab