Abstract

To construct a prediction model for all-cause mortality in elderly diabetic nephropathy (DN) patients, in this cohort study, the data of 511 DN patients aged ≥65 years were collected and the participants were divided into the training set (n = 358) and the testing set (n = 153). The median survival time of all participants was 2 years. The data in the training set were grouped into the survival group (n = 203) or the death group (n = 155). Variables with P ≤ 0.1 between the two groups were selected as preliminary predictors and involved into the multivariable logistic regression model and the covariables were gradually adjusted. The receiver operator characteristic (ROC), Kolmogorov-Smirnov (KS), and calibration curves were plotted for evaluating the predictive performance of the model. Internal validation of the performance of the model was verified in the testing set. The predictive values of the model were also conducted in terms of people with different genders and ages or accompanied with chronic kidney disease (CKD) or cardiovascular diseases (CVD), respectively. In total, 216 (42.27%) elderly DN patients were dead within 2 years. The prediction model for the 2-year mortality of elderly patients with DN was established based on length of stay (LOS), temperature, heart rate, peripheral oxygen saturation (SpO2), serum creatinine (Scr), red cell distribution width (RDW), the simplified acute physiology score-II (SAPS-II), hyperlipidemia, and the Chronic Kidney Disease Epidemiology Collaboration equation for estimated glomerular filtration rate (eGFR-CKD-EPI). The AUC of the model was 0.78 (95% CI: 0.73–0.83) in the training set and 0.72 (95% CI: 0.63–0.80) in the testing set. The AUC of the model was 0.78 (95% CI: 0.65–0.91) in females and 0.78 (95%CI: 0.68–0.88) in patients ≤75 years. The AUC of the model was 0.74 (95% CI: 0.64–0.84) in patients accompanied with CKD. The model had good predictive value for the mortality of elderly patients with DN within 2 years. In addition, the model showed good predictive values for female DN patients, DN patients ≤75 years, and DN patients accompanied with CKD.

1. Introduction

Diabetic nephropathy (DN) is a common microvascular complication of diabetes mellitus (DM) [1]. Approximately 30% of DM patients are diagnosed with renal complications including DN [2]. DN in patients can lead to end‐stage renal failure and disability, which is associated with high mortality all over the world [3]. DN patients tend to be elderly and may be associated with various complications, such as cerebrovascular, cardiovascular, peripheral vascular, connective tissue, liver, and chronic pulmonary diseases and tumors [4, 5]. DN is associated with higher mortality rates and worse clinical outcomes, which were largely due to the serious complications [6]. Therefore, predicting the all-cause mortality in DN patients was of great value for providing timely interventions in these patients and improving the outcomes of these patients.

Previously, various studies have explored the risk factors for the mortality in DN patients [79], but the risk of mortality could not be estimated based on the findings of these studies, as they did not form a prediction model. Currently, the model for predicting the mortality of DN patients was rare. In 2017, Sato et al. [10] established a prediction model for all-cause mortality in DN patients [10]. The model had an area under the curve (AUC) of 0.791, which had good predictive ability for the mortality of DN patients. Previously, multiple studies have indicated that prediction model based on combined variables might be better than those including only one variable [11]. The prediction model by Sato et al. [10] was focused on predialysis neutrophil-lymphocyte ratio, and validation was not performed to verify the performance of the model. Due to the poor prognosis of DN patients at old age [12], a suitable prediction model was required for the all-cause mortality in elderly DN patients to quickly identify those at high risk of mortality and provide timely treatments for these patients.

In this study, the purpose was to construct a prediction model for all-cause mortality in elderly DN patients. The predictors were screened out and included in the model. The internal validation was performed to evaluate the predictive value of the model. Subgroup analysis was also conducted in terms of gender and being complicated with chronic kidney disease (CKD) or cardiovascular diseases (CVD).

2. Methods

2.1. Study Population

In this cohort study, the data of 522 DN patients aged ≥65 years were derived from Medical Information Mart for Intensive Care (MIMIC-III) database. MIMIC-III database is an extensive and single-center database, constructed by Institutional Review Boards (IRB) of the Massachusetts Institute of Technology (Cambridge, MA, USA) and Beth Israel Deaconess Medical Center. It contained the data of over 50000 hospital patients admitted to intensive care units (ICUs) between 2001 and 2012 including the demographic details, admission and discharge times, dates of death, procedures such as dialysis, imaging studies, blood chemistry, hematology, urine analysis, microbiology test results, administration records of intravenous medications, medication orders, free text notes such as provider progress notes and hospital discharge summaries, and nurse-verified vital signs [13]. After excluding participants without the data on Sequential Organ Failure Assessment (SOFA) score, the simplified acute physiology score-II (SAPS-II), and temperature, 511 patients were finally involved in our study.

