Next Article in Journal
Enhancing Breast Cancer Detection through Advanced AI-Driven Ultrasound Technology: A Comprehensive Evaluation of Vis-BUS
Next Article in Special Issue
Novel Treatment for Pre-XDR Tuberculosis Linked to a Lethal Case of Acute Myocarditis
Previous Article in Journal
Contemporary Diagnostics of Cardiac Sarcoidosis: The Importance of Multimodality Imaging
Previous Article in Special Issue
A Review of Current Evidence for the Use of Steroids in the Medical Intensive Care Unit
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Feature Identification Using Interpretability Machine Learning Predicting Risk Factors for Disease Severity of In-Patients with COVID-19 in South Florida

1
Christine E. Lynn College of Nursing, Florida Atlantic University, Boca Raton, FL 33431, USA
2
College of Engineering & Computer Science, Florida Atlantic University, Boca Raton, FL 33431, USA
3
Memorial Healthcare System, Hollywood, FL 33021, USA
*
Author to whom correspondence should be addressed.
Submission received: 12 July 2024 / Revised: 16 August 2024 / Accepted: 21 August 2024 / Published: 26 August 2024
(This article belongs to the Special Issue Pulmonary Disease: Diagnosis and Management)

Abstract

:
Objective: The objective of the study was to establish an AI-driven decision support system by identifying the most important features in the severity of disease for Intensive Care Unit (ICU) with Mechanical Ventilation (MV) requirement, ICU, and InterMediate Care Unit (IMCU) admission for hospitalized patients with COVID-19 in South Florida. The features implicated in the risk factors identified by the model interpretability can be used to forecast treatment plans faster before critical conditions exacerbate. Methods: We analyzed eHR data from 5371 patients diagnosed with COVID-19 from South Florida Memorial Healthcare Systems admitted between March 2020 and January 2021 to predict the need for ICU with MV, ICU, and IMCU admission. A Random Forest classifier was trained on patients’ data augmented by SMOTE, collected at hospital admission. We then compared the importance of features utilizing different model interpretability analyses, such as SHAP, MDI, and Permutation Importance. Results: The models for ICU with MV, ICU, and IMCU admission identified the following factors overlapping as the most important predictors among the three outcomes: age, race, sex, BMI, diarrhea, diabetes, hypertension, early stages of kidney disease, and pneumonia. It was observed that individuals over 65 years (‘older adults’), males, current smokers, and BMI classified as ‘overweight’ and ‘obese’ were at greater risk of severity of illness. The severity was intensified by the co-occurrence of two interacting features (e.g., diarrhea and diabetes). Conclusions: The top features identified by the models’ interpretability were from the ‘sociodemographic characteristics’, ‘pre-hospital comorbidities’, and ‘medications’ categories. However, ‘pre-hospital comorbidities’ played a vital role in different critical conditions. In addition to individual feature importance, the feature interactions also provide crucial information for predicting the most likely outcome of patients’ conditions when urgent treatment plans are needed during the surge of patients during the pandemic.

1. Introduction

The COVID-19 pandemic has significantly affected global health and economies, overwhelming hospitals and straining limited resources [1]. As of April 2024, there have been over 700 million confirmed cases of COVID-19 and over 7 million deaths worldwide, with Florida alone accounting for almost 1.14% of all confirmed cases and approximately 1.36% of worldwide deaths [2]. The immense strain on healthcare systems, particularly in South Florida, necessitated the development of a predictive model for Intensive Care Unit (ICU) with Mechanical Ventilation (MV) requirement, ICU, or InterMediate Care Unit (IMCU) admission to effectively allocate resources, optimize patient care, and improve outcomes for patients with COVID-19. All patients who required MVs were admitted to the ICU; however, not all patients in the ICU necessarily needed MV. For ease of presentation, we shall refer to ICU patients with MV as only MV for brevity.
The COVID-19 pandemic has emphasized the critical need for optimizing resource allocation in healthcare systems. The demand for ICU beds surged due to COVID-19, with estimates of ICU admission rates ranging from 5% to 30% [3,4]. This increased the demand for critical care resources, including ventilators, posed challenges, and led to shortages during the peak of the pandemic [5]. Accurately predicting patients at a higher risk of requiring intensive care interventions allows healthcare providers to proactively allocate resources such as ICU beds, ventilators, and specialized staff to the patients who need them the most [6]. This targeted resource allocation ensures that critical care resources are utilized effectively and judiciously, maximizing their impact on patient outcomes.
In addition, predicting IMCU admission allows for proper triage and prevents the underutilization or overburdening of ICU resources. The early identification of patients requiring IMCU admission enables healthcare providers to intervene promptly and provide appropriate care. Patients in the IMCU may require specialized monitoring, non-invasive ventilation, or other interventions to manage their respiratory status or other critical care needs. Not all patients require the care provided in an ICU, but they may still benefit from closer monitoring and specialized interventions available in the IMCU. This improves care coordination, patient flow, and cost-effective resource management. It allows hospitals to anticipate the number of patients requiring intensive care, plan staffing and equipment needs, as well as coordinate care across different healthcare facilities. By adopting a proactive approach, healthcare systems can better manage an influx of patients, maintain quality care, and ensure that critical care is provided to those in need [7]. Timely, individualized interventions can also help prevent disease progression, reduce complications, and improve patient outcomes [8].
Given the demands for critical decisions in the healthcare system, clinicians can also benefit from AI-driven decision support systems for deciding optimal treatment plans for hospitalized patients with COVID-19 when there is an urgent need for decision making. AI-driven decision support systems can provide insights into disease severity, prognosis, and enable healthcare providers to better communicate disease diagnoses to patients and their families, make informed decisions about end-of-life care, and allocate resources appropriately based on the patient’s likelihood of recovery [9]. AI has been successfully deployed in clinical settings to aid in many healthcare decisions [10], such as detecting diseases [11], risk assessments [12], and personalized outcomes [13,14].

1.1. An Overview of Machine Learning Studies on COVID-19 Disease Severity

The literature on COVID-19 disease severity has undergone thorough analyses using various predictive methods. Machine Learning (ML) has been used to develop scoring tools based on essential features to measure COVID-19 disease severity [15,16,17] and to predict the disease severity scores of patients presented by a standardized severity scale, such as the National Early Warning Score 2 (NEWS2) [18]. Thus, ML tools can learn to predict outcomes based on a known severity metric or establish a new severity scale that improves the known risk assessment methods by incorporating novel features not included in the standardized scales [17]. Furthermore, past studies have focused on predicting the severity conditions of in-patients, such as MV requirements [18,19,20,21], days spent in ICU [20,21,22,23,24], whether rehospitalization is necessary for recurring health problems [15,25], and forecasting the need for other therapeutic interventions [24].
The studies have used various datasets to train their models; some relied solely on lab biomarkers [24,26], while others utilized medical histories like electronic Health Records (eHRs) and demographic information [21], and a few combined both data [19,23,25]. Additionally, some studies incorporated features from medical imaging, eHRs, and laboratory tests [22,27,28]. It is worth noting that the performance of the ML models reported in the studies varied considerably. Some achieved a high performance with an Area Under the Curve (AUC) exceeding 0.90, specifically in studies with smaller cohorts [27,29]. This could be potentially due to overfitting the specific population tested. For instance, Hong et al. [29] reported an accuracy of over 0.90 for ML prediction when trained on only 63 patients with COVID-19 to predict the severity of illness. Similarly, Liu et al. [27] achieved very high accuracy (AUC > 0.96) on a small patient cohort of approximately 158 individuals. Nevertheless, studies that involve larger cohorts of over 5000 patients [15,21,25] exhibit slightly lower accuracy (AUC < 0.90), but offer generalizability over a larger population.
Furthermore, numerous studies on lab biomarkers have yielded notably robust accuracy. The same performance can also be achieved when the model has been trained on demographics and eHR data, suggesting that a viable alternative during a pandemic crisis when hospital resources are limited and obtaining lab test results become challenging. Multiple ML approaches have been employed, including boosting-based classifiers [18,24,25,26,27], regression methods [15,17,25,30], Support Vector Machine [18,24,27], Random Forest (RF) classifier [18,24,27,30], and Naïve Bayes classifier [17,26]. In comparative studies assessing the accuracy of various ML models, it is worth noting that boosting methods, such as Catboost, and XGBoost, have been reported to deliver the highest accuracy, as documented by Noy et al. [18], Liu et al. [27] and Hong et al. [29]. Nevertheless, several other comparative studies [24,30] have indicated that an RF classifier outperforms other ML models. Consequently, our selection of classification methods remains competitive with previous studies.
Additionally, investigations into Deep Learning (DL) neural network-based approaches [21,23] produced similar accuracy results to their ML counterparts and have also been shown to outperform in certain cases [23]. Furthermore, a hybrid approach has been reported combining DL and ML into a single predictive architecture, where the DL extracted features for chest CT scan images utilized by a CatBoost model downstream along with lab biomarkers and eHR data to predict COVID-19 disease severity with high accuracy were employed [22]. This allows for integrating radiological information with other laboratory and clinical reports, further enabling predictive models to make critical decisions [22].
Many studies have tried to discover the main contributing features to risk prediction utilizing statistical measures [17,20], coefficient scores of regression models [15,21,25,27,30], and other interpretability approaches, such as the Gini index of feature importance [15,23], permutation-based approaches, and SHAP (SHapley Additive exPlanations)-based feature analysis [22,26]. The frequently emphasized demographic factors in predicting COVID-19 disease severity are age, hypertension, diabetes, heart disease, gender, smoking status, and Chronic Kidney Disease (CKD) [15,17,20,25].
It has been shown that the studies incorporating eHRs and lab tests tend to assign higher importance to lab biomarkers in their feature ranking as they are direct measures of current health conditions pertaining to disease severity [18,22,24]. However, studies that relied exclusively on eHR data also performed comparably, identifying comorbidities and other demographic variables as top-ranking features [20]. SHAP analysis can provide additional insights by emphasizing individual features with a more pronounced impact on severity [18]. For instance, it has been found that older age, lower platelet counts, and lower lymphocyte levels are closely associated with COVID-19 disease severity [18]. It is worth noting that the exploration of individual class-based SHAP scores for demographic data has been limited and mainly exists in studies focusing on lab biomarkers. Our current study presents SHAP scores for each feature that contributes to predicting MV requirement, ICU, or IMCU admission.

1.2. Contributions of the Current Study

Our current study focuses on a large South Florida cohort, a previously explored dataset for COVID-19 disease severity predictive analysis. Previously, Datta et al. [31] utilized the same dataset to find the most critical features underlying mortality risk, utilizing an RF classifier and SHAP-based interpretability. A comparative analysis utilizing the performance of traditional ML, DL, and the fine-tuned Large Language Model (LLM) in our previous study revealed that the RF model has the highest precision and accuracy compared to other traditional ML (e.g., XGBoost and KNN) or DL (e.g., MLP) models [32].
Here, we extended the work to understand the essential features indicative of COVID-19 disease severity for in-patients from South Florida for predicting therapeutic interventions, such as ICU with MV requirement, ICU, and IMCU admission. A similar severity prediction was performed using DL on the same dataset [33]; however, the severity scale was not well-defined. Furthermore, the analysis lacked the ability to interpret the model’s decision. The current study explores a better metric for patients, caregivers, and clinicians, which is a direct prediction of therapeutic interventions. Additionally, numerous studies focused solely on comparing ML models and did not examine different interpretability approaches for feature analysis.
Thus, our contributions to this work are two-fold: we determine the severity of COVID-19 disease by utilizing the three distinct RF classifiers, and compare the importance of features across the cases susceptible to the conditions. This comparison shows how different features are similar across different severity cases and vary across the conditions.
To enhance the reliability of the findings, we utilized multiple interpretability methods to better understand the essential features, as their importance varies across many studies. In addition, we analyzed the interactions between SHAP features to understand how combinations of features impact the severity of COVID-19 disease.

2. Materials and Methods

2.1. Dataset Collection and Subject Information

