In the survival analysis, when there is a cure fraction and the occurrence times of events are correlated, the cure frailty model is utilized. The main objective is to propose a method of analysis for two types of correlated frailty in the non-mixture cured model in order to separate the individual and shared heterogeneity between subjects. The cure models with correlated frailty and promotion time are considered. In both models, the likelihood function are based on piecewise exponential distribution for hazard function. To estimate the parameters, hierarchical Bayesian modeling is employed. Due to non-closed forms of the posteriors, they are estimated by MCMC algorithms. The Cox correlated frailty model is used as a benchmark and models are compared by DIC Criterion . The results show the superiority of cure models with correlated frailty.

The survival analysis methods are usually conducted based on assumption that the population is homogeneity. However, generally, this assumption in most cases is unrealistic, because of unobserved risk factors or subject specific random effect. Disregarding the heterogeneity leads to unbiased results. So frailty model as a mixed model was used to adjust for uncertainty that cannot be explained by observed factors in survival analysis. In this paper, family of power variance function distributions that includes gamma and inverse Gaussian distribution were introduced and evaluated for frailty effects. Finally the proportional hazard frailty models with Weibull baseline hazard as a parametric model used for analyzing survival data of the colorectal cancer patients.