Nomogram model predicts the risk of visual impairment in diabetic retinopathy: a retrospective study
BMC Ophthalmology volume 22, Article number: 478 (2022)
To develop a model for predicting the risk of visual impairment in diabetic retinopathy (DR) by a nomogram.
Patients with DR who underwent both optical coherence tomography angiography (OCTA) and fundus fluorescein angiography (FFA) were retrospectively enrolled. FFA was conducted for DR staging, swept-source optical coherence tomography (SS-OCT) of the macula and 3*3-mm blood flow imaging by OCTA to observe retinal structure and blood flow parameters. We defined a logarithm of the minimum angle of resolution visual acuity (LogMAR VA) ≥0.5 as visual impairment, and the characteristics correlated with VA were screened using binary logistic regression. The selected factors were then entered into a multivariate binary stepwise regression, and a nomogram was developed to predict visual impairment risk. Finally, the model was validated using the area under the receiver operating characteristic (ROC) curve (AUC), calibration plots, decision curve analysis (DCA), and clinical impact curve (CIC).
A total of 29 parameters were included in the analysis, and 13 characteristics were used to develop a nomogram model. Finally, diabetic macular ischaemia (DMI) grading, disorganization of the retinal inner layers (DRIL), outer layer disruption, and the vessel density of choriocapillaris layer inferior (SubVD) were found to be statistically significant (P < 0.05). The model was found to have good accuracy based on the ROC (AUC = 0.931) and calibration curves (C-index = 0.930). The DCA showed that risk threshold probabilities in the (3–91%) interval models can be used to guide clinical practice, and the proportion of people at risk at each threshold probability is illustrated by the CIC.
The nomogram model for predicting visual impairment in DR patients demonstrated good accuracy and utility, and it can be used to guide clinical practice.
Chinese Clinical Trial Registry, ChiCTR2200059835. Registered 12 May 2022,
Diabetes is one of the fastest-growing chronic diseases in the world , and diabetic retinopathy (DR) is one of the major causes of vision loss in diabetic patients . Studies have been performed on early risk screening for type 1 diabetes , type 2 diabetes  and DR , and the initial results will be progressively applied to clinical practice . Currently, fundus fluorescein angiography (FFA) is the gold standard for the diagnostic grading of DR , and swept-source optical coherence tomography (SS-OCT) can construct clear two-dimensional cross-sectional images of the retina . The novel recently developed optical coherence tomography angiography (OCTA) procedure is very sensitive to retinal microvascular changes, and it can yield clear images among some patients with poorly defined refractive media (cataract, vitreous opacities) .
It has been shown that DR staging  and retinal structural abnormalities (disorganization of the retinal inner layers (DRIL) , diabetic macular ischaemia (DMI) , diabetic macular oedema (DME) , outer layer disruption , structural changes in the foveal avascular zone (FAZ) , etc.) are correlated with visual impairment in DR patients, but no studies have revealed the magnitude of the influence of various correlated factors in the formation of visual impairment in DR patients. Furthermore, no studies have integrated these risk factors to construct an efficient, accurate, simple and intuitive model for predicting the risk of visual impairment. A nomogram would be useful to achieve this goal. Based on the reported risk factors for visual acuity (VA) in DR patients and some new retinal structure indicators, this study screened out risk factors that were strongly correlated with VA through regression analysis and creatively integrated and quantified the impact of the risk factors for visual impairment in DR patients using a nomogram . By scoring each factor, a convenient and practical risk prediction model for visual impairment was constructed to guide clinical decision-making.
This research was implemented in accordance with the requirements of the Declaration of Helsinki, and the protocol was approved by the Ethics Committee of the Second Affiliated Hospital of Army Medical University. This study was registered in the Clinical Trials Registry (ChiCTR2200059835). A total of 252 eyes in 133 patients with DR who attended the ophthalmology department of the Second Affiliated Hospital of the Army Medical University were examined from August 2020 to January 2022. The inclusion criteria were as follows: 1. patients of either sex, ages ≥18 years; 2. patients diagnosed with DR according to the International Clinical Classification Criteria for Diabetic Retinopathy (2002) ; and 3. clear images were obtained for patients who had undergone OCTA and FFA within 1 week. The exclusion criteria were as follows: 1. patients with other fundus retinopathies in addition to DR (age-related macular degeneration, polypoid chorioretinopathy, optic nerve atrophy, etc.); 2. patients with refractive media clouding (cataracts affecting vision, vitreous haemorrhage) for which clear images cannot be obtained; and 3. patients with other ophthalmic diseases affecting vision. The exclusion process and reasons for all exclusions are shown in Online Additional file 1.
The patient’s best-corrected visual acuity (BCVA) was assessed using the International Standard Visual Acuity Scale and converted to logarithm of the minimum angle of resolution (LogMAR) VA. Slit lamp, intraocular pressure (IOP) and fundus photograph were routinely examined, and the patient’s duration of diabetes at this visit, comorbidities, and examination test results were recorded.
OCTA and OCT image acquisition and processing
The OCTA instrument (DRI OCT Triton; Topcon Inc., Tokyo, Japan) we utilize used the latest technology. The device type was SS-OCT, which was based on the optical coherence tomography angiography ratio analysis (OCTARA) algorithm, the wavelength was 1050 nm, the acquisition speed was 100 kHz, and the axial and lateral resolution of the tissue was 7 μm and 20 μm, respectively. The capillary plexus was automatically segmented as follows (IMAGEnet6): the superficial capillary plexus (SCP) was located 3 mm below the internal limiting membrane (ILM) to 15 mm below the junction of the inner plexiform layer (IPL) and inner nuclear layer (INL); the deep capillary plexus (DCP) was located 15 mm below the IPL/INL to 70 mm below the INL, and the choriocapillaris layer was located Bruch Membrane (BM) to 20 mm below the BM . The following images were excluded: (1) quality scores < 45; (2) blurred images (inability to distinguish capillaries from background signal); (3) artefacts (white lines, vessel displacement) due to movement and blinking; and (4) images that did not distinguish between DCP and SCP. SS-OCT of the macula was performed on the patient and analysed for DRIL, DME occurrence, and outer layer disruption of the retina (Online Additional file 2). A 3*3-mm (each b-scan includes 320 A-scans for a total of 320 b-scans) flow imaging was performed from SCP, which was used to evaluate DMI grading in artificial means (Online Additional file 3) , and foveal avascular zone (FAZ) (d = 1 mm) and paracentral fovea (d = 0.75 mm) vessel density (VD) origin in DCP and SCP were analysed using IMAGEnet6 software. FAZ area, perimeter , circularity  (Online Additional file 4), and nonperfusion area (NPA) within 3 mm*3 mm  of DCP and SCP were analysed using ImageJ software (Online Additional file 5) (Table 1).
FFA was performed on the patient using a Heidelberg (SPECTRALIS HRA), arm retinal time (ART) was recorded and DR staging was performed on the affected eye .
Manual outlining of the images was performed using ImageJ (version 1.51a), and the grade was assessed by two blinded specialists (Chen and Zou). The specialists showed good agreement for manual measurement indicators (Online Additional file 6), and disagreements were resolved via discussion.
SPSS 26.0 (IBM Corporation, New York, USA) and the ‘rms’ and ‘rmda’ packages of RStudio (version 4.1.3; https://www.R-project.org) software were used for statistical analysis. The normality of the measurement data was determined using the Kolmogorov–Smirnov test, and those that did not satisfy the normal distribution were expressed as medians (quartiles). The enumeration data were presented as frequencies (percentages). LogMAR VA was transformed into binary outcomes with or without impairment according to the cut-off value criterion (Standard: Blindness and Visual Impairment Criteria (International Classification of Diseases, WHO, 2009) defining VA ≥ 0.5 as having impaired vision and the opposite is not) , and each parameter was subjected to univariate logistic regression analysis and validated using the area under the receiver operating characteristic (ROC) curve (AUC). The parameters that were strongly correlated with VA were subsequently entered into a multivariate logistic regression (LR) in a stepwise manner to construct the nomogram. The model accuracy and fit were assessed using ROC and calibration curves, and decision curve analysis (DCA)  was conducted to assess the rate of the benefit of the model to patients. Clinical impact curve (CIC) was used to stratify risk proportions for each threshold probability in the model . P < 0.05 was considered to indicate a statistically significant difference.
Of the 133 patients, the median age was 56 years (50–61 years), the median duration of diabetes was 10 years (7–16 years), the median HbA1c value was 8.6% (7.3–10.5%), hypertension was observed in a total of 76 patients (57.14%), and diabetic nephropathy was observed in a total of 77 patients (59.69%). Table 2 showed the underlying clinical characteristics of the included patients. There were 252 eyes with nonproliferative diabetic retinopathy (NPDR) lesions, 49 (19.44%) eyes with mild DR, 89 (35.32%) eyes with moderate DR, 53 (21.03%) eyes with severe DR, and 61 (24.21%) eyes with proliferative diabetic retinopathy PDR.
Development of the nomogram models
A total of 44 eyes with impairment and 206 eyes without impairment were examined, indicating an impairment rate of 17.5%. First, each risk factor was analysed by univariate binary regression (Table 3), and risk factors with statistical significance (P < 0.05) and AUC > 0.650 in both binary regression and ROC curve analysis were screened. The following factors were included in the regression model: DR, DMI, DME, DRIL, outer layer disruption, NPA (SCP, SCP + DCP), VD (superior of SCP, centre and superior of DCP, inferior of choriocapillaris), FAZ (area and perimeter of SCP), and systemic factors such as sex, age, HbA1c, EGFR, TG, TCH, duration of diabetes, and the presence of hypertension and diabetic nephropathy. The model was statistically significant (P < 0.001) (Table 4).
Validation of nomogram models
The accuracy of the model in determining the presence of visual impairment was 91.2% at the optimal cut-off value (P = 0.187, P denotes the model-derived probability), with a sensitivity of 93.2% and a specificity of 86.5%. The Youden index was maximized, and DMI, Sub VD, DRIL, and Outer Layer Disruption were statistically significant (P < 0.05) factors in the model. Thus, these independent risk factors were included in the nomogram (Fig. 1). The ROC indicated that the nomogram had high discrimination (AUC = 0.931 (0.889–0.973), P < 0.05) (Fig. 2). Calibration curves (C-index = 0.930, P < 0.05) (Fig. 3) revealed good model fit consistency, and DCA illustrated that if the visual impairment risk threshold probabilities ranged from 3 to 91%, patients would significantly benefit from using this study’s model to aid decision-making. Finally, the CIC displayed the proportional number of people at risk of visual impairment at each threshold probability (Fig. 4).
Nomograms are considered reliable and practical predictive tools, capable of generating individual probabilities of clinical events by integrating multiple prognostic influences , and quantifying risk , thereby satisfying our desire for integrated biological and clinical models, enabling the need for personalized medicine, aiding better clinical decision-making through a user-friendly digital interface and more generalized conclusions . This study concluded that DRIL and outer layer disruption (provided by OCT images of the macula), DMI and SubVD (provided by 3*3-mm blood flow images) were strongly correlated with visual impairment in patients with DR. The AUC showed that the nomogram shows better discrimination than the use of individual risk factors (Fig. 2). The visualization of influencing factors using the nomogram reveals the magnitude of influence for each risk factor for the development of visual impairment, thereby facilitating clinical discrimination of indicators that require more attention. The DCA showed that patients with risk threshold probabilities in the range of 3–91% had a higher net benefit than those in the intervention for the all patients’ scenario or the ‘no intervention at all’ scenario . These findings showed that the model had a wide range of clinical applications. The CIC showed the number of true positives at risk of visual impairment at different threshold probabilities and the number predicted by the model, providing a reference for clinicians to determine the condition.
As an example of our findings, an affected eye with DRIL = 1 (21 points), DMI = 3 (31 points), Outer Layer Disruption =1 (15 points), subVD = 50% (20 points) would have a total score of 87 points, corresponding to a risk probability of substandard VA (LogMAR VA ≥ 0.5) of approximately 90% (Fig. 1). Even if patients currently had a fair vision on examination, there was a high risk that their vision would decline later; assuming a threshold probability of 20%, the DCA showed a net benefit rate of 13% in the vertical coordinate at this point, indicating that the model screened 13 additional high-risk patients for noncompliance per 100 DR patients tested for vision ((net benefit of the model - net benefit of all treatments)/(threshold probability/(1-threshold probability)) × 100), without increasing the number of false-positives . CIC showed that the true number of DR patients with substandard vision was approximately 170 (out of 1000 patients) compared to the model’s predicted number of at-risk patients (260), thus indicating that the clear results are useful in clinical practice (Fig. 4). Based on the model, clinicians can implement early medical interventions, such as changing treatment regimens or increasing the frequency of follow-up visits, to reduce the risk of vision loss in DR patients, which is important for disease management in high-risk DR patients. This model can also be used for other clinical applications, such as among DR patients who need to be treated with cataract or vitrectomy surgery. If OCT/OCTA images of the macula are available, using this study’s prediction model, a general judgement can be made about the patient’s postoperative VA risk profile, which will assist in surgical planning and preoperative conversations.
The DRIL in this study was heavily represented in the model and may represent a disruption of the visual conduction pathway in the inner retina that greatly affects visual acuity in DR patients . Its baseline length and degree of change over time have been shown to be correlated with VA . However, the severity of DRIL has not yet been graded . The relationship with visual acuity needs to continue to be explored, and the DRIL grading in the model appeared to be grade 1. The higher score than grade 2 in our study may be related to the fact that fewer eyes (1.59%) were grade 2, for reasons that will have to be explained by further research. Second, DaCosta indicated that visual acuity in DR patients gradually decreased with increasing severity of DMI  and that the incidence of DMI increased with increasing severity of DR , which was consistent with our study (Online Additional file 7). However, research on DMI is still in its infancy, and there is no definition of DMI characteristics with or without visual threat . Outer layer disruption, which can be used to assess retinal photoreceptor function, was closely associated with VA  and was shown in this study to be an important risk factor for VA in patients with DR. Interestingly, a study suggested that microvascular changes in DR are associated with early VA loss ,and our research also showed SubVD (provided by 3*3-mm blood flow images) were strongly correlated with visual impairment in patients with DR.
DR is still one of the leading causes of visual impairment, and the development of irreversible damage to vision should not be the ideal endpoint for all predictors but should be preceded by identification of the disease and timely intervention. The relationship between retinal changes and VA in DR patients needs to be further explored, and indicators that can more sensitively reflect visual impairment in DR patients are important for early screening and prognosis of the disease. In this study, we used multivariate stepwise regression to construct a pooled nomogram model for predicting the risk of visual impairment in DR patients based on reported risk factors for VA and new retinal structural parameters. To broaden the application of the model, we validated the feasibility and superiority of this integrated risk factor approach using VA impairment grading criteria and a multifaceted validation model using ROC, calibration plots, DCA and CIC, and the results showed that both models performed well. The entire construction and validation process was rigorous and complete, thereby providing an accurate and intuitive nomogram model for predicting the risk of visual impairment in patients with DR. This nomogram can be used to guide clinical practice.
There were also some limitations to this study. First, this was a retrospective study, and the sample size was limited. Second, DRIL was graded as occurring on one side of the central recess and both sides, and the actual length was not measured, which may also explain why the grade 2 score in the model was less than the grade 1 score. Third, this study only used data from a single visit and did not follow up on changes in patient VA; thus, our model can only predict the current risk of visual impairment and not changes in VA at a future time. Fourth, the model was not validated using an external case group. Therefore, prospective studies with larger samples and external validation are necessary to improve the model in future studies.
In summary, the nomogram prediction model has demonstrated greater accuracy and utility in integrating risk factors related to VA in patients with DR. The DMI, DRIL, outer layer disruption, and SubVD are independent risk factors in the nomogram model for predicting visual impairment risk in patients with DR.
Availability of data and materials
The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request.
optical coherence tomography angiography
fundus fluorescein angiography
swept-source optical coherence tomography
area under the curve
receiver operating characteristic
decision curve analysis
clinical impact curves
- LogMAR VA:
logarithm of the minimum angle of resolution visual acuity
diabetic macular ischemia
diabetic macular edema
disorder retinal inner layer
- Outer layer disruption:
retinal pigment epithelium or ellipsoid zone rupture
superficial capillary plexus
deep capillary plexus
no perfusion area
foveal avascular zone
vessel density of choriocapillaris layer inferior
arm retinal time
ultrasound of cervical vascular
high density lipoprotein cholesterol
low-density lipoprotein cholesterol
estimated glomerular filtration rate
Bhanushali D, Anegondi N, Gadde SG, Srinivasan P, Chidambara L, Yadav NK, et al. Linking retinal microvasculature features with severity of diabetic retinopathy using optical coherence tomography angiography. Invest Ophthalmol Vis Sci. 2016;57(9):Oct519–25.
Cheung N, Mitchell P, Wong TY. Diabetic retinopathy. Lancet. 2010;376(9735):124–36.
Kang EY, Lo FS, Wang JP, Yeh LK, Wu AL, Tseng YJ, et al. Nomogram for prediction of non-proliferative diabetic retinopathy in juvenile-onset type 1 diabetes: a cohort study in an Asian population. Sci Rep. 2018;8(1):12164.
Liu Y, Yang J, Tao L, Lv H, Jiang X, Zhang M, et al. Risk factors of diabetic retinopathy and sight-threatening diabetic retinopathy: a cross-sectional study of 13 473 patients with type 2 diabetes mellitus in mainland China. BMJ Open. 2017;7(9):e016280.
Mo R, Shi R, Hu Y, Hu F. Nomogram-based prediction of the risk of diabetic retinopathy: a retrospective study. J Diabetes Res. 2020;2020:7261047.
Babenko B, Mitani A, Traynis I, Kitade N, Singh P, Maa AY, et al. Detection of signs of disease in external photographs of the eyes via deep learning. Nat Biomed Eng. 2022:1–14. https://doi.org/10.1038/s41551-022-00867-5.
The effect of intensive diabetes treatment on the progression of diabetic retinopathy in insulin-dependent diabetes mellitus. The Diabetes Control and Complications Trial. Arch Ophthalmol. 1995;113(1):36–51. https://doi.org/10.1001/archopht.1995.01100010038019.
Youngquist RC, Carr S, Davies DE. Optical coherence-domain reflectometry: a new optical evaluation technique. Opt Lett. 1987;12(3):158–60.
Schwartz DM, Fingler J, Kim DY, Zawadzki RJ, Morse LS, Park SS, et al. Phase-variance optical coherence tomography: a technique for noninvasive angiography. Ophthalmology. 2014;121(1):180–7.
Sun JK, Radwan SH, Soliman AZ, Lammer J, Lin MM, Prager SG, et al. Neural retinal disorganization as a robust marker of visual acuity in current and resolved diabetic macular edema. Diabetes. 2015;64(7):2560–70.
Cheung CMG, Fawzi A, Teo KY, Fukuyama H, Sen S, Tsai WS, Sivaprasad S. Diabetic macular ischaemia- a new therapeutic target? Prog Retin Eye Res. 2022;89:101033. https://doi.org/10.1016/j.preteyeres.2021.101033. Epub 2021 Dec 11.
Campochiaro PA, Wykoff CC, Shapiro H, Rubio RG, Ehrlich JS. Neutralization of vascular endothelial growth factor slows progression of retinal nonperfusion in patients with diabetic macular edema. Ophthalmology. 2014;121(9):1783–9.
Das R, Spence G, Hogg RE, Stevenson M, Chakravarthy U. Disorganization of inner retina and outer retinal morphology in diabetic macular edema. JAMA Ophthalmol. 2018;136(2):202–8.
Balaratnasingam C, Inoue M, Ahn S, McCann J, Dhrami-Gavazi E, Yannuzzi LA, et al. Visual acuity is correlated with the area of the Foveal avascular zone in diabetic retinopathy and retinal vein occlusion. Ophthalmology. 2016;123(11):2352–67.
Balachandran VP, Gonen M, Smith JJ, DeMatteo RP. Nomograms in oncology: more than meets the eye. Lancet Oncol. 2015;16(4):e173–80.
Wu L, Fernandez-Loaiza P, Sauma J, Hernandez-Bogantes E, Masis M. Classification of diabetic retinopathy and diabetic macular edema. World J Diabetes. 2013;4(6):290–4.
Stanga PE, Tsamis E, Papayannis A, Stringa F, Cole T, Jalil A. Swept-source optical coherence tomography Angio™ (Topcon Corp, Japan): technology review. Dev Ophthalmol. 2016;56:13–7.
Soares M, Neves C, Marques IP, Pires I, Schwartz C, Costa M, et al. Comparison of diabetic retinopathy classification using fluorescein angiography and optical coherence tomography angiography. Br J Ophthalmol. 2017;101(1):62–8.
La Mantia A, Kurt RA, Mejor S, Egan CA, Tufail A, Keane PA, et al. Comparing fundus fluorescein angiography and swept-source optical coherence tomography angiography in the evaluation of diabetic macular perfusion. Retina. 2019;39(5):926–37.
Migicovsky Z, Harris ZN, Klein LL, Li M, McDermaid A, Chitwood DH, et al. Rootstock effects on scion phenotypes in a 'Chambourcin' experimental vineyard. Hortic Res. 2019;6:64.
Sonoda S, Sakamoto T, Yamashita T, Shirasawa M, Uchino E, Terasaki H, et al. Choroidal structure in normal eyes and after photodynamic therapy determined by binarization of optical coherence tomographic images. Invest Ophthalmol Vis Sci. 2014;55(6):3893–9.
ETDRS Research Group Investigators. Classification of diabetic retinopathy from fluorescein angiograms. ETDRS report number 11. Early treatment diabetic retinopathy study research group. Ophthalmology. 1991;98(5 Suppl):807–22.
Pascolini D, Mariotti SP. Global estimates of visual impairment: 2010. Br J Ophthalmol. 2012;96(5):614–8.
Fitzgerald M, Saville BR, Lewis RJ. Decision curve analysis. Jama. 2015;313(4):409–10.
Chen M, Li Z, Yan Z, Ge S, Zhang Y, Yang H, et al. Predicting neurological deterioration after moderate traumatic brain injury: development and validation of a prediction model based on data collected on admission. J Neurotrauma. 2022;39(5–6):371–8.
Liang G, Chen X, Zha X, Zhang F. A Nomogram to improve predictability of small-incision Lenticule extraction surgery. Med Sci Monit. 2017;23:5168–75.
Gold JS, Gönen M, Gutiérrez A, Broto JM, García-del-Muro X, Smyrk TC, et al. Development and validation of a prognostic nomogram for recurrence-free survival after complete surgical resection of localised primary gastrointestinal stromal tumour: a retrospective analysis. Lancet Oncol. 2009;10(11):1045–52.
Vickers AJ, Elkin EB. Decision curve analysis: a novel method for evaluating prediction models. Med Decis Mak. 2006;26(6):565–74.
Radwan SH, Soliman AZ, Tokarev J, Zhang L, van Kuijk FJ, Koozekanani DD. Association of Disorganization of retinal inner layers with vision after resolution of center-involved diabetic macular edema. JAMA Ophthalmol. 2015;133(7):820–5.
Sun JK, Lin MM, Lammer J, Prager S, Sarangi R, Silva PS, et al. Disorganization of the retinal inner layers as a predictor of visual acuity in eyes with center-involved diabetic macular edema. JAMA Ophthalmol. 2014;132(11):1309–16.
DaCosta J, Bhatia D, Talks J. The use of optical coherence tomography angiography and optical coherence tomography to predict visual acuity in diabetic retinopathy. Eye (Lond). 2020;34(5):942–7.
Sun Z, Tang F, Wong R, Lok J, Szeto SKH, Chan JCK, et al. OCT angiography metrics predict progression of diabetic retinopathy and development of diabetic macular edema: a prospective study. Ophthalmology. 2019;126(12):1675–84.
Thanks to all the patients for their participation and the researchers from the Second Affiliated Hospital of the Army Medical University for their contribution to this study.
This study was supported by scientific research funds from the Department of Ophthalmology, Second Affiliated Hospital of Army Medical University.
Ethics approval and consent to participate
The study is following the principles of the Declaration of Helsinki, was approved by the Medical Ethics Committee of the Second Affiliated Hospital of the Army Medical University and was audited by the China Clinical Trials Registry (ChiCTR2200059835). Since this project is a retrospective study, written informed consent was waived with the consent of the ethics committee.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Zhao, Y., Yu, R., Sun, C. et al. Nomogram model predicts the risk of visual impairment in diabetic retinopathy: a retrospective study. BMC Ophthalmol 22, 478 (2022). https://doi.org/10.1186/s12886-022-02710-6
- Visual impairment
- Decision curve analysis
- Clinical impact curve
- Optical coherence tomography angiography
- Diabetic retinopathy