An Accelerated Failure Time Cure Model with Shifted Gamma Frailty and Its Application to Epidemiological Research
Abstract
:1. Introduction
2. Regression Models for Survival Time Response
2.1. Literature Review
2.2. Problem Formulation
2.3. Mixture Cure Models
2.4. Accelerated Failure Time Models
2.5. Frailty Models
3. A Novel Accelerated Failure Time Frailty Mixture Cure Model
3.1. Proposed Model
3.2. Estimation Method and Its Algorithm
- Step 1
- Set the initial values .
- Step 2
- Calculate the sample version of Equation (22) for . That is,
- Step 3
- Find the updated value for each parameter
- Step 4
- If the convergence condition is satisfied, terminate the algorithm and set the estimated values to . Otherwise, increase the value of k by 1 and return to step 2.
4. Numerical Examples
4.1. Simulations
4.1.1. Setting
- (i)
- ,
- (ii)
- ,
4.1.2. Results
4.2. Real Data Example
4.2.1. Dataset and Previous Study
- Proportional hazard model;
- AFT model;
- Mixture cure model with the proportional hazard model;
- Mixture cure model with the AFT model;
- Proposed model.
4.2.2. Results
5. Discussion
Author Contributions
Funding
Institutional Review Board Statement
Informed Consent Statement
Data Availability Statement
Acknowledgments
Conflicts of Interest
References
- Cox, D.R. Regression Models and Life-Tables. J. R. Stat. Soc. Ser. B Stat. Methodol. 1972, 34, 87–220. [Google Scholar] [CrossRef]
- Orbe, J.; Ferreira, E.; Núñez-Antón, V. Comparing proportional hazards and accelerated failure time models for survival analysis. Stat. Med. 2002, 21, 3493–3510. [Google Scholar] [CrossRef] [PubMed]
- Berkson, J.; Gage, R. Survival Curve for Cancer Patients Following Treatment. J. Am. Stat. Assoc. 1952, 47, 501–515. [Google Scholar] [CrossRef]
- Chen, M.; Ibrahim, J.; Sinha, D. A New Bayesian Model for Survival Data with a Surviving Fraction. J. Am. Stat. Assoc. 1999, 94, 909–919. [Google Scholar] [CrossRef]
- Sy, J.P.; Taylor, J. Estimation in a Cox Proportional Hazards Cure Model. Biometrics 2000, 56, 227–236. [Google Scholar] [CrossRef]
- Vaupel, J.W.; Manton, K.G.; Stallard, E. The Impact of Heterogeneity in Individual Frailty on the Dynamics of Mortality. Demography 1979, 16, 439–454. [Google Scholar] [CrossRef]
- Yamaguchi, K. Accelerated Failure-Time Regression Models with a Regression Model of Surviving Fraction:An Application to the Analysis of “Permanent Employment” in Japan. J. Am. Stat. Assoc. 1992, 87, 284–292. [Google Scholar]
- Li, C.; Taylor, J. A semi-parametric accelerated failure time cure model. Stat. Med. 2002, 21, 3235–3247. [Google Scholar] [CrossRef]
- Yu, B. A frailty mixture cure model with application to hospital readmission cata. Biom. J. 2008, 50, 386–394. [Google Scholar] [CrossRef]
- Hutton, J.L.; Monaghan, P.F. Choice of Parametric Accelerated Life and Proportional Hazards Models for Survival Data: Asymptotic Results. Lifetime Data Anal. 2002, 8, 375–393. [Google Scholar] [CrossRef]
- Patel, K.; Kay, R.; Rowell, L. Comparing proportional hazards and accelerated failure time models: An application in influenza. Pharmaceut. Statist. 2006, 5, 213–224. [Google Scholar] [CrossRef] [PubMed]
- Aalen, O.O. Heterogeneity in survival analysis. Statist. Med. 1988, 7, 1121–1137. [Google Scholar] [CrossRef] [PubMed]
- Pan, W. Using Frailties in the Accelerated Failure Time Model. Lifetime Data Anal. 2001, 7, 55–64. [Google Scholar] [CrossRef] [PubMed]
- Price, D.; Manatunga, A. Modelling survival data with a cured fraction using frailty models. Stat. Med. 2001, 20, 1515–1527. [Google Scholar] [CrossRef]
- Chen, P.; Zhang, J.; Zhang, R. Estimation of the accelerated failure time frailty model under generalized gamma frailty. Comput. Stat. Data Anal. 2013, 62, 171–180. [Google Scholar] [CrossRef]
- He, M. Some Flexible Families of Mixture Cure Frailty Models and Associated Inference. Ph.D. Thesis, McMaster University, Hamilton, ON, Canada, 2021. Available online: http://hdl.handle.net/11375/26258 (accessed on 20 July 2022).
- O’Quingley, J. Proportional Hazards Regeression; Springer: New York, NY, USA, 2008; 542p. [Google Scholar]
- Dempster, A.P.; Laird, N.M.; Rubin, D.B. Maximum Likelihood from Incomplete Data via the EM algorithm. J. R. Stat. Soc. Ser. B Stat. Methodol. 1977, 39, 1–38. [Google Scholar]
- Scolas, S.; Ghouch, A.; Legrabd, C.; Oulhaj, A. Variable selection in a flexible parametric mixture cure model with interval-censored data. Stat. Med. 2016, 35, 1210–1225. [Google Scholar] [CrossRef] [Green Version]
- Kalbfreisch, J.; Prentice, R. The Statistical Analysis of Failure Time Data, 2nd ed.; Wiley Series in Probability and Statistics; John Wiley & Sons, Inc.: Hoboken, NJ, USA, 2002; 462p. [Google Scholar]
- Wei, L.J. The accelerated failure time model: A useful alternative to the cox regression model in survival analysis. Stat. Med. 1992, 11, 1871–1879. [Google Scholar] [CrossRef]
- Aalen, O.O. Two Examples of Modelling Heterogeneity in Survival Analysis. Scand. Stat. Theory Appl. 1987, 14, 19–25. [Google Scholar]
- Peng, Y.; Zhang, J. Estimation method of the semiparametric mixture cure gamma frailty model. Stat. Med. 2008, 27, 5177–5194. [Google Scholar] [CrossRef]
- Lambert, P.; Collett, D.; Kimber, A.; Johnson, R. Parametric accelerated failure time models with ramdom effects and an application to kidney transplant survival. Stat. Med. 2004, 23, 3177–3192. [Google Scholar] [CrossRef] [PubMed] [Green Version]
- Elbers, C.; Ridder, G. True and spurious duration dependence: The identifiability of the proportional hazards model. Rev. Econ. Stud. 1982, 49, 403–409. [Google Scholar] [CrossRef]
- Stacy, E.W. A generalization of the gamma distribution. Ann. Math. Statist. 1962, 33, 1187–1192. [Google Scholar] [CrossRef]
- Prentice, R.L. A log gamma model and its maximum likelihood estimation. Biometrika 1974, 61, 539–544. [Google Scholar] [CrossRef]
- Kim, S.; Lee, J.Y.; Sung, D.K. A Shifted Gamma Distribution Model for Long-Range Dependent Internet Traffic. IEEE Commun. Lett. 2003, 7, 124–126. [Google Scholar]
- You, X. Approximation of the median of the gamma distribution. J. Number Theory 2017, 174, 487–493. [Google Scholar] [CrossRef]
- Casella, G.; Berger, R.L. Statistical Inference, 2nd ed.; Thomson Learning: Pacific Grove, CA, USA, 2002; 660p. [Google Scholar]
- Tsutatani, H.; Funamoto, M.; Sugiyama, D.; Kuwabara, K.; Miyamatsu, N.; Watanabe, K.; Okamura, T. Association between lifestyle factors assessed by standard question items of specific health checkup and the incidence of metabolic syndrome and hypertension in community dwellers: A five-year cohort study of National Health Insurance beneficiaries in Habikino City. Nihon Koshu Eisei Zasshi 2017, 64, 258–269. [Google Scholar]
- Grambsch, P.M.; Therneau, T.M. Proportional Hazards Tests and Diagnostics Based on Weighted Residuals. Biometrika 1994, 81, 515–526. [Google Scholar] [CrossRef]
- Heinze, G.; Schemper, M. A solution to the problem of separation in logistic regression. Stat. Med. 2002, 21, 2409–2419. [Google Scholar] [CrossRef]
- Firth, D. Bias reduction of Maximum Likelihood Estimates. Biometrika 1993, 80, 27–38. [Google Scholar] [CrossRef]
- Akaike, H. Information theory and an extension of the maximum likelihood principle. In Proceedings of the 2nd International Symposium on Information Theory, Tsahkadsor, Armenia, 2–8 September 1971; Akademiai Kiado: Budapest, Hungary, 1973; pp. 267–281. [Google Scholar]
- Hjort, N.L.; Claeskens, G. Focused Information Criteria and Model Averaging for the Cox Hazard Regression Model. J. Am. Stat. Assoc. 2006, 101, 1449–1464. [Google Scholar] [CrossRef]
- Jullum, M.; Hjort, N.L. What price semiparametric Cox regression? Lifetime Data Anal. 2019, 25, 405–438. [Google Scholar] [CrossRef] [PubMed] [Green Version]
- Xu, L.; Zhang, J. An EM-like algorithm for the semiparametric accelerated failure time gamma frailty model. Comput. Stat. Data. Anal. 2010, 54, 1467–1474. [Google Scholar] [CrossRef]
- Sun, J. The Statistical Analysis of Interval-Censored Failure Time Data; Springer: New York, NY, USA, 2006; 304p. [Google Scholar]
- Karthik, S.; Bhadoria, R.S.; Lee, J.G.; Sivaraman, A.K.; Samanta, S.; Balasundaram, A.; Chaurasia, B.K.; Ashokkumar, S. Prognostic Kalman Filter Based Bayesian Learning Model for Data Accuracy Prediction. Comput. Mater. Contin. 2022, 72, 243–259. [Google Scholar] [CrossRef]
- Singh, L.K.; Garg, H.; Khanna, M.; Bhadoria, R.S. An enhanced deep image model for glaucoma diagnosis using feature-based detection in retinal fundus. Med. Biol. Eng. Comput. 2021, 59, 333–353. [Google Scholar] [CrossRef]
Literature | Cured Patients | Uncured Model | Frailty |
---|---|---|---|
Sy and Taylor [5] | ✓ | PH | – |
Vaupel [6], Aalen [12] | – | PH | gamma |
Pan [13] | – | AFT | gamma, log-normal |
Chen et al. [15] | – | AFT | generalized gamma |
Yu [9], Price and Manatunga [14] | ✓ | PH | gamma |
He [16] | ✓ | AFT | generalized gamma |
Present study | ✓ | AFT | shifted gamma |
Parameter | True Value | Mean (SD) | |||||
---|---|---|---|---|---|---|---|
n = 100 | n = 500 | n = 1000 | |||||
−0.5 | |||||||
0.5 | |||||||
−0.8 | |||||||
3 | |||||||
1 | |||||||
q | 2 | ||||||
0.1 | |||||||
0.5 | |||||||
−1 | |||||||
0.6 |
Parameter | True Value | Mean (SD) | |||||
---|---|---|---|---|---|---|---|
n = 100 | n = 500 | n = 1000 | |||||
−0.5 | |||||||
0.5 | |||||||
−0.8 | |||||||
3 | |||||||
1 | |||||||
q | 2 | ||||||
3 | |||||||
0.1 | |||||||
0.5 | |||||||
−1 | |||||||
0.6 |
Model | Distribution | Number of Parameters | AIC |
---|---|---|---|
Proportional hazard (PH) | Exponential | 19 | 7120.407 |
Weibull | 20 | 7030.615 | |
AFT model | Log-normal | 20 | 7002.352 |
Generalized gamma | 21 | 6999.232 | |
Mixture cure + PH | Exponential | 38 | 7103.040 |
Weibull | 39 | 7047.378 | |
Mixture cure + AFT | Log-normal | 39 | 7005.633 |
Generalized gamma* | 40 | 7004.674 | |
Generalized gamma * | 29 | 6992.934 | |
Mixtrue cure + AFT frailty | Generalized gamma Shifted gamma* | 41 | 7003.011 |
Generalized gamma Shifted gamma * | 30 | 6987.012 |
Model | Distribution | Number of Parameters | AIC |
---|---|---|---|
Proportional hazard (PH) | Exponential | 19 | 11,804.01 |
Weibull | 20 | 11,644.16 | |
AFT | Log-normal | 20 | 11,596.45 |
Generalized gamma | 21 | 11,586.00 | |
Mixture cure + PH | Exponential | 38 | 11,798.48 |
Weibull | 39 | 11,669.49 | |
Mixture cure + AFT | Log-normal | 39 | 11,609.63 |
Generalized gamma * | 40 | 11,600.23 | |
Generalized gamma * | 27 | 11,579.07 | |
Mixture cure + AFT frailty | Generalized gamma Shifted gamma * | 41 | 11,596.26 |
Generalized gamma Shifted gamma * | 26 | 11,575.05 |
Covariate | Inference of (Regression Coefficients for the Uncured Group) | ||
---|---|---|---|
Estimates | 95% CI | p-Value | |
age | |||
waist | |||
exe1h_day | |||
exe30_2_week | |||
sleep_good | |||
walk_speed | |||
eat_speed_n | |||
eat_speed_f | |||
eat_b_sleep | |||
snacking | |||
breakfast | |||
weight_move | |||
plus10kg | |||
smoking | |||
drink_amount2 | |||
drink_amount3 | 0.054 * | ||
drink_amount4 | |||
drink_amount5 | |||
Covariate | Inference of (Regression Coefficients for the Cured Group) | ||
Point Estimates | 95% CI | p-Value | |
Intercept | |||
age | |||
eat_speed_n | |||
eat_speed_f | 0.0664 * | ||
kanshyoku | |||
plus10kg | |||
drink_amount2 | 0.0605 * | ||
drink_amount4 |
Parameter | Distribution | |||
---|---|---|---|---|
Generalized Gamma | Shifted Gamma | |||
μ | σ | q | θ | |
Estimates | ||||
Standard error |
Covariate | Inference of (Regression Coefficients for the Uncured Group) | ||
---|---|---|---|
Estimates | 95% CI | p-Value | |
age | * | ||
waist | |||
exe1h_day | |||
exe30_2_week | |||
sleep_good | |||
walk_speed | |||
eat_speed_n | |||
eat_speed_f | |||
eat_b_sleep | |||
snacking | |||
breakfast | |||
weight_move | |||
plus10kg | |||
smoking | |||
drink_amount2 | |||
drink_amount3 | |||
drink_amount4 | |||
drink_amount5 | |||
Covariate | Inference of (Regression Coefficients for the Cured Group) | ||
Point Estimates | 95% CI | p-Value | |
Intercept | <0.001 | ||
age | <0.001 | ||
eat_speed_f | |||
plus10kg |
Parameter | Distribution | |||
---|---|---|---|---|
Generalized Gamma | Shifted Gamma | |||
μ | σ | q | θ | |
Estimates | ||||
Standard error |
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations. |
© 2022 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).
Share and Cite
Aida, H.; Hayashi, K.; Takeuchi, A.; Sugiyama, D.; Okamura, T. An Accelerated Failure Time Cure Model with Shifted Gamma Frailty and Its Application to Epidemiological Research. Healthcare 2022, 10, 1383. https://doi.org/10.3390/healthcare10081383
Aida H, Hayashi K, Takeuchi A, Sugiyama D, Okamura T. An Accelerated Failure Time Cure Model with Shifted Gamma Frailty and Its Application to Epidemiological Research. Healthcare. 2022; 10(8):1383. https://doi.org/10.3390/healthcare10081383
Chicago/Turabian StyleAida, Haro, Kenichi Hayashi, Ayano Takeuchi, Daisuke Sugiyama, and Tomonori Okamura. 2022. "An Accelerated Failure Time Cure Model with Shifted Gamma Frailty and Its Application to Epidemiological Research" Healthcare 10, no. 8: 1383. https://doi.org/10.3390/healthcare10081383
APA StyleAida, H., Hayashi, K., Takeuchi, A., Sugiyama, D., & Okamura, T. (2022). An Accelerated Failure Time Cure Model with Shifted Gamma Frailty and Its Application to Epidemiological Research. Healthcare, 10(8), 1383. https://doi.org/10.3390/healthcare10081383