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.

Relation between three classes of structural models for the effect of a time-varying exposure on survival

Adjusting for time‐varying confounders in survival analysis using structural nested cumulative survival time models

Instrumental variables estimation of exposure effects on a time-to-event response using structural cumulative survival models

Estimating effects of time-varying exposures on mortality risk

Estimation of the causal effects of time-varying treatments in nested case-control studies using marginal structural Cox models

A multi-state model incorporating estimation of excess hazards and multiple time scales

Semiparametric Additive Modeling of the Restricted Mean Survival Time

A Tutorial on Multilevel Survival Analysis: Methods, Models and Applications

Multiple imputation in Cox regression when there are time-varying effects of exposures

Comparison of methods for analyzing environmental mixtures effects on survival outcomes and application to a population-based cohort study

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

Causal models for estimating the effects of weight gain on mortality

Methods to Analyze Time-to-Event Data: The Cox Regression Analysis

Statistical Considerations for Evaluating Treatment Effect under Various Non-proportional Hazard Scenarios

A penalized Cox proportional hazards model with multiple time-varying exposures

Modeling of cumulative effects of time-varying drug exposures on within-subject changes in a continuous outcome

Comparison of nested case-control and survival analysis methodologies for analysis of time-dependent exposure

Modified Cox Models: A Simulation Study on Different Survival Distributions, Censoring Rates, and Sample Sizes

Causal effect estimation in survival analysis with high dimensional confounders

Matching methods for obtaining survival functions to estimate the effect of a time-dependent treatment