Abstract
Survival data analysis often uses the Cox proportional hazards (PH) model. This model is widely applied due to its straightforward interpretation of the hazard ratio under the assumption that the hazard rates for two subjects remain constant over time. However, in several randomized clinical trials with long-term survival data comparing two new treatments, it is frequently observed that Kaplan-Meier plots exhibit crossing survival curves. This violation of the PH assumption of the Cox PH model can not be applied to evaluate the treatment's effect on survival. This paper introduces a novel long-term survival model with non-PH that incorporates a frailty term into the hazard function. This model allows us to examine the effect of prognostic factors on survival and quantify the degree of unobservable heterogeneity. The model parameters are estimated using the maximum likelihood estimation procedure, and we evaluate the performance of the proposed models through simulation studies. Additionally, we demonstrate the applicability of our approach by fitting the models to a real skin cancer dataset.