Descriptive Statistics
From 20 February to 18 April 2020, 4791 adult patients with definitive diagnosis of COVID-19 were admitted to 15 university hospitals in Tehran, Iran, 249 patients were excluded due to their missing data. Of the 4542 hospitalized patients, who were between 18 and 97 years of age and had severe symptoms, 822 (18.09%) died during hospitalization, and 3720 (81.90%) patients recovered and discharged. The mean±SD and median [interquartile range] age of patients were, respectively, 55.55±16.84 and 56.00 [27.00] for the survived group, and 68.72±14.93 and 71.00 [21.25] for the deceased. 2473 (54.44%) patients were male. The demographic characteristics of the patients are displayed in Table 1.
The recruited patients had severe symptoms such as fever above 380 C, severe cough, shortness of breath, respiratory rate > 30 time/minute, hypoxia (PO2 < 93%), chest pain, and loss of consciousness. Criteria for discharge from the hospital was 72 hours without a fever and no need for an antipyretic drug, PO2 > 93%, and improvement in the clinical, and respiratory symptoms.
The mean±SD of hospitalization days for the survived and deceased groups were 4.76±4.77 and 6.45±7.21 days, respectively (Table 1). The treatment program was identical in all hospitals, and according to the COVD-19 treatment guidelines of Iran’s Ministry of Health, unless there were concurrent diseases requiring specific treatments.
Table 1: Demographic and hospitalization characteristics of the COVID-19 patients
As mentioned in the method section, before entering all the independent variables in the MLR model, a feature selection algorithm was used to define the most important predictors associated with the target variable. According to the feature selection results, apart from RBC, platelets, hemoglobin, hematocrit, PDW, and PCO2, all laboratory tests, age, gender, and the need for special care such as oxygen therapy and ventilation for survived patients entered the model as potential variables correlated with disease prognosis (Table 2).
Table 2: Predictors statistical characteristics in patients with COVID-19
MLR results
In this step, the SPSS modeler ignored predictors with a lot of missing value. Thus the amount of the missing value in the datasets reached under 15%. Accuracy and AUC index of the testing dataset as fitness criteria was 85.87% and 85.00, respectively. So the model has been fitted well, without overfitting problems [14].
As shown in Table 3 and based on the MLR model, age, ventilation, oxygen therapy, CRP, AST, ALT, ALP, CPK, serum calcium, serum sodium, serum potassium, and creatinine, have a significant relationship with the survival of COVID-19 patients.
According to the results, the group in which ventilation was required had a higher mortality rate (OR=.45, 95% CI: .38-.52). From the 222 patients that received this intervention, only 16.67% (37 people) survived. Regarding oxygen therapy, the results are the same (OR=.93, 95% CI: .88-.98), with a weaker correlation (Beta = -.08). Regarding laboratory test predictors, the results are as follows:
Negative association of levels of CRP (OR=.77, 95% CI: .71-.84), AST (OR=.54, 95% CI: .45-.64), ALP (OR=.79, 95% CI: .71-.87), CPK (OR=.78, 95% CI: .71-.86), sodium (OR=.82, 95% CI: .74-.91), potassium (OR=.87, 95% CI: .78-.98) and creatinine (OR=.92, 95% CI: .87-.99) as well as a positive relationship with ALT (OR=1.23, 95% CI: 1.07-1.42), and calcium (OR=1.38, 95% CI: 1.24-1.54), have been shown in the MLR model for surviving patients.
Table 3: Multiple logistic regression results (target variable is the final outcome (survived or deceased))
ULR results
In the following stage, significant predictors were altered into categorical variables in order to enter them in the ULR model and calculate the Relative Risk. Thresholds used for categorization of these predictors were their standard rages. This transformation enabled us to identify the critical values of laboratory tests, attributed to the patients' survival status (See appendix 1).
According to the ULR results, the age-related RR increased in older ages. The results showed that ages 40-49, 50-59, 60-69, 70-79 and above 80 are associated with elevated RR levels of 1.80 (CI: 1.13-2.87), 2.63 (CI: 1.71-4.02), 4.40 (CI: 2.92-6.63), 7.49 (CI: 5.01-11.19), and 13.85 (CI: 9.23-2.77) respectively. So, age is an essential COVID-19 risk factor that is not to be overlooked. Other significant risk factors are related to liver injury. These include ALT higher than the reference range 55-100 (RR=2.20, CI: 1.69-2.86) and above 100 (RR=2.67, CI = 2.86-5.34), AST range 40-100 IU/L (RR=2.51, CI: 2.04-3.10) and above 100 IU/L (RR=5.93, CI: 4.75-7.39), and ALP higher than 200 IU/L (RR=2.46, CI: 1.80-3.37). Serum electrolytes included sodium under 135 mEq/l (RR=1.69, CI: 1.35-2.11), and over 145 mEq/l (RR=7.24, CI: 5.07-1.33), potassium above 5.50 mEq/l (RR=7.53, CI: 4.15-13.64), and calcium under 8.50 mEq/l (RR=3.39, CI: 2.81-4.09) are related to increased mortality rate due to COVID-19. Also, elevated CPK range 307-600 IU/L (male) and 192-400 IU/L (female) (RR=2.73, CI: 2.12-3.53) and above 600 IU/L (male), and above 400 (female) (RR=4.41, CI: 3.40-5.71), CRP between 8-100 mg/l (RR=1.40, CI: .94-2.08) and more than 100 mg/l are (RR=3.22, CI: 1.99-5.20), and creatinine more than 1.5 (RR=6.37, CI: 5.30-7.66) were associated with COVID-19 mortality rate (MR) (Table 4).
Table 4: RR related to different ranges of the studied predictors
CHAID model: Identifying high risk patients
Based on the MLR analysis, of the 35 laboratory test parameters that were analyzed, 9 had a significant relationship with patients’ survival rate. These 9 parameters were then converted into categorical variables and entered into the CHAID node to develop a decision tree model.
According to the CHAID model, the patients were classified into five categories based their risk factors: low risk (under 20% MR), moderate risk 1 (20-30% MR), moderate risk 2 (30-40% MR), high risk (40-50% MR), and very high risk (more than 50% MR). These risk factors were separately defined for each age group (Figure 1).
In the ages under 40, sodium, CPK, AST, and creatinine (accuracy=96.25%, AUC=.63), in the age group of 40-59, AST, calcium, creatinine, and sodium (accuracy=88.57%, AUC=.75), within 60-69 year-old patients, AST, creatinine, CPK, ALP, and CRP (accuracy=82.70%, AUC=.71), in 70-79 age group, AST, creatinine, calcium, ALP, and CPK (accuracy=8.00%, AUC=.83), and for those above 80, creatinine, AST and sodium (accuracy=62.16%, AUC=.63) were associated with high-risk COVID-19 patients.
In the age under 40 (4.64% MR), with normal creatinine, an increase of AST ≥ 100 IU/L, with CPK ≥ 307 IU/L in male patients and CPK ≥ 192 IU/L in female patients were associated with 41.67% MR. So, this group is considered high-risk patients. Also, in this age group, we see 75.00% MR correlated with elevated creatinine levels ≥ 1.5 mg/l and an abnormal sodium level.
In the age group 40-59 (11.84% MR), high-risk group patients have been seen:
Patients with creatinine ≥ 1.5 mg/l that have 41.52% MR, patients with hypocalcemia, and AST ≥ 100 IU/L that have 56.08% MR, patients with creatinine levels ≥ 1.5 mg/l and hypocalcemia that have 46.67% MR, and patients with creatinine levels of ≥ 1.5 mg/l and AST ≥ 40 IU/L that have 59.55% MR, hence regarded as high-risk COVID-19 patients.
In the age group 60-69 (18.69% MR), elevated AST ≥ 100 IU/L with CPK ≥ 307 IU/L in male patients and CPK ≥ 192 IU/L in female patients were associated with 67.57% MR.
In the age group 70-79 (31.47% MR), six groups of high-risk patients have been seen:
Male patients with CPK ≥ 307 IU/L and female patients with CPK ≥ 192 IU/L have 5.00% MR.
In patients with AST level 40-99 IU/L and elevated ALP ≥ 200 IU/L, MR is 48.75%, and with elevated CPK ≥ 600 IU/L (male) and CPK ≥ 400 IU/L (female), is 61.54%. Also, in this age group, there is 52.63% MR when AST ≥ 100 IU/L, 67.19% MR when AST ≥ 100 IU/L and creatinine ≥ 1.5 mg/l, and 66.67% MR when AST ≥ 100 IU/L and there is hypocalcemia too.
In the age above 80 years (4.40% MR), MR, due to COVID-19, is 49.09% when AST ≥ 100 IU/L, 59.90% when creatinine ≥ 1.5 mg/l (43.59% when AST < 40 IU/L, and 71.56% when AST ≥ 40 IU/L), and 79.59% when creatinine ≥ 1.5 mg/l and there are abnormal sodium levels too. More information on this is presented in Figure 1.