With the exemption of informed consent and the HIPAA waiver, the Institutional Review Board (IRB) approved the study and is exempt from further review. From 14 March 2020 to 16 January 2021, data were obtained from the Memorial Healthcare System (MHS), Hollywood, FL, USA and investigated by the co-authors from the Christine E. Lynn College of Nursing and College of Engineering and Computer Science at Florida Atlantic University.
For this project, retrospective data for 5371 hospitalized patients with confirmed cases of COVID-19, as defined by the RT-PCR sampling of nasal and pharyngeal swabs, were retrieved from a comprehensive healthcare system in South Florida. Initially, 203 (independent patient variables) contained data on ‘patients’ sociodemographic characteristics’ (e.g., age, sex, BMI, and smoking status), ‘pre-hospital comorbidities’ (e.g., diarrhea, diabetes, and pneumonia), and ‘medications’ (e.g., Angiotensin Receptor Blockers (ARBs) and Angiotensin Converting Enzyme (ACE) inhibitors) [31].

2.2. Study Design Considerations

Figure 1 presents the visual abstract of a complete data science cycle. The input dataset for our analysis consisted of 5594 patients admitted to the hospital with COVID-19-related symptoms. Among the 5371 patients with COVID-19, 4296 (80%) were in the training dataset, while the remaining 1075 (20%) were in the test dataset. Preprocessing steps ensured data quality by eliminating repeating variables and removing the features with over 10% missing values [34] from the patients’ eHRs. In the current study, only BMI had 6.5% missing values among the remaining independent variables, for which we utilized the Bayesian Ridge Imputation method [30]. Three separate models were trained utilizing 25 variables, comprising 24 independent variables for each dependent variable: MV, ICU, or IMCU.

2.3. Data Classification

The current research is based on binary classification problems with one of three dependent variables (MV, ICU, or IMCU). Each dependent variable was used in a different model and segmented into binary consequences, classified either as ‘MV requirement’ (‘1’) vs. ‘no MV requirement’ (‘0’) or ‘ICU admission’ (‘1’) vs. ‘no ICU admission’ (‘0’), or ‘IMCU admission’ (‘1’) vs. ‘no IMCU admission’ (‘0’). The labels were then converted to binary numerical values before training the models. A classification model was trained on the observed values (independent variables, see Table 1) based on the inputs being binary (e.g., sex and diarrhea) or multiclass (e.g., age and BMI), and the model predicted the outputs (dependent variable) of the binary class.
Of the 24 independent variables used in this research, age and BMI were converted from continuous to categorical variables. Age was categorized based on age matrices [35], where ‘younger adults’ ranged between 20 and 34, ‘middle adults’ ranged between 35 and 64, and ‘older adults’ ranged between 65 and 90 years. Similarly, BMI was categorized based on the BMI metrics [36], with ‘underweight’ defined as below 18.50, ‘normal weight’ between 18.5 and 24.9, ‘overweight’ between 25 and 29.9, and ‘obese’ defined as 30 or greater (see Table 1).
In addition, an improved alternative for dummy coding demonstrated in Table 1 was adapted from our previous study on the same dataset [31] in place of ‘one-hot-encoding’ [37]. This approach harnessed the capabilities and efficiency of AI modeling while leveraging interpretability and domain knowledge. Thus, the dummy coding facilitated the effective comprehension of the feature analysis. It is important to acknowledge that the model’s optimal performance may be affected due to biases introduced in this approach.
However, it simultaneously equips healthcare providers with a well-defined understanding of the health status based on the features analyzed. For example, this research delved into inquiries, such as which age group or presence of diabetes had a more pronounced impact, whether patients taking medication for hypertension (ARBs and ACEIs) yielded benefits, whether individuals with diarrhea exhibited a higher predictability of requiring ICU admission, or if the combination of two variables (interaction effect) exerted a more substantial influence than a single variable.

2.4. Correlation Check

Tetrachoric correlation analysis was employed to acquire a more in-depth understanding of the effectiveness and suitability of 24 independent variables. We chose the tetrachoric correlation due to the categorical nature of our variables, which cannot be analyzed using the Pearson correlation analysis. In addition, the tetrachoric correlation can also handle skewness and outliers in the dataset, capture non-linear correlations between variables, and deal with ordinal data [38].
The analysis revealed that the correlation coefficients among all pairs of variables were generally low (<0.50), with one exception being the positive correlation observed between CKD stage 5 and the dependence on renal dialysis, presenting a correlation coefficient of 0.66. The correlation analysis was conducted solely for exploratory data analysis and was not used as a guiding principle for feature selection. The model’s accuracy can potentially be enhanced when two correlated features are part of the same dataset, as discussed by Deb et al. [39].

2.5. Data Splitting

The Scikit-learn library [40] randomly split the dataset into two: training (N = 4296) and test (N = 1075). The datasets (train and test) contained 8%, 11%, and 19% of the data from ‘MV’, ‘ICU’, and ‘IMCU’, respectively; similarly, 92%, 89%, and 81% of the data were from the ‘no MV’, ‘no ICU’, and ‘no IMCU’ classes, respectively.

2.6. Resampling Data

The data were found to be imbalanced based on the information above. Oversampling was utilized to balance uneven datasets. Synthesizing new examples instead of duplicating was performed using the Synthetic Minority Over-sampling TechniquE (SMOTE) to balance the data [41]. SMOTE was applied to the training dataset, but the test dataset did not undergo any modifications to prevent data leakage issues [39].

3. Results

3.1. Cohort Description

The table below describes the categorical variables and corresponding dummy coding values.
We present the rest of the results in the following order: MV requirement analysis is reported first, followed by ICU admission and IMCU admission. Section 3 includes statistical analysis, model prediction, feature interpretability, and interactions.

3.2. Statistical Analysis

Three individual chi-squares were used to estimate the predictive value of 24 independent variables in identifying individuals likely to require MV, or be admitted to ICU or IMCU [42].
As shown in Table 2, 12 of the 24 variables were statistically significant in predicting the likelihood of requiring MV. These variables included age, sex, diabetes, hypertension, CKD stages 1–4, CKD stage 5, heart failure, coronary artery disease, liver disease, pneumonia, diarrhea, and dependence on renal dialysis. The highest risk factors associated with MV were diarrhea (OR = 7.2), CKD stages 1–4 (OR = 2.89), and hypertension (OR = 2.75).
As shown in Table 3, among the 24 variables, 15 were statistically significant in predicting the likelihood of admitting to the ICU. These variables included age, BMI, sex, race, diabetes, hypertension, Chronic Obstructive Pulmonary Disease (COPD), CKD stages 1–4, CKD stage 5, heart failure, coronary artery disease, liver disease, pneumonia, diarrhea, and dependence on renal dialysis. The variables associated with the highest risk factors were diarrhea (OR = 8.86) and age (OR = 3.43), with older adults being 3.43-times more likely to be admitted to the ICU than younger adults.
As shown in Table 4, 15 of the 24 variables were statistically significant in predicting the likelihood of admitting to the IMCU. These variables included age, sex, smoking status, diabetes, hypertension, COPD, CKD stages 1–4, CKD stage 5, heart failure, cardiac arrhythmias, coronary artery disease, pneumonia, ARBs, ACEIs, and diarrhea. The highest risk factors were diarrhea (OR = 2.86) and age (OR = 2.74), with older adults being 2.74-times more likely to be admitted to the IMCU than younger adults.
Using traditional statistical approaches, multiple options exist for selecting the best key features. One of the most popular approaches involves both forward and backward stepwise approaches. We chose the backward Wald stepwise binary logistic regression method due to its conservativeness and low likelihood of introducing false positives to the model [43,44,45].
Several problems have been noted with this traditional approach. The biggest issue is that these approaches can overfit the model specifically for that sample, and therefore, selecting the key features can vary from sample to sample [46,47]. To compensate for this potential problem, a 10-fold cross-validation approach was used with Wald’s backward binary logistic regression. This validation routine randomly parsed the data into 9 separate training datasets of 573 patients and a 10th testing dataset with 574 patients. Mean Squared Errors (MSEs) and R2 were compared in each fold. Lastly, a final imputed dataset was constructed from the features identified in the training dataset and compared with the values in the test dataset to obtain the final variable selection and R2. Table 5, Table 6 and Table 7 report the three separated Wald’s backward binary logistic regression results for MV, ICU, and IMCU.
As seen in Table 5, of the initial 24 features used to predict the likelihood to require MV, the following 14 are retained by the model: age, BMI, sex, race, diabetes, hypertension, CKD stages 1–4, cardiac arrhythmias, cerebrovascular disease, pneumonia, ARBs, ACEIs, diarrhea, and dependence on renal dialysis. These features were statistically significant in predicting MV requirement and accounted for approximately 20% of the variability in the model [χ2(8) =373.96, p < 0.001, Nagelkerke R2 = 0.20] with a 92.5% overall correct classification rate. Of these retained variables, diarrhea had the largest odds ratio in the multivariate binary logistic regression model; those diagnosed with diarrhea were 6.31-times more likely to require MV than those who did not.
As seen in Table 6, of the initial 24 features used to predict patients admitted to the ICU, the following 15 are retained by the model: age, BMI, sex, race, diabetes, hypertension, asthma, CKD stages 1–4, heart failure, cardiac arrhythmias, cerebrovascular disease, pneumonia, ARBs, ACEIs, and diarrhea. These features were statistically significant in predicting ICU admissions and accounted for approximately 26% of the variability in the model [χ2(8) = 374.96, p < 0.001, Nagelkerke R2 = 0.26] with an 89.4% overall correct classification rate. Of these retained variables, diarrhea had the highest odds ratio in the multivariate binary logistic regression model; those diagnosed with diarrhea were 8.52-times more likely to be admitted to the ICU than those who did not.
Lastly, Table 7 presents the results for predicting the likelihood of being admitted to the IMCU; of the initial 24 features, the following 14 were retained by the model: age, BMI, sex, race, ethnicity, diabetes, hypertension, COPD, CKD stages 1–4, heart failure, cardiac arrhythmias, pneumonia, ACEIs, and diarrhea. These features were statistically significant in predicting IMCU admissions and accounted for approximately 11% of the variability in the model [χ2(15) = 374.96, p < 0.001, Nagelkerke R2 = 0.11] with an 81.5% overall correct classification rate. Of these retained variables, diarrhea had the largest odds ratio in the multivariate binary logistic regression model; those diagnosed with diarrhea were 2.56-times more likely to be admitted to IMCU than those who did not.
To further assess the potential of invariances across gender, age, and ethnicity, sets of backward binary logistic regressions were conducted and cross-validated for each group. Table 8, Table 9 and Table 10 report all the variables we identified as important for predicting MV requirement and ICU and IMCU admissions across groups.
The classification accuracy for MV by sex ranged from 91.5% for males to 93.6% for females. There were more significant differences in model accuracy between the age categories, with older adults having the lowest accuracy (89.7%) compared to younger adults (97.3%) and ‘middle’ adults (94.1%). The accuracy for the non-Hispanic and Hispanic groups was approximately equal for MV, with accuracy values of 92.4% and 92.7%, respectively.
The highest overall accuracy across all groups was observed for predicting ICU admissions. All groups reported an accuracy value above 90% for predicting ICU admission, except for ethnicity, approaching 90% accuracy (non-Hispanic = 89.6% and Hispanic = 89.8%).
When assessing IMCU admissions, the overall accuracy decreased, ranging from a low of 76.4% for older adults to a high of 90.5% for younger adults. As can be seen, some variabilities across the features are important for predicting MV, ICU, and IMCU outcomes among the demographic groups of sex, age, and ethnicity. This suggests the necessity of retaining demographic variables in all future analyses.

3.3. Predictive Analysis

3.3.1. Model Performance

