- 1Scientific Research Department, First People’s Hospital of Zigong City, Zigong, China
- 2Graduate School, Guangxi University of Chinese Medicine, Nanning, China
- 3Department of Cardiology, The Second Affiliated Hospital of Chongqing Medical University, Chongqing, China
- 4Department of Neurosurgery, The First Affiliated Hospital of Chongqing Medical University, Chongqing, China
- 5AMITA Health Saint Joseph Hospital Chicago, Chicago, IL, United States
- 6Department of Cardiology, The First Affiliated Hospital of Chongqing Medical University, Chongqing, China
- 7Department of Cardiology, First People’s Hospital of Zigong City, Zigong, China
- 8Information Department, First People’s Hospital of Zigong City, Zigong, China
- 9Department of Infectious Diseases, The First People’s Hospital of Changde City, Changde, China
- 10General Affairs Section, The People’s Hospital of Tongnan District, Chongqing, China
Background: Heart failure (HF) combined with hypertension is an extremely important cause of in-hospital mortality, especially for the intensive care unit (ICU) patients. However, under intense working pressure, the medical staff are easily overwhelmed by the large number of clinical signals generated in the ICU, which may lead to treatment delay, sub-optimal care, or even wrong clinical decisions. Individual risk stratification is an essential strategy for managing ICU patients with HF combined with hypertension. Artificial intelligence, especially machine learning (ML), can develop superior models to predict the prognosis of these patients. This study aimed to develop a machine learning method to predict the 28-day mortality for ICU patients with HF combined with hypertension.
Methods: We enrolled all critically ill patients with HF combined with hypertension in the Medical Information Mart for IntensiveCare Database-IV (MIMIC-IV, v.1.4) and the eICU Collaborative Research Database (eICU-CRD) from 2008 to 2019. Subsequently, MIMIC-IV was divided into training cohort and testing cohort in an 8:2 ratio, and eICU-CRD was designated as the external validation cohort. The least absolute shrinkage and selection operator (LASSO) Cox regression with internal tenfold cross-validation was used for data dimension reduction and identifying the most valuable predictive features for 28-day mortality. Based on its accuracy and area under the curve (AUC), the best model in the validation cohort was selected. In addition, we utilized the Shapley Additive Explanations (SHAP) method to highlight the importance of model features, analyze the impact of individual features on model output, and visualize an individual’s Shapley values.
Results: A total of 3,458 and 6582 patients with HF combined with hypertension in MIMIC-IV and eICU-CRD were included. The patients, including 1,756 males, had a median (Q1, Q3) age of 75 (65, 84) years. After selection, 22 out of a total of 58 clinical parameters were extracted to develop the machine-learning models. Among four constructed models, the Neural Networks (NN) model performed the best predictive performance with an AUC of 0.764 and 0.674 in the test cohort and external validation cohort, respectively. In addition, a simplified model including seven variables was built based on NN, which also had good predictive performance (AUC: 0.741). Feature importance analysis showed that age, mechanical ventilation (MECHVENT), chloride, bun, anion gap, paraplegia, rdw (RDW), hyperlipidemia, peripheral capillary oxygen saturation (SpO2), respiratory rate, cerebrovascular disease, heart rate, white blood cell (WBC), international normalized ratio (INR), mean corpuscular hemoglobin concentration (MCHC), glucose, AIDS, mean corpuscular volume (MCV), N-terminal pro-brain natriuretic peptide (Npro. BNP), calcium, renal replacement therapy (RRT), and partial thromboplastin time (PTT) were the top 22 features of the NN model with the greatest impact. Finally, after hyperparameter optimization, SHAP plots were employed to make the NN-based model interpretable with an analytical description of how the constructed model visualizes the prediction of death.
Conclusion: We developed a predictive model to predict the 28-day mortality for ICU patients with HF combined with hypertension, which proved superior to the traditional logistic regression analysis. The SHAP method enables machine learning models to be more interpretable, thereby helping clinicians to better understand the reasoning behind the outcome and assess in-hospital outcomes for critically ill patients.
Introduction
Cardiac diseases are among the leading causes of overall mortality and hospitalization globally. Amongst them, heart failure (HF) is of the highest socio-economic relevance, and it is a global epidemic with high morbidity, mortality, and readmission rates, affecting more than 64 million people worldwide (1–4). In the United States, the estimated prevalence of HF is expected to increase by 24% to approximately 8.5 million in 2030 (5, 6). Hypertension, the most frequent comorbidity of HF, promotes the development of the disease and contributes to its progression and poor outcome (7). In the United States, around 10–51% of hospitalized patients with HF have been documented with ICU admission, and ICU-admitted patients have significantly higher adjusted in-hospital mortality compared with those admitted to the general medical floor (8–10). In addition, the in-hospital mortality rate for patients treated in the ICU was 10.6%, compared with 4.0% for all HF patients (11). Because ICU physicians receive large amounts of data from many patients stored in electronic patient-data management systems (PDMS) surpassing the amount limits of the human brain to process information, it is often difficult for physicians to extract the most important information in a short period to make the best decisions for patient care. In addition, the limited ability of humans to process this vast amount of data makes them prone to data overload, change blindness, and task fixation (12), which also increases the risk of clinicians failing to identify and interpret relevant information and act accordingly (13, 14). Low nurse-to-patient ratios in the ICU and insufficient numbers of ICU physicians are associated with higher ICU mortality in patients whose conditions deteriorate and who do not receive timely and appropriate treatment (15–18). Despite recent advances in diagnosis and treatment as well as evidence-based management, the results regarding HF remain unsatisfactory (19).
Risk stratification as a common method for risk classification, deciding the duration of intervention, and assessing the mortality in patients with HF combined with hypertension provide not only a fundamental strategy for clinical decision-making but also practical information for health policy and insurance services (20). For this reason, several in-hospital mortality prediction models have been developed and evaluated for risk stratification and mortality prediction of HF patients in the ICU (21–30). However, an interpretable machine learning (ML) model has not been established to predict 28-day in-hospital mortality for ICU patients with HF combined with hypertension.
Artificial intelligence, such as ML techniques, is excellent at analyzing complex signals in data-rich environments (31). The large amount of data collected in ICU and the public availability of datasets such as MIMIC-III (32) and emergency intensive care unit (eICU) (33) are critical to developing ML in this context. ML is the use of computational algorithms that can learn to identify underlying patterns and classes from large amounts of data. It is an alternative method of using previous or existing data to train computer models to make predictions about the outcome. In this study, we aimed to develop and validate an interpretable prediction model to predict 28-day in-hospital mortality in patients admitted to the ICU with HF combined with hypertension using ML algorithms and leveraging data from the Medical Information Mart for Intensive Care (MIMIC-IV) and eICU database.
Materials and methods
Data source and outcome
This study was a retrospective cohort study based on cohort data extracted from the MIMIC-IV (v.1.4) database, which contains over 70,000 ICU admissions across the United States collected from 2008 to 2019. The database is a large, single-center, publicly available, and de-identified patient database containing comprehensive patient information [e.g., demographics, admission records, International Classification of Diseases-9th and Diseases-10th (ICD-9 and ICD-10) revision diagnoses, laboratory tests, medications, procedures, fluid balance, discharge summaries, vital sign measurements undertaken at the bedside, caregivers notes, radiology reports, and survival data] (34, 35). In addition, the eICU-CRD database (version 2.0), a multicenter database of more than 200,000 ICU admissions in the United States, was used as an independent external validation set. We studied these courses in depth and obtained permission to use the database (record ID:42039823). The requirement for individual patient consent and an ethical approval statement was waived as the program does not affect clinical practice and all patient privacy information in the database was de-identified. The selected primary outcome for this study was the all-cause mortality within 28 days of patients with HF combined with hypertension who were admitted to ICU.
Study patients and definitions
All adult patients in the MIMIC-IV database with a diagnosis of HF who were admitted to the ICU were recruited (only the first admission was included for analysis). The diagnosis was identified by a manual review of ICD-9 and ICD-10 codes. The exclusion criteria for participation in the study were as follows: (1) patients without hypertension, (2) patients with ICU length of stay less than 24 h or more than 28 days, (3) patients with severe liver disease, (4) patients with malignant cancer, and (5) patients with more than 30% missing data. In this study, the patient’s first hospitalization time was taken as the starting point for statistics on whether he died, and the patient’s death or whether the patient did not die within the period recorded in the database was taken as the statistical endpoint. The primary outcome of this study was in-hospital mortality, defined as the survival status at the time of hospital discharge.
Data collection and variable extraction
Following the variable selection method of Deshmukh et al. (35), 58 candidate variables that were associated with the results were selected. The extracted variables included the general demographic variables of patients and other important variables, as follows: gender, age, ethnicity, body weight, comorbidities, vital signs, laboratory findings, medical treatments, and first care unit. The severity of the disease was assessed using Sequential Organ Failure Assessment (SOAF), Simplified Acute Physiology Score II (SAPS-II), Oxford Acute Severity of Illness Score (OASIS), and Logistic Organ Dysfunction System scoring system.
Charlson comorbidity index was used, and the comorbidities included hyperlipidemia, atrial fibrillation, paraplegia, renal disease, aids, dementia, diabetes without, diabetes with, peripheral vascular disease, cerebrovascular disease, dementia, chronic pulmonary disease, rheumatic disease, peptic ulcer disease, myocardial infarction, and congestive heart failure. For vital signs, the mean values in the ICU for the following variables were selected: respiratory rate, heart rate (HR), body temperature, mean blood pressure, diastolic blood pressure, systolic blood pressure, creatine kinase MB isoenzyme, creatinine phosphokinase, Troponin, and N-terminal prohormone of B-type natriuretic peptide (NT-proBNP). For the results of the first laboratory examination after admission to the ICU, the mean value for the following variables were selected: RBC, hematocrit, hemoglobin, platelets, MCV, WBC, MCHC, MCH, RDW, anion gap, bicarbonate, BUN, creatinine, calcium, chloride, sodium, potassium, PT, PTT, and INR. Medical treatments included MECHVENT and RRT. Finally, the cumulative urine output within the first 24 h. To reduce the impact of missing data on classification, a modified KNN-based (K-nearest neighbor) classification algorithm to fill in the missing values was proposed. Considering a large number of features still presented in the cohort, the least absolute shrinkage and selection operator (LASSO) regression in the variable selection was utilized to effectively prevent overfitting.
Missing data handling
Missing data with < 30% in each feature was processed by KNN-based classification algorithm using the “DMwR2” package in R. KNN-based classification algorithm is used as a missing value estimation method, which is a non-parametric ML algorithm based on neighbors. The imputed value is the average of the neighbor’s measurements or multiple neighbors’ measurements. The estimation of missing values is obtained by using the average of the non-missing values of its neighbors, and in addition, the case where all neighbors in a given set are missing can be circumvented by averaging the overall column for that particular feature. Therefore, this technique is advantageous for dealing with data sets that contain a large number of variables with missing values (36).
Machine learning model building
The original raw data included gender, age, vital signs, laboratory tests, and comorbidities. The processed data included 58 characteristics. In our study, five common algorithms, i.e., logistic regression (LR), Neural Networks (NN), Multi-Layer perceptron (MLP), Naive Bayes (NB), and Random Forest (RF), were applied to build models for predicting 28-day in-hospital mortality of ICU patients with HF combined with hypertension. To improve the stability of the prediction models, all continuous variables were rescaled to a distribution with a mean of 0 and a standard deviation of 1 with scale transformation. Threefold cross-validation of the ML models to be tuned (LR, NN, MLP, NB, and RF) was performed to select the best prediction model for each algorithm with different tuning parameters. During the search process, the accuracy or ROC was set as the metric. The test set was not used during model tuning and was used only for model evaluation after the entire model selection and training process.
Model assessment
The final models were evaluated using the confusion matrix metrics such as sensitivity, specificity, positive prediction value (PPV), negative prediction value (NPV), accuracy, and the area under the curve (AUC) of the receiver operating characteristic (ROC). ROC curves were constructed based on the prediction probabilities and the area under the curve (AUC) values of the models in the testing dataset were compared to identify the model with the best predictive performance.
Features important
Feature ranking evaluation is a measure that evaluates the importance of each feature in a feature set through its impact on the final classification result. We analyzed the importance of features using the DALEX package, which explains the predictions of any classifier in an interpretable and faithful manner by learning an interpretable model locally around the prediction. By calculating the relative importance of variables, the impact of features on the prediction model was plotted.
Statistical analysis
Patients were divided into two groups based on whether they died or survived during their 28-day stay in the ICU. Then, categorical variables were presented as a percentage of the total and continuous variables as mean ± SD or median and IQR, according to the normality of the distribution. For categorical and continuous variables, between-group differences were compared by using a two-sided Pearson’s χ2 test or Fisher’s and two-sided one-way ANOVA or Wilcoxon rank sum test, respectively. Logistic regression with the LASSO penalization method was performed for predictor selection, which helped to reduce the dimensionality of the prediction model. LASSO regression shrinked the coefficient estimates toward zero, with the degree of shrinkage dependent on an additional parameter, lambda. To determine the penalty factor (lambda), we constructed a tenfold cross-validated error plot for the LASSO model. After that, the patients were randomly divided into two groups, of which 80% were used as the training cohort and the remaining 20% as the test cohort. Five common ML methods were applied to develop the models in the test cohort. The quantitative performance of the models was assessed by comparing the AUC and accuracy in the test cohort. The optimized model with the best mortality prediction performance in the test cohort (i.e., the Neural Network-NN) was defined as the final model. Then the top 7 most important clinical features were screened out of the 22 most influential features in the NN model, and they were used to build the best NN model again and finally used for external validation. The sensitivity performance analysis of the NN model was compared with the seven most important clinical features in the LR model. In addition, the Shapley additive explanations (SHAP) method was adopted to improve the interpretability of the final model. The SHAP values of features were evaluated by the lime package. We selected four cases for the feature’s SHAP value evaluation.
All statistical analyses were carried out using R software (v. 3.6.3, R Foundation for Statistical Computing), and statistical significance was set at p < 0.05.
Results
Baseline characteristics
As shown in Figure 1, data of 15,354 critically ill patients were downloaded from the MIMIC-IV database. Among them, 3,458 HF patients with hypertension were included in our study. Among the included patients, 459 patients passed away and 2,999 survived within 28 days, respectively. Table 1 summarizes the comparison of baseline characteristics, vital signs, and laboratory parameters within 28 days between non-survivors and survivors. In the non-survivor group, gender, age, ethnicity, body weight, SOFA, SAPS II, OASIS, LODS, Charlson comorbidity index, hyperlipidemia, atrial fibrillation, paraplegia, cerebrovascular disease, congestive heart failure, respiratory rate, HR, creatinine phosphokinase, Troponin, NT-proBNP, platelets, MCV, WBC, MCHC, hematocrit, RDW, anion gap, bicarbonate, bun, creatinine, calcium, chloride, sodium, potassium, PT, INR, MECHVENT, RRT, hematocrit, and the cumulative urine output within 24 h differ significantly compared to those who survived. However, dementia, diabetes and non-diabetes, peripheral vascular disease, chronic pulmonary disease, rheumatic disease, peptic ulcer disease, body temperature, mean blood pressure, diastolic blood pressure, systolic blood pressure, creatine kinase MB isoenzyme, RBC, hemoglobin, and MCH showed no significant difference between the two groups. Figure 1 is the flow chart describing the procedure for subject selection.
Table 1. Baseline characteristics, vital signs, laboratory parameters and statistic results of mimic-IIV patients with HF combined with hypertension.
Features selected in models
The LASSO regularization process resulted in 22 potential predictors based on 2,766 patients in the training cohort (Figures 2A,B). Using MLP, NN, RF, NB, and LR, the 58 selected variables were used to identify patients who died during their hospital stay in the training cohort. We show the proportional importance of the top 22-ranked input variables in the NN model and the LR model, respectively. Figure 2B shows the LASSO-selected predictors (shrinkage parameter, λ = 0.01914052).
Figure 2. Demographic and clinical feature selection using the least absolute shrinkage and selection operator (LASSO) binary logistic regression model. (A) Tuning parameter (λ) selection in the LASSO model used 10- fold cross-validation via minimum criteria. The partial likelihood deviance (binomial deviance) curve was plotted versus log(λ). Dotted vertical lines were drawn at the optimal values by using the minimum criteria and the one SE of the minimum criteria (the 1-SE criteria). λ value of 0.01914052, with log(λ), −3.9545 was chosen (1- SE criteria) according to 10-fold cross-validation. (B) LASSO coefficient profiles of the 58 features. A coefficient profile plot was produced against the log(λ) sequence. The vertical line was drawn at the value selected using 10-fold cross-validation, where optimal resulted in 20 features with non-zero coefficients.
Development and comparison of machine learning models
A total of 58 clinical features were collected during the first 24 h after ICU admission. KNN was used to impute missing data. LASSO regression was employed to identify signature variables for hospital mortality in patients with HF combined with hypertension. Ultimately, 22 out of 58 clinical features were associated with prognosis, and these results are presented in Table 2. In addition, we have constructed five ML binary classifiers, namely MLP, NN, RF, NB, and LR, to predict the risk of death in HF patients with hypertension (Figure 1). Then the obtained hyperparameters were used to train the ML model with the entire training data, and the performance of the model was evaluated using the testing cohort. Figure 4 and Table 3 describe the performance of these predictive models, showing that the NN model with all available variables relatively outperformed the other four models or predictive factors in testing cohorts with an AUC of 0.764 and an accuracy of 0.8731 in the testing cohort, compared with the other ML models (AUC: LR, 0.640; RF, 0.748; NB, 0.751; MLP, 0.730). Therefore, we selected the NN as the most promising approach among the five ML algorithms for further prediction in this study.
Figure 3. SHAP summary plot for the top 22 clinical features contributing to the NN (A) and LR (B) model. SHAP feature importance is measured as one minus AUC loss after permutations. This matrix plot depicts the importance of each covariate in the development of the final predictive model.
Figure 4. The receiver operator characteristic (ROC) curves for the ML models predict in-hospital mortality (The training cohort and testing cohort). ROC curves for five ML models predicting in-hospital mortality in the training (A) and testing cohort (B), respectively; (C) ROC curves for in-hospital mortality in the test set predicted by the un-simplified NN model and the ROC curves for in-hospital mortality in the test set and external validation set predicted by the simplified NN model, respectively.
Table 3. Predictive performances of the five machine learning models for predicting in-hospital mortality.
Significant predictors and development of the simplified model
We identified the 22 most significant predictors by permuting feature importance techniques, transforming the NN model and LR model into a universally applicable prediction model (Figure 3). The features specific to death included age, MECHVENT, chloride, bun, anion gap, paraplegia, RDW, hyperlipidemia, SpO2, respiratory rate, cerebrovascular disease, HR, WBC, INR, MCHC, glucose, aids, MCV, NT-proBNP, calcium, RRT, and PTT. For the convenience of clinical providers and patients, the order of importance of these features is not the same in the two models. Furthermore, we also assembled a simplified ML model for HF risk stratification by artificial intelligence with the top seven most important high-ranking and readily available variables, namely bun, paraplegia, anion gap, RDW, MECHVENT, chloride, and age. The AUC of SMART-HF reached 0.741 and 0.674 in the test cohort and external validation cohort, respectively (Figure 4C).
Shapley additive explanations values depending on variables
The impact of the top seven factors on the NN model’s mortality risk prediction was further explored using the SHAP dependency plot. The probability of death in ICU patients with HF combined with hypertension increases with the following indicators: age, increased anion gap, elevated BUN, elevated RDW, increased serum chloride level, and paraplegia (Figure 5). In contrast, the mortality risk decreases with the increase of the MECHVENT index.
Discussion
Given that HF patients with hypertension have a relatively high chance of ICU admission, prioritization of patients who require higher levels of care or immediate medical attention is critical. Accurately predicting prognosis is the foundation of both patient-centered care and shared decision-making, such as selecting treatment strategies and informing patients. The present study has shown the potential for the NN model to assist physicians with predicting 28-day all-cause in-hospital mortality. Using data derived from the MIMIC-IV database, this study explored and validated five real-time diagnostic and prognostic prediction models based on a ML algorithm for 28-day in-hospital mortality in ICU patients with HF combined with hypertension. These learning models incorporated static and dynamic variables. Ultimately, it was found that an NN model best-stratified patients’ risks with good external validation. The algorithm showed that age, MECHVENT, chloride, BUN, anion gap, paraplegia, RDW, hyperlipidemia, SpO2, respiratory rate, cerebrovascular disease, HR, WBC, INR, MCHC, glucose, AIDS, MCV, NT-proBNP, calcium, RRT, and PTT were associated with an increased risk of death. The NN prediction model may facilitate clinical decision-making for advanced management of ICU-admitted HF patients with hypertension. The proposed model has several advantages over traditional clinical risk models. First, a highly specific cohort of ICU patients with HF combined with hypertension we used, rather than a more generalized cohort such as patients admitted to the general medical floor. Second, a simplified ML model for HF risk stratification was also assembled by artificial intelligence with the top seven most important high-ranking variables to avoid the collection of a large number of variables for a prognostic model in real clinical settings. In addition, improving the model predictive performance can directly improve risk assessment even before patients receive a more comprehensive diagnostic evaluation in the ICU.
Hospitalized patients with HF often require admission to the ICU, especially when their condition is complicated by various comorbidities such as hypertension. Data from 341 hospitals in the USA showed a median ICU admission rate of 10% (IQR, 6–16%) for hospitalized HF patients (8). Numerous studies have shown that in-hospital mortality in patients with advanced HF admitted to the ICU is significantly higher than that of HF patients admitted to hospital wards only. All-cause in-hospital mortality was 10.6% for HF patients admitted to the ICU, compared with 4.0% for all HF patients in the ADHERE study (11). Meanwhile, 17.3% mortality among ICU patients versus 6.5% among all hospitalized HF patients was reported in the RO-AHFS study, and 17.8% death rate among ICU patients versus 4.5% among all hospitalized HF patients was observed in the ALARM-HF study (37, 38). It is noteworthy that a decision on whether a patient with HF requires intensive care depends on both clinical judgment and resource availability, which adds unmeasured differences to outcome studies.
Our study population comprised ICU-admitted HF patients with hypertension, and the in-hospital mortality rate was 15.3% (n = 459 patients). This rate was substantially higher than other prediction models for in-hospital mortality based on all HF patients, regardless of ICU admission. However, in the ADHERE in-hospital mortality risk stratification model, the in-hospital mortality rate of their study population was only 4.2% (27). In the optimized heart failure prediction model, the GWTG-HF risk scoring model, and a single-center elderly Chinese patient-based model established by Jia et al., the rates were 3.8, 2.86, and 5.58%, respectively (28, 29, 39).
Our model contains only variables that are easily accessible: its simplicity makes risk prediction applicable for different purposes during the hospitalization of HF patients in the ICU. For example, when the calculated risk of death for an individual is high, it indicates the need for more aggressive monitoring or resource allocation, which can help assign patients to different levels of care. This is particularly useful when healthcare resources are limited. The discriminative performance of the model is very high, and its validation was confirmed by testing the model in a cohort of HF patients with hypertension in the eICU database.
Although several published studies provide a wealth of computational tools or predictive models that can be easily used in a variety of settings to assess risk in patients with HF (21, 24, 40–45). Such calculators unfortunately require tedious data entry. Real-time processing of the predictive model directly from the Electronic Health Record (EHR) provides immediate and seamless calculation, and the score from this calculation is well suited to support clinical decision-making and prioritization when the healthcare system is overloaded. An accurate prognosis is a basis for many clinical decisions regarding patients admitted to the ICU with HF (24). To avoid the shortcomings of using traditional LR analyses such as overfitting and predictor variables with skewed distributions, data on demographic characteristics, vital signs, comorbidities, and laboratory variables were used in the present study for LASSO regression analysis to screen for independent risk factors for in-hospital mortality.
Other predictive models have been published previously, and many variables have been reported to correlate with mortality in HF patients. Variables, such as gender, BUN level, BMI, age, sodium levels, health status, systolic blood pressure, diabetes mellitus, serum creatinine levels, low SBP, chronic obstructive lung disease, NYHA (New York Heart Association) classification, left ventricular ejection fraction (LVEF), smoking, not receiving ACEIs/ARBs (Angiotensin-Converting Enzyme Inhibitors/Angiotensin II Receptor Blockers), and not receiving beta-blockers, have been reported to explain the predictive model (25, 44). Consistent with previous studies, our study identified age as a strong prognostic predictor. When HF worsens, especially in the elderly, it can lead to severe ischemia and hypoxia, respiratory failure, and ultimately death. In our model, BUN levels also substantially contributed to the predicted probabilities, and elevated BUN levels substantially contributed to increased in-hospital mortality, which is consistent with previously published studies (27–29). Elevated BUN levels suggest a high probability of prerenal injury, which may be related to reduced renal blood supply due to insufficient effective blood volume or a decrease in cardiac output after the onset of HF, as well as fluid management (urine output). Some HF mortality prediction models believe that HR affects prognosis strongly (28, 30), while some other models disagree (27, 46, 47). In our study, HR was included in the final model, in contrast to BMI, which was not included in the final model because our study failed to prove that BMI was a predictor of in-hospital mortality of ICU-admitted HF patients with hypertension. This may be related to differences in study populations and our relatively small sample size. Whether the “obesity paradox” biases our results is uncertain, as the “obesity paradox” presented for both critical care-related outcomes (48) and HF (49). Therefore, further studies are still needed to elaborate on the effect of obesity or BMI.
Our study showed that chloride and hyperlipidemia correlated with an unfavorable outcome. Consistent with previous findings, chloride has a more prominent contribution to the pathophysiology and affects the prognosis of HR, which may be related to hyperchloremia, acidosis, inflammation, and renal injury secondary to hyper chlorination. In clinical practice, it is uncommon to focus on this one marker alone in patients with HF, but it is often integrated with blood gas analysis and evaluation of anion gap (50). Although several studies in acute and chronic HF populations have demonstrated the prognostic value of hypochlorhydria, interventional clinical trials exploring serum chloride as a therapeutic target have been inconclusive to date. Ongoing prospective randomized controlled studies may shed light on the role of serum chloride as a therapeutic target to improve outcomes in patients with HF. These studies should clarify whether serum chloride should be included in current models for predicting prognosis in HF (51). In the general population and patients with atherosclerotic cardiovascular disease, hypercholesterolemia has consistently been shown to be associated with poor outcomes, including mortality, cardiovascular events, and the development of HF (52, 53). Conversely, in patients with established HF, several analyses have now demonstrated an inverse relationship between cholesterol levels and outcomes. That is, low cholesterol levels have been shown to be independently associated with increased mortality, while higher cholesterol levels have been associated with improved survival. It is unclear whether low cholesterol levels play a pathogenic role in adverse outcomes in patients with HF or whether low cholesterol levels simply reflect advanced disease status. Given the observed inverse relationship between cholesterol levels and mortality in patients with HF, the applicability of cholesterol treatment goals recommended for the general population and patients with atherosclerotic CVD to patients with HF is unclear and remains to be determined (54). Consistent with previous studies (39, 55), our study also identified the anion gap as a strong prognostic predictor and found that this factor was independently associated with an increased risk of death.
Our study also showed that the anion gap was also a very important prognostic feature. The anion gap formula (AnionGap=SNa+SK−SHCO3_−SCl–) itself already shows an interdependent relationship between serum chloride, sodium, potassium, and bicarbonate (50), so it is not surprising why hypochlorhydria often occurs in conjunction with metabolic alkalosis. It is known as “chloride-depleted alkalosis,” a state of extracellular fluid volume constriction caused primarily by diuretic-induced diuresis. Although the exact role of pH as a prognostic indicator of HF has not been fully explored, pH is influenced by chloride levels in the form of chloride-depleted alkalosis or hyperchloremic metabolic acidosis (56). Chloride-depleted alkalosis is an independent predictor of in-hospital mortality in patients with decompensated HF (57). In HF, electrolyte depletion is primarily the result of salt restriction and cyclic and thiazide diuretic therapy, whereas metabolic alkalosis often occurs as a result of diuretic usage (58).
Notably, our study found that the use of MECHVENT improved the prognosis of patients with HF (Figure 5). The use of MECHVENT often indicates that patients are in serious conditions, such as the occurrence of acute HF and respiratory failure, but the use of MECHVENT as treatment improves the patient’s prognosis by rapidly improving respiratory symptoms and ventilatory function compared to conventional drug therapy. The previous study has shown that in-hospital mortality in HF patients receiving non-invasive ventilation (NIV) or non-invasive ventilation (NIV) + invasive mechanical ventilation (IMV) decreases significantly over time, even if the clinical profile is worsening (59).
In our model, paraplegia and high RDW were also high-risk factors for ICU-admitted patients with HR combined with hypertension. When paralysis, prolonged bed rest, and reduced activity occur, the likelihood of thrombosis, crushing pneumonia, and infection increases. Consistent with the previous study (60), RDW is a powerful predictor of poor long-term outcomes in HF patients with acute exacerbation (AHF), and its prognostic value outperforms that of other well-established risk variables or biomarkers.
In addition, we further enhanced the readability of the model by using the SHAP framework, making the individual variables that contribute to the overall prediction easily available and understandable to physicians in real-time, along with the model’s risk score.
In comparison with the reported GWTG-HF risk score (29), a well-validated tool for predicting in-hospital mortality in HF patients, both the NN model and the LR model showed superiority in predictive power in our study population. Both models showed good discrimination and calibration power in both the derivation and validation sets. To obtain a more concise and broader range of net benefit threshold probabilities, we chose the NN model to develop our simplified ML model.
Neural networks are constructed from basic units called neurons, which can be easily arranged into layers. The layers are easily connected, and the entire network can be trained end-to-end using a stochastic gradient descent algorithm (61). While single-layer networks can approximate any function to arbitrary accuracy [as implied by the general approximation theorem (62)], the real power of these models lies in providing automatic abstraction by stacking multiple layers into deep neural networks (DNNs). Each layer abstracts its input, providing the next layer with a representation of the data that is more likely to work within the scope of the task being solved. The most advanced models of NNs are DNNs models and have been shown to provide superhuman performance (63, 64) on many tasks involving difficult to abstract data, such as those involving image and audio processing.
In recent literature (65, 66), shallow NNs have been used to predict mortality in heart failure despite training on unbalanced datasets, showing superior performance to other learning methods. DNNs have been used to predict mortality (67, 68) or the risk of heart failure and acute heart failure (69). The authors of these two works compared DNNs with other ML techniques and showed improved performance. Another interesting recent application of NNs in this field is exploiting their ability to process very complex and related data. This is the case with the Deep Cox Mixture Model (70), where NNs assist the Cox Regression Model to fit the risk ratio of the regression. This work is based on a sound statistical and ML background, is fully disclosed, and provides state-of-the-art performance when working with diverse groups of individuals.
Compared to traditional ML algorithms such as logistic regression, neural networks typically require more data, at least thousands or even millions of labeled samples. This is not an easy problem to solve, but if other algorithms are used, the related ML problem can be solved with less data. At the same time, compared to traditional algorithms, neural networks are computationally more expensive than traditional algorithms. Advanced deep learning algorithms can take weeks to train successfully. Whereas most traditional ML takes less than a few minutes, hours, or days. Of course, the computational power required by NN depends heavily on the size of the data and also on the depth and complexity of the network. The smaller the data set as well as the smaller the depth and complexity of the network, the less computational power is required by the NN.
This study has several limitations: First, as a retrospective study, selection bias was hardly avoidable. However, the inclusion criteria were set strictly so that the cases included in the study reflected the actual conditions as accurately as possible. Second, data were collected from patient medical records and the final performance of our predictive model was strictly dependent on the accuracy of the records. Third, as a single-center study, the scope and number of study populations involved were relatively small. Fourth, the impact of this predictive model on routine patient care in different clinical settings has not been well investigated. Therefore, data with mortality outcomes from other independent healthcare systems will be required to fully assess its generalizability. Finally, our model may only facilitate the rapid identification of critical clinical situations at the bedside but does not provide additional information about the underlying life-threatening pathophysiological mechanisms.
Conclusion
Using ML techniques, we developed a predictive model to predict the 28-day mortality for ICU-admitted patients with HF combined with hypertension based on NN. It was proved with a better predictive value than the traditional logistic regression analysis. With a high AUC of 0.764 and an accuracy of 0.8731 in the testing cohort, this model is promising for routine use in the ICUs to automatically warn the staff at any stage of the disease. The SHAP method enables ML models to be more interpretable, thereby helping clinicians to better understand the reasoning behind the outcome and evaluate in-hospital outcomes for critically ill patients, especially those with uncertain survival outcomes. Also, the model involves a small number of routinely collected variables that can be easily used at the bedside.
Data availability statement
The original contributions presented in this study are included in the article/supplementary material, further inquiries can be directed to the corresponding author.
Ethics statement
The studies involving human participants were reviewed and approved by the Medical Ethics Committee at Zi Gong First People’s Hospital. Written informed consent for participation was not required for this study in accordance with the national legislation and the institutional requirements.
Author contributions
SgP, JH, XzL, JD, and SaP contributed to the conception and design. SgP, JH, JD, and HC collected and analyzed the data. JH drew the figures and tables. SgP wrote the draft. CS, JT, WC, WW, and XgL contributed to the manuscript writing and revision. All authors approved the final manuscript.
Funding
This work was supported by the Zigong Science and Technology Project (No. 2015ZC02) and the Medical and Health Own Fund Guiding Project of the Zigong Academy of Medical Sciences (No. 2020YXY03).
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
References
1. GBD 2017 Disease and Injury Incidence and Prevalence Collaborators. Global, regional, and national incidence, prevalence, and years lived with disability for 354 diseases and injuries for 195 countries and territories, 1990-2017: a systematic analysis for the global burden of disease study 2017. Lancet. (2018) 392:1789–858. doi: 10.1016/S0140-6736(18)32279-7
2. Conrad N, Judge A, Canoy D, Tran J, Pinho-Gomes AC, Millett ERC, et al. Temporal trends and patterns in mortality after incident heart failure: a longitudinal analysis of 86 000 individuals. JAMA Cardiol. (2019) 4:1102–11. doi: 10.1001/jamacardio.2019.3593
3. Taylor CJ, Ordóñez-Mena JM, Roalfe AK, Lay-Flurrie S, Jones NR, Marshall T, et al. Trends in survival after a diagnosis of heart failure in the United Kingdom 2000-2017: population based cohort study. BMJ. (2019) 364:l223. doi: 10.1136/bmj
4. Tromp J, Ferreira JP, Janwanishstaporn S, Shah M, Greenberg B, Zannad F, et al. Heart failure around the world. Eur J Heart Fail. (2019) 21:1187–96. doi: 10.1002/ejhf.1585
5. Heidenreich PA, Albert NM, Allen LA, Bluemke DA, Butler J, Fonarow GC, et al. Forecasting the impact of heart failure in the united states: a policy statement from the American heart association. Circ Heart Fail. (2013) 6:606–19. doi: 10.1161/HHF.0b013e318291329a
6. Benjamin EJ, Virani SS, Callaway CW, Chamberlain AM, Chang AR, Cheng S, et al. Heart Disease and stroke statistics-2018 update: a report from the American heart association. Circulation. (2018) 137:e67–492. doi: 10.1161/CIR.0000000000000558
7. Krum H, Gilbert RE. Demographics and concomitant disorders in heart failure. Lancet. (2003) 362:147–58. doi: 10.1016/S0140-6736(03)13869-X
8. Safavi KC, Dharmarajan K, Kim N, Strait KM, Li SX, Chen SI, et al. Variation exists in rates of admission to intensive care units for heart failure patients across hospitals in the United States. Circulation. (2013) 127:923–9. doi: 10.1161/CIRCULATIONAHA.112.001088
9. van Diepen S, Bakal JA, Lin M, Kaul P, McAlister FA, Ezekowitz JA. Variation in critical care unit admission rates and outcomes for patients with acute coronary syndromes or heart failure among high- and low-volume cardiac hospitals. J Am Heart Assoc. (2015) 4:e001708. doi: 10.1161/JAHA.114.001708
10. Wunsch H, Angus DC, Harrison DA, Collange O, Fowler R, Hoste EA, et al. Variation in critical care services across North America and Western Europe. Crit Care Med. (2008) 36:2787–e9. doi: 10.1097/CCM.0b013e318186aec8
11. Adams KF Jr., Fonarow GC, Emerman CL, LeJemtel TH, Costanzo MR, Abraham WT, et al. Characteristics and outcomes of patients hospitalized for heart failure in the United States: rationale, design, and preliminary observations from the first 100,000 cases in the acute decompensated heart failure national registry (ADHERE). Am Heart J. (2005) 149:209–16. doi: 10.1016/j.ahj.2004.08.005
12. Goehner NS, Roman PEF. Monitoring technologies in acute care environments: a comprehensive guide to patient monitoring technology. Can J Anesth. (2014) 61:891–2. doi: 10.1007/s12630-014-0177-y
13. Fackler JC, Watts C, Grome A, Miller T, Crandall B, Pronovost P. Critical care physician cognitive task analysis: an exploratory study. Crit Care. (2009) 13:R33. doi: 10.1186/cc7740
14. Wright MC, Dunbar S, Macpherson BC, Moretti EW, Del Fiol G, Bolte J, et al. Toward designing information display to support critical care. A qualitative contextual evaluation and visioning effort. Appl Clin Inform. (2016) 7:912–29. doi: 10.4338/ACI-2016-03-RA-0033
15. Duke G, Green J, Briedis J. Survival of critically ill medical patients is time-critical. Crit Care Resusc. (2004) 6:261–7.
16. Numata Y, Schulzer M, van der Wal R, Globerman J, Semeniuk P, Balka E, et al. Nurse staffing levels and hospital mortality in critical care settings: literature review and meta-analysis. J Adv Nurs. (2006) 55:435–48. doi: 10.1111/j.1365-2648.2006.03941.x
17. Falk AC, Wallin EM. Quality of patient care in the critical care unit in relation to the nurse-patient ratio: a descriptive study. Intensive Crit Care Nurs. (2016) 35:74–9. doi: 10.1016/j.iccn.2016.01.002
18. Wallace DJ, Angus DC, Barnato AE, Kramer AA, Kahn JM. Nighttime intensivist staffing and mortality among critically ill patients. N Engl J Med. (2012) 366:2093–101. doi: 10.1056/NEJMsa1201918
19. Ponikowski P, Voors AA, Anker SD, Bueno H, Cleland JG, Coats AJ, et al. 2016 ESC Guidelines for the diagnosis and treatment of acute and chronic heart failure: the task force for the diagnosis and treatment of acute and chronic heart failure of the European society of cardiology (ESC). Developed with the special contribution of the heart failure association (HFA) of the ESC. Eur J Heart Fail. (2016) 18:891–975. doi: 10.1002/ejhf.592
20. Yancy CW, Jessup M, Bozkurt B, Butler J, Casey DE Jr., Drazner MH, et al. 2013 ACCF/AHA guideline for the management of heart failure: a report of the American college of cardiology foundation/American heart association task force on practice guidelines. Circulation. (2013) 128:e240–327. doi: 10.1161/CIR.0b013e31829e8776
21. Frizzell JD, Liang L, Schulte PJ, Yancy CW, Heidenreich PA, Hernandez AF, et al. Prediction of 30-Day All-cause readmissions in patients hospitalized for heart failure: comparison of machine learning and other statistical approaches. JAMA Cardiol. (2017) 2:204–9. doi: 10.1001/jamacardio.2016.3956
22. Wang W, Cai D. Complement components sC5b-9 and CH50 predict prognosis in heart failure patients combined with hypertension. Am J Hypertens. (2020) 33:53–60. doi: 10.1093/ajh/hpz140
23. Tohyama T, Ide T, Ikeda M, Kaku H, Enzan N, Matsushima S, et al. Machine learning-based model for predicting 1 year mortality of hospitalized patients with heart failure. ESC Heart Fail. (2021) 8:4077–85. doi: 10.1002/ehf2.13556
24. Li F, Xin H, Zhang J, Fu M, Zhou J, Lian Z. Prediction model of in-hospital mortality in intensive care unit patients with heart failure: machine learning-based, retrospective analysis of the MIMIC-III database. BMJ Open. (2021) 11:e044779. doi: 10.1136/bmjopen-2020-044779
25. Pocock SJ, Ariti CA, McMurray JJ, Maggioni A, Køber L, Squire IB, et al. Predicting survival in heart failure: a risk score based on 39 372 patients from 30 studies. Eur Heart J. (2013) 34:1404–13. doi: 10.1093/eurheartj/ehs337
26. Collier TJ, Pocock SJ, McMurray JJ, Zannad F, Krum H, van Veldhuisen DJ. The impact of eplerenone at different levels of risk in patients with systolic heart failure and mild symptoms: insight from a novel risk score for prognosis derived from the EMPHASIS-HF trial. Eur Heart J. (2013) 34:2823–9. doi: 10.1093/eurheartj/eht247
27. Fonarow GC, Adams KF Jr., Abraham WT, Yancy CW, Boscardin WJ, Adhere Scientific Advisory Committee, Study Group, and Investigators. Risk stratification for in-hospital mortality in acutely decompensated heart failure: classification and regression tree analysis. JAMA. (2005) 293:572–80. doi: 10.1001/jama.293.5.572
28. Abraham WT, Fonarow GC, Albert NM, Stough WG, Gheorghiade M, Greenberg BH, et al. Predictors of in-hospital mortality in patients hospitalized for heart failure: insights from the organized program to initiate lifesaving treatment in hospitalized patients with heart failure (OPTIMIZE-HF). J Am Coll Cardiol. (2008) 52:347–56. doi: 10.1016/j.jacc.2008.04.028
29. Peterson PN, Rumsfeld JS, Liang L, Albert NM, Hernandez AF, Peterson ED, et al. A validated risk score for in-hospital mortality in patients with heart failure from the American heart association get with the guidelines program. Circ Cardiovasc Qual Outcomes. (2010) 3:25–32. doi: 10.1161/CIRCOUTCOMES.109.854877
30. Yagyu T, Kumada M, Nakagawa T. Novel risk stratification with time course assessment of in-hospital mortality in patients with acute heart failure. PLoS One. (2017) 12:e0187410. doi: 10.1371/journal.pone.0187410
31. Krumholz HM. Big data and new knowledge in medicine: the thinking, training, and tools needed for a learning health system. Health Aff (Millwood). (2014) 33:1163–70. doi: 10.1377/hlthaff.2014.0053
32. Johnson AE, Pollard TJ, Shen L, Lehman LW, Feng M, Ghassemi M, et al. MIMIC-III, a freely accessible critical care database. Sci Data. (2016) 3:160035. doi: 10.1038/sdata.2016.35
33. Pollard TJ, Johnson AEW, Raffa JD, Celi LA, Mark RG, Badawi O. The eICU Collaborative Research Database, a freely available multi-center database for critical care research. Sci Data. (2018) 5:180178. doi: 10.1038/sdata.2018.178
34. Zhou S, Zeng Z, Wei H, Sha T, An S. Early combination of albumin with crystalloids administration might be beneficial for the survival of septic patients: a retrospective analysis from MIMIC-IV database. Ann Intensive Care. (2021) 11:42. doi: 10.1186/s13613-021-00830-8
35. Deshmukh F, Merchant SS. Explainable machine learning model for predicting GI bleed mortality in the intensive care unit. Am J Gastroenterol. (2020) 115:1657–68. doi: 10.14309/ajg.0000000000000632
36. Gromski PS, Xu Y, Kotze HL, Correa E, Ellis DI, Armitage EG, et al. Influence of missing values substitutes on multivariate analysis of metabolomics data. Metabolites. (2014) 4:433–52. doi: 10.3390/metabo4020433
37. Chioncel O, Ambrosy AP, Filipescu D, Bubenek S, Vinereanu D, Petris A, et al. Patterns of intensive care unit admissions in patients hospitalized for heart failure: insights from the RO-AHFS registry. J Cardiovasc Med (Hagerstown). (2015) 16:331–40. doi: 10.2459/JCM.0000000000000030
38. Follath F, Yilmaz MB, Delgado JF, Parissis JT, Porcher R, Gayat E, et al. Clinical presentation, management and outcomes in the acute heart failure global survey of standard treatment (ALARM-HF). Intensive Care Med. (2011) 37:619–26. doi: 10.1007/s00134-010-2113-0
39. Jia Q, Wang YR, He P, Huang XL, Yan W, Mu Y, et al. Prediction model of in-hospital mortality in elderly patients with acute heart failure based on retrospective study. J Geriatr Cardiol. (2017) 14:669–78. doi: 10.11909/j.issn.1671-5411.2017.11.002
40. Levy WC, Mozaffarian D, Linker DT, Sutradhar SC, Anker SD, Cropp AB, et al. The Seattle heart failure model: prediction of survival in heart failure. Circulation. (2006) 113:1424–33. doi: 10.1161/CIRCULATIONAHA.105.584102
41. Rich JD, Burns J, Freed BH, Maurer MS, Burkhoff D, Shah SJ. Meta-analysis global group in chronic (MAGGIC) heart failure risk score: validation of a simple tool for the prediction of morbidity and mortality in heart failure with preserved ejection fraction. J Am Heart Assoc. (2018) 7:e009594. doi: 10.1161/JAHA.118.009594
42. Komajda M, Carson PE, Hetzel S, McKelvie R, McMurray J, Ptaszynska A, et al. Factors associated with outcome in heart failure with preserved ejection fraction: findings from the irbesartan in heart failure with preserved ejection fraction study (I-PRESERVE). Circ Heart Fail. (2011) 4:27–35. doi: 10.1161/CIRCHEARTFAILURE.109.932996
43. Mortazavi BJ, Downing NS, Bucholz EM, Dharmarajan K, Manhapra A, Li SX, et al. Analysis of machine learning techniques for heart failure readmissions. Circ Cardiovasc Qual Outcomes. (2016) 9:629–40. doi: 10.1161/CIRCOUTCOMES.116.003039
44. Angraal S, Mortazavi BJ, Gupta A, Khera R, Ahmad T, Desai NR, et al. Machine learning prediction of mortality and hospitalization in heart failure with preserved ejection fraction. JACC Heart Fail. (2020) 8:12–21. doi: 10.1016/j.jchf.2019.06.013
45. Cui X, Zhou J, Pivodic A, Dahlström U, Ge J, Fu M. Temporal trends in cause-specific readmissions and their risk factors in heart failure patients in Sweden. Int J Cardiol. (2020) 306:116–22. doi: 10.1016/j.ijcard.2020.02.048
46. Voors AA, Ouwerkerk W, Zannad F, van Veldhuisen DJ, Samani NJ, Ponikowski P, et al. Development and validation of multivariable models to predict mortality and hospitalization in patients with heart failure. Eur J Heart Fail. (2017) 19:627–34. doi: 10.1002/ejhf.785
47. Lupón J, de Antonio M, Vila J, Peñafiel J, Galán A, Zamora E, et al. Development of a novel heart failure risk tool: the Barcelona bio-heart failure risk calculator (BCN bio-HF calculator). PLoS One. (2014) 9:e85466. doi: 10.1371/journal.pone.0085466
48. Patel JJ, Rosenthal MD, Miller KR, Codner P, Kiraly L, Martindale RG. The critical care obesity paradox and implications for nutrition support. Curr Gastroenterol Rep. (2016) 18:45. doi: 10.1007/s11894-016-0519-8
49. Horwich TB, Fonarow GC, Clark AL. Obesity and the obesity paradox in heart failure. Prog Cardiovasc Dis. (2018) 61:151–6. doi: 10.1016/j.pcad.2018.05.005
50. Zandijk AJL, van Norel MR, Julius FEC, Sepehrvand N, Pannu N, McAlister FA, et al. Chloride in heart failure: the neglected electrolyte. JACC Heart Fail. (2021) 9:904–15. doi: 10.1016/j.jchf.2021.07.006
51. Rivera FB, Alfonso P, Golbin JM, Lo K, Lerma E, Volgman AS, et al. The role of serum chloride in acute and chronic heart failure: a narrative review. Cardiorenal Med. (2021) 11:87–98. doi: 10.1159/000515604
52. Pekkanen J, Linn S, Heiss G, Suchindran CM, Leon A, Rifkind BM, et al. Ten-year mortality from cardiovascular disease in relation to cholesterol level among men with and without preexisting cardiovascular disease. N Engl J Med. (1990) 322:1700–7. doi: 10.1056/NEJM199006143222403
53. Orso F, Fabbri G, Maggioni AP. Epidemiology of heart failure. Handb Exp Pharmacol. (2017) 243:15–33. doi: 10.1007/164_2016_74
54. Bozkurt B, Aguilar D, Deswal A, Dunbar SB, Francis GS, Horwich T, et al. Contributory risk and management of comorbidities of hypertension, obesity, diabetes mellitus, hyperlipidemia, and metabolic syndrome in chronic heart failure: a scientific statement from the American heart association. Circulation. (2016) 134:e535–78. doi: 10.1161/CIR.0000000000000450
55. Park JJ, Choi DJ, Yoon CH, Oh IY, Lee JH, Ahn S, et al. The prognostic value of arterial blood gas analysis in high-risk acute heart failure patients: an analysis of the Korean heart failure (KorHF) registry. Eur J Heart Fail. (2015) 17:601–11. doi: 10.1002/ejhf.276
56. Kazory A, Ronco C. Emergence of chloride as an overlooked cardiorenal connector in heart failure. Blood Purif. (2020) 49:219–21. doi: 10.1159/000503774
57. Khan NN, Nabeel M, Nan B, Ghali JK. Chloride depletion alkalosis as a predictor of in-hospital mortality in patients with decompensated heart failure. Cardiology. (2015) 131:151–9. doi: 10.1159/000377669
58. Elfar A, Sambandam KK. The basic metabolic profile in heart failure-marker and modifier. Curr Heart Fail Rep. (2017) 14:311–20. doi: 10.1007/s11897-017-0344-x
59. de Miguel-Díez J, Jiménez-García R, Méndez-Bailón M, Muñoz-Rivas N, Hernández-Barrera V, Puente-Maestu L, et al. National trends in mechanical ventilation among patients hospitalized with heart failure: a population-based study in Spain (2001-2017). Eur J Intern Med. (2020) 78:76–81. doi: 10.1016/j.ejim.2020.04.016
60. Melchio R, Rinaldi G, Testa E, Giraudo A, Serraino C, Bracco C, et al. Red cell distribution width predicts mid-term prognosis in patients hospitalized with acute heart failure: the RDW in acute heart failure (RE-AHF) study. Intern Emerg Med. (2019) 14:239–47. doi: 10.1007/s11739-018-1958-z
61. Robbins H, Monro S. A stochastic approximation method. Ann Math Stat. (1951) 22:400–7. doi: 10.1007/978-1-4612-5110-1_9
62. Baker MR, Patil RB. Universal approximation theorem for interval neural networks. Reliab Comput. (1998) 4:235–9. doi: 10.1023/A:1009951412412
63. Silver D, Huang A, Maddison CJ, Guez A, Sifre L, van den Driessche G, et al. Mastering the game of go with deep neural networks and tree search. Nature. (2016) 529:484–9. doi: 10.1038/nature16961
64. He K, Zhang X, Ren S, Sun J. Delving deep into rectifiers: surpassing human-level performance on imagenet classification. In: Proceedings of the 2015 IEEE International Conference on Computer Vision (ICCV). Santiago: (2015). p. 1026–34. doi: 10.1109/ICCV.2015.123
65. Awan SE, Bennamoun M, Sohel F, Sanfilippo FM, Dwivedi G. Machine learning-based prediction of heart failure readmission or death: implications of choosing the right model and the right metrics. ESC Heart Fail. (2019) 6:428–35. doi: 10.1002/ehf2.12419
66. Awan SE, Bennamoun M, Sohel F, Sanfilippo FM, Chow BJ, Dwivedi G. Feature selection and transformation by machine learning reduce variable numbers and improve prediction for heart failure readmission or death. PLoS One. (2019) 14:e0218760. doi: 10.1371/journal.pone.0218760
67. Mancia G, Facchetti R, Seravalle G, Cuspidi C, Grassi G. Adding home and/or ambulatory blood pressure to office blood pressure for cardiovascular risk prediction. Hypertension. (2021) 77:640–9. doi: 10.1161/HYPERTENSIONAHA.120.16303
68. Jallepalli S, Pathak P, Gupta P, Kar S, Gupta M. Development and Validation of Artificial Intelligence – Based Cardiovascular Disease (ai-cvd) Risk Score. Rochester, NY: Social Science Electronic Publishing (2019). doi: 10.2139/ssrn.3444410
69. Kwon JM, Kim KH, Jeon KH, Lee SE, Lee HY, Cho HJ. Artificial intelligence algorithm for predicting mortality of patients with acute heart failure. PLoS One. (2019) 14:e0219302. doi: 10.1371/journal.pone.0219302
Keywords: MIMIC-IV, interpretable machine learning, neural networks, heart failure, hypertension
Citation: Peng S, Huang J, Liu X, Deng J, Sun C, Tang J, Chen H, Cao W, Wang W, Duan X, Luo X and Peng S (2022) Interpretable machine learning for 28-day all-cause in-hospital mortality prediction in critically ill patients with heart failure combined with hypertension: A retrospective cohort study based on medical information mart for intensive care database-IV and eICU databases. Front. Cardiovasc. Med. 9:994359. doi: 10.3389/fcvm.2022.994359
Received: 14 July 2022; Accepted: 15 September 2022;
Published: 12 October 2022.
Edited by:
Roney Orismar Sampaio, Heart Institute, University of São Paulo, BrazilReviewed by:
Zeyi Cheng, Sichuan University, ChinaAn-Tian Chen, Chinese Academy of Medical Sciences and Peking Union Medical College, China
Copyright © 2022 Peng, Huang, Liu, Deng, Sun, Tang, Chen, Cao, Wang, Duan, Luo and Peng. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Shuang Peng, MTQ1MTg2ODA5OUBxcS5jb20=
†These authors have contributed equally to this work