On the predictive performance of two Bayesian joint models: a simulation study
- 2 November 2022
- journal article
- research article
- Published by Informa UK Limited in Communications in Statistics - Simulation and Computation
- Vol. 51 (11), 6388-6397
- https://doi.org/10.1080/03610918.2020.1804579
Abstract
Joint modeling of longitudinal and survival data is becoming increasingly popular. It can properly handle multiple issues commonly encountered in longitudinal studies, such as endogenous time-dependent covariates and informative missingness. There are several statistical packages for conducting Bayesian joint modeling analysis, two of which are JMbayes, which has been applied in biomedical research of various fields, and rstanarm, which has been developed recently. Both packages are very flexible in specifying different association structures, and are capable of handling multivariate outcomes. However, no studies have ever been conducted to compare their performance. In this study, we conducted simulation studies to compare the performance of the two packages, with a focus on their prediction. We found that rstanarm often had better predictive performance than JMbayes, but its computation was more intensive and increased dramatically with larger sample sizes. In contrast, JMbayes was fast in model fitting and prediction, and in some cases, its performance could be improved by using larger sample sizes or longer periods of longitudinal data. Model mis-specification seemed to have a greater influence on the predictive performance of rstanarm than JMbayes.Keywords
Funding Information
- national natural science foundation of Hunan (2018JJ2265)
- National Natural Science Foundation of China (81771493)
- NIH/NIA (R01AG036042)
- the Illinois Department of Public Health
This publication has 30 references indexed in Scilit:
- A Bayesian semiparametric multivariate joint model for multiple longitudinal outcomes and a time‐to‐eventStatistics in Medicine, 2011
- Basic Concepts and Methods for Joint Models of Longitudinal and Survival DataJournal of Clinical Oncology, 2010
- Boosting for high-dimensional time-to-event data with competing risksBioinformatics, 2009
- A Joint Model for Longitudinal Measurements and Survival Data in the Presence of Multiple Failure TypesBiometrics, 2007
- A Bayesian model for joint segmentation and registrationNeuroImage, 2006
- Jointly Modeling Longitudinal and Event Time Data With Application to Acquired Immunodeficiency SyndromeJournal of the American Statistical Association, 2001
- Joint modelling of longitudinal measurements and event time dataBiostatistics, 2000
- TIME-DEPENDENT COVARIATES IN THE COX PROPORTIONAL-HAZARDS REGRESSION MODELAnnual Review of Public Health, 1999
- Hybrid Monte CarloPhysics Letters B, 1987
- Monte Carlo sampling methods using Markov chains and their applicationsBiometrika, 1970