The models’ performances on the imbalanced dataset were assessed with three metrics: precision, recall, F1-score (weighted), and AUC. Precision measures how many of the positive predictions made by the classifier are correct, i.e., ‘of all the predicted severe cases, how many are actually severe?’ Recall quantifies how well the classifier identifies all positive cases, i.e., ‘what proportion of actual severe cases did the model correctly predict as severe?’
F1 score balances precision and recall, considering both false positives and negatives. The harmonic mean of precision and recall is used to estimate the F1 score, where the best value is 1.0 and the worst value is 0.0 [48]. Upon training the RF classifier, we obtained F1-scores (Figure 2a) of 89% (precision: 88% and recall: 90%), 87% (precision: 87% and recall: 88%), and 75% (precision: 73% and recall: 78%), respectively, for MV, ICU, and IMCU.
Model performance was relatively poor in predicting IMCU admissions because, in this case, the model displayed an inadequate performance in the majority class, which resulted in a poor F1 score. The above results are reported after the models have been optimized for hyperparameter tuning using the grid search K-fold cross-validation (10-fold) method [49,50].
Confusion matrices (Figure 2b) were used to assess the model’s classification, which indicated that the model accurately classified 972 (TP: 960 and TN: 12), 942 (TP: 896 and TN: 46), and 842 (TP: 810 and TN: 32) instances, and misclassified 103 (FP: 34 and FN: 69), 133 (FP: 60 and FN: 73), and 233 (FP: 63 and FN: 170) instances for MV, ICU, and IMCU, respectively. The model predicted that the IMCU had the lowest TP and highest FP, leading to low specificity and sensitivity.
The models’ performances were also evaluated using the AUC of the Receiver Operating Characteristic (ROC) analysis. A perfect model would provide an AUC of 1, and an uninformed model would provide an AUC of 0.5. A probability curve plots the TP versus the FP rate at different threshold values and distinguishes between ‘severity’ and ‘no-severity’ cases.
The reported ROC-AUC (Figure 3a) values for the predictions of MV, ICU, and IMCU were 73% (95% CI: 0.66–0.79), 80% (95% CI: 0.75–0.85), and 64% (95% CI: 0.60–0.69), respectively. In the current study, AUC performance for MV and IMCU was poorer than ICU.
In the test data, we observed an inadequate model performance in predicting minority classes. The imbalanced nature of the dataset caused this discrepancy. The model is biased toward the data’s predominant nature to predict the majority class; for instance, predicting ‘ICU’ is lower in accuracy than ‘no ICU’.

3.3.2. Model Interpretability

Three interpretability (post hoc) methods were used to explain the models’ performance. These included MDI (Mean Decrease in Impurity), Permutation Importance, and SHAP for being the most common interpretable techniques based on our background literature reviews [15,22,23,26,31].
The first interpretability method was MDI, based on the Gini impurity values of the RF, where splitting rules maximize impurity reduction [51,52,53,54]. Impurity indicates how well a node can classify a dataset. The Gini index determines the impurity by calculating the probability of misclassifying a randomly selected element from the dataset [55]. The term ‘impurity’ determines how homogeneous or mixed the classes are in a node, where ‘zero’ impurity has only one class, and a node with maximum impurity has an equal mix of all classes [52]. MDI measures how much each feature reduces the impurity of the nodes where it is used to split the data, and it calculates feature importance as the sum of impurity decreases across all the splits that include the feature, averaging over all the trees in the ensemble [52]. The higher the MDI (Figure 4a, Figure 5a and Figure 6a), the more influential the feature is for the model. A split with a significant decrease in impurity is important for the RF classifier; consequently, the feature and the corresponding split are essential for the model’s decision.
The second method was model-agnostic interpretability based on permutation-based analysis [56,57]. In this method, the values of an input feature on the dataset are shuffled, and the model’s change in prediction accuracy is recorded for the shuffled feature. The exact process is repeated for other features, keeping the rest of the dataset unshuffled. The features are then ranked based on the variability in model performance [56,57]. This permutation-based feature ranking (Figure 4b, Figure 5b and Figure 6b) provides the most important values for the model’s performance.
The third (another model agnostic) interpretability method is based on the cooperative game theory that computes SHAP values for each player in a multiplayer game to understand the outcome [58,59,60,61]. Shapley values are measured for each feature to identify its contribution to changing the model’s decision. The SHAP value for a particular feature is calculated as the change in the difference in prediction from randomly sampling the feature value from a distribution compared to the average prediction across all data instances (see Figure 4c,d, Figure 5c,d and Figure 6c,d).
SHAP and MDI have a similar performance when identifying important features [51]; however, SHAP considers the local and global distributions of features, whereas MDI only provides global importance. SHAP overestimates irrelevant features over relevant features under binary networks [62] and incorrectly identifies important features in out-of-distribution scenarios [63]. However, the permutation-based approach is unique as it can capture other important features undetected by SHAP [64].
The ‘bee swarm plot’ (Figure 4c, Figure 5c and Figure 6c) depicts the SHAP values of each feature distributed for all populations, where the sub-categories of each feature are presented as a ‘heatmap’ as per the dummy coding (see Table 1). ‘Cool color’ (blue) represents a sub-category coded by a lower number (i.e., ‘0’ denotes ‘no-diabetes’), and ‘warm color’ (red) represents a higher number (i.e., ‘1’ denotes ‘diabetes’). On the X-axis, the ‘0’ in the middle of the SHAP scale indicates a ‘Neutral Point’—a value toward the right supports a positive decision (predicting ‘MV’, ‘ICU’, or ‘IMCU’); consequently, a value toward the left supports a negative decision (predicting ‘no-MV’, ‘no-ICU’, or ‘no-IMCU’). Features are arranged along the Y-axis based on their importance in ranking order; the more important features are placed at the top, which is assigned using their absolute Shapley values.
SHAP analysis can also be comprehended through the ‘waterfall plot’ (Figure 4d, Figure 5d and Figure 6d), where each feature’s contribution can be analyzed to understand model explainability. The features on the Y-axis are ranked based on the compositional score (X-axis, top) estimated as the mean of SHAP values across the population for the presented feature. The cumulative score (X-axis, bottom) shows the additive values of features on the model that explain the model’s interpretation.
In Figure 4a, the top five features that determine a model’s prediction for MV include diarrhea, age, BMI, race, and diabetes, in ranking order, using MDI. The Permutation Importance (Figure 4b) indicates similar top five features; however, sex (ranked 5th) is more important than race (ranked 7th). The SHAP analysis (Figure 4c,d) agrees with the other two methods over the top four features: diarrhea, age, diabetes, and BMI. It ranked hypertension 5th among the most important features. Despite the slight difference, all interpretability methods tend to agree on the most important features contributing to the model’s decision making.
Figure 4d shows that the top 13 features contribute 90% to the model’s interpretation. Among them, five are from sociodemographic characteristics (age, race, sex, smoking status, and ethnicity), six are from comorbidities (diarrhea, diabetes, BMI, hypertension, pneumonia, and CKD stages 1–4), and two are from the medication category (ACEIs and ARBs). Thus, the findings suggest that the trained models can use the information from each feature category to make a successful prediction.
The features were analyzed similarly for ICU, where MDI (Figure 5a), Permutation Importance (Figure 5b), and SHAP (Figure 5c,d) show that there is consistency across the feature importance for both models: ICU and MV. SHAP showed that sex was among the top 5 features instead of race, as indicated by MDI and Permutation Importance. Figure 5d indicates that the top 12 features contributed to the model’s 90% interpretation. Among them, four were from sociodemographic characteristics (age, sex, race, and smoking status), six were from comorbidities (diarrhea, diabetes, BMI, hypertension, pneumonia, and CKD stages 1–4), and two were from the medication category (ARBs, and ACEIs). We found one less socio-demographic feature (ethnicity) than in MV, contributing to the 90% decision.
In the IMCU, MDI (Figure 6a) and Permutation Importance (Figure 6b) demonstrated that BMI, age, diarrhea, and race are the top 4 most important features. However, SHAP (Figure 6c,d) indicated that the top 3 features are identical to the other two methods. In SHAP, BMI is the most important feature, followed by age, diarrhea, hypertension, and sex.
Figure 6d shows that the top 13 features contribute to the model’s 90% interpretation. Among them, four were from sociodemographic characteristics (age, sex, race, and smoking status), eight were from comorbidities (BMI, diarrhea, hypertension, diabetes, pneumonia, CKD stages 1–4, COPD, and heart failure), and one was from the medication category (ACEIs). Compared to MV, we found that one sociodemographic (ethnicity) and one medication (ARBs) feature were replaced by COPD and heart failure. Similarly, COPD and heart failure were more important features than ARBs compared to ICU, contributing to the 90% interpretation.
It is worth noting that the results show similar features between MV and ICU, compared to IMCU, as patients requiring MV were also admitted to the ICU.

3.3.3. SHAP Dependence Plot

A ‘SHAP dependence plot’ demonstrates how the model output differs depending on the interaction of two features [16,65]. By examining the scatterplot, we can analyze the pattern and trend of the relationship between a variable and the model’s output while considering the other variable [66,67]. If there is an interaction effect, it will be evident through distinct patterns along the Y-axis [16,66].
Figure 7 below shows that the sub-groups of one of the categorical variables are on the X-axis, whereas the sub-groups of other categorical variables are on the Y-axis (right). The Y-axis (left) represents the SHAP value that interacts between the two categorical variables associated with each patient, represented as dots.

4. Discussion

The success of Random Forest classifiers in previous studies [23,24,27,30,31,32] led us to select the same approach for our current study, and the performance of our models after training is comparable with previous studies. This suggests that ML can be a valuable tool for predicting the severity of COVID-19 disease for rapid therapeutic interventions, like MV requirement, ICU, or IMCU admission. Under high-demand conditions, healthcare systems can rely on an alternate, fast intelligence system for therapeutic decision making under critical care. Our study contributes toward understanding COVID-19 disease severity in a large sample from South Florida and identifies essential features utilizing interpretability approaches in AI/ML techniques. This study identified the primary aspects of eHR data that play a crucial role in COVID-19 disease severity.
We extended the feature analysis by comparing multiple interpretability techniques. Our analysis explored how each feature impacts the critical condition and how the interplay between features contributes to the severity of COVID-19 disease.
In this study, we trained three classifiers (RF) to predict each of the three different conditions: (1) MV requirement, (2) ICU, and (3) IMCU admission. We trained the models using 24 independent variables containing information on patients’ sociodemographic characteristics, comorbidities, and medications. The F1-scores (weighted) across our three prediction models were 0.89 for MV, 0.87 for ICU, and 0.75 for IMCU, with an AUC of 0.73 for MV, 0.80 for ICU, and 0.64 for IMCU. This performance was comparable to previous findings [15,19,21,25]. It is worth noting that the model’s accuracy was higher for MV and ICU compared to IMCU.
There is evidence of a type I error in predicting IMCU, possibly due to misclassifying patients’ most important features (rank-wise) as false positives, which may not be true in actual cases and are better predictors for MV or ICU instead. In addition, IMCU patients may comprise two different populations: one directly admitted to the IMCU during their initial hospitalization, and another group stepped down from initial hospitalization to the ICU. Thus, the patients in the IMCU have a large diversity of risk factors, and it might be hard for the classifier to learn a better decision boundary. Also, there is evidence of type II errors in predicting IMCU admissions from ‘other COVID-19 patients’ who have not experienced severe events (MV, ICU, and IMCU). In those cases, presumably, the feature distribution of the ‘other COVID-19 patients’ class (‘no-MV’, ‘no-ICU’, and ‘no-IMCU’) may not be distinct from IMCU as the illness of the patient admitted to the IMCU is less severe compared to those in the MV or ICU groups. Thus, the model failed to distinguish features between the two classes (IMCU vs. ‘no-MV’, ‘no-ICU’, and ‘no-IMCU’, when taken together), resulting in more false negatives.
We utilized SHAP analysis to better understand the model’s decision and compared it with other interpretability methods: MDI and Permutation Importance. The SHAP analysis showed that the top five features (diarrhea, age, diabetes, BMI, and hypertension) for predicting MV were similar across other interpretive methods but in different ranking orders (MDI: diarrhea, age, BMI, race, diabetes; Permutation Importance: diarrhea, diabetes, age, BMI, and sex).
In the ICU, the SHAP analysis showed the top five features (diarrhea, age, diabetes, BMI, and sex), which were almost consistent with the other two interpretability methods (MDI: diarrhea, age, BMI, race, and diabetes; Permutation Importance: diarrhea, diabetes, age, BMI, and race). This shows that the other two methods placed race in the top five features instead of sex.
Similarly, in the IMCU, the SHAP analysis showed the top five features (BMI, age, diarrhea, hypertension, and sex) that are only matched for the top three features with the other two interpretability methods (MDI: BMI, age, diarrhea, race, and smoking status; Permutation Importance: BMI, diarrhea, age, race, and pneumonia). As discussed, IMCU is less severe than the other two disease severities (MV and ICU); hence, the model interprets the features differently than the other two disease severities. These findings were consistent with previous research on the clinical attributes and the prevalence of coexisting medical conditions in COVID-19 patients [15,17,20,25]. Even though limited studies compare feature importance across severe conditions, it is important to highlight that only one such study compared feature importance across different age groups, and the risk factors associated with COVID-19 vary across ages [68].
This research investigated how patients with compromised health conditions contribute to worsening the severity of COVID-19 disease. These include diarrhea, diabetes, hypertension, pneumonia, and CKD stages 1–4. This research revealed that patients on medications such as ARBs and ACEs, commonly used to manage high blood pressure and heart failure, decreased the likelihood of COVID-19 disease severity. Other researchers have also reported the protective nature of the medications (ARBs and ACEs) for hypertension [69,70,71,72] and agree with our findings.
Finally, we also found some crucial interactions across features. Patients with pneumonia and diabetes, as well as diabetes and diarrhea, are more likely to require MV. In the ICU, diarrhea affects middle-aged adults the most and interacts strongly with diabetes. In the IMCU, the interaction between diarrhea and pneumonia, as well as hypertension and middle-aged or older adults, increase the risk of severity.
We found some common features for those reported as significant in the chi-squared test and features retained by binary logistic regression models, reported in the statistical results, with the top features reported in the three interpretability approaches. The common features reported in MV include age, sex, diarrhea, diabetes, hypertension, CKD stages 1–4, and pneumonia. Race and BMI stand out as top features across interpretability approaches, but they were not significant in the chi-squared test.
In the ICU, age, sex, race, BMI, diarrhea, diabetes, hypertension, and pneumonia were important features in all statistical (chi-square and binary logistic regression) and interpretability methods. Lastly, in the IMCU, age, sex, diarrhea, diabetes, hypertension, CKD stages 1–4, and pneumonia were important features. BMI, race, and smoking status were important across all three interpretability approaches, but not in the statistical methods.
It is important to note that heart-related comorbidities significant in the statistical methods for all three severity conditions were not among the top 10 important features across the three interpretability techniques.
Not many studies have analyzed features using traditional statistics and ML interpretable techniques, as explored in the current study. Wu et al. [73] used seven different interpretability techniques to identify important features from lab biomarker data for predicting COVID-19 disease severity. However, the several biomarkers analyzed in this study are not commonly administered in regular checkups nor frequently tested by clinicians. Thus, our study uniquely contributes to identifying important features over easily accessible eHR data across different COVID-19 severities of disease.

