Joint Modeling of Multivariate Longitudinal and Survival Outcomes with the R package INLAjoint

Denis Rustand,Janet van Niekerk,Elias Teixeira Krainski,Håvard Rue
2024-02-13
Abstract:This paper introduces the R package INLAjoint, designed as a toolbox for fitting a diverse range of regression models addressing both longitudinal and survival outcomes. INLAjoint relies on the computational efficiency of the integrated nested Laplace approximations methodology, an efficient alternative to Markov chain Monte Carlo for Bayesian inference, ensuring both speed and accuracy in parameter estimation and uncertainty quantification. The package facilitates the construction of complex joint models by treating individual regression models as building blocks, which can be assembled to address specific research questions. Joint models are relevant in biomedical studies where the collection of longitudinal markers alongside censored survival times is common. They have gained significant interest in recent literature, demonstrating the ability to rectify biases present in separate modeling approaches such as informative censoring by a survival event or confusion bias due to population heterogeneity. We provide a comprehensive overview of the joint modeling framework embedded in INLAjoint with illustrative examples. Through these examples, we demonstrate the practical utility of INLAjoint in handling complex data scenarios encountered in biomedical research.
Computation,Methodology
What problem does this paper attempt to address?
The paper aims to address the issue of jointly modeling multivariate longitudinal data and survival outcomes in biomedical research. Specifically: 1. **Method Introduction**: The paper introduces the R package `INLAjoint`, which is a toolbox for fitting various regression models that can simultaneously handle longitudinal data and survival outcomes. `INLAjoint` is based on the Integrated Nested Laplace Approximations (INLA) method, an efficient Bayesian inference alternative that offers faster speed and higher accuracy compared to the Markov Chain Monte Carlo (MCMC) method. 2. **Problem Background**: In biomedical research, longitudinal markers and survival times are often associated. Modeling them separately may introduce information bias or confounding bias, thus requiring joint models to comprehensively consider both factors. For example, the CD4 lymphocyte count and survival time in AIDS patients, and the dynamic changes of prostate-specific antigen (PSA) and cancer recurrence in prostate cancer patients. 3. **Limitations of Existing Methods**: Traditional two-stage methods and other software such as SAS's NLMIXED, Stata's merlin command, and several R packages (e.g., JM, joineRML, and frailtypack) can handle simple joint models but face high computational complexity and convergence difficulties when dealing with multivariate longitudinal and survival outcomes. 4. **Advantages of INLAjoint**: By treating each sub-model as a modular building block, `INLAjoint` can easily assemble complex joint models and leverage the efficient computational capabilities of the INLA method. Additionally, it provides user-friendly syntax and rich plotting and prediction tools, making it more convenient for users to apply the INLA method for joint modeling. 5. **Application Examples**: The paper demonstrates the practicality of `INLAjoint` in handling complex data scenarios through specific examples, particularly in biomedical research applications. In summary, the paper primarily addresses how to provide an efficient, accurate, and user-friendly joint modeling method for handling multivariate longitudinal data and survival outcomes.