We read with great interest the recent article by Sung et al. [1] titled “Association of uterine leiomyoma with type 2 diabetes mellitus in young women: a population-based cohort study” in Diabetes Metabolism Journal. The study found that young women with uterine leiomyoma had a higher risk of developing type 2 diabetes mellitus (T2DM), and undergoing myomectomy appeared to reduce this risk. However, we note several biases in the use of the cox proportional hazards (CoxPH) model that the authors did not address.
The established criteria may result in mixed censoring outcomes, i.e., right-censoring and interval-censoring events [2,3]. Events of T2DM diagnosed through medical records could result in interval-censoring if they occurred between follow-up visits, and right-censoring for diagnosed between the end of follow-up and the time of data analysis. The CoxPH model primarily handles right-censored data. In contrast, the accelerated failure time (AFT) model is often preferred for scenarios involving various types of censored data [4]. The AFT model can effectively handle left-censored, right-censored, and interval-censored data by appropriately adjusting the likelihood function [5]. By using the ‘survival’ and ‘icenReg’ packages, mixed censored data can be fitted and analyzed, and event times can be estimated [6].
Moreover, the CoxPH model requires the proportional hazards assumption, meaning that covariate effects are constant over time [7]. If this assumption is violated, the model may not provide unbiased estimates of the coefficients, and the predictions may not be reliable. The authors should utilize Schoenfeld residuals or alternative methods to evaluate the proportional hazards assumption for the association between covariates and the risk of T2DM. Schoenfeld residuals are calculated as the differences between the observed and expected values of covariates at each failure time [8]. If the residuals exhibit a systematic change over time, it suggests that the effect of the covariate may be time-dependent. When the proportional hazards assumption does not hold, authors should use a stratified cox model, a cox model with time-varying effects, or an AFT model instead of the standard CoxPH model [4,9].
In conclusion, we believe that a re-evaluation considering the potential impact of censoring events and the proportional hazards assumption is necessary. Further research is anticipated to provide more empirical data and clearer insights into this field.
NOTES
-
CONFLICTS OF INTEREST
No potential conflict of interest relevant to this article was reported.
-
FUNDING
This work was supported by the Taishan Scholar Program of Shandong Province (grant number tstp20230665); the National Natural Science Foundation of China (grant numbers 82270724); the Qingdao Key Health Discipline Development Fund; and the Qingdao Key Clinical Specialty Elite Discipline.
REFERENCES
- 1. Sung JH, Kim KS, Han K, Park CY. Association of uterine leiomyoma with type 2 diabetes mellitus in young women: a population-based cohort study. Diabetes Metab J 2024;48:1105-13.PubMed
- 2. Thomas LE, Turakhia MP, Pencina MJ. Competing risks, treatment switching, and informative censoring. JAMA Cardiol 2021;6:871-3.ArticlePubMed
- 3. Leung KM, Elashoff RM, Afifi AA. Censoring issues in survival analysis. Annu Rev Public Health 1997;18:83-104.ArticlePubMed
- 4. Lee KH, Rondeau V, Haneuse S. Accelerated failure time models for semi-competing risks data in the presence of complex censoring. Biometrics 2017;73:1401-12.ArticlePubMedPMCPDF
- 5. Ning J, Qin J, Shen Y. Score estimating equations from embedded likelihood functions under accelerated failure time model. J Am Stat Assoc 2014;109:1625-35.ArticlePubMedPMC
- 6. Crowther MJ, Royston P, Clements M. A flexible parametric accelerated failure time model and the extension to time-dependent acceleration factors. Biostatistics 2023;24:811-31.ArticlePubMedPMCPDF
- 7. Sheng A, Ghosh SK. Effects of proportional hazard assumption on variable selection methods for censored data. Stat Biopharm Res 2020;12:199-209.ArticlePubMedPMC
- 8. Xue X, Xie X, Gunter M, Rohan TE, Wassertheil-Smoller S, Ho GY, et al. Testing the proportional hazards assumption in casecohort analysis. BMC Med Res Methodol 2013;13:88.ArticlePubMedPMCPDF
- 9. Blythe R, Parsons R, Barnett AG, Cook D, McPhail SM, White NM. Prioritising deteriorating patients using time-to-event analysis: prediction model development and internal-external validation. Crit Care 2024;28:247.ArticlePubMedPMCPDF
Citations
Citations to this article as recorded by