2.2. Potential Predictors

Potential predictors were analyzed in this study including gender, marital status (divorced, married, separated, single, widowed, or unknown), ethnicity (Asian, Black, Hispanic or Latino, White, others, or unknown), length of stay (LOS, day), age (years), respiratory rate (times/min), temperature (°C), heart rate (times/min), systolic blood pressure (SBP, mmHg), diastolic blood pressure (DBP, mmHg), mean arterial pressure (MAP, mmHg), peripheral oxygen saturation (SpO2, %), white blood cells (WBC, 103/μL), red blood cells (RBC, 103/μL), sodium (mEq/L), potassium (mEq/L), phosphate (mEq/L), calcium (mEq/L), magnesium (mEq/L), platelets (PLT, k/μL), lactate, international normalized ratio (INR), mean corpuscular volume (MCV, fl), glucose (mg/dL), serum creatinine (Scr, mg/dL), blood urea nitrogen (BUN, mg/dL), bicarbonate, hematocrit, hemoglobin, mean corpuscular hemoglobin concentration (MCHC, 10 g/L), red cell distribution width (RDW, %), chronic obstructive pulmonary disease (COPD, no or yes), atrial fibrillation (AF, no or yes), liver cirrhosis (no or yes), respiratory failure (no or yes), hyperlipidemia (no or yes), malignant cancer (no or yes), SAPS-II, SOFA score, insulin (no or yes), metformin (no or yes), survival time, the Chronic Kidney Disease Epidemiology Collaboration equation for estimated glomerular filtration rate (eGFR-CKD-EPI, mL/min/m2), the Modification of Diet in Renal Disease equation for estimated glomerular filtration rate (eGFR-MDRD, mL/min/m2), CVD (no or yes), CKD (no or yes), myocardial infarction (no or yes), hypertension (no or yes), and peripheral vascular disease (no or yes).

2.3. Outcome Variables

The outcome variable was the death of elderly DN patients within 2 years. The follow-up time was 10 years and the median survival time was 2 years.

2.4. Definitions of Variables

eGFR-MDRD = 175.0 × Scr  −1.154 × age−0.203 × 0.742 (if female) × 1.212 (if black); eGFR-CKD-EPI = 141 × min (Scr/κ, 1) α × max (Scr/κ, 1) − 1.029 × 0.993 age × 1.108 (if female) × 1.159 (if black). κ is 0.7 for females and 0.9 for males, α is −0.329 for females and −0.411 for males, min indicates the minimum of Scr/κ or 1, and max indicates the maximum of Scr/κ or 1. LOS is the length of stay in the ICUs.

2.5. Logistic Regression Model

Logistic regression is a classification method applied for binary or classification method generalizing logistic regression to multiclass problems multinomial outcome variables. It evaluates the associations between a dependent categorical outcome and one or more independent predictor variables, which provides predicted probabilities for each category [14] (1). The detailed formula of the logistic regression model is as follows:

2.6. Statistical Analysis

The normal distributed measurement data were expressed as mean ± standard deviation (mean ± SD), and comparisons between groups were subjected to independent-sample t-test. Nonnormal distributed data were described as M (Q1, Q3), and the Mann-Whitney U rank-sum test was used for comparing differences between groups. The enumeration data were displayed as n (%), and comparisons between groups were performed by χ2 test or Fisher’s exact probability method [15]. All the data were divided into the training set (n = 358) and the testing set (n = 153) at a ratio of 7 : 316. The prediction model was constructed in the training set and verified in the testing set. The data in the training set were grouped into the survival group (n = 203) or the death group (n = 155), and comparisons between the two groups were performed. Variables with P ≤ 0.1 were selected as preliminary predictors. The preliminarily screened predictors were then involved in the multivariable logistic regression model and the covariables were gradually adjusted. Subgroup analysis was conducted in male group and female group, CKD group and non-CKD group, CVD group and non-CVD group, age ≤75 years group, and age >75 years group, respectively. The area under the curve (AUC), Kolmogorov-Smirnov (KS), calibration curve, sensitivity, specificity, negative predictive value (NPV), positive predictive value (PPV), and accuracy were employed for evaluating the predictive performance of the model. A nomogram was also plotted to evaluate the possibility of mortality of elderly patients with DN. The confidence level was 0.05 and Python 3 was used for statistical analysis.

