Home Choice of baseline hazards in joint modeling of longitudinal and time-to-event cancer survival data
Article
Licensed
Unlicensed Requires Authentication

Choice of baseline hazards in joint modeling of longitudinal and time-to-event cancer survival data

  • Anand Hari , Edakkalathoor George Jinto , Divya Dennis , Kumarapillai Mohanan Nair Jagathnath Krishna EMAIL logo , Preethi S. George , Sivasevan Roshni and Aleyamma Mathew
Published/Copyright: May 13, 2024

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.


Corresponding author: Kumarapillai Mohanan Nair Jagathnath Krishna, Division of Cancer Epidemiology and Biostatistics, Regional Cancer Centre, Thiruvananthapuram, Kerala, India, E-mail:

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.

  1. Research ethics: The local Institutional Scientific Review committee with IRB No.12/2018/13, deemed the study exempted from human ethics committee review.

  2. Author contributions: All authors have accepted responsibility for the entire content of this manuscript and approved its submission.

  3. Competing interests: Authors state no conflict of interest.

  4. 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.

  5. 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

Received: 2023-10-21
Accepted: 2024-04-23
Published Online: 2024-05-13

© 2024 Walter de Gruyter GmbH, Berlin/Boston

Downloaded on 17.9.2025 from https://www.degruyterbrill.com/document/doi/10.1515/sagmb-2023-0038/html
Scroll to top button