nodesplit_plot {rnmamod} | R Documentation |
End-user-ready results for the node-splitting approach
Description
nodesplit_plot
hosts a toolkit of functions that
facilitates the comparison of the consistency model
(via run_model
) with the node-splitting approach
(via run_nodesplit
) regarding the posterior summaries of the
direct and indirect effects and inconsistency factor of the split
nodes, the between-trial standard deviation and model assessment
parameters (Spiegelhalter et al., 2002) after each split node in the
network.
Usage
nodesplit_plot(full, node, drug_names, save_xls)
Arguments
full |
|
node |
An object of S3 class |
drug_names |
A vector of labels with the name of the interventions in
the order they appear in the argument |
save_xls |
Logical to indicate whether to export the tabulated results
to an 'xlsx' file (via the |
Details
intervalplot_inconsistency_factor
includes as many interval
plots as the number of split nodes in the network. Each interval plot
illustrates the posterior median and 95% credible interval of the direct and
indirect effect of the split nodes and the corresponding inconsistency
factor.
The line that corresponds to the inconsistency factor is highlighted with
green, when it does not cross the vertical line of no difference (between
the direct and indirect effect), and red otherwise. If there are more than
30 split nodes, the function presents the interval plots on split nodes
with conclusive inconsistency factor (green intervals) or those with
an opposite sign in the direct and indirect effects.
intervalplot_tau
is an interval plot on the median and 95% credible
interval of tau after each split node. The lines that correspond to
the split nodes are sorted in ascending order of the deviance information
criterion (DIC) which appears at the top of each line.
The estimated median and 95% credible intervals of tau under the
consistency model appear in the interval plot as a solid and two dotted
parallel blue lines, respectively. The different levels of heterogeneity
appear as green, yellow, orange, and red rectangulars to indicate a low,
reasonable, fairly high, and fairly extreme heterogeneity, respectively,
following the classification of Spiegelhalter et al. (2004).
When a fixed-effect model has been performed, nodesplit_plot
does
not return the intervalplot_tau
.
table_model_assessment
also includes the column
DIC-based better fit that indicates the preferred model in terms of
parsimony for each split node. Therefore, the DIC of the model after each
split node is compared with the DIC of the consistency model
(Dias et al., 2010). If the difference in DIC exceeds 5, the consistency
model is preferred; if the difference in DIC is less than -5, the model
after the split node is preferred; otherwise, there is little to choose
between the compared models.
For a binary outcome, when measure
is "RR" (relative risk) or "RD"
(risk difference) in run_model
, nodesplit_plot
currently presents the results in the odds ratio scale. This is because,
the odds ratio is used as the 'best-case' effect measure in
run_model
. Then, relative risk, and risk difference are
obtained as a function of the odds ratio and the selected baseline risk
(See 'Details' in run_model
).
The split nodes have been automatically selected via the
mtc.nodesplit.comparisons
function of the R-package
gemtc.
See 'Details' in run_nodesplit
.
Furthermore, nodesplit_plot
exports both data-frames to separate
'xlsx' files (via the write_xlsx
function
of the R-package
writexl) to the working
directory of the user.
nodesplit_plot
can be used only for a network of interventions and
when there is at least one split node. Otherwise, the execution of the
function will be stopped and an error message will be printed on the R
console.
Value
nodesplit_plot
returns the following list of elements:
table_effect_size |
A data-frame with the posterior median, posterior standard deviation and 95% credible interval of the direct and indirect effect and the inconsistency factor of each split node. |
table_model_assessment |
A data-frame with the model assessment parameters (DIC, posterior mean of total residual deviance, and number of effective parameters), the posterior median, posterior standard deviation and 95% credible interval of tau under the consistency model and after each split node. See 'Details'. |
intervalplot_inconsistency_factor |
A panel of interval plots on the direct and indirect effect of the split nodes and the corresponding inconsistency factor. See 'Details'. |
intervalplot_tau |
An interval plot on tau after each split node. See 'Details'. |
Author(s)
Loukia M. Spineli
References
Dias S, Welton NJ, Caldwell DM, Ades AE. Checking consistency in mixed treatment comparison meta-analysis. Stat Med 2010;29(7-8):932–44. doi: 10.1002/sim.3767
Spiegelhalter DJ, Abrams KR, Myles JP. Bayesian approaches to clinical trials and health-care evaluation. John Wiley and Sons, Chichester, 2004.
Spiegelhalter DJ, Best NG, Carlin BP, van der Linde A. Bayesian measures of model complexity and fit. J R Stat Soc B 2002;64(4):583–616. doi: 10.1111/1467-9868.00353
See Also
mtc.nodesplit.comparisons
,
run_model
, run_nodesplit
,
write_xlsx
Examples
data("nma.baker2009")
# Read results from 'run_model' (using the default arguments)
res <- readRDS(system.file('extdata/res_baker.rds', package = 'rnmamod'))
# Read results from 'run_nodesplit' (using the default arguments)
node <- readRDS(system.file('extdata/node_baker.rds', package = 'rnmamod'))
# The names of the interventions in the order they appear in the dataset
interv_names <- c("placebo", "budesonide", "budesonide plus formoterol",
"fluticasone", "fluticasone plus salmeterol",
"formoterol", "salmeterol", "tiotropium")
# Plot the results from both models
nodesplit_plot(full = res,
node = node,
drug_names = interv_names)