# NOT RUN {
## Smoking cessation NMA
# Set up network of smoking cessation data
head(smoking)
smk_net <- set_agd_arm(smoking,
study = studyn,
trt = trtc,
r = r,
n = n,
trt_ref = "No intervention")
# Print details
smk_net
# }
# NOT RUN {
# Fitting a random effects model
smk_fit_RE <- nma(smk_net,
trt_effects = "random",
prior_intercept = normal(scale = 100),
prior_trt = normal(scale = 100),
prior_het = normal(scale = 5))
smk_fit_RE
# }
# NOT RUN {
# }
# NOT RUN {
# Plot prior vs. posterior, by default all parameters are plotted
plot_prior_posterior(smk_fit_RE)
# Plot prior vs. posterior for heterogeneity SD only
plot_prior_posterior(smk_fit_RE, prior = "het")
# Customise plot
plot_prior_posterior(smk_fit_RE, prior = "het",
prior_args = list(colour = "darkred", size = 2),
post_args = list(alpha = 0.6))
# }
# NOT RUN {
# }
Run the code above in your browser using DataCamp Workspace