3. Results

3.1. Missing Value Manipulation and Sensitivity Analysis

The missing values of variables are shown in Supplementary Table 1. The missing data were manipulated via multiple interpolation using R mice. Sensitivity analysis was performed in the data before and after the manipulation. The results delineated that there was no statistical difference between the data before and after the manipulation, indicating that the data after manipulation could be used for further analysis.

3.2. Baseline Characteristics of Participants

In total, 522 DN patients aged ≥65 years from MIMIC-III were involved in our study. Participants without the data on SOFA score and SAPS-II (n = 9) and those without the data on temperature (n = 2) were excluded, and 511 patients were finally included. The detailed screen process is shown in Figure 1. Among them, 292 people were males, accounting for 57.14%. The median LOS was 2.6 days. The median age of all participants was 74.39 years. The median glucose level was 166 mg/dL. The median Scr level was 2.7 mg/dL. The median BUN was 45 mg/dL. The median survival time of all patients was 652.00 days. The median eGFR-CKD-EPI was 21.44 mL/min/m2 and the median eGFR-MDRD was 21.8 mL/min/m2. There were 389 patients accompanied with CVD, accounting for 76.13%, and 333 patients accompanied with CKD, accounting for 65.17%. The median survival time of all participants was 730 days and 216 people died within 2 years, accounting for 42.27%. The LOS in the survival group was shorter than that in the death group (2.15 days versus 3.24 days). The median survival time of the participants in the survival group was longer than that in the death group (730.00 days versus 61.50 days) (Table 1). The equilibrium test revealed that there was no significant difference between the data of participants in the training set and the testing set (Table 2).

LOS: length of stay, SBP: systolic blood pressure, DBP: diastolic blood pressure, MAP: mean arterial pressure, SpO2: peripheral oxygen saturation, WBC: white blood cells, RBC: red blood cells, INR: international normalized ratio, MCV: mean corpuscular volume, MCHC: mean corpuscular hemoglobin concentration, RDW: red cell distribution width, COPD: chronic obstructive pulmonary disease, AF: atrial fibrillation, eGFR-CKD-EPI: the Chronic Kidney Disease Epidemiology Collaboration equation for estimated glomerular filtration rate, eGFR-MDRD: the Modification of Diet in Renal Disease equation for estimated glomerular filtration rate, CKD: chronic kidney disease, CVD: cardiovascular diseases, SOFA: Sequential Organ Failure Assessment, SAPS-II: the simplified acute physiology score-II.

3.3. Comparisons between the Characteristics of Patients in the Survival Group and Death Group in the Training Set

The median LOS (2.15 days versus 3.01 days, Z = 3.734), age (73.59 years versus 76.03 years, Z = 1.770), INR (1.20 versus 1.30, Z = 2.767), Scr (2.30 mg/dL versus 2.90 mg/dL, Z = 2.100), BUN (43.00 mg/dL versus 50.00 mg/dL, Z = 2.447), SOFA score (5.00 versus 6.00, Z = 4.397), the average heart rate (80.12 times/min versus 85.81 times/min, t = −2.95), SpO2 (96.52 versus 97.42, t = −1.77), RBC (3.72 103/μL versus 3.58 103/μL, t = 1.85), SAPS-II (40.28 versus 45.77, t = −4.62), and the proportion of patients with respiratory failure (23.65% versus 33.55%, χ2 = 4.282) were lower in the survival group than in the death group. The median eGFR-MDRD (25.42 mL/min/m2 versus 20.41 mL/min/m2, Z = −2.266), eGFR-CKD-EPI (25.60 mL/min/m2 versus 19.68 mL/min/m2, Z = −2.705), the average temperature (36.59°C versus 36.41°C, t = 1.75), calcium (8.82 mEq/L versus 8.57 mEq/L, t = 2.48), hemoglobin (11.02 versus 10.61 t = 2.03), MCHC (33.05 10 g/L versus 32.59 10 g/L, t = 2.78), and the proration of patients with hyperlipidemia (60.59% versus 34.19%, χ2 = 4.282), CKD (70.44% versus 60.65%, χ2 = 3.771), diabetic retinopathy (21.18% versus 12.26%, χ2 = 4.888), and insulin use (94.58% versus 89.68%, χ2 = 3.031) in the survival group were higher than those in the death group. The proportion of patients with different marital status was statistically different between the survival group and the death group (χ2 = 10.722) (Table 3).