5. Limitations

We had no opportunity to intervene in the clinicians’ data collection process in the current study. Re-designing the prospective study may help eliminate ‘patient selection’, ‘patient interaction’, and ‘clinician reporting’ biases [74]. Data were collected when there was a surge of in-patients along with a high rate of death and severity of illness due to limited access to COVID-19 treatments (medications and vaccinations), leading to instances of incomplete, missing, or inaccurate clinical reporting.
Some of these variables had to be removed from the dataset due to unavailability for many patients, such as ‘sputum at admission’ and ‘fever’. This information could have improved our models’ predictions if reported properly. Furthermore, the sociodemographic distribution of COVID-19 patients in South Florida could influence the model’s performance [74,75,76], leading to a better performance in certain sub-groups (e.g., sex, age, and ethnicity) with a greater population density.
Similarly, the imbalanced test dataset posed a challenge for the model in predicting the minority class (severe cases) accurately. It is important to reduce type II errors where the model fails to detect the minority class, as those are crucial for patient recovery from the severity of disease. Thus, clinicians should not solely rely on the model’s predictions, but instead use current clinical trends and judgments.
Models can still perform poorly under outliers and unseen cases due to biased estimations and learning about spuriously correlated features, despite the great strides made to improve the models’ performance [74,75,76,77].
Thus, to aid clinicians in stratifying the most important features, we kept all independent variables and avoided feature selection based on statistical analysis, which could lead to underperformance.

6. Future Direction

In our future work, we will explore the case-specific underperformance of our models, as we found that IMCU prediction results in lower accuracy compared to the other two models. Data augmentation and feature engineering [10,75] based on the significant features reported in the statistics can also help optimize the model’s performance. Specifically, feature selection methods, such as Uniform Manifold Approximation and Projection (UMAP), have been shown to significantly improve a classifier’s performance when predicting COVID-19 severity [78]. Hence, we will implement UMAP feature selection in our future works for a more robust performance.
Second, we will explore a multiclass classifier that simultaneously predicts all three severity cases rather than training separate models for each condition. A comprehensive understanding of disease progression can enhance the transparent and trustworthy knowledge of end users and aid them in making better-informed decisions with the model.
In addition to multiclass classifications, we will utilize SHAP scores to cluster different patient groups, as Khadem et al. [79] reported, identifying the most susceptible cluster based on mortality counts. Cluster analysis may be necessary in future studies to identify the clusters representing a higher severity risk of COVID-19.

7. Conclusions

Developing an AI-driven decision support system for predicting the critical clinical events of in-patients with COVID-19 disease not only addresses the immediate needs of the pandemic, but also advances the field of AI/ML in healthcare. By leveraging cutting-edge technologies and algorithms, such as ML, researchers and healthcare professionals can unlock the potential of data-driven insights to transform patient care. The application of AI/ML in healthcare extends beyond the COVID-19 disease, holding promise for improving diagnosis, treatment selection, disease surveillance, and patient outcomes across various medical specialties and healthcare settings. By pushing the boundaries of AI/ML in healthcare, we can foster innovation, enhance decision making, and ultimately improve health outcomes for individuals and populations. This knowledge empowers public health authorities to proactively plan and implement targeted interventions, mitigating the impact of disease outbreaks and optimizing healthcare delivery.

Author Contributions

Conceptualization, D.D.; methodology, D.D.; software, D.D.; validation, D.D., S.R. and D.N.; formal analysis, D.D.; investigation, P.E.; resources, S.G.D., P.E. and C.S.; data curation, P.E. and C.S.; writing—original draft preparation, S.R., D.D., D.N. and L.M.; writing—review and editing, D.D., S.R., J.H., C.S., P.E., L.M., D.N. and S.G; visualization, D.D.; supervision, S.G.D. and D.N.; project administration, S.G.D.; funding acquisition, L.M., D.D., D.N. and S.G.D. All authors have read and agreed to the published version of the manuscript.

Funding

Internal grant funding from Florida Atlantic University Center for SMART Health.

Institutional Review Board Statement

The study was conducted in accordance with the Declaration of Helsinki and approved by the Institutional Review Board (or Ethics Committee) of FLORIDA ATLANTIC UNIVERSITY (protocol code 1946692-1, 12 September 2022).

Informed Consent Statement

Patient consent was waived due to unidentifiable data.

Data Availability Statement

