Abstract
This work aims at jointly modelling longitudinal and survival HIV data by considering the sharing of a set of parameters of interest. For the CD4 longitudinal stochastic process we propose a regression model where individual heterogeneity is allowed to vary in terms of the mean and the variance, relaxing the usual assumption of a common variance for the longitudinal residuals. Along, we will be considering a hazard regression model to analyse the time between HIV/AIDS diagnostic and death. For introducing enough flexibility in the structure linking the longitudinal and survival processes, we consider time-varying coefficients. That is achieved using Penalized Splines and allows the relationship to vary in time. The CD4 residuals standard deviation is considered as a covariate in the hazard model, thus enabling to study the effect of the CD4 counts’ stability on the survival. The proposed framework surpasses the performance of the most “traditional” joint models, which generally consider a common variance and a time-invariant link.
Similar content being viewed by others
References
Baghfalaki T, Ganjali M (2015) A Bayesian approach for joint modeling of skew-normal longitudinal measurements and time to event data. REVSTAT-Stat J 13(2):169–191
Brezger A, Lang S (2006) Generalized structured additive regression based on Bayesian p-splines. Comput Stat Data Anal 50:967–991
Carlin B, Louis T (2001) Bayes and empirical Bayes methods for data analysis, 2nd edn. Chapman & Hall, London
Chen J, Huang Y (2015) A Bayesian mixture of semiparametric mixed-effects joint models for skewed-longitudinal and time-to-event data. Stat Med 34(20):2820–2843
Christensen R, Johnson W, Branscum A, Hanson TE (2011) Bayesian ideas and data analysis: an introduction for scientists and statisticians. Chapman & Hall/CRC Press, London
Costa MJ, Shaw JEH (2009) Parametrization and penalties in spline models with an application to survival analysis. Comput Stat Data Anal 53(3):657–670
Cox D (1972) Regression models and life tables. J R Stat Soc B 34:187–220
Eilers PHC, Marx BD (2010) Splines, knots, and penalties. WIREs Comput Stat 2(6):637–653
Farcomeni A, Viviani S (2014) Longitudinal quantile regression in the presence of informative dropout through longitudinal-survival joint modeling. Stat Med 34(7):1199–1213. https://doi.org/10.1002/sim.6393
Faucett CL, Thomas DC (1996) Simultaneously modelling censored survival data and repeatedly measured covariates: a Gibbs sampling approach. Stat Med 15(15):1663–1685
Gao F, Miller J, Xiong C, Beiser J, Gordon M (2011) A joint-modeling approach to assess the impact of biomarker variability on the risk of developing clinical outcome. Stat Method Appl 20:83–100
Gelman A (2006) Prior distributions for variance parameters in hierarchical models (comment on article by browne and draper). Bayesian Anal 1(3):515–534
Gelman A, Hwang J, Vehtari A (2014) Understanding predictive information criteria for Bayesian models. Stat Comput 24(6):997–1016
Guo X, Carlin P (2004) Separate and joint modelling of longitudinal and event time data using standard computer packages. Am Stat 58(1):16–24
Hanson T, Branscum A, Johnson W (2011) Predictive comparison of joint longitudinal-survival modeling: a case study illustrating competing approaches. Lifetime Data Anal 17:3–28. https://doi.org/10.1007/s10985-010-9162-0
Henderson R, Diggle P, Dobson A (2000) Joint modelling of longitudinal measurements and event time data. Biostatistics 1(4):465–480
Hennerfeind A, Brezger A, Fahrmeir L (2006) Geoadditive survival models. JASA 101:1065–1075
Hofner B, Kneib T, Hartl W, Kuchenhoff H (2011) Building Cox-type structured hazard regression models with time-varying effects. Stat Model 11(1):3–24
Hsieh F, Tseng Y, Wang J (2006) Joint modeling of survival and longitudinal data: likelihood approach revisited. Biometrics 62:1037–1043
Huang Y, Hu XJ, Dagne GA (2014) Jointly modeling time-to-event and longitudinal data: a Bayesian approach. Stat Method Appl 23(1):95–121
Ibrahim J, Chen MH, Sinha D (2004) Bayesian methods for joint modeling of longitudinal and survival data with applications to cancer vaccine studies. Stat Sin 14(3):863–883
Jiang B, Elliott MR, Sammel MD, Wang N (2015) Joint modeling of cross-sectional health outcomes and longitudinal predictors via mixtures of means and variances. Biometrics 71(2):487–97
Kalbfleisch JD, Prentice RL (2002) The statistical analysis of failure time data. John Wiley, New York
Kneib T, Fahrmeir L (2007) A mixed model approach for geoadditive hazard regression. Scand J Stat 34(1):207–228
Lang S, Brezger A (2004) Bayesian p-splines. J Comput Graph Stat 13(1):183–212
Lange N, Carlin BP, Gelfand AE (1992) Hierarchical Bayes models for the progression of HIV infection using longitudinal CD4 T-cell numbers. J Am Stat Assoc 87(419):615–626
Lin DY, Fischl MA, Schoenfeld DA (1993) Evaluating the role of CD4-lymphocyte counts as surrogate endpoints in human immunodeficiency virus clinical trials. Stat Med 12(9):835–842
Lin X, Raz J, Harlow S (1997) Linear mixed models with heterogeneous within-cluster variances. Biometrics 53(3):910–923
Lunn D, Thomas A, Best N, Spiegelhalter D (2000) Winbugs—a bayesian modelling framework: concepts, structure, and extensibility. Stat Comput 10:325–337
Lunn D, Jackson C, Best N, Thomas A, Spiegelhalter D (2013) The BUGS book: a practical introduction to Bayesian analysis. Chapman & Hall/CRC, London
Lyles RH, Munõz A, Xu J, Taylor JMG, Chmiel JS (1999) Adjusting for measurement error to assess health effects of variability in biomarkers. Stat Med 18:1069–1086
Martins R, Silva GL, Andreozzi V (2016) Bayesian joint modeling of longitudinal and spatial survival aids data. Stat Med 35(19):3368–3384. https://doi.org/10.1002/sim.6937
Martins R, Silva G, Andreozzi V (2017) Joint analysis of longitudinal and survival aids data with a spatial fraction of long-term survivors: a Bayesian approach. Biom J 59(6):1166–1183
May MT, Gompels M, Delpech V, Porter K, Orkin C, Kegg S, Sabin C (2014) Impact on life expectancy of hiv-1 positive individuals of cd4+ cell count and viral load response to antiretroviral therapy. AIDS 28(8):1193–1202
McLain AC, Lum KJ, Sundaram R (2012) A joint mixed effects dispersion model for menstrual cycle length and time-to-pregnancy. Biometrics 68(2):648–656
Rizopoulos D (2012) Joint models for longitudinal and time-to-event data with applications in R. Chapman and Hall/CRC, London
Rizopoulos D (2016) The R package jmbayes for fitting joint models for longitudinal and time-to-event data using MCMC. J Stat Softw 72(7):1–46
Rizopoulos D, Ghosh P (2011) A Bayesian semiparametric multivariate joint model for multiple longitudinal outcomes and a time-to-event. Stat Med 30(12):1366–1380
Rizopoulos D, Verbeke G, Lesaffre E (2009) Fully exponential Laplace approximations for the joint modelling of survival and longitudinal data. J R Stat Soc B 71:637–654
Rubin DB (1976) Inference and missing data. Biometrika 63(3):581–592
Song X, Wang C (2008) Semiparametric approaches for joint modeling of longitudinal and survival data with time-varying coefficients. Biometrics 64(2):557–566
Souza-Jr P, Szwarcwald CL, Castilho EA (2007) Delay in introducing antiretroviral therapy in patients infected by HIV in Brazil, 2003–2006. Clin Sci 62(5):579–584
Spiegelhalter D, Best N, Carlin B, Van der Linde A (2002) Bayesian measures of model complexity and fit (with discussion). J R Stat Soc B 64(4):583–639
Tang NS, Tang AM, Pan DD (2014) Semiparametric Bayesian joint models of multivariate longitudinal and survival data. Comput Stat Data Anal 77:113–129
Tseng YK, Hsieh FS, Wang JL (2005) Joint modelling of accelerated failure time and longitudinal data. Biometrika 92(3):587–603
Watanabe S (2010) Asymptotic equivalence of Bayes cross validation and widely applicable information criterion in singular learning theory. J Mach Learn Res 11(455):3571–3591
Wulfsohn M, Tsiatis A (1997) A joint model for survival and longitudinal data measured with error. Biometrics 53:330–339
Yang M, Luo S, DeSantis S (2018) Bayesian quantile regression joint models: Inference and dynamic predictions. Stat Methods Med Res 28(8):2524–2537. https://doi.org/10.1177/0962280218784757
Yu Z, Liu L, Bravata DM, Williams LS (2014) Joint model of recurrent events and a terminal event with time-varying coefficients. Biom J 56(2):183–197
Zhang H, Huang Y (2019) Quantile regression-based Bayesian joint modeling analysis of longitudinal-survival data, with application to an AIDS cohort study. Lifetime Data Anal 26(2):339–368. https://doi.org/10.1007/s10985-019-09478-w
Zhu H, Ibrahim J, Chi YY, Tang N (2012) Bayesian influence measures for joint models for longitudinal and survival data. Biometrics 68(3):954–964
Zhuang R, Chen YQ (2019) Measuring surrogacy in clinical research. Stat Biosci pp 1–29
Acknowledgements
The author wish to thank Maria Goretti Fonseca and Cláudia Coeli for the database. This work was partially funded by Fundação para a Ciência e a Tecnologia projects UID/BIM/04585/2016 and UIDB/00006/2020.
Author information
Authors and Affiliations
Corresponding author
Additional information
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary Information
Appendix
Appendix
Following Gelman et al. (2014) we define WAIC as
where \(V_{s=1}^S [a_s]\), represents the posterior sample variance and the over line stands for the posterior sample mean. Like the DIC models with smaller WAIC values should be preferred.
Rights and permissions
About this article
Cite this article
Martins, R. A flexible link for joint modelling longitudinal and survival data accounting for individual longitudinal heterogeneity. Stat Methods Appl 31, 41–61 (2022). https://doi.org/10.1007/s10260-021-00566-6
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s10260-021-00566-6