3.4. Predictors for Mortality of Elderly Patients with DN

Variables with P ≤ 0.1 in the survival group and the death group were included in the multivariable logistical analysis. Stepwise regression was applied to identify the predictors for mortality of elderly patients with DN within 2 years. As depicted in Table 4, LOS (OR = 1.10, 95% CI: 1.03–1.17), temperature (OR = 0.74, 95% CI: 0.63–0.88), heart rate (OR = 1.03, 95% CI: 1.01–1.04), SpO2 (OR = 1.06, 95% CI: 1.01–1.11), Scr (OR = 0.83, 95% CI: 0.69–0.98), RDW (OR = 1.25, 95% CI: 1.10–1.42), SAPS-II (OR = 1.02, 95% CI: 1.01–1.05), hyperlipidemia (OR = 0.43, 95% CI: 0.27–0.70), and eGFR-CKD-EPI (OR = 0.97, 95% CI: 0.94–0.99) were predictors associated with the risk of mortality in elderly patients with DN within 2 years. The final model was Log (p/1 − p) = 0.09 × LOS −  0.29 × temperature − 0.19 ×  creatinine + 0.03 × heart rate + 0.05 × SpO2 + 0.22 × RDW +  0.02 × SAPS-II-0.84 × hyperlipidemia − 0.03 × eGFR-CKD-EPI.

3.5. Predictive Value of the Model

According to the data in Table 5, for the model in the training set, the sensitivity was 0.85 (95% CI: 0.80–0.91), the specificity was 0.59 (95% CI: 0.52–0.65), the PPV was 0.61 (95% CI: 0.55–0.68), the NPV was 0.84 (95% CI: 0.78–0.90), the AUC was 0.78 (95% CI: 0.73–0.83), and the accuracy was 0.70 (95% CI: 0.65–0.75). The ROC, KS, and calibration curves in the training set are shown in Figure 2. For the model in the testing set, the sensitivity was 0.90 (95% CI: 0.83–0.98), the specificity was 0.47 (95% CI: 0.37–0.57), the PPV was 0.53 (95% CI: 0.43–0.62), the NPV was 0.88 (95% CI: 0.79–0.97), the AUC was 0.72 (95% CI: 0.63–0.80), the accuracy was 0.64 (95% CI: 0.56–0.72). The ROC, KS, and calibration curves in the testing set are exhibited in Figure 3. The nomogram was plotted and a sample was selected, which showed that the total score of the patient was 284, and the predicted mortality probability was 0.155, which was lower than the cut-off, 0.33 (Figure 4). The predicted outcome of the patient was survival, which was consistent with the actual outcome.

3.6. The Predictive Value of the Model concerning Different Subgroups
3.6.1. Gender

In the male group, the sensitivity was 0.90 (95% CI: 0.80–0.99), the specificity was 0.39 (95% CI: 0.25–0.52), the PPV was 0.54 (95% CI: 0.42–0.66), the NPV was 0.83 (95% CI: 0.67–0.98), the AUC was 0.66 (95% CI: 0.55–0.78), and the accuracy was 0.61 (95% CI: 0.51–0.72). In the female group, the sensitivity was 0.91 (95% CI: 0.79–1.00), the specificity was 0.56 (95% CI: 0.41–0.71), the PPV was 0.51 (95% CI: 0.36–0.67), the NPV was 0.92 (95% CI: 0.82–1.00), the AUC was 0.78 (95% CI: 0.65–0.91), and the accuracy was 0.68 (95% CI: 0.56–0.79) (Table 6).

3.6.2. Age

In patients >75 years group, the sensitivity was 0.88 (95% CI: 0.75–1.00), the specificity was 0.36 (95% CI: 0.22–0.50), the PPV was 0.43 (95% CI: 0.30–0.57), the NPV was 0.84 (95% CI: 0.68–1.00), the AUC was 0.65 (95% CI: 0.52–0.78), and the accuracy was 0.54 (95% CI: 0.43–0.66). In patients ≤75 years group, the sensitivity was 0.92 (95% CI: 0.83–1.00), the specificity was 0.57 (95% CI: 0.43–0.72), the PPV was 0.62 (95% CI: 0.49–0.75), the NPV was 0.90 (95% CI: 0.79–1.00), the AUC was 0.78 (95% CI: 0.68–0.88), and the accuracy was 0.72 (95% CI: 0.63–0.82) (Table 6).