This dataset is provided by the South Florida Memorial Healthcare System for analysis. We cannot make this publicly available without their permission. You can contact the corresponding authors to request permission to access the data from the contact at South Florida Memorial Healthcare System.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Miller, I.F.; Becker, A.D.; Grenfell, B.T.; Metcalf, C.J.E. Disease and healthcare burden of COVID-19 in the United States. Nat. Med. 2020, 26, 1212–1217. Available online: https://www.nature.com/articles/s41591-020-0952-y (accessed on 20 August 2024). [CrossRef] [PubMed]
  2. Worldometer. COVID-19 Coronavirus Pandemic, 2023. World Health Organization (WHO). COVID-19 Weekly Epidemiological Update 2023. Available online: https://www.worldometers.info/coronavirus/ (accessed on 12 July 2024).
  3. Richardson, S.; Hirsch, J.S.; Narasimhan, M.; Crawford, J.M.; McGinn, T.; Davidson, K.W.; the Northwell COVID-19 Research Consortium. Presenting Characteristics, Comorbidities, and Outcomes among 5700 Patients H Hospitalized with COVID-19 in the New York City Area. JAMA 2020, 323, 2052–2059. [Google Scholar] [CrossRef]
  4. Meille, G.; Decker, S.L.; Owens, P.L.; Selden, T.M. COVID-19 admission rates and changes in US hospital inpatient and intensive care unit occupancy. In JAMA Health Forum; American Medical Association: Chicago, IL, USA, 2023; Volume 4, p. e234206. [Google Scholar]
  5. Ranney, M.L.; Griffeth, V.; Jha, A.K. Critical supply shortages—The need for ventilators and personal protective equipment during the COVID-19 pandemic. N. Engl. J. Med. 2020, 382, e41. [Google Scholar] [CrossRef] [PubMed]
  6. Schwab, P.; DuMont Schütte, A.; Dietz, B.; Bauer, S. Clinical predictive models for COVID-19: Systematic study. J. Med. Internet Res. 2020, 22, e21439. [Google Scholar] [CrossRef] [PubMed]
  7. Bhatraju, P.K.; Ghassemieh, B.J.; Nichols, M.; Kim, R.; Jerome, K.R.; Nalla, A.K.; Greninger, A.L.; Pipavath, S.; Wurfel, M.M.; Evans, L.; et al. COVID-19 in critically ill patients in the Seattle region—Case series. N. Engl. J. Med. 2020, 382, 2012–2022. [Google Scholar] [CrossRef]
  8. Cummings, M.J.; Baldwin, M.R.; Abrams, D.; Jacobson, S.D.; Meyer, B.J.; Balough, E.M.; Aaron, J.G.; Claassen, J.; Rabbani, L.E.; Hastie, J.; et al. Epidemiology, clinical course, and outcomes of critically ill adults with COVID-19 in New York City: A prospective cohort study. Lancet 2020, 395, 1763–1770. [Google Scholar] [CrossRef]
  9. Gupta, S.; Hayek, S.S.; Wang, W.; Chan, L.; Mathews, K.S.; Melamed, M.L.; Brenner, S.K.; Leonberg-Yoo, A.; Schenck, E.J.; Radbel, J.; et al. Factors associated with death in critically ill patients with coronavirus disease 2019 in the US. JAMA Intern. Med. 2020, 180, 1436–1447. [Google Scholar] [CrossRef]
  10. Zhang, A.; Xing, L.; Zou, J.; Wu, J.C. Shifting machine learning for healthcare from development to deployment and from models to data. Nat. Biomed. Eng. 2022, 6, 1330–1345. [Google Scholar] [CrossRef]
  11. Abràmoff, M.D.; Lavin, P.T.; Birch, M.; Shah, N.; Folk, J.C. Pivotal trial of an autonomous AI-based diagnostic system for detection of diabetic retinopathy in primary care offices. NPJ Digit. Med. 2018, 1, 39. [Google Scholar] [CrossRef]
  12. Johnson, K.B.; Wei, W.; Weeraratne, D.; Frisse, M.E.; Misulis, K.; Rhee, K.; Zhao, J.; Snowdon, J.L. Precision medicine, AI, and the future of personalized health care. Clin. Transl. Sci. 2021, 14, 86–93. [Google Scholar] [CrossRef] [PubMed]
  13. Dixon, D.; Sattar, H.; Moros, N.; Kesireddy, S.R.; Ahsan, H.; Lakkimsetti, M.; Fatima, M.; Doshi, D.; Sadhu, K.; Hassan, M.J. Unveiling the Influence of AI Predictive Analytics on Patient Outcomes: A Comprehensive Narrative Review. Cureus 2024, 16, e59954. [Google Scholar] [CrossRef] [PubMed]
  14. Hilton, C.B.; Milinovich, A.; Felix, C.; Vakharia, N.; Crone, T.; Donovan, C.; Proctor, A.; Nazha, A. Personalized predictions of patient outcomes during and after hospitalization using artificial intelligence. NPJ Digit. Med. 2020, 3, 51. [Google Scholar] [CrossRef] [PubMed]
  15. Chen, Z.; Russo, N.W.; Miller, M.M.; Murphy, R.X.; Burmeister, D.B. An observational study to develop a scoring system and model to detect risk of hospital admission due to COVID-19. J. Am. Coll. Emerg. Physicians Open 2021, 2, e12406. [Google Scholar] [CrossRef]
  16. Liang, W.; Liang, H.; Ou, L.; Chen, B.; Chen, A.; Li, C.; Li, Y.; Guan, W.; Sang, L.; Lu, J.; et al. Development and validation of a clinical risk score to predict the occurrence of critical illness in hospitalized patients with COVID-19. JAMA Intern. Med. 2020, 180, 1081–1089. [Google Scholar] [CrossRef]
  17. Zhao, Z.; Chen, A.; Hou, W.; Graham, J.M.; Li, H.; Richman, P.S.; Thode, H.C.; Singer, A.J.; Duong, T.Q. Prediction model and risk scores of ICU admission and mortality in COVID-19. PLoS ONE 2020, 15, e0236618. [Google Scholar] [CrossRef]
  18. Noy, O.; Coster, D.; Metzger, M.; Atar, I.; Shenhar-Tsarfaty, S.; Berliner, S.; Rahav, G.; Rogowski, O.; Shamir, R. A machine learning model for predicting deterioration of COVID-19 inpatients. Sci. Rep. 2022, 12, 2630. [Google Scholar] [CrossRef]
  19. Ferrari, D.; Milic, J.; Tonelli, R.; Ghinelli, F.; Meschiari, M.; Volpi, S.; Faltoni, M.; Franceschi, G.; Iadisernia, V.; Yaacoub, D.; et al. Machine learning in predicting respiratory failure in patients with COVID-19 pneumonia—Challenges, strengths, and opportunities in a global health emergency. PLoS ONE 2020, 15, e0239172. [Google Scholar] [CrossRef] [PubMed]
  20. Ryan, C.; Minc, A.; Caceres, J.; Balsalobre, A.; Dixit, A.; Ng, B.K.; Schmitzberger, F.; Syed-Abdul, S.; Fung, C. Predicting severe outcomes in COVID-19 related illness using only patient demographics, comorbidities and symptoms. Am. J. Emerg. Med. 2021, 45, 378–384. [Google Scholar] [CrossRef]
  21. Singh, V.; Kamaleswaran, R.; Chalfin, D.; Buño-Soto, A.; Roman, J.S.; Rojas-Kenney, E.; Molinaro, R.; von Sengbusch, S.; Hodjat, P.; Comaniciu, D.; et al. A deep learning approach for predicting severity of COVID-19 patients using a parsimonious set of laboratory markers. iScience 2021, 24, 103523. [Google Scholar] [CrossRef]
  22. Chieregato, M.; Frangiamore, F.; Morassi, M.; Baresi, C.; Nici, S.; Bassetti, C.; Bnà, C.; Galelli, M. A hybrid machine learning/deep learning COVID-19 severity predictive model from CT images and clinical data. Sci. Rep. 2022, 12, 4329. [Google Scholar] [CrossRef]
  23. Li, X.; Ge, P.; Zhu, J.; Li, H.; Graham, J.; Singer, A.; Richman, P.S.; Duong, T.Q. Deep learning prediction of likelihood of ICU admission and mortality in COVID-19 patients using clinical variables. PeerJ 2020, 8, e10337. [Google Scholar] [CrossRef] [PubMed]
  24. Magunia, H.; Lederer, S.; Verbuecheln, R.; Gilot, B.J.; Koeppen, M.; Haeberle, H.A.; Mirakaj, V.; Hofmann, P.; Marx, G.; Bickenbach, J.; et al. Machine learning identifies ICU outcome predictors in a multicenter COVID-19 cohort. Crit. Care 2021, 25, 295. [Google Scholar] [CrossRef]
  25. Beiser, D.G.; Jarou, Z.J.; Kassir, A.A.; Puskarich, M.A.; Vrablik, M.C.; Rosenman, E.D.; McDonald, S.A.; Meltzer, A.C.; Courtney, D.M.; Kabrhel, C.; et al. Predicting 30-day return hospital admissions in patients with COVID-19 discharged from the emergency department: A national retrospective cohort study. J. Am. Coll. Emerg. Physicians Open 2021, 2, e12595. [Google Scholar] [CrossRef]
  26. Garcia-Gutiérrez, S.; Esteban-Aizpiri, C.; Lafuente, I.; Barrio, I.; Quiros, R.; Quintana, J.M.; Uranga, A. Machine learning-based model for prediction of clinical deterioration in hospitalized patients by COVID 19. Sci. Rep. 2022, 12, 7097. [Google Scholar]
  27. Liu, Q.; Pang, B.; Li, H.; Zhang, B.; Liu, Y.; Lai, L.; Le, W.; Li, J.; Xia, T.; Zhang, X.; et al. Machine learning models for predicting critical illness risk in hospitalized patients with COVID-19 pneumonia. J. Thorac. Dis. 2021, 13, 1215–1229. [Google Scholar] [CrossRef]
  28. Purkayastha, S.; Xiao, Y.; Jiao, Z.; Thepumnoeysuk, R.; Halsey, K.; Wu, J.; Tran, T.M.L.; Hsieh, B.; Choi, J.W.; Wang, D.; et al. Machine learning-based prediction of COVID-19 severity and progression to critical illness using CT imaging and clinical data. Korean J. Radiol. 2021, 22, 1213. [Google Scholar] [CrossRef]
  29. Hong, W.; Zhou, X.; Jin, S.; Lu, Y.; Pan, J.; Lin, Q.; Yang, S.; Xu, T.; Basharat, Z.; Zippi, M.; et al. A comparison of XGBoost, random forest, and nomograph for the prediction of disease severity in patients with COVID-19 pneumonia: Implications of cytokine and immune cell profile. Front. Cell. Infect. Microbiol. 2022, 12, 819267. [Google Scholar] [CrossRef] [PubMed]
  30. Patel, D.; Kher, V.; Desai, B.; Lei, X.; Cen, S.; Nanda, N.; Gholamrezanezhad, A.; Duddalwar, V.; Varghese, B.; AOberai, A. Machine learning based predictors for COVID-19 disease severity. Sci. Rep. 2021, 11, 4673. [Google Scholar] [CrossRef]
  31. Datta, D.; Dalmida, S.G.; Martinez, L.; Newman, D.; Hashemi, J.; Khoshgoftaar, T.M.; Shorten, C.; Sareli, C.; Eckardt, P. Using machine learning to identify patient characteristics to predict mortality of in-patients with COVID-19 in south Florida. Front. Digit. Health 2023, 5, 1193467. [Google Scholar] [CrossRef]
  32. Shorten, C.; Cardenas, E.; Khoshgoftaar, T.M.; Hashemi, J.; Dalmida, S.G.; Newman, D.; Datta, D.; Martinez, L.; Sareli, C.; Eckard, P. Exploring Language-Interfaced Fine-Tuning for COVID-19 Patient Survival Classification. In Proceedings of the 2022 IEEE 34th International Conference on Tools with Artificial Intelligence (ICTAI), Macao, China, 31 October–2 November 2022; pp. 1449–1454. [Google Scholar]
  33. Shorten, C.; Khoshgoftaar, T.M.; Hashemi, J.; Dalmida, S.G.; Newman, D.; Datta, D.; Martinez, L.; Sareli, C.; Eckard, P. Predicting the Severity of COVID-19 Respiratory Illness with Deep Learning. In Proceedings of the International FLAIRS Conference Proceedings, Jensen Beach, FL, USA, 15–18 May 2022; Volume 35. [Google Scholar]
  34. Bennett, D.A. How can I deal with missing data in my study? Aust. N. Z. J. Public Health 2001, 25, 464–469. [Google Scholar] [CrossRef]
  35. Statsenko, Y.; Al Zahmi, F.; Habuza, T.; Almansoori, T.M.; Smetanina, D.; Simiyu, G.L.; Gorkom, K.N.-V.; Ljubisavljevic, M.; Awawdeh, R.; Elshekhali, H.; et al. Impact of Age and Sex on COVID-19 Severity Assessed From Radiologic and Clinical Findings. Front. Cell. Infect. Microbiol. 2022, 11, 777070. [Google Scholar] [CrossRef] [PubMed]
  36. Romaine, D.S.; Randall, O.S. The Encyclopedia of the Heart and Heart Disease; Facts on File: New York, NY, USA, 2005. [Google Scholar]
  37. Hancock, J.T.; Khoshgoftaar, T.M. Survey on categorical data for neural networks. J. Big Data 2020, 7, 28. [Google Scholar] [CrossRef]
  38. Kubinger, K.D. On artificial results due to using factor analysis for dichotomous variables. Psychol. Sci. 2003, 45, 106–110. [Google Scholar]
  39. Deb, D.; Smith, R.M. Application of Random Forest and SHAP Tree Explainer in Exploring Spatial (In) Justice to Aid Urban Planning. ISPRS Int. J. Geo Inf. 2021, 10, 629. [Google Scholar] [CrossRef]
  40. Pedregosa, F.; Varoquaux, G.; Gramfort, A.; Michel, V.; Thirion, B.; Grisel, O.; Blondel, M.; Prettenhofer, P.; Weiss, R.; Dubourg, V.; et al. Scikit-learn: Machine learning in Python. J. Mach. Learn. Res. 2011, 12, 2825–2830. [Google Scholar]
  41. Chawla, N.V.; Bowyer, K.W.; Hall, L.O.; Kegelmeyer, W.P. SMOTE: Synthetic minority over-sampling technique. J. Artif. Intell. Res. 2002, 16, 321–357. [Google Scholar] [CrossRef]
  42. Abd El-Raheem, G.O.H.; Mohamed, D.S.I.; Yousif, M.A.A.; Elamin, H.E.S. Characteristics and severity of COVID-19 among Sudanese patients during the waves of the pandemic. Sci. Afr. 2021, 14, e01033. [Google Scholar] [CrossRef]
  43. Hosmer, D.W.; Lemeshow, S.; Sturdivant, R.X. Applied Logistic Regression, 3rd ed.; John Wiley & Sons: Hoboken, NJ, USA, 2013. [Google Scholar]
  44. Steyerberg, E.W.; Eijkemans, M.J.; Habbema, J.F. Stepwise selection in small data sets: A simulation study of bias in logistic regression analysis. J. Clin. Epidemiol. 2018, 101, 76–83. [Google Scholar] [CrossRef]
  45. Harrell, F.E., Jr. Regression Modeling Strategies: With Applications to Linear Models, Logistic Regression, and Survival Analysis; Springer: Berlin/Heidelberg, Germany, 2015. [Google Scholar]
  46. Ambler, G.; Brady, A.R.; Royston, P. Simplifying a prognostic model: A simulation study based on clinical data. Stat. Med. 2002, 21, 3803–3822. [Google Scholar] [CrossRef]
  47. Wang, H.; Li, R.; Tsai, C.L. Tuning parameter selectors for the smoothly clipped absolute deviation method. Biometrika 2007, 94, 553–568. [Google Scholar] [CrossRef]
  48. Hicks, S.A.; Strümke, I.; Thambawita, V.; Hammou, M.; Riegler, M.A.; Halvorsen, P.; Parasa, S. On evaluation metrics for medical applications of artificial intelligence. Sci. Rep. 2022, 12, 5979. [Google Scholar] [CrossRef] [PubMed]
  49. Hamida, S.; El Gannour, O.; Cherradi, B.; Ouajji, H.; Raihani, A. Optimization of Machine Learning Algorithms Hyper-Parameters for Improving the Prediction of Patients Infected with COVID-19. In Proceedings of the 2020 IEEE 2nd international conference on electronics, control, optimization and computer science (ICECOCS), Kenitra, Morocco, 2–3 December 2020; pp. 1–6. [Google Scholar]
  50. Sah, S.; Surendiran, B.; Dhanalakshmi, R.; Yamin, M. COVID-19 cases prediction using SARIMAX Model by tuning hyperparameter through grid search cross-validation approach. Expert Syst. 2023, 40, e13086. [Google Scholar] [CrossRef] [PubMed]
  51. Liu, B.; Udell, M. Impact of accuracy on model interpretations. arXiv 2020, arXiv:2011.09903. [Google Scholar]
  52. Ishwaran, H. The effect of splitting on random forests. Mach. Learn. 2015, 99, 75–118. [Google Scholar] [CrossRef] [PubMed]
  53. Kim, Y.; Kim, Y. Explainable heat-related mortality with random forest and SHapley additive exPlanations (SHAP) models. Sustain. Cities Soc. 2022, 79, 103677. [Google Scholar] [CrossRef]
  54. Zhai, B.; Perez-Pozuelo, I.; Clifton, E.A.; Palotti, J.; Guan, Y. Making sense of sleep: Multimodal sleep stage classification in a large, diverse population using movement and cardiac sensing. ACM Interact. Mob. Wearable Ubiquitous Technol. 2020, 4, 1–33. [Google Scholar] [CrossRef]
  55. Rodríguez-Pérez, R.; Bajorath, J. Interpretation of compound activity predictions from complex machine learning models using local approximations and shapley values. J. Med. Chem. 2019, 63, 8761–8777. [Google Scholar] [CrossRef]
  56. Altmann, A.; Toloşi, L.; Sander, O.; Lengauer, T. Permutation importance: A corrected feature importance measure. Bioinformatics 2010, 26, 1340–1347. [Google Scholar] [CrossRef]
  57. Van Lissa, C.J.; Stroebe, W.; Leander, N.P.; Agostini, M.; Draws, T.; Grygoryshyn, A.; Gützgow, B.; Kreienkamp, J.; Vetter, C.S.; Abakoumkin, G.; et al. Using machine learning to identify important predictors of COVID-19 infection prevention behaviors during the early phase of the pandemic. Patterns 2022, 3. [Google Scholar] [CrossRef]
  58. Moncada-Torres, A.; van Maaren, M.C.; Hendriks, M.P.; Siesling, S.; Geleijnse, G. Explainable machine learning can outperform cox regression predictions and provide insights in breast cancer survival. Sci. Rep. 2021, 11, 6968. [Google Scholar] [CrossRef]
  59. Passarelli-Araujo, H.; Passarelli-Araujo, H.; Urbano, M.R.; Pescim, R.R. Machine learning and comorbidity network analysis for hospitalized patients with COVID-19 in a city in southern Brazil. Smart Health 2022, 26, 100323. [Google Scholar] [CrossRef] [PubMed]
  60. Batunacun; Wieland, R.; Lakes, T.; Nendel, C. Using SHAP to interpret XGBoost predictions of grassland degradation in Xilingol, China. Geosci. Mod. Dev. Discuss. 2020, 2020, 1–28. [Google Scholar] [CrossRef]
  61. Gómez-Ramírez, J.; Ávila-Villanueva, M.; Fernández-Blázquez, M.Á. Selecting the most important self-assessed features for predicting conversion to mild cognitive impairment with random forest and permutation-based methods. Sci. Rep. 2020, 10, 20630. [Google Scholar] [CrossRef] [PubMed]
  62. Huang, X.; Marques-Silva, J. On the failings of Shapley values for explainability. Int. J. Approx. Reason. 2024, 171, 109112. [Google Scholar] [CrossRef]
  63. Kumar, I.E.; Venkatasubramanian, S.; Scheidegger, C.; Friedler, S. Problems with Shapley-Value-Based Explanations as Feature Importance Measures. In Proceedings of the International Conference on Machine Learning, Virtual Site, 12–18 July 2020; pp. 5491–5500. [Google Scholar]
  64. Molnar, C.; Freiesleben, T.; König, G.; Herbinger, J.; Reisinger, T.; Casalicchio, G.; Wright, M.N.; Bischl, B. Relating the Partial Dependence Plot and Permutation Feature Importance to the Data Generating Process. In Proceedings of the World Conference on Explainable Artificial Intelligence, Lisbon, Portugal, 26–28 July 2023; Springer Nature: Cham, Switzerland, 2023; pp. 456–479. [Google Scholar]
  65. Nohara, Y.; Inoguchi, T.; Nojiri, C.; Nakashima, N. Explanation of Machine Learning Models of Colon Cancer Using SHAP Considering Interaction Effects. arXiv 2022, arXiv:2208.03112. [Google Scholar]
  66. Lundberg, S.M.; Erion, G.; Chen, H.; DeGrave, A.; Prutkin, J.M.; Nair, B.; Katz, R.; Himmelfarb, J.; Bansal, N.; Lee, S.-I. From local explanations to global understanding with explainable AI for trees. Nat. Mach. Intell. 2020, 2, 56–67. [Google Scholar] [CrossRef]
  67. Qiu, W.; Chen, H.; Dincer, A.B.; Lundberg, S.; Kaeberlein, M.; Lee, S.I. Interpretable machine learning prediction of all-cause mortality. Commun. Med. 2022, 2, 125. [Google Scholar] [CrossRef]
  68. Molani, S.; Hernandez, P.V.; Roper, R.T.; Duvvuri, V.R.; Baumgartner, A.M.; Goldman, J.D.; Ertekin-Taner, N.; Funk, C.C.; Price, N.D.; Rappaport, N.; et al. Risk factors for severe COVID-19 differ by age for hospitalized adults. Sci. Rep. 2022, 12, 6568. [Google Scholar] [CrossRef]
  69. Ebinger, J.E.; Achamallah, N.; Ji, H.; Claggett, B.L.; Sun, N.; Botting, P.; Nguyen, T.-T.; Luong, E.; Kim, E.H.; Park, E.; et al. Pre-existing traits associated with COVID-19 illness severity. PLoS ONE 2020, 15, e0236240. [Google Scholar] [CrossRef]
  70. Şenkal, N.; Meral, R.; Medetalibeyoğlu, A.; Konyaoğlu, H.; Kose, M.; Tukek, T. Association between chronic ACE inhibitor exposure and decreased odds of severe disease in patients with COVID-19. Anatol. J. Cardiol. 2020, 24, 21–29. [Google Scholar]
  71. Zhang, X.; Cai, H.; Hu, J.; Lian, J.; Gu, J.; Zhang, S.; Ye, C.; Lu, Y.; Jin, C.; Yu, G.; et al. Epidemiological, clinical characteristics of cases of SARS-CoV-2 infection with abnormal imaging findings. Int. J. Infect. Dis. 2020, 94, 81–87. [Google Scholar] [CrossRef] [PubMed]
  72. Ge, L.; Meng, Y.; Ma, W.; Mu, J. A retrospective prognostic evaluation using unsupervised learning in the treatment of COVID-19 patients with hypertension treated with ACEI/ARB drugs. PeerJ 2024, 12, e17340. [Google Scholar] [CrossRef] [PubMed]
  73. Wu, H.; Ruan, W.; Wang, J.; Zheng, D.; Liu, B.; Geng, Y.; Chai, X.; Chen, J.; Li, K.; Li, S.; et al. Interpretable machine learning for covid-19: An empirical study on severity prediction task. IEEE Trans. Artif. Intell. 2021, 4, 764–777. [Google Scholar] [CrossRef]
  74. Ueda, D.; Kakinuma, T.; Fujita, S.; Kamagata, K.; Fushimi, Y.; Ito, R.; Matsui, Y.; Nozaki, T.; Nakaura, T.; Fujima, N.; et al. Fairness of artificial intelligence in healthcare: Review and recommendations. Jpn. J. Radiol. 2024, 42, 3–15. [Google Scholar] [CrossRef] [PubMed]
  75. Ghassemi, M.; Naumann, T.; Schulam, P.; Beam, A.L.; Chen, I.Y.; Ranganath, R. A review of challenges and opportunities in machine learning for health. AMIA Summits Transl. Sci. Proc. 2020, 2020, 191. [Google Scholar]
  76. Kelly, C.J.; Karthikesalingam, A.; Suleyman, M.; Corrado, G.; King, D. Key challenges for delivering clinical impact with artificial intelligence. BMC Med. 2019, 17, 195. [Google Scholar] [CrossRef] [PubMed]
  77. Cohen, J.P.; Cao, T.; Viviano, J.D.; Huang, C.-W.; Fralick, M.; Ghassemi, M.; Mamdani, M.; Greiner, R.; Bengio, Y. Problems in the deployment of machine-learned models in health care. Can. Med. Assoc. J. 2021, 193, E1391–E1394. [Google Scholar] [CrossRef]
  78. Laatifi, M.; Douzi, S.; Bouklouz, A.; Ezzine, H.; Jaafari, J.; Zaid, Y.; El Ouahidi, B.; Naciri, M. Machine learning approaches in COVID-19 severity risk prediction in Morocco. J. Big Data 2022, 9, 5. [Google Scholar] [CrossRef]
  79. Khadem, H.; Nemat, H.; Elliott, J.; Benaissa, M. Interpretable machine learning for inpatient COVID-19 mortality risk assessments: Diabetes mellitus exclusive interplay. Sensors 2022, 22, 8757. [Google Scholar] [CrossRef]
