oldsets <- settings.meta(digits = 2, digits.tau2 = 2, digits.tau = 2)
data(smokingcessation)
# Transform data from arm-based format to contrast-based format
#
pw1 <- pairwise(list(treat1, treat2, treat3),
event = list(event1, event2, event3), n = list(n1, n2, n3),
data = smokingcessation, sm = "OR")
# Conduct random effects network meta-analysis
#
nma1 <- netmeta(pw1, common = FALSE)
# Calculate and print concise results for all pairwise
# meta-analyses
#
np1 <- netpairwise(nma1)
np1
print(np1, details.method = FALSE)
# \donttest{
# Random effects model
#
pw2 <- pairwise(studlab = study, treat = treatment,
n = n, mean = mean, sd = sd, data = Senn2013,
varnames = c("MD", "seMD"))
nma2 <- netmeta(pw2, common = FALSE, reference = "plac")
# Calculate and print concise results for all pairwise
# meta-analyses
#
np2 <- netpairwise(nma2)
np2
print(np2, details.method = FALSE)
forest(np2)
# Print detailed information for each pairwise comparison
#
np3 <- netpairwise(nma2, separate = TRUE)
forest(np3)
funnel(np3)
radial(np3)
funnel(np3, k.min = 1)
# }
settings.meta(oldsets)
Run the code above in your browser using DataLab