3.6.3. Accompanied with CKD or Not

In patients accompanied with CKD group, the sensitivity was 0.90 (95% CI: 0.80–1.00), the specificity was 0.51 (95% CI: 0.39–0.63), the PPV was 0.47 (95% CI: 0.34–0.59), the NPV was 0.92 (95% CI: 0.83–1.00), the AUC was 0.74 (95% CI: 0.64–0.84), and the accuracy was 0.64 (95% CI: 0.54–0.73). In patients not complicated with CKD group, the sensitivity was 0.90 (95% CI: 0.79–1.00), the specificity was 0.37 (95% CI: 0.19–0.55), the PPV was 0.61 (95% CI: 0.47–0.76), the NPV was 0.77 (95% CI: 0.54–1.00), the AUC was 0.67 (95% CI: 0.52–0.82), and the accuracy was 0.65 (95% CI: 0.53–0.77) (Table 6).

3.6.4. Accompanied with CVD or Not

In patients accompanied with CVD group, the sensitivity was 0.90 (95% CI: 0.82–0.98), the specificity was 0.46 (95% CI: 0.34–0.58), the PPV was 0.57 (95% CI: 0.46–0.67), the NPV was 0.86 (95% CI: 0.75–0.97), the AUC was 0.71 (95% CI: 0.61–0.80), and the accuracy was 0.66 (95% CI: 0.57–0.74). In patients not accompanied with CVD group, the sensitivity was 0.89 (95% CI: 0.68–1.00), the specificity was 0.48 (95% CI: 0.28–0.68), the PPV was 0.38 (95% CI: 0.17–0.59), the NPV was 0.92 (95% CI: 0.78–1.00), the AUC was 0.71 (95% CI: 0.50–0.92), and the accuracy was 0.59 (95% CI: 0.42–0.75) (Table 6).

The comparisons of the AUCs of different subgroups delineated that the model had good predictive values for female DN patients, DN patients ≤75 years, and DN patients accompanied with CKD. The predictive values of the model for DN patients accompanied with CVD and DN patients not accompanied with CVD were similar (Figure 5).

4. Discussion

This study extracted the data of 511 DN patients aged ≥65 years and screened the predictors to establish a prediction model for the mortality of DN patients within 2 years. The results revealed that the model had good predictive ability for the mortality of DN patients within 2 years. Additionally, the predictive values of female DN patients, DN patients ≤75 years, DN patients accompanied with CKD, and patients with or without CVD were also good. The findings of our study might offer a tool for identifying DN patients with high risk of death within 2 years and the clinicians should provide timely interventions to those patients to improve their outcomes.

This study established a prediction model for the mortality of elderly DN patients within 2 years. In previous prediction models for the mortality of DN patents, many studies were focused on evaluating the risk of renal survival in DN patients [9, 16].Our study constructed a model and evaluated its predictive value for all-cause mortality in DN patients. DN patients were associated with various complications and the all-cause mortality of DN patients was high and should be brought to attention [17]. Sato et al. [10] established a prediction model for all-cause mortality in DN patients, but this model was based on only one laboratory index (predialysis neutrophil-lymphocyte ratio) and the sample size was small (n = 78). In addition, internal validation was also not performed to verify the performance of the model [10]. In our study, the prediction model was constructed based on the predictors including LOS, temperature, heart rate, SpO2, Scr, RDW, the simplified acute physiology score-II (SAPS-II), hyperlipidemia, and eGFR-CKD-EPI, which presented a better predictive ability compared to the model involving one predictor. The sample size in this study was larger than that in the previous study. Additionally, internal validation was performed and it was found that the predictive value of the model for the mortality of DN patents within 2 years was good. The prediction model in our study might provide a tool for the clinicians for quickly identifying DN patients with high risk of death and timely interventions should be provided in those patients for improving their outcomes. We also plotted a nomogram of the prediction model based on the results from the logistic regression. The nomogram can quickly and intuitively obtain the probability of mortality of each patient. Meanwhile, subgroup analysis was also conducted to evaluate the predictive values for patients with different gender, age, being accompanied with CKD or not, and being accompanied with CVD or not. The results revealed that the model had better predictive values for female DN patients, DN patients ≤75 years, and DN patients accompanied with CKD. The predictive values of the model for DN patients accompanied with CVD and DN patients not accompanied with CVD were similar. This indicated that the model might be more suitable for female DN patients, DN patients ≤75 years, and DN patients accompanied with CKD. These results suggested that the model could benefit specific patients with DN.

