Journal of Multidisciplinary Applied Natural Science
##plugins.themes.gdThemes.journalSlogan##
Scopus CiteScore 2025
2.1
Calculated on 05 May, 2025
SJR 2025
0.25
Powered by scimagojr.com
##plugins.themes.gdThemes.journalSlogan##
2.1
Calculated on 05 May, 2025
0.25
Powered by scimagojr.com
https://doi.org/10.47352/jmans.2774-3047.448
Jerry Dwi Trijoyo Purnomo
Purhadi Purhadi
Shofi Andari
Wahyu Dwi Rahmawati
Sri Harini
About 51% of all occurrences of blindness worldwide are caused by cataracts, which are the most common cause of blindness. According to the Rapid Assessment of Avoidable Blindness (RAAB) survey conducted in Indonesia between 2014 and 2016, East Java had the highest rate of blindness (4.4%), with cataracts accounting for 81.1% of cases. Approximately 10,000 cataract cases were reported at the Undaan Eye Hospital in Surabaya, one of the major eye hospitals in East Java, in 2023. This study aims to analyze the recovery time of patients undergoing cataract surgery on both eyes, considering the influencing factors. Data from bilateral cataract patients at the Undaan Eye Hospital in Surabaya between January 2023 and December 2024 were utilized. Gender, age, history of hypertension, heart disease, diabetes, gastric illness, stroke and cholesterol were among the characteristics that were examined. To determine the relationship between the recovery durations of both eyes following surgery, a bivariate survival analysis with Clayton and Gumbel copula functions was used for the analysis. The Gumbel copula was determined to be the better model for describing the recovery of bilateral cataract patients, with a QIC value of 1529.078 and a Kendall's Tau of 0.6636, indicating a moderate and positive dependency between the recovery periods of the right and left eye. The model estimate results showed a substantial correlation between recovery time and a few covariates, including age (hazard ratio of 0.516), history of hypertension (hazard ratio of 0.357), and history of diabetes (hazard ratio of 0.615). Compared to younger individuals, older patients recover more slowly. Furthermore, people with a history of diabetes and hypertension typically recover at a lower rate than those without these conditions.
[1] L. Lin, Y. Liang, G. Jiang, Q. Gan, T. Yang, P. Liao, and H. Liang. (2025). "Global, regional, and national burden of cataract: A comprehensive analysis and projections from 1990 to 2021". PLOS ONE. 20 (6): e0326263. 10.1371/journal.pone.0326263.
DOI: https://doi.org/10.1371/journal.pone.0326263[2] T. Sun and Y. Ding. (2020). "CopulaCenR: Copula Based Regression Models for Bivariate Censored Data in R". The R Journal. 12 (1): 266-282. 10.32614/RJ-2020-025.
DOI: https://doi.org/10.32614/RJ-2020-025[3] G. Cheng, L. Zhou, X. Chen, and J. Z. Huang. (2014). "Efficient Estimation of Semiparametric Copula Models for Bivariate Survival Data". Journal of Multivariate Analysis. 123 : 330-344. 10.1016/j.jmva.2013.10.008.
DOI: https://doi.org/10.1016/j.jmva.2013.10.008[4] Y. Zhang, Y. Chen, and H. Li. (2020). "Copula-Based Analysis of Bilateral Breast Cancer Data". Statistical Methods in Medical Research. 29 (4): 1023-1037.
[5] K. Suresh, J. M. G. Taylor, and A. Tsodikov. (2021). "A Gaussian Copula Approach for Dynamic Prediction of Survival with a Longitudinal Biomarker". Biostatistics. 22 (3): 504-521. 10.1093/biostatistics/kxz049.
DOI: https://doi.org/10.1093/biostatistics/kxz049[6] J. Huang, H. Zhou, and N. Ebrahimi. (2022). "Bayesian Bivariate Cure Rate Models Using Copula Functions". International Journal of Statistics and Probability. 11 (3): 9-21. 10.5539/ijsp.v11n3p9.
DOI: https://doi.org/10.5539/ijsp.v11n3p9[7] D. Petti, A. Eletti, G. Marra, and R. Radice. (2022). "Copula Link-Based Additive Models for Bivariate Time-to-Event Outcomes with General Censoring Scheme". Computational Statistics & Data Analysis. 175 : 107550. 10.1016/j.csda.2022.107550.
DOI: https://doi.org/10.1016/j.csda.2022.107550[8] W. He, G. Y. Yi, and A. Yuan. (2024). "Analysis of Multivariate Survival Data under Semiparametric Copula Models". Canadian Journal of Statistics. 52 (2): 380-413. 10.1002/cjs.11776.
DOI: https://doi.org/10.1002/cjs.11776[9] W. dos Reis Miranda Filho and F. N. Demarqui. (2025). "A Class of Semiparametric Models for Bivariate Survival Data". Lifetime Data Analysis. 31 : 102-125. 10.1007/s10985-024-09642-x.
DOI: https://doi.org/10.1007/s10985-024-09642-x[10] S. Hong, W. Park, Y. Eom, H. M. Kim, and J. S. Song. (2022). "Comparisons of Outcomes and Complications of Immediate Sequential Bilateral Cataract Surgery and Unilateral Cataract Surgery in a Tertiary Hospital in South Korea". Scientific Reports. 12 : 22382. 10.1038/s41598-022-26851-2.
DOI: https://doi.org/10.1038/s41598-022-26851-2[11] D. R. Cox and D. O. Oakes. (1984). "Analysis of Survival Data". Chapman & Hall/CRC, New York.
[12] T. H. Scheike, K. K. Holst, and J. B. Hjelmborg. (2015). "Measuring Early or Late Dependence for Bivariate Lifetimes of Twins". Lifetime Data Analysis. 21 (2): 280-299. 10.1007/s10985-014-9309-5.
DOI: https://doi.org/10.1007/s10985-014-9309-5[13] A. Sklar. (1959). "Fonctions de Répartition à n Dimensions et Leurs Marges". Publications de l'Institut de Statistique de l'Université de Paris. 8 (3): 229-231.
[14] T. Emura and Y.-H. Chen. (2018). "Analysis of Survival Data with Dependent Censoring: Copula-Based Approaches". Springer, Singapore. 10.1007/978-981-10-7164-5.
DOI: https://doi.org/10.1007/978-981-10-7164-5[15] P. A. Shaw and M. P. Fay. (2016). "A Rank Test for Bivariate Time-to-Event Outcomes When One Event Is a Surrogate". Statistics in Medicine. 35 (19): 3413-3423. 10.1002/sim.6950.
DOI: https://doi.org/10.1002/sim.6950[16] E. J. Gumbel. (1960). "Bivariate Exponential Distributions". Journal of the American Statistical Association. 55 (292): 698-707. 10.1080/01621459.1960.10483368.
DOI: https://doi.org/10.1080/01621459.1960.10483368[17] D. G. Clayton. (1978). "A Model for Association in Bivariate Life Tables and Its Application in Epidemiological Studies of Familial Tendency in Chronic Disease Incidence". Biometrika. 65 (1): 141-151. 10.1093/biomet/65.1.141.
DOI: https://doi.org/10.1093/biomet/65.1.141[18] P. B. Mpofu, G. Bakoyannis, C. T. Yiannoutsos, A. W. Mwangi, and M. Mburu. (2020). "A Pseudo-Likelihood Method for Estimating Misclassification Probabilities in Competing-Risks Settings When True-Event Data Are Partially Observed". Biometrical Journal. 62 (7): 1747-1768. 10.1002/bimj.201900198.
DOI: https://doi.org/10.1002/bimj.201900198[19] M. Mei. (2016). "A Goodness-of-Fit Test for Semi-Parametric Copula Models of Right-Censored Bivariate Survival Times". Simon Fraser University.
[20] T. Sun and Y. Ding. (2021). "Copula-Based Semiparametric Regression Method for Bivariate Data under General Interval Censoring". Biostatistics. 22 (2): 315-330. 10.1093/biostatistics/kxz032.
DOI: https://doi.org/10.1093/biostatistics/kxz032[21] D. W. Hosmer, S. Lemeshow, and S. May. (2008). "Applied Survival Analysis: Regression Modeling of Time-to-Event Data". John Wiley & Sons, New Jersey. 10.1002/9780470258019.
DOI: https://doi.org/10.1002/9780470258019[22] L. Y. Hin and T.-G. Wang. (2009). "Working-Correlation-Structure Identification in Generalized Estimating Equations". Statistics in Medicine. 28 (4): 642-658. 10.1002/sim.3489.
DOI: https://doi.org/10.1002/sim.3489