Figure 1. Work flowchart. The flowchart was modified from our previous study on the same dataset [31]. The data inclusion strategy is shown in the flowchart. Out of 5594 hospitalized cases, 5371 are confirmed cases with COVID-19. The inclusion criteria for 25 variables (24 independent and one of three dependent variables (ICU with MV, ICU, or IMCU)) are included in the study after data preprocessing.
Figure 1. Work flowchart. The flowchart was modified from our previous study on the same dataset [31]. The data inclusion strategy is shown in the flowchart. Out of 5594 hospitalized cases, 5371 are confirmed cases with COVID-19. The inclusion criteria for 25 variables (24 independent and one of three dependent variables (ICU with MV, ICU, or IMCU)) are included in the study after data preprocessing.
Diagnostics 14 01866 g001
Figure 2. (a) Performance matrices show that the models’ performances in the test dataset yield an F1-score of 0.89 for MV, 0.87 for ICU admission, and 0.75 for IMCU. (b) Confusion matrices demonstrate that the models correctly classify 972 instances for MV, 942 instances for ICU, and 842 instances for IMCU.
Figure 2. (a) Performance matrices show that the models’ performances in the test dataset yield an F1-score of 0.89 for MV, 0.87 for ICU admission, and 0.75 for IMCU. (b) Confusion matrices demonstrate that the models correctly classify 972 instances for MV, 942 instances for ICU, and 842 instances for IMCU.
Diagnostics 14 01866 g002
Figure 3. The figures above depict the ROC analysis to evaluate the models’ performances in predicting MV, ICU, and IMCU for training and test cohorts. (a) This shows the ROC curve for the training dataset (b). This shows ROC curve for the test data, the prediction of ICU achieves the best performance. It should be noted that FP is 1 minus the specificity (TN), which means that the closer FP is to 0, the higher the sensitivity (TN). Therefore, the point should be in the top-left corner of the ROC curve to obtain the optimal values for specificity and sensitivity.
Figure 3. The figures above depict the ROC analysis to evaluate the models’ performances in predicting MV, ICU, and IMCU for training and test cohorts. (a) This shows the ROC curve for the training dataset (b). This shows ROC curve for the test data, the prediction of ICU achieves the best performance. It should be noted that FP is 1 minus the specificity (TN), which means that the closer FP is to 0, the higher the sensitivity (TN). Therefore, the point should be in the top-left corner of the ROC curve to obtain the optimal values for specificity and sensitivity.
Diagnostics 14 01866 g003
Figure 4. Feature interpretation. (a) Feature importance using MDI in MV; the most important features are depicted in the higher bars. (b) Features ranked (Y-axis) based on the decrease in accuracy (X-axis) under the effect of permutations. The ‘bee swarm plot’ (c) depicts the graphical representation of feature values based on dummy coding (Table 1), where red indicates higher values and blue indicates lower values. Each dot is representative of a patient and highlights which feature values correspond to SHAP values (X-axis). (c) The most important feature is diarrhea, followed by age, diabetes, BMI, and hypertension, the top 5 most important features in ranking order for predicting MV use. The ‘Waterfall plot’ (d) indicates that diarrhea results in a composition ratio of >25% (X-axis, top). Furthermore, the top 13 features accumulate 90% of the model’s cumulative ratio (X-axis, bottom).
Figure 4. Feature interpretation. (a) Feature importance using MDI in MV; the most important features are depicted in the higher bars. (b) Features ranked (Y-axis) based on the decrease in accuracy (X-axis) under the effect of permutations. The ‘bee swarm plot’ (c) depicts the graphical representation of feature values based on dummy coding (Table 1), where red indicates higher values and blue indicates lower values. Each dot is representative of a patient and highlights which feature values correspond to SHAP values (X-axis). (c) The most important feature is diarrhea, followed by age, diabetes, BMI, and hypertension, the top 5 most important features in ranking order for predicting MV use. The ‘Waterfall plot’ (d) indicates that diarrhea results in a composition ratio of >25% (X-axis, top). Furthermore, the top 13 features accumulate 90% of the model’s cumulative ratio (X-axis, bottom).
Diagnostics 14 01866 g004
Figure 5. Feature interpretation. (a) Feature importance using MDI in ICU; the most important features are depicted in the higher bars. (b) Features ranked (Y-axis) based on the decrease in accuracy (X-axis) under the effect of permutations. The ‘bee swarm plot’ (c) shows that the most important feature is diarrhea, followed by age, diabetes, BMI, and sex, the top 5 most important features in ranking order for predicting ICU. The ‘Waterfall plot’ (d) indicates that diarrhea results in a composition ratio of >25% (X-axis, top). Furthermore, the top 12 features accumulate 90% of the model’s cumulative ratio (X-axis, bottom).
Figure 5. Feature interpretation. (a) Feature importance using MDI in ICU; the most important features are depicted in the higher bars. (b) Features ranked (Y-axis) based on the decrease in accuracy (X-axis) under the effect of permutations. The ‘bee swarm plot’ (c) shows that the most important feature is diarrhea, followed by age, diabetes, BMI, and sex, the top 5 most important features in ranking order for predicting ICU. The ‘Waterfall plot’ (d) indicates that diarrhea results in a composition ratio of >25% (X-axis, top). Furthermore, the top 12 features accumulate 90% of the model’s cumulative ratio (X-axis, bottom).
Diagnostics 14 01866 g005
Figure 6. Feature interpretation. (a) Feature importance using MDI in IMCU; the most important features are depicted in the higher bars. (b) Features ranked (Y-axis) based on the decrease in accuracy (X-axis) under the effect of permutations. The ‘bee swarm plot’ (c) shows that the most important feature is BMI, followed by age, diarrhea, hypertension, and sex, the top 5 most important features in ranking order. The ‘Waterfall plot’ (d) indicates that BMI results in a composition ratio of ~20% (X-axis, top). Furthermore, the top 13 features accumulate 90% of the model’s cumulative ratio (X-axis, bottom).
Figure 6. Feature interpretation. (a) Feature importance using MDI in IMCU; the most important features are depicted in the higher bars. (b) Features ranked (Y-axis) based on the decrease in accuracy (X-axis) under the effect of permutations. The ‘bee swarm plot’ (c) shows that the most important feature is BMI, followed by age, diarrhea, hypertension, and sex, the top 5 most important features in ranking order. The ‘Waterfall plot’ (d) indicates that BMI results in a composition ratio of ~20% (X-axis, top). Furthermore, the top 13 features accumulate 90% of the model’s cumulative ratio (X-axis, bottom).
Diagnostics 14 01866 g006
Figure 7. SHAP dependence plot. (a,b) Interactions between pneumonia–diarrhea and diarrhea–diabetes have a higher impact on predicting MV. An increase in the interaction between diarrhea (red dots) and pneumonia (scattered dots on the right) impacts the model’s output (higher SHAP values). (a) Higher density of red dots in the presence of pneumonia (indicated as categorical values on the x−axis), resulting in an increased likelihood of MV (shown as greater SHAP values) during the presence of two variables. (c,d) Interactions between age–diarrhea and diabetes–diarrhea have a higher impact on predicting ICU admissions. In (c), younger and middle-aged adults, diarrhea shows a sharp rise in prediction (higher SHAP values with red dots); however, diarrhea does not affect the older population much (less density of red dots). The X-axis of (c) represents the three different age categories (‘younger, ‘middle’, and ‘older’ adults). (e,f) Interactions between pneumonia–diarrhea and age–hypertension have higher impacts on predicting IMCU admissions.
Figure 7. SHAP dependence plot. (a,b) Interactions between pneumonia–diarrhea and diarrhea–diabetes have a higher impact on predicting MV. An increase in the interaction between diarrhea (red dots) and pneumonia (scattered dots on the right) impacts the model’s output (higher SHAP values). (a) Higher density of red dots in the presence of pneumonia (indicated as categorical values on the x−axis), resulting in an increased likelihood of MV (shown as greater SHAP values) during the presence of two variables. (c,d) Interactions between age–diarrhea and diabetes–diarrhea have a higher impact on predicting ICU admissions. In (c), younger and middle-aged adults, diarrhea shows a sharp rise in prediction (higher SHAP values with red dots); however, diarrhea does not affect the older population much (less density of red dots). The X-axis of (c) represents the three different age categories (‘younger, ‘middle’, and ‘older’ adults). (e,f) Interactions between pneumonia–diarrhea and age–hypertension have higher impacts on predicting IMCU admissions.
Diagnostics 14 01866 g007
Table 1. Parameters and characteristics.
Table 1. Parameters and characteristics.
Dummy Coding
Patients’ CharacteristicsAge‘Young adults’: 0, ‘Middle adults’: 1, ‘Older adults’: 2
Sex‘Female’: 0, ‘Male’: 1
Race‘Black’: 0, ‘Others’: 1, ‘White’: 2
Ethnicity‘Hispanic’: 0, ‘Not Hispanic’: 1
Smoking Status‘Never’: 0, ‘Former’: 1,‘Current’: 2
Pre-hospital ComorbiditiesCOPD‘False’: 0, ‘True’: 1
Kidney Disease (stg1_4)‘False’: 0, ‘True’: 1
Kidney Disease (stg5)‘False’: 0, ‘True’: 1
Diarrhea‘No Diarrhea’: 0, ‘Diarrhea’: 1
Hypertension‘No Hypertension’: 0, ‘Hypertension‘: 1
Diabetes ‘No Diabetes’: 0, ‘Diabetes‘: 1
Pneumonia‘No Pneumonia‘: 0, ‘Pneumonia‘: 1
Heart Failure‘False’: 0, ‘True’: 1
Cardiac Arrhythmias‘False’: 0, ‘True’: 1
Coronary Artery Disease‘False’: 0, ‘True’: 1
Dependence on Renal Dialysis‘No’: 0, ‘Yes’: 1
Cerebrovascular Disease ‘No’: 0, ‘Yes’: 1
BMI‘Underweight’: 0, ‘Normal Weight’: 1, ‘Overweight’: 2, ‘Obesity’: 3
Liver Disease‘No’: 0, ‘Yes’: 1
Asthma‘No’: 0, ‘Yes’: 1
HIV‘No’: 0, ‘Yes’: 1
Cancer‘No’: 0, ‘Yes’: 1
MedicationsARBs‘No’: 0, ‘Yes’: 1
ACEIs‘No’: 0, ‘Yes’: 1
Dummy coding is adopted from our previous study [31].
Table 2. Significant features across patients’ likelihood to require MV.
Table 2. Significant features across patients’ likelihood to require MV.
FeaturesNo MV
(N = 4964)
MV
(N = 407)
N%N%χ2dfpOR
Age
Young adult60912.3225.451.432<0.0012.74
Middle Adult235347.415036.9
Older Adult200240.323557.7
BMI
Underweight671.382.04.9830.1731.77
Normal82116.55413.3
Overweight170334.313432.9
Obese237347.821151.8
Sex (Male)247849.923758.210.4010.0011.4
Race
Black155331.312330.25.5820.0611.09
Other254051.222956.3
White87117.55513.5
Ethnicity (Not Hispanic)331766.826264.41.0110.3140.9
Smoking Status
Never412183.033081.11.0120.6031.29
Former71614.46516.0
Current1272.6122.9
Diabetes193439.024059.062.501<0.0012.25
Hypertension321864.834083.558.901<0.0012.75
COPD4258.64611.33.5310.0601.36
Asthma1513.0143.40.2010.6551.135
CKD Stages 1 to 474415.011728.752.901<0.0012.89
CKD Stage 51563.1256.110.4010.0012.02
Heart Failure66313.48019.712.521<0.0011.59
Cancer2845.7297.11.3510.2451.26
Cardiac Arrhythmias62612.64811.80.2310.6320.93
Cerebrovascular Disease3186.4194.71.9310.1650.72
Coronary Artery Disease77015.59022.112.191<0.0011.55
Liver Disease1022.1163.96.1610.0131.95
HIV531.130.70.4010.5280.69
Pneumonia199240.121452.624.091<0.0011.65
ARBs128325.811027.00.2710.6011.06
ACEIs171034.415036.90.9610.3271.11
Diarrhea61312.320550.4421.161<0.0017.2
Dependence on Renal Dialysis484697.638293.920.581<0.0012.69
Individual chi-square results of the 24 demographic features predicting MV requirement.
Table 3. Significant features for patients’ likelihood of being admitted to the ICU.
Table 3. Significant features for patients’ likelihood of being admitted to the ICU.
FeaturesNo ICU
(N = 4775)
ICU
(N = 596)
N%N%χ2dfpOR
Age
Young adult60095.1314.970.192<0.0013.43
Middle Adult227590.92289.1----
Older Adult190084.933715.1----
BMI 0.0
Underweight6688912.010.5130.0150.667
Normal80191.5748.5----
Overweight164289.419510.6----
Obese226687.731812.3----
Sex (Male)236887.234712.815.791<0.0011.42
Race 0.0
Black151290.21649.816.762<0.0010.86
Other241687.335312.7----
White84791.5798.5----
Ethnicity (Not Hispanic)157888.121411.91.9510.1630.881
Smoking Status 0.0
Never396989.248210.81.9420.3791.23
Former68587.79612.3----
Current12187.11812.9----
Diabetes181883.635616.4103.161<0.0012.41
Hypertension306686.249213.879.711<0.0012.64
COPD39583.97616.113.291<0.0011.62
Asthma14084.82515.22.8410.0921.45
CKD Stages 1 to 469981.216218.861.921<0.0012.78
CKD Stage 515082.93117.16.9110.0091.69
Heart Failure62383.812016.222.331<0.0011.68
Cancer27186.64213.41.8210.1781.26
Cardiac Arrhythmias59988.97511.10.0010.9781
Cerebrovascular Disease30791.1308.91.7510.1851.23
Coronary Artery Disease73285.112814.914.891<0.0011.51
Liver Disease9782.22117.85.4910.0191.76
HIV5394.635.41.8910.1691.51
Pneumonia190786.429913.622.911<0.0011.51
ARBs123288.416111.60.4110.5241.06
ACEIs16378822312.02.3010.1301.15
Diarrhea51162.530737.5683.481<0.0018.86
Dependence on Renal Dialysis11177.63222.418.951<0.0012.38
Individual chi-square results of the 24 demographic features predicting ICU admission.
Table 4. Significant features for patients’ likelihood of being admitted to IMCU.
Table 4. Significant features for patients’ likelihood of being admitted to IMCU.
FeaturesNo IMCU
(N = 4361)
IMCU
(N = 1010)
N%N%χ2dfpOR
Age
Young adult56689.76510.378.792<0.0012.74
Middle Adult209483.740916.3----
Older Adult170176.053624.0----
BMI 0.0
Underweight6688.0912.05.7930.1221.78
Normal72983.314616.7----
Overweight148680.935119.1----
Obese208080.550419.5----
Sex (Male)222183.643516.420.271<0.0011.37
Race 0.0
Black137281.930418.11.3020.5221.09
Other223280.653719.4----
White75781.716918.3----
Ethnicity (Not Hispanic)145881.433418.60.0510.8251.02
Smoking Status 0.0
Never364181.881018.26.8420.0331.28
Former60877.817322.2----
Current11280.62719.4----
Diabetes164875.852624.269.501<0.0011.79
Hypertension274877.281022.8108.311<0.0012.38
COPD33170.314029.740.321<0.0011.96
Asthma12978.23621.81.0110.3141.21
CKD Stages 1 to 462272.223927.853.841<0.0011.86
CKD Stage 513071.85128.210.7810.0011.73
Heart Failure54072.720327.340.971<0.0011.78
Cancer24277.37122.73.2810.0701.29
Cardiac Arrhythmias52778.214721.84.5610.0331.24
Cerebrovascular Disease27280.76519.30.0610.8151.03
Coronary Artery Disease64274.721825.328.721<0.0011.59
Liver Disease9278.02622.00.8210.3641.23
HIV4885.7814.30.7610.3840.72
Pneumonia170677.350022.736.551<0.0011.5
ARBs107677.231722.819.241<0.0011.4
ACEIs148079.638020.44.9210.0261.17
Diarrhea53164.928735.1167.521<0.0012.86
Dependence on Renal Dialysis10875.53524.53.0910.0791.14
Individual chi-square results of the 24 demographic features predicting IMCU admission.
Table 5. Backward binary logistic regression for predicting MV requirement.
Table 5. Backward binary logistic regression for predicting MV requirement.
95% CI for OR
FeaturesBSEWalddfpORLowerUpper
Age0.430.1115.271<0.0011.541.241.91
BMI0.140.083.3110.0691.150.991.34
Sex0.320.117.9410.0051.381.101.72
Black 9.4820.009
Other0.120.130.8110.3691.120.871.44
White−0.410.194.8110.0280.670.460.96
Diabetes0.490.1216.711<0.0011.631.292.05
Hypertension0.730.1718.981<0.0012.081.502.89
CKD Stages 1 to 40.570.1417.471<0.0011.761.352.30
Cardiac Arrhythmias−0.380.184.3410.0370.690.480.98
Cerebrovascular Disease−0.520.263.9010.0480.600.361.00
Pneumonia0.350.1110.2410.0011.431.151.77
ARBs−0.440.1311.111<0.0010.650.500.84
ACEIs−0.280.124.9710.0260.760.600.97
Diarrhea1.840.11268.201<0.0016.315.067.87
Dependence on Renal Dialysis0.690.267.0110.0081.991.203.30
Constant−4.940.30267.941<0.0010.01
Note. Nagelkerke R square = 0.202; classification = 92.5%.
Table 6. Backward binary logistic regression for predicting ICU admission.
Table 6. Backward binary logistic regression for predicting ICU admission.
95% CI for OR
FeaturesBSEWalddfpORLowerUpper
Age0.460.1022.701<0.0011.581.311.91
BMI0.220.0710.6610.0011.251.091.42
Sex0.410.1016.891<0.0011.511.241.83
Black 23.872<0.001
Other0.320.117.9710.0051.381.101.72
White−0.360.164.7810.0290.700.510.96
Diabetes0.610.1034.671<0.0011.841.502.26
Hypertension0.720.1424.921<0.0012.041.542.71
Asthma0.570.255.0810.0241.761.082.88
CKD Stages 1 to 40.430.1212.431<0.0011.541.211.95
Heart Failure0.410.148.6310.0031.511.151.99
Cardiac Arrhythmias−0.370.165.2710.0220.690.500.95
Cerebrovascular Disease−0.440.223.9310.0470.650.421.00
Pneumonia0.230.105.5710.0181.261.041.52
ARBs−0.510.1219.551<0.0010.600.480.75
ACEIs−0.290.117.2010.0070.750.600.92
Diarrhea2.140.10463.521<0.0018.527.0110.36
Constant−4.980.27346.661<0.0010.01
Note. Nagelkerke R square = 0.26; classification = 89.4%.
Table 7. Backward binary logistic regression for predicting IMCU admission.
Table 7. Backward binary logistic regression for predicting IMCU admission.
95% CI for OR
FeaturesBSEWalddfpORLowerUpper
Age0.270.0714.121<0.0011.301.141.50
BMI0.150.059.1810.0021.161.061.28
Sex0.320.0718.761<0.0011.381.191.59
Black 7.0420.030
Other0.220.114.2210.0401.241.011.53
White−0.080.120.5110.4760.920.741.15
Ethnicity0.170.102.8810.0901.190.971.45
Diabetes0.270.0811.751<0.0011.311.121.53
Hypertension0.590.1032.141<0.0011.801.472.20
COPD0.350.128.6610.0031.421.121.78
CKD Stage 1 to 40.220.105.1710.0231.251.031.51
Heart Failure0.290.117.0710.0081.341.081.66
Cardiac Arrhythmias−0.210.123.1010.0780.810.641.02
Pneumonia0.290.0715.541<0.0011.341.161.54
ACEIs−0.240.088.4010.0040.790.670.93
Diarrhea0.940.09119.111<0.0012.562.173.04
Constant−3.420.22249.201<0.0010.03
Note. Nagelkerke R square = 0.11; classification = 81.5%.
Table 8. Backward binary logistic regressions cross-validated by demographic groups for predicting MV requirement.
Table 8. Backward binary logistic regressions cross-validated by demographic groups for predicting MV requirement.
SexAgeEthnicity
FeaturesFemaleMaleYoungMiddleOlderNon-HispanicHispanic
Male------1.89 **--1.47 *1.38 *
Age
  Young Adult--------------
  Middle0.881.24------2.55 **--
  Older1.722.02------8.43 **--