The impaired glomerular filtration rate (GFR) was regarded as a marker of DN in DM patients [18]. A previous meta-analysis revealed that the impaired GFR was an independent risk factor for progressive CKD, end-stage renal failure, and all-cause mortality in general population [19]. The eGFR-CKD-EPI is an extensively used equation for estimating GFR [20]. The decline of eGFR-CKD-EPI was associated with renal hyperfiltration and impaired GFR in DM patients [21]. These supported the results in our study, which revealed that the eGFR-CKD-EPI was a predictor for the mortality of DN patients within 2 years. Patients with rapid decline of eGFR-CKD-EPI should be brought to the forefront and special treatments should be provided to prevent the mortality of DN patients. DN was associated with higher Scr levels in patients, and high Scr levels indicated a declining renal function [22, 23]. This allied with the results in this study, which indicated that the Scr level was an important predictor for the mortality of elderly DN patients within 2 years. Clinicians should pay special attention to DN patients with high level of Scr. SpO2 is an index for oxygenation status of people and tissue hypoxia is an important contributor to diabetic complications [24]. Frequent abnormal blood oxygen in patients was reported to be associated with elevated inflammation in patients [25]. Herein, SpO2 was a predictor for the mortality of elderly DN patients within 2 years. In this study, RDW was another predictor for the mortality of elderly DN patients within 2 years. This was supported by several previous studies. Zhang et al. [26] identified that patients with DN were found to be with high level of RDW and RDW was associated with increased risk of progression to ESRD in patients with DN [26]. Another study also demonstrated that high level of RDW was an indicator of prognosis in DN patients and high level of RDW in T2D patients indicated a poor prognosis for DN [27]. SAPS-II is an indicator evaluating the outcomes of patients in ICUs and estimating their risk of mortality [28]. SAPS-II has good power to predict the deaths in ICU, which has been recommended for the identification and mortality prognostication of patients in ICUs [29]. In our study, SAPS-II was found to be a predictor for the mortality in ICU patients with DN. High-risk patients were associated with longer LOS in ICUs and with higher hospital mortality [30]. The prolonged LOS in ICUs has been reported to be a risk factor for infections, which might also increase the risk of death in patients [31]. These gave evidence to the findings in this study, showing that LOS in ICUs was a predictor for the mortality of DN patients in ICUs.

Several limitations existed in our study. Firstly, this study extracted the data from MIMIC-III database, which lacked several important variables including the medications of DN patients, as well as the control of blood glucose of the subjects, and these were closely associated with the outcomes of these patients. Secondly, external validation of the predictive value of the model was not performed. In the future, studies with large scale of sample size were required to validate the findings in our study. Currently, there were numerous machine learning algorithms that can be used for predicting the mortality of elderly patients with DN. Some recent studies have also used principal component analysis- (PCA-) firefly based deep learning model for predicting the occurrence or the detection of diabetic retinopathy [3234]. The predictive accuracy was evidently improved using these methods. Diabetic nephropathy and DN are common microvascular complications of diabetes mellitus. In our study, we only used logistic regression model, and, in the future, PCA-firefly based deep learning model might be applied in our further studies to improve the predictive ability for the mortality of DN patients and achieve a better tool for the clinicians to quickly and accurately identify those with high risk of death.

5. Conclusion

This study established a prediction model for the mortality of DN patients within 2 years based on LOS, temperature, heart rate, SpO2, Scr, RDW, SAPS-II, hyperlipidemia, and eGFR-CKD-EPI. The model had good predictive value for the mortality of elderly patients with DN within 2 years. In addition, the model showed good predictive values for female DN patients, DN patients ≤75 years, and DN patients accompanied with CKD.

Data Availability

The data used to support the findings of this study can be obtained from the corresponding author upon request.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Supplementary Materials

Supplementary Table 1: sensitivity analysis of the data before and after the manipulation of the missing value. (Supplementary Materials)