Choice of baseline hazards in joint modeling of longitudinal and time-to-event cancer survival data
-
Anand Hari
, Preethi S. George
Abstract
Longitudinal time-to-event analysis is a statistical method to analyze data where covariates are measured repeatedly. In survival studies, the risk for an event is estimated using Cox-proportional hazard model or extended Cox-model for exogenous time-dependent covariates. However, these models are inappropriate for endogenous time-dependent covariates like longitudinally measured biomarkers, Carcinoembryonic Antigen (CEA). Joint models that can simultaneously model the longitudinal covariates and time-to-event data have been proposed as an alternative. The present study highlights the importance of choosing the baseline hazards to get more accurate risk estimation. The study used colon cancer patient data to illustrate and compare four different joint models which differs based on the choice of baseline hazards [piecewise-constant Gauss–Hermite (GH), piecewise-constant pseudo-adaptive GH, Weibull Accelerated Failure time model with GH & B-spline GH]. We conducted simulation study to assess the model consistency with varying sample size (N = 100, 250, 500) and censoring (20 %, 50 %, 70 %) proportions. In colon cancer patient data, based on Akaike information criteria (AIC) and Bayesian information criteria (BIC), piecewise-constant pseudo-adaptive GH was found to be the best fitted model. Despite differences in model fit, the hazards obtained from the four models were similar. The study identified composite stage as a prognostic factor for time-to-event and the longitudinal outcome, CEA as a dynamic predictor for overall survival in colon cancer patients. Based on the simulation study Piecewise-PH-aGH was found to be the best model with least AIC and BIC values, and highest coverage probability(CP). While the Bias, and RMSE for all the models showed a competitive performance. However, Piecewise-PH-aGH has shown least bias and RMSE in most of the combinations and has taken the shortest computation time, which shows its computational efficiency. This study is the first of its kind to discuss on the choice of baseline hazards.
Funding source: Indian council for Medical research-Department of Health Research, Government of India
Award Identifier / Grant number: R.11012/03/2021-GIA/HR
Acknowledgment
The authors are grateful to the parent institution for the support given to carry out this work. The authors are also thankful to the unknown referees for the critical comments which helped in improving the manuscript.
-
Research ethics: The local Institutional Scientific Review committee with IRB No.12/2018/13, deemed the study exempted from human ethics committee review.
-
Author contributions: All authors have accepted responsibility for the entire content of this manuscript and approved its submission.
-
Competing interests: Authors state no conflict of interest.
-
Research funding: The authors acknowledge the funding received from Indian Council for Medical Research – Department of Health Research, Government of India with award no. R.11012/03/2021-GIA/HR.
-
Data availability: Not applicable.
References
Brown, E.R. and Ibrahim, J.G. (2003). A Bayesian semiparametric joint hierarchical model for longitudinal and survival data. Biometrics 59: 221–228. https://doi.org/10.1111/1541-0420.00028.Search in Google Scholar PubMed
Cox, D.R. (1972). Regression models and life-tables. J. R. Stat. Soc. Series B Methodol. 34: 187–202. https://doi.org/10.1111/j.2517-6161.1972.tb00899.x.Search in Google Scholar
Crowther, M.J., Abrams, K.R., and Lambert, P.C. (2012). Flexible parametric joint modelling of longitudinal and survival data. Stat. Med. 31: 4456–4471. https://doi.org/10.1002/sim.5644.Search in Google Scholar PubMed
Crowther, M.J., Abrams, K.R., Lambert, P.C., and Proust-Lima, C. (2013). Joint modeling of longitudinal and survival data. Stat. Med. 13: 165–184. https://doi.org/10.1177/1536867x1301300112.Search in Google Scholar
Furgal, A.K.C., Sen, A., and Taylor, J.M. (2019). Review and Comparison of computational approaches for joint longitudinal and time-to-event models. Int. Stat. Rev. 87: 393–418. https://doi.org/10.1111/insr.12322.Search in Google Scholar PubMed PubMed Central
Harris, E.I., Lewin, D.N., Wang, H.L., Lauwers, G.Y., Srivastava, A., Shyr, Y., Shakhtour, B., Revetta, F., and Washington, M.K. (2008). Lymphovascular invasion in colorectal cancer: an interobserver variability study. Am. J. Surg. Pathol. 32: 1816–1821. https://doi.org/10.1097/pas.0b013e3181816083.Search in Google Scholar PubMed PubMed Central
Henderson, R., Diggle, P., and Dobson, A. (2000). Joint modelling of longitudinal measurements and event time data. Biostatistics 1: 465–480. https://doi.org/10.1093/biostatistics/1.4.465.Search in Google Scholar PubMed
Hickey, G.L., Philipson, P., Jorgensen, A., and Kolamunnage-Dona, R. (2016). Joint modelling of time-to-event and multivariate longitudinal outcomes: recent developments and issues. BMC Med. Res. Methodol. 16: 1–5. https://doi.org/10.1186/s12874-016-0212-5.Search in Google Scholar PubMed PubMed Central
Hickey, G.L., Philipson, P., Jorgensen, A., and Kolamunnage-Dona, R. (2018). Joint models of longitudinal and time-to-event data with more than one event time outcome: a review. Int. J. Biostat. 14: 20170047, https://doi.org/10.1515/ijb-2017-0047.Search in Google Scholar PubMed
Hsieh, F., Tseng, Y.K., and Wang, J.L. (2006). Joint modeling of survival and longitudinal data: likelihood approach revisited. Biometrics 62: 1037–1043. https://doi.org/10.1111/j.1541-0420.2006.00570.x.Search in Google Scholar PubMed
Lawrence, G.A., Boye, M.E., Crowther, M.J., Ibrahim, J.G., Quartey, G., Micallef, S., and Bois, F.Y. (2015). Joint modeling of survival and longitudinal non-survival data: current methods and issues. Report of the DIA Bayesian joint modeling working group. Stat. Med. 34: 2181–2195. https://doi.org/10.1002/sim.6141.Search in Google Scholar PubMed PubMed Central
Liebig, C., Ayala, G., Wilks, J., Verstovsek, G., Liu, H., Agarwal, N., Berger, D.H., and Albo, D. (2009). Perineural invasion is an independent predictor of outcome in colorectal cancer. J. Clin. Oncol 27: 5131–5137, https://doi.org/10.1200/JCO.2009.22.4949.Search in Google Scholar PubMed PubMed Central
Mchunu, N.N., Mwambi, H.G., Reddy, T., Yende-Zuma, N., and Naidoo, K. (2020). Joint modelling of longitudinal and time-to-event data: an illustration using CD4 count and mortality in a cohort of patients initiated on antiretroviral therapy. BMC Infect. Dis. 20: 1–9. https://doi.org/10.1186/s12879-020-04962-3.Search in Google Scholar PubMed PubMed Central
Mohd Suan, M.A., Tan, W.L., Soelar, S.A., Ismail, I., and Abu Hassan, M.R. (2015). Intestinal obstruction: predictor of poor prognosis in colorectal carcinoma? Epidemiol. Health 37: e2015017. https://doi.org/10.4178/epih/e2015017.Search in Google Scholar PubMed PubMed Central
Mwanyekange, J., Mwalili, S., Ngesa, O., and Proust-Lima, C. (2018). Bayesian inference in a joint model for longitudinal and time to event data with Gompertz baseline hazards. Mod. Appl. Sci. 12: 159–172. https://doi.org/10.5539/mas.v12n9p159.Search in Google Scholar
Rizopoulos, D. (2010). JM: an R package for the joint modelling of longitudinal and time-to-event data. J. Stat. Softw. 35: 1–33. https://doi.org/10.18637/jss.v035.i09.Search in Google Scholar
Rizopoulos, D. (2012a). Joint models for longitudinal and time-to-event data: With applications in R. CRC Press: New York.10.1201/b12208Search in Google Scholar
Rizopoulos, D. (2012b). Fast fitting of joint models for longitudinal and event time data using a pseudo-adaptive Gaussian quadrature rule. Comput. Stat. Data Anal. 56: 491–501. https://doi.org/10.1016/j.csda.2011.09.007.Search in Google Scholar
Rustand, D., Briollais, L., Tournigand, C., and Rondeau, V. (2022). Two-part joint model for a longitudinal semicontinuous marker and a terminal event with application to metastatic colorectal cancer data. Biostatistics 23: 50–68. https://doi.org/10.1093/biostatistics/kxaa012.Search in Google Scholar PubMed PubMed Central
Saulnier, T., Philipps, V., Meissner, W.G., Rascol, O., Pavy-Le Traon, A., Foubert-Samier, A., and Proust-Lima, C. (2022). Joint models for the longitudinal analysis of measurement scales in the presence of informative dropout. Methods 203: 142–151. https://doi.org/10.1016/j.ymeth.2022.03.003.Search in Google Scholar PubMed
Shizgal, B. (1981). A Gaussian quadrature procedure for use in the solution of the Boltzmann equation and related problems. J. Comput. Phys. 41: 309–328. https://doi.org/10.1016/0021-9991(81)90099-1.Search in Google Scholar
Sudell, M., Kolamunnage-Dona, R., and Tudur-Smith, C. (2016). Joint models for longitudinal and time-to-event data: a review of reporting quality with a view to meta-analysis. BMC Med. Res. Methodol. 16: 1. https://doi.org/10.1186/s12874-016-0272-6.Search in Google Scholar PubMed PubMed Central
Tang, A.M. and Tang, N.S. (2015). Semiparametric Bayesian inference on skew–normal joint modeling of multivariate longitudinal and survival data. Stat. Med. 34: 824–843. https://doi.org/10.1002/sim.6373.Search in Google Scholar PubMed
Tseng, Y.K., Hsieh, F., and Wang, J.L. (2005). Joint modelling of accelerated failure time and longitudinal data. Biometrika 92: 587–603. https://doi.org/10.1093/biomet/92.3.587.Search in Google Scholar
Tsiatis, A.A. and Davidian, M. (2004). Joint modeling of longitudinal and time-to-event data: an overview. Stat. Sin. 14: 809–834.Search in Google Scholar
Wolbers, M., Babiker, A., Sabin, C., Young, J., Dorrucci, M., Chene, G., Mussini, C., Porter, K., and Bucher, H.C. (2010). CASCADE Collaboration, Pretreatment CD4 cell slope and progression to AIDS or death in HIV-infected patients initiating antiretroviral therapy—the CASCADE collaboration: a collaboration of 23 cohort studies. PLoS Med. 7: e1000239. https://doi.org/10.1371/journal.pmed.1000239.Search in Google Scholar PubMed PubMed Central
Wu, L., Liu, W., Yi, G.Y., and Huang, Y. (2012). Analysis of longitudinal and survival data: joint modeling, inference methods, and issues. J. Probab. Stat. 2012: 640153, https://doi.org/10.1155/2012/640153.Search in Google Scholar
Wulfsohn, M.S. and Tsiatis, A.A. (1997). A joint model for survival and longitudinal data measured with error. Biometrics 53: 330–339, https://doi.org/10.2307/2533118.Search in Google Scholar
© 2024 Walter de Gruyter GmbH, Berlin/Boston
Articles in the same Issue
- Frontmatter
- Research Articles
- Empirically adjusted fixed-effects meta-analysis methods in genomic studies
- A CNN-CBAM-BIGRU model for protein function prediction
- A heavy-tailed model for analyzing miRNA-seq raw read counts
- Flexible model-based non-negative matrix factorization with application to mutational signatures
- Choice of baseline hazards in joint modeling of longitudinal and time-to-event cancer survival data
- Assessing the feasibility of statistical inference using synthetic antibody-antigen datasets
- A global test of hybrid ancestry from genome-scale data
- Integrative pathway analysis with gene expression, miRNA, methylation and copy number variation for breast cancer subtypes
- Bayesian LASSO for population stratification correction in rare haplotype association studies
Articles in the same Issue
- Frontmatter
- Research Articles
- Empirically adjusted fixed-effects meta-analysis methods in genomic studies
- A CNN-CBAM-BIGRU model for protein function prediction
- A heavy-tailed model for analyzing miRNA-seq raw read counts
- Flexible model-based non-negative matrix factorization with application to mutational signatures
- Choice of baseline hazards in joint modeling of longitudinal and time-to-event cancer survival data
- Assessing the feasibility of statistical inference using synthetic antibody-antigen datasets
- A global test of hybrid ancestry from genome-scale data
- Integrative pathway analysis with gene expression, miRNA, methylation and copy number variation for breast cancer subtypes
- Bayesian LASSO for population stratification correction in rare haplotype association studies