BMI--------------
  Normal--0.340.03 **--1.381.38--
  Overweight--0.360.01 **--1.971.86--
  Obese--0.580.01 **--3.06 *2.74 *--
Race
  Black-------- ----
  Other1.191.06--0.801.62 *--1.19
  White0.76 *0.57 *--0.51 *0.82--0.70 *
Hispanic----3.36 **--------
Smoke
  Never --------------
  Former--------------
  Current--------------
Diabetes1.88 **1.48 * 2.16 * 1.49 *1.69 **
Hypertension1.91 **2.35 ** 2.52 *1.661.99 *2.21 **
COPD--------------
Asthma--------------
CKD Stages 1–41.61 **1.68 **----1.95 **1.56 *1.82 **
CKD Stage 5----------0.151.69
Heart Failure 1.49
Cancer--------------
Cardiac Arrhythmias 0.39 0.52 **
Cerebrovascular Disease0.35----------0.57
Coronary Artery Disease1.43--7.95 **1.59------
Liver Disease----8.74 **--------
HIV--------------
Pneumonia------1.301.56 **----
ARBs0.630.62--0.48 **0.707 **0.51 **1.47 **
ACEIs0.740.73--0.59 * --0.67 *
Diarrhea5.547.9043.17 **7.72 **5.18 **5.59 **0.75
Dependence on Renal Dialysis2.58 ----3.14 **7.61 **6.88 **
R20.180.240.430.220.1824.60.20
Correct Classification %93.6%91.5%97.3%94.1%89.7%92.4%92.7%
Note. Only variables that were identified as key predictors are reported for each model. Unique variance; * p < 0.05; ** p < 0.01. ‘--' in the table indicates the features not retained by the model for the sub-class.
Table 9. Backward binary logistic regressions cross-validated by demographic groups for predicting ICU admissions.
Table 9. Backward binary logistic regressions cross-validated by demographic groups for predicting ICU admissions.
SexAge Ethnicity
Features FemaleMaleYoungMiddle Older Non-HispanicHispanic
Male------1.87 *1.91 * 1.75 **1.39 **
Age
  Young Adult---------- ----
  Middle0.891.26------ 2.080.92
  Older1.712.11------ 5.6 **1.29
