This network meta-analysis is an update to Fu et al, 2014, J Cardiovasc Med (Hagerstown). 2014 Dec; 15(12): 873-879.
forest-#.png files, are in this folderforest-#.png files, are in this folderforest-#.png files, are in this folderforest-#.png files, are in this folderUse the gemtc package. Check the link for references.
library(openxlsx)
library(data.table)
library(gemtc)
library(xtable)
Set effect to be used for the linearModel parameter in the mtc.model function.
effect <- "fixed"
Set sampler parameters.
nAdapt <- 5000 # Burn-in
# nChain <- 4 # gemtc defaults to 4 MCMC chains; no need to set
nIter <- 25000
thin <- 10
message(sprintf("Each MCMC chain will have %d samples.\nInferences will be based on a total of %d samples.",
nIter / thin,
4 * (nIter / thin)))
## Each MCMC chain will have 2500 samples.
## Inferences will be based on a total of 10000 samples.
Software version information.
sessionInfo()
## R version 3.2.2 (2015-08-14)
## Platform: x86_64-w64-mingw32/x64 (64-bit)
## Running under: Windows 8 x64 (build 9200)
##
## locale:
## [1] LC_COLLATE=English_United States.1252
## [2] LC_CTYPE=English_United States.1252
## [3] LC_MONETARY=English_United States.1252
## [4] LC_NUMERIC=C
## [5] LC_TIME=English_United States.1252
##
## attached base packages:
## [1] stats graphics grDevices utils datasets methods base
##
## other attached packages:
## [1] xtable_1.7-4 gemtc_0.7-1 coda_0.17-1 data.table_1.9.6
## [5] openxlsx_3.0.0
##
## loaded via a namespace (and not attached):
## [1] Rcpp_0.11.6 lattice_0.20-33 digest_0.6.8 truncnorm_1.0-7
## [5] chron_2.3-47 grid_3.2.2 plyr_1.8.3 meta_4.3-2
## [9] formatR_1.2 magrittr_1.5 evaluate_0.8 stringi_0.4-1
## [13] rjags_4-5 rmarkdown_0.8 tools_3.2.2 stringr_1.0.0
## [17] igraph_1.0.1 yaml_2.1.13 htmltools_0.2.6 knitr_1.11