BMI
  Normal--0.29 *0.02 **---- ----
  Overweight--0.310.01 **---- ----
  Obese--0.500.01 **---- ----
Race
  Black---------- ----
  Other--1.07--1.321.32 --1.69 *
  White--0.57 *--0.650.65 --0.69 *
Hispanic----2.89 *1.231.23 * ----
Smoke
  Never ---------- ----
  Former---------- ----
  Current---------- ----
Diabetes1.91 *1.47 *--2.42 *2.41 **2.08 *1.76 **
Hypertension1.84 *2.37 *3.01 *2.32 *2.32 **--2.39 **
COPD--0.65---- 1.76*--
Asthma---------- --1.78 *
CKD Stages 1–41.60 *1.71 *10.42 *1.70 *1.70 * 1.221.56 *
CKD Stage 5----6.26 *---- 0.15 *--
Heart Failure--1.62 *3.09 *2.24 *-- --1.39 *
Cancer---------- ----
Cardiac Arrhythmias--0.39------ 0.53 *--
Cerebrovascular Disease0.39 *----0.47 *0.47 0.52--
Coronary Artery Disease1.40-------- 1.31--
Liver Disease----3.09*---- ----
HIV---------- ----
Pneumonia--1.731.69---- 1.261.24 *
ARBs0.65 *0.62 *0.440.52 *0.51 * --0.66 *
ACEIs0.73 *0.730.230.760.76 0.46 **0.69 *
Diarrhea5.56 *7.93 **41.46 **10.13 **10.13 **7.86 **9.28 **
Dependence on Renal Dialysis2.55 *-------- 6.39 **--
R20.180.240.460.480.27 0.290.26
Correct Classification%93.6%91.5%95.9%95.7%97.1% 89.6%89.8%
Note. Only variables that were identified as key predictors are reported for each model. Unique variance; * p < 0.05; ** p < 0.01. ‘--' in the table indicates the features not retained by the model for the sub-class.
Table 10. Backward binary logistic regressions cross-validated by demographic groups for predicting IMCU admissions.
Table 10. Backward binary logistic regressions cross-validated by demographic groups for predicting IMCU admissions.
SexAgeEthnicity
Features FemaleMaleYoungMiddle OlderNon-HispanicHispanic
Male------1.64 **1.25 *1.40 *1.38 *
Age
  Young Adult--------------
  Middle0.86--------0.851.16
  Older1.44--------1.321.46
BMI--------------
  Normal1.52------1.583.14--
  Overweight2.08------2.07 *5.21--
  Obese2.32 *------2.16 *5.39--
Race
  Black--------------
  Other--------1.24 *--1.12
  White--------0.89--0.91
Hispanic--------------
Smoke
  Never --------------
  Former----0.00--------
  Current----2.49 *--------
Diabetes1.35 **1.36 **4.46 **1.45 **--1.63 **1.25 *
Hypertension1.53 **2.14 **0.81 *1.94 **1.62 **1.89 **1.94 *
COPD1.34 *1.51 **--1.59 *1.36 *1.66 **1.29 *
Asthma----2.79--------
CKD Stages 1–4--1.29 *3.48 *1.56 *--1.66 **--
CKD Stage 5----0.42--------
Heart Failure1.36 *1.3 *0.521.44 *1.36 *--1.41 *
Cancer------1.50----0.78
Cardiac Arrhythmias--0.72 *--0.68------
Cerebrovascular Disease--------------
Coronary Artery Disease--------------
Liver Disease--------------
HIV--------------
Pneumonia1.42 **1.28 *2.25 *1.2 *1.35 **1.401.32 *
ARBs----4.21 *--0.81 *0.79--
ACEIs--0.72 *--0.79 **0.74 *0.760.79 *
Diarrhea2.88 **2.42 *2.93 *2.81 **2.30 **2.82 **2.52 *
Dependence on Renal Dialysis0.51 *----------0.76
R20.110.090.230.110.070.150.09
Correct Classification %83.7%78.9%90.5%83.8%76.4%82.1081.1
Note. Only variables that were identified as key predictors are reported for each model. Unique variance; * p < 0.05; ** p < 0.01. ‘--' in the table indicates the features not retained by the model for the sub-class.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Datta, D.; Ray, S.; Martinez, L.; Newman, D.; Dalmida, S.G.; Hashemi, J.; Sareli, C.; Eckardt, P. Feature Identification Using Interpretability Machine Learning Predicting Risk Factors for Disease Severity of In-Patients with COVID-19 in South Florida. Diagnostics 2024, 14, 1866. https://doi.org/10.3390/diagnostics14171866

AMA Style

Datta D, Ray S, Martinez L, Newman D, Dalmida SG, Hashemi J, Sareli C, Eckardt P. Feature Identification Using Interpretability Machine Learning Predicting Risk Factors for Disease Severity of In-Patients with COVID-19 in South Florida. Diagnostics. 2024; 14(17):1866. https://doi.org/10.3390/diagnostics14171866

Chicago/Turabian Style

Datta, Debarshi, Subhosit Ray, Laurie Martinez, David Newman, Safiya George Dalmida, Javad Hashemi, Candice Sareli, and Paula Eckardt. 2024. "Feature Identification Using Interpretability Machine Learning Predicting Risk Factors for Disease Severity of In-Patients with COVID-19 in South Florida" Diagnostics 14, no. 17: 1866. https://doi.org/10.3390/diagnostics14171866

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop