Skip to main content

ORIGINAL RESEARCH article

Front. Cell. Infect. Microbiol., 19 April 2022
Sec. Virus and Host

Statistical Analysis and Machine Learning Prediction of Disease Outcomes for COVID-19 and Pneumonia Patients

Yu Zhao,&#x;Yu Zhao1,2†Rusen Zhang&#x;Rusen Zhang3†Yi Zhong,Yi Zhong1,2Jingjing WangJingjing Wang4Zuquan Weng,*Zuquan Weng2,5*Heng Luo,,*Heng Luo1,2,6*Cunrong Chen*Cunrong Chen4*
  • 1College of Computer and Data Science, Fuzhou University, Fuzhou, China
  • 2Centre for Big Data Research in Burns and Trauma, Fuzhou University, Fuzhou, China
  • 3Department of Cardiovascular Medicine, Affiliated Fuzhou First Hospital of Fujian Medical University, Fuzhou, China
  • 4Department of Critical Care Medicine, Union Hospital of Fujian Medical University, Fuzhou, China
  • 5College of Biological Science and Engineering, Fuzhou University, Fuzhou, China
  • 6MetaNovas Biotech Inc., Foster City, CA, United States

The Coronavirus Disease 2019 (COVID-19) has spread all over the world and impacted many people’s lives. The characteristics of COVID-19 and other types of pneumonia have both similarities and differences, which confused doctors initially to separate and understand them. Here we presented a retrospective analysis for both COVID-19 and other types of pneumonia by combining the COVID-19 clinical data, eICU and MIMIC-III databases. Machine learning models, including logistic regression, random forest, XGBoost and deep learning neural networks, were developed to predict the severity of COVID-19 infections as well as the mortality of pneumonia patients in intensive care units (ICU). Statistical analysis and feature interpretation, including the analysis of two-level attention mechanisms on both temporal and non-temporal features, were utilized to understand the associations between different clinical variables and disease outcomes. For the COVID-19 data, the XGBoost model obtained the best performance on the test set (AUROC = 1.000 and AUPRC = 0.833). On the MIMIC-III and eICU pneumonia datasets, our deep learning model (Bi-LSTM_Attn) was able to identify clinical variables associated with death of pneumonia patients (AUROC = 0.924 and AUPRC = 0.802 for 24-hour observation window and 12-hour prediction window). The results highlighted clinical indicators, such as the lymphocyte counts, that may help the doctors to predict the disease progression and outcomes for both COVID-19 and other types of pneumonia.

1 Highlights

1. In this study, multiple approaches including statistical analysis, machine learning and deep learning were utilized to understand the relationship between clinical variables and disease outcomes of COVID-19 and pneumonia.

2. The best models obtained good performance on the test sets when predicting severity of COVID-19 and the mortality of pneumonia patients in ICUs with AUROC larger than 0.9 and AUPRC greater than 0.8.

3. Feature importance and two-level attention mechanisms were utilized to interpret the models and highlight clinical variables associated with disease outcomes. Lymphocyte counts were found to be an important biomarker shared by COVID-19 and other types of pneumonia.

2 Introduction

In December 2019, an outbreak of highly infectious respiratory disease was reported in Wuhan, China, which was later identified as Coronavirus Disease 2019 (COVID-19) (Ward et al., 2020; Yamasaki et al., 2020). Ever since, COVID-19 has spread all over the world and impacted many people’s lives (Velavan and Meyer, 2020). In order to defeat COVID-19 pandemic, researchers developed machine learning models to help predict the outbreak so that the medical system can get ready and allocate resources ahead of time (Ardabili et al., 2020a; Ardabili et al., 2020b; Niazkar and Niazkar, 2020). Additionally, machine learning has also been utilized to help with the diagnosis and understand the clinical indicators of COVID-19 (Hussain et al., 2020; Li et al., 2020; Podder and Mondal, 2020; Alballa and Al-Turaiki, 2021; Podder et al., 2021). One of the major symptoms caused by COVID-19 is pneumonia, which is also a common infectious diseases threatening both human health and medical resources (Zimlichman et al., 2013). Patients with mild pneumonia only need outpatient treatment, while severe patients have to be hospitalized or even admitted to the intensive care unit (ICU) for rescue (Wei, 2020). Timely treatment for severe patients is significantly associated with reduced mortality (Attaway et al., 2021). Therefore, it is important to have early assess of the disease stage and make effective treatment plans (Watkins, 2020; González-Nóvoa et al., 2021).

In the early days, to classify the severity of respiratory and infectious diseases, the disease research societies have proposed different scoring standards (Liapikou et al., 2009). Among them, CURB-65, a combination of five endpoints, confusion, uremia, respiratory rate, blood pressure and age (65 years as a cutoff), is widely assessed in patients with community-acquired pneumonia (CAP) (Jones et al., 2011; Satici et al., 2020). For the COVID-19 patients, researchers also used CURB-65 to predict the risk of death (Nguyen et al., 2020). In addition to CURB-65, Zhang et al. (2019) developed a classification and regression tree (CART) scoring system by machine learning based on demographic data and clinical characteristics to evaluate the survival of CAP patients in ICUs. The results of Zhang et al. (2019) showed that the CART obtained better evaluation performance compared to CURB-65. Goodwin and Demner-Fushman, (2019) further proposed a deep learning system, PRONTO, to predict the risk of patient death 1-2 days ahead. Based on extracorporeal membranous oxygenation (ECMO), Zhou et al. (2020) designed the pneumonia ECMO-eligible risk (PEER) system, which can predict the risk of death in a nomogram for both low-risk and high-risk pneumonia patients in the future. Yan et al. (2020) developed a machine learning model to propose a clinical decision rule composed of three biomarkers for COVID-19 patients, which can predict patient mortality 10 days in advance. Although many models have been developed to predict the disease outcomes, fewer studies focused on feature interpretability or clinical variable analysis across disease types. Even though some researchers compared differences between COVID-19 and regular pneumonia (Zhao et al., 2020), due to heterogenicity of different data sources, more analysis is helpful for a deeper understanding and better treatment approaches. Therefore, we compared and analyzed the clinical features of patients with COVID-19 and other types of pneumonia and assessed their performance of predicting disease outcomes.

In this study, a combination of machine learning models, including XGBoost and deep learning, was developed to predict the severity of COVID-19 and the patient mortality of other types of pneumonia in different time windows. Feature importance of both temporal and discrete clinical variables was analyzed to understand the association factors for different disease outcomes. Our goal is to provide early warnings for patients with dire outcomes so that doctors could have time to come up with appropriate monitoring and intervention procedures to prevent a worse situation.

3 Materials and Methods

3.1 The Overall Framework

The overall workflow of this study is summarized in Figure 1, which includes the following steps: (1) data collection of COVID-19 and other pneumonia data, (2) statistical analysis and development, training and evaluation of machine learning and deep learning models, and (3) result analysis including model evaluation and feature interpretation. The details were introduced in the following sections.

FIGURE 1
www.frontiersin.org

Figure 1 The overall workflow of this study. EHR, electronic health record; LR, logistic regression; RF, random forest; SBP, systolic blood pressure; WBCs, white blood cell counts.

3.2 Collection and Analysis of COVID-19 Clinical Data

The demographic and clinical data of 80 COVID-19 patients were collected during the diagnosis and treatment of patients in Hubei province from February 1, 2020 to March 30, 2020, assisted by a medical team from Fujian Medical University Union Hospital. The study followed the principles of good clinical practice (GCP) and relevant national legal requirements with consent and approval from the patients’ families. It was also approved by the Ethics Committee of Fujian Medical University Union Hospital (No. 2021KY052).

Patients from this dataset were divided into severe COVID-19 (SC) and non-severe COVID-19 (non-SC) groups. Patients without any severity information were excluded. The missing clinical measurements were filled with average values.

To analyze the data, mean (along with standard deviation) or median (along with interquartile range (IQR)) values were calculated as numeric measurements. The Kolmogorov-Smirnov normality test was used to examine if the sample data were normally distributed. Additional statistical analysis methods including T-test, Kruskal Wallis rank test and Fisher’s exact test were utilized to summarize and analyze the clinical features of COVID-19 patients (Zhang et al., 2018). A p value less than 0.05 was considered statistically significant.

3.3 Collection and Processing of MIMIC-III and eICU Data

The Medical Information Mart for Intensive Care (MIMIC-III) is a database that collected de-identified health records of over 50,000 patients who stayed in critical care units of Beth Israel Deaconess Medical Center between 2001 and 2012 (Johnson et al., 2016). Another database, eICU, is a combined clinical database from many ICUs in the United States. It contains clinical data of more than 200,000 hospitalized patients from 2014 to 2015 (Pollard et al., 2019). Since both the MIMIC-III and eICU data were harvested much earlier than the COVID-19 pandemic, it is safe to assume there were no COVID-19 patients in the records. The clinical data included demographics, vital signs, laboratory test results, diagnoses and other information of the patients. In this study, more than 3,000 patients with pneumonia by the International Classification of Diseases, Ninth Revision (ICD-9) codes of 481, 482 and 486 (Table 1) were identified in both MIMIC-III and eICU. Multiple filters were used to filter the patients (Figure 2). Adolescents and patients who did not stay in an ICU were excluded, and only patients with lengths of ICU stay no less than 30 hours were extracted.

TABLE 1
www.frontiersin.org

Table 1 ICD-9 codes used in this study for pneumonia.

FIGURE 2
www.frontiersin.org

Figure 2 The flowchart of pneumonia patient selection from the MIMIC-III database.

The pneumonia patients from MIMIC-III and eICU were grouped into two groups, survival or death, based on the records of ICU admission, discharge and death information. We referred Ng et al. (2016) for data extraction. Given an outcome of ICU discharge (survival) or death, the period right before the outcome was considered as the prediction window (6, 12, 24, 36, 48, 60 or 72 hours), and the time before the prediction window was used as the observation window (24, 48 or 72 hours). The clinical records in the observation window were used as input to develop models to predict outcomes (Figure 3). The abnormal values in the observation windows were deleted according to the criteria of Harutyunyan et al. (2019).

FIGURE 3
www.frontiersin.org

Figure 3 Data extraction framework for pneumonia patients in MIMIC-III. DBP, diastolic blood pressure; HR, heart rate; RR, respiratory rate; SBP, systolic blood pressure.

3.4 Model Development

3.4.1 Traditional Machine Learning Models

We used a series of machine learning models in the following sections to predict disease outcomes based on the clinical features for both COVID-19 and other types of pneumonia.

3.4.1.1 Logistic Regression (LR)

Logistic regression (Bouyer, 1991) is a generalized linear model used to solve classification problems. It is linear regression with a layer of the sigmoid function to map features towards labels. By introducing regularization, it can reduce the influence of multicollinearity. Given n data samples with m features as xj(j∈ (1,2, … ,m)), the predicted probability yi^ of sample i can be calculated as follows:

y^i=11+eθTx(i)
θTX(i)=θ0+j=1mθjxj(i)

where θ represents coefficients.

3.4.1.2 Random Forest (RF)

Random forest (Biau and Scornet, 2016) is composed of multiple decision trees (classified regression trees, or CARTs). The model randomly drafts N training subsets M = {M1, M2,…,Mn} based on bootstraps. The probability P of each sample not being drawn is calculated as follows:

P=(11N)N

The N decision trees T = {T1, T2,…,Tn} are developed based on their corresponding training subsets. The Gini index is used for CARTs to select tree nodes. Gini index is calculated as follows:

G(M)=k=1Kpk(1pk)=1k=1Kpk2

Here M is the independent training subset, and pk indicates the probability that the sample belongs to the k-th category.

3.4.1.3 XGBoost

XGBoost (Chen et al., 2015) is a gradient lifting tree model composed of multiple CARTs. The algorithm divides samples into two branches according to the feature thresholds. After multiple groupings, the end of each CART (leaf node) contains samples with the same label. Given n samples with m features as xi(i∈ (1,2, … , m)), the predicted probability y^i of sample i can be calculated as follows:

y^i=K=1Kfk(xi),fkF

where fk is the prediction score of a single decision tree, and ϝ is the space of all trees.

In order to get the optimal solution, the following loss function with regularization is optimized:

(ϕ)=il(y^i,yi)+kΩ(fk)
Ω(f)=γT+12λωi2

Where ℒ is the loss function, Ω is the penalty term, T is the number of leaves, ωi is the score of leaf node i, and γ and λ are the coefficient parameters.

3.4.1.4 Data Processing

For the COVID-19 data, the latest clinical features measured during the hospitalization of the COVID-19 patients were used to train the models to classify the outcomes of severe COVID-19 (SC) and non-severe COVID-19 (non-SC).

For the MIMIC-III and eICU data, the distribution parameters (maximum, minimum, median, mean and standard deviation) of vital signs and laboratory measurements in the observation window, as well as the drug prescriptions were prepared as features. For the oxygenation indexes (SpO2 and SaO2), the maximum values were not used since most of them were close to 100%. The models were trained to predict the outcome of survival for each patient.

For logistic regression, both L1 and L2 regularizations were applied. The top 10 most important features in XGBoost were analyzed to understand the clinical indicators for disease outcomes.

3.4.2 Deep Learning Models

3.4.2.1 Model Design

Since the MIMIC-III and eICU data contain temporal clinical features, two time-series deep learning models, including bi-directional long short-term memory (LSTM) models without and with attention (namely Bi-LSTM and Bi-LSTM_Attn, respectively), were developed to predict survival of the pneumonia patients. However, due to the limited samples and discrete features of COVID-19 data, we were not able to apply these models and had to rely on the traditional models.

For MIMIC-III and eICU data, the clinical features include vital signs, laboratory measurements and drug usages (Supplementary Table S1). The clinical data of one-hour intervals from the observation windows were extracted from both databases. If there were multiple measurements within a given hour, their mean values were taken as features. If no data within the one-hour internal were observed, the mean values of the entire temporal series were used.

The architecture of the deep learning model (Figure 4) was developed based on Zhang et al. (Kaji et al., 2019) and the RETAIN (Choi et al., 2016) model. Taking the bi-directional long short-term memory (LSTM) model with attention (Bi-LSTM_Attn) as an example, the hourly data of vital signs and laboratory measurements in the observation windows were used as the inputs. Self-attention mechanisms were utilized to calculate the attention weights of both time-series clinical data and the discrete drug prescription features. Then both sides were concatenated by the context vector and a final sigmoid function was used to make a classification prediction. The model included three parts: 1) the vital sign and laboratory measurement (VL) encoder, 2) the drug (DG) encoder and 3) the classification module, which were introduced in the sections below.

FIGURE 4
www.frontiersin.org

Figure 4 The deep learning model architecture of bi-directional long short-term memory (LSTM) with attention (Bi-LSTM_Attn). It consists of three parts: vital sign and laboratory measurement (VL) encoder, drug (DG) encoder and the classification module.

The LSTM model without the attention layers (Bi-LSTM) had a similar but simpler architecture. Since no attention module was presented, the outputs from the head and tail LSTM modules were merged with the drug feature vectors and sent to the linear layers with the sigmoid function for classification.

3.4.2.2 The VL Encoder

To obtain the embeddings of vital signs and laboratory measurements, we used the bidirectional LSTM (Bi-LSTM) networks, which combined both the forward and backward information within a sequence. Given a temporal feature vector at a specific time (t), Xit(i{1,2,3,,N},t{1,2,3,,T}) of patient i, the embedding of the input vector Xi(t) was calculated by the Bi-LSTM model as follows:

hi(t)=BiLSTM(hi(t1),Xi(t))

The bidirectional output vector hi(t) was defined as hi(t)=[hi(t),hi(t)]. Since hi(t) contained the sequential information around Xi(t) , in order to let the model know the importance of each time point in the observation window, the attention weight of each moment was calculated by self-attention mechanism. The temporal context vector cα was calculated as follows:

αi(t)=Softmax(WαThi(t)+bα)
cα=i=1Tαi(t)hi(t)

Where WαT and bα are learnable parameters and αi(t) is the attention weight of each time point.

3.4.2.3 The DG Encoder

Given a list of drug prescriptions for each patient as Di= [d0, d1, d2, …,dM] (di∈{0, 1}, i ∈{1, 2, 3, … , M}), multi-layer artificial neural networks were used to get the embedding vector of the prescription vi as vi = [v0,v1,v2,…,vM] follows:

Vi=MLP(Di)

Here, the Rectified Linear Unit (ReLU) was utilized as the activation function. To help the model understand the importance of the drug prescription features, MLP layers were used to calculate attention weight βi for the embedding feature vector Vi of each drug prescription. Then a sum of dot product of all the embedding vectors Vi and attention vectors βi were calculated as follows:

βi=Softmax(WβTVi+bβ)
cβ=i=1TβiVi

Where WβT and bβ are learnable parameters, ⊙ is the element-wise multiplication and cβ is the drug context vector.

Because the sample data we selected is unbalanced, so we use the class-balanced (CB) loss (Cui et al., 2019) to calculate the classification loss as follows:

Loss(p,y)=1NiN(yilogy^i+(1yi)log(1y^i))
Bloss(p,y)=1EnyLoss(p,y)=1β1βnyLoss(p,y)

where N is total sample, β ∈[0,1) is a hyperparameter and ny is the number of y labels in the sample.

3.4.3 Model Prediction and Evaluation

Since the COVID-19 dataset is a relatively small dataset, the COVID-19 patients were randomly divided into training (including validation) and test sets using 4:1 ratio in order to have enough patients in the test set. For the MIMIC-III and eICU data, pneumonia patients were randomly divided using 8:1:1 ratio into the training, validation and test sets. The training and validation sets were combined for 5-fold cross-validations during the development of traditional machine learning models. For deep learning models, the training set was used to train the models; the validation set was utilized to learn and optimize the hyperparameters and the test set was used to evaluate the models. The deep learning models contained a single-layer Bi-LSTM with 256 hidden layers to generate embeddings for temporal features, and two MLP layers (128 and 103 neurons) to generate embeddings for the drug prescription features. Both lists of embeddings were merged and sent to a linear layer of 128 neurons with a sigmoid activation function as the classification module. The dropout rate of 0.5 was used before the linear layer. During training, the Adam optimizer was used with a learning rate of 0.001, the batch size was set to 32, and the models were trained by 150 epochs. The area under the receiver operating characteristic curve (AUROC) and area under the precision-recall curve (AUPRC) were calculated to evaluate the model performance. Additional details regarding the evaluation metrics were attached in Supplementary Materials.

4 Results

4.1 Statistical Analysis of the COVID-19 Data

The statistical analysis for the COVID-19 data was shown in Table 2. Of the 80 COVID-19 patients, male accounted for 47.50%. The average age was 64.21 ± 20 years with Kolmogorov-Smirnov normality test p > 0.05, indicating a normal distribution. Cough (71.25%) and fever (63.75%) were the most frequent symptoms, while a few patients had other symptoms, such as shortness of breath (17.50%), diarrhea (5.00%) and chest pain (3.75%). Some patients had other diseases, including diabetes (23.75%), heart disease (5.00%), respiratory system disease (5.00%), other symptoms (28.75%). The median of eastern cooperative oncology group (ECOG) score and acute physiology and chronic health evaluation II (APACHE II) score were 1.00 [IQR, 1.00-2.00] and 5.00 [IQR, 3.00-7.25]. Among the drugs used in the treatment of COVID-19, arbidol was the most common one (83.75%), followed by oseltamivir (20.00%) and ribavirin (4.00%). For blood test results, the medians of albumin, lymphocyte counts and prothrombin time (PT) were 36.80 g/L [IQR, 33.30-40.65], 1.39×109/L [IQR, 0.97-1.72] and 13.40 seconds (s) [IQR, 13.00-14.25], respectively.

TABLE 2
www.frontiersin.org

Table 2 Statistics of the clinical variables for COVID-19 data.

Of the 80 patients with COVID-19, 13 were severe COVID-19 (SC) patients and 67 were non-severe COVID-19 (non-SC) patients. The severe patients had an average age of 73.15 ± 13.64 years, while non-severe patients had a younger average age of 62.48 ± 12.48 years (t-test p < 0.05). Here, age may be a factor that is associated with the severeness of COVID-19. In addition, the prognosis scores, including ECOG scores and APACHE II scores of severe patients, were higher than those of the non-severe patients both with Kruskal Wallis rank test p < 0.001, as larger prognosis scores indicate more severe conditions. For blood test results, albumin and lymphocyte counts were lower in severe patients than non-severe patients, while the prothrombin time (PT) was longer in severe COVID-19 patients than non-severe COVID-19 patients (all Kruskal Wallis rank p values < 0.001).

4.2 Model Performance and Interpretation for the COVID-19 Data

Figure 5 showed the distribution of AUROC and AUPRC values of different models after 100 rounds of 5-fold cross-validations on the COVID-19 data. The hyperparameters of the models were attached in Supplementary Table S2. All AUROC values of logistic regression with L1 regularization (LR(L1)), logistic regression with L2 regularization (LR(L2)) and XGBoost models in predicting the severity of patients with COVID-19 were greater than 0.5, and the AUPRC values of LR(L1), LR(L2) and XGBoost were all greater than 0.25. While random forest (RF) had slightly worse performance, LR(L1) (AUROC = 0.986 ± 0.033 and AUPRC = 0.948 ± 0.121), LR(L2) (AUROC = 0.983 ± 0.032 and AUPRC = 0.931 ± 0.135) and XGBoost (AUROC = 0.979 ± 0.048 and AUPRC = 0.927 ± 0.154) were more effective to predict the severity of COVID-19 patients.

FIGURE 5
www.frontiersin.org

Figure 5 Performance comparison of multiple models during 100 rounds of 5-fold cross-validations on the COVID-19 data. (A) AUROC and (B) AUPRC curve showing the performance of the classifier using different model using cross-validation on the COVID-19 data.

For the XGBoost model, we calculated the contributions of each clinical variables towards the severity of COVID-19 by information gain, and showed the top 10 most important features in Figure 6. Results indicate that the APACHE II has the largest importance, followed by neutrophil counts, lymphocyte counts and AST (aspartate aminotransferase). In order to get an interpretable decision tree, we set the parameter n_estimators to 1 in XGBoost and fed the top 10 features into XGBoost one by one. Finally, we got the best decision tree from the APACHE II feature (AUROC = 1.00, AUPRC = 0.833) to predict the severity of COVID-19 patients (Supplementary Figure S1). Additionally, the top 10 most important features from the LR(L1), LR(L2) and RF models were included in Supplementary Table S3, which were found similar to the top features from the XGBoost model.

FIGURE 6
www.frontiersin.org

Figure 6 Top 10 most important features from the XGBoost model on the COVID-19 data. APACHE II, acute physiology and chronic health evaluation I; AST, aspartate aminotransferase; BUN, blood urea nitrogen; PT, prothrombin time; WBC, white blood cell count; TT, thrombin time.

4.3 Model Performance and Interpretation of the MIMIC-III and eICU Data

For MIMIC-III data, the performance of models was dependent on the selection of the observation and prediction windows (Figure 7). From the data point of view, shorter prediction window or observation window corresponded to more patients in the test set, especially the survival patients. For model performance, when the observation window was fixed, the performance of all models generally decreased as the length of the prediction window increased. This makes sense since longer prediction window presents more challenges for the models to predict the future.

FIGURE 7
www.frontiersin.org

Figure 7 Performance of different models to predict pneumonia patient survival using different prediction and observation windows based on the test set of MIMIC-III. The left column showed AUROC values and the right column showed AUPRC values. (A, B) are AUROC and AUPRC of different models under different 24-hour observation windows and different prediction windows, respectively. (C, D) are AUROC and AUPRC of different models under different 48-hour observation windows and different prediction windows, respectively. (E, F) are AUROC and AUPRC of different models under different 72-hour observation windows and different prediction windows, respectively.

The hyperparameters of LR(L1) and XGBoost models were attached in Supplementary Table S4 and 10-fold and 15-fold cross-validation results of these models were included in Supplementary Table S5. Comparing across different models, the AUROC values of the traditional machine learning models (LR(L1) and XGBoost) were higher than the deep learning models when both the observation and prediction windows were within 48 hours. However, when both the observation and prediction window were large (no less than 60 hours), the Bi-LSTM_Attn model outperformed others in terms of AUROC [Figure 7(e)]. Since the deep learning models were designed to learn and infer based on temporal features, it was expected to obtain a better performance when time windows were longer.

Like before, we calculated the feature importance from the XGBoost model via information gain and showed the top 10 most important features in Figure 8. Among the 10 features, norepinephrine, vasopressin, morphine and midazolam were known drugs for the treatment of patients with pneumonia. Additionally, we found that the minimum respiratory rates (RR) value, the minimum pH, the average platelet count, the maximum and average blood urea nitrogen (BUN), and the lymphocytes minimum were important clinical features affecting the patient survival with pneumonia in ICUs.

FIGURE 8
www.frontiersin.org

Figure 8 Top 10 most important features from the MIMIC-III data based on analysis of the XGBoost model. BUN, blood urea nitrogen; RR, respiratory rate.

For eICU data, the performance evaluation of the deep learning models (Bi-LSTM and Bi-LSTM_Attn) under different time windows was shown in Supplementary Figure S2. We found that for the eICU data, the best time window is consistent with MIMIC-III data for predicting the mortality of pneumonia patients during ICU stay (AUROC = 0.856 and AUPRC = 0.566), which is a 24-hour observation window and a 12-hour prediction window.

4.4 Visualization of the Attention Mechanism in the Deep Learning Model for Predicting MIMIC-III Pneumonia

Since the deep learning model Bi-LSTM_Attn included two attention layers, we used an example patient to visualize the attention weights of temporal clinical features from the best predictive windows (24-hour observation window and 12-hour prediction window) in Figure 9. For non-temporal drug prescription features, the top 10 most important features were ranked by average attention weights and shown in Supplementary Table S6.

FIGURE 9
www.frontiersin.org

Figure 9 Clinical features and attention weights of one patient in the 24-hour observation window who passed away after 12 hours (prediction window, not shown in figure) from MIMIC-III. BMP, beats per minute; DBP, diastolic blood pressure; HR, heart rate; MBP, mean blood pressure; RR, respiratory rate; SBP, systolic blood pressure.

Figure 9 showed the average vital signs and attention weights of the patient who suffered from severe pneumonia in MIMIC-III data, with complications of acute respiratory failure, congestive heart failure, acute renal failure and other diseases, and eventually died in the ICU (ICUSTAY_ID = 268985). We observed that during the first 23 hours of the observation window, the patient’s heart rate was too high compared to the normal range (60-100 BMP), possibly due to congestive heart failure. Since the pneumonia patient suffered acute respiratory failure, the respiratory rate (RR) increased beyond the normal range (12-20 BMP), especially in the later part of the observation window. Though the vital signs fluctuated and changed over time, the attention weights showed in the bottom row of Figure 9 increased to a relatively larger value after the 17th hour. It reached its maximum at the 23rd hour, when all the vital signs of the patients showed a sharp drop afterwards. Though it was still 12 hours before the patient passed away, the model already picked up some unusual signs at this time point.

5 Discussion

In this study, we collected data of both COVID-19 and other types of pneumonia and analyzed the clinical variables that may affect the severity of COVID-19 and other types of pneumonia. Though the quick detection technology of COVID-19 has been widely adopted for diagnosis, it is still useful to do a retrospective analysis on these clinical variables to better understand the disease characteristics for better prevention towards fatal outcomes.

APACHE II score system is the most commonly used comprehensive index to assess the severity of pneumonia patients in ICUs. A higher APACHE II score indicates a worse condition and a greater risk of death for the patient. Wang et al. (2020) determined that the median of APACHE II in severe patients with COVID-19 was 17 (IQR: 10-22) and its lower quartile was slightly greater than the critical value of APACHE II in Supplementary Figure S1, which indicated that APACHE II < 8.5 identified by our model could be a useful cutoff to classify the severity of patients with COVID-19. Studies (Narasaraju et al., 2011; Liu et al., 2016; Lefrançais et al., 2018) have shown that excessive activation of neutrophils can release neutrophil extracellular traps (NETs), and the excessive formation of NETs can lead to a series of inflammatory reactions. These reactions may cause permanent damage to the lungs as well as the cardiovascular and renal systems in COVID-19 patients. In our results, lymphocytes are the third important factor affecting the severity of patients with COVID-19. During the clinical observation, the laboratory examination results showed that the number of lymphocytes in severe patients with COVID-19 decreased significantly, and the numbers of CD3, CD8, and CD4 T cells also continued to decline, which led to the aggravation of inflammatory reactions, increased cytokine levels and seriously damaged lung function (Chen et al., 2020; Tatum et al., 2020). Neutrophil-to-lymphocyte ratio (NLR) is an effective and convenient inflammatory marker for predicting systemic inflammation (de Jager et al., 2012; Cataudella et al., 2017; Tatum et al., 2020). However, NLR can also be used to determine the severity of the disease in patients with COVID-19. Liu et al. (2020) found that the NLRs of severe patients with COVID-19 were higher than non-severe patients. Severe patients with COVID-19 had symptoms of decreased lymphopenia and increased neutropenia, which caused the increase of proinflammatory cytokine levels in patients. The cytokine storm may make the immune system to lose control and damage the host cells, leading to multiple organ failure and patient fatality (Bermejo-Martin et al., 2020; Tan et al., 2020). Therefore, clinicians can identify high-risk COVID-19 patients using NLR. In addition to the clinical indicators discussed above, other features highlighted in Figure 6 were also found to have associations with COVID-19. Higuera-De-La-Tijera et al. (2021) showed the aspartate aminotransferase (AST) levels of severe COVID-19 patients were significantly higher than non-severe COVID-19 patients since AST is related to liver function. Additionally, Bastug et al. (2020) identified total protein and AST as important clinical indicators of COVID-19 patients admitted to ICUs. Ou et al. (2020) found that the white blood cell counts (WBCs) and procalcitonin in severe COVID-19 patients increased significantly, indicating their clinical predictability of COVID-19 severity. The increase of blood urea nitrogen (BUN) can lead to renal function injury (STARK, 1980), which may deteriorate the severity of COVID-19 (Kang et al., 2021). Lippi and Plebani (2020) reported the procalcitonin levels of severe COVID-19 patients were twice as much as non-severe COVID-19 patients, and Oblokulov et al. (2021) found this indicator helpful to come up with the treatment plan. Zhang et al. (2020) showed that the prothrombin time (PT) and thrombin time (TT) of severe COVID-19 patients were significantly higher than those of non-severe patients, as the coagulation function were affected in most severe COVID-19 patients.

Among the top 10 most important drugs for pneumonia, ipratropium bromide, aspirin and linezolid have treatment effects. Ipratropium bromide is a safe and effective drug in the treatment of chronic obstructive pulmonary disease (COPD) (Taylor et al., 2001). When ipratropium bromide was used in the treatment of acute exacerbation of COPD, blood oxygenation did not decrease initially (Ogale et al., 2010), which helped to improve lung function and reduce the risk of pneumonia. Aspirin has antiplatelet and anti-inflammatory effects, which can reduce the risk of cardiovascular disease (Gasparyan et al., 2008). A large primary care institution in the UK has shown that aspirin can prevent pneumonia patients from cardiovascular disease, thus improving the survival rate (Hamilton et al., 2021). Meanwhile, aspirin and macrolide drugs have a synergistic effect. For ICU patients with severe pneumonia, taking aspirin and macrolide drugs at the same time can produce different complementary effects and moderate the inflammatory reaction caused by pneumonia (Falcone et al., 2019). The clinical and microbial experiments of Wunderink et al. (2012) showed that linezolid was more effective than vancomycin in the treatment of pneumonia. On the contrary, fluticasone propionate can aggravate the condition of pneumonia. Fluticasone propionate is a kind of inhaled corticosteroids (ICS) for the treatment of chronic obstructive pulmonary disease (COPD). However, it may cause bacterial infection of the respiratory system due to poor dissolution and fluid lining in the airway and increase the risk of pneumonia (Janson et al., 2017).

COVID-19 is an acute respiratory infectious disease caused by viral infection (Fauci et al., 2020), while a regular pneumonia is a respiratory infectious disease caused by bacterial or viral infection (Bartlett and Mundy, 1995). They may have different pathogens and treatment strategies. Despite the differences, there are similarities among the top contributing clinical features of the two. It is observed that lymphocytes are the third most important feature to predict COVID-19 (Figure 6) as well as the top 10th most important features for pneumonia (Figure 8). It was reported that both COVID-19 and regular pneumonia patients may experience inflammatory reactions caused by infection and result in lymphopenia (Yamasaki et al., 2020; Zhao et al., 2020). In terms of treatment, COVID-19 is mainly treated with antiviral symptomatic support approaches (Kim et al., 2020), while other types of pneumonia were usually treated with sensitive antibiotics according to the pathogen (Lutfiyya et al., 2006). A few drugs may be useful to treat both conditions. Aspirin, for example, can reduce inflammation of COVID-19 and other pneumonia by inhibiting the synthesis of prostaglandins or other substances that may cause inflammatory reactions (Hamilton et al., 2021). Among the top 10 important features of both datasets shown in Figure 6 and Figure 8, lymphocytes and blood urea nitrogen (BUN) are shared indicators to predict the patient survival in MIMIC-III as well as the severity of COVID-19. For the COVID-19 data, except APACHE II, the other nine are laboratory measurements, while for the MIMIC-III data, there are six laboratory measurements and four drug treatment features (norepinephrine, vasopressin, morphine and midazolam). While the laboratory measurements are very important clinical features to predict disease outcomes in both cases, compared to COVID-19 data, the drug treatment features for MIMIC-III pneumonia patients gained more importance and may be effective in alleviating the condition of pneumonia. Since COVID-19 is a relatively new disease and drug treatments are still in active development, and the severity levels of COVID-19 and pneumonia data are different, it makes sense that more drug features were found in MIMIC-III pneumonia data. The understanding of both the similarities and differences of COVID-19 and other pneumonia may help us to monitor the critical clinical variables, make effective treatment plans and minimize fatalities. In the future work, we hope to extend the models to other disease areas and aid doctors to better understand diseases and make timely and effective clinical decisions.

The current research has a few limitations. The COVID-19 data from Wuhan, Hubei Province has a small sample size, and they only represent a very local patient pool and may not be applicable for patient groups in other regions or other COVID-19 variants. Similar limitation may also exist in the MIMIC-III and eICU database as their data are mainly from the United States. Thus, the conclusions of this study need to be tested in other large databases, cohorts and broader regions, to better identify the similarity and differences between COVID-19 and other types of pneumonia and the clinical indicators associated with them.

6 Conclusion

In this study, we collected pneumonia data from three sources, the COVID-19 clinical records, eICU and MIMIC-III. Statistical analysis was carried out to analyze the distributions and differences of clinical variables based on disease severeness. We developed multiple machine learning models to predict the severity of COVID-19 patients and the mortality of patients with other types of pneumonia in ICUs. The models were evaluated and compared by multiple metrics and different settings. We further analyzed the important features and interpreted the models. For the pneumonia patients, we found the same optimal time window to predict the death of ICU patients from both MIMIC-III and eICU (24-hour observation window and 12-hour prediction window). In addition, we further introduced deep learning models with two levels of attention mechanisms to visualize and interpret the temporal clinical variables via attention weights for pneumonia. It is observed while differences remain between COVID-19 and other pneumonia, lymphocytes are important biomarkers in both cases. This study provided a comprehensive analysis and comparison between COVID-19 and other types of pneumonia, which may have implications to aid clinicians for better disease understanding and clinical decision making.

Data Availability Statement

The data analyzed in this study is subject to the following licenses/restrictions: Only credentialed users who sign the specified DUA can access the files. Requests to access these datasets should be directed to MIMIC-III (https://physionet.org/content/mimiciii/1.4/) and eICU (https://physionet.org/content/eicu-crd/2.0/).

Ethics Statement

The studies involving human participants were reviewed and approved by Institutional Review Boards of Beth Israel Deaconess Medical Center (Boston, MA) and the Massachusetts Institute of Technology (Cambridge, MA), Ethics Committee of Fujian Medical University Union Hospital. The patients/participants provided their written informed consent to participate in this study. Written informed consent was obtained from the individual(s) for the publication of any potentially identifiable images or data included in this article.

Author Contributions

Concept and design, YiZ, JW, RZ and CC. Procedures, YuZ and CC. Writing of article, YuZ, RZ, ZW, and HL. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by National Natural Science Foundation of China (81971837), Natural Science Foundation of Fujian Province, China (2020J05109), and Funds of Scientific Research-Support Project, Fujian Provincial Department of Finance (2021XH018).

Conflict of Interest

Author HL was employed by the company MetaNovas Biotech Inc. at present. Related work in this study was done in Fuzhou University.

The remaining 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.

Acknowledgments

Thank the Centre for Big Data Research in Burns and Trauma for providing the research environment, resources and technology, and colleagues for their help and support.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fcimb.2022.838749/full#supplementary-material

References

Alballa, N., Al-Turaiki, I. (2021). Machine Learning Approaches in COVID-19 Diagnosis, Mortality, and Severity Risk Prediction: A Review. Inf. Med. Unlocked 24, 100564. doi: 10.1016/j.imu.2021.100564

CrossRef Full Text | Google Scholar

Ardabili, S., Mosavi, A., Band, S. S., Varkonyi-Koczy, A. R. (2020a). “Coronavirus Disease (COVID-19) Global Prediction Using Hybrid Artificial Intelligence Method of ANN Trained With Grey Wolf Optimizer”, In 2020 IEEE 3rd International Conference and Workshop in Óbuda on Electrical and Power Engineering (CANDO-EPE). (Budapest, Hungary: IEEE), 000251–000254.

Google Scholar

Ardabili, S. F., Mosavi, A., Ghamisi, P., Ferdinand, F., Varkonyi-Koczy, A. R., Reuter, U., et al. (2020b). Covid-19 Outbreak Prediction With Machine Learning. Algorithms 13 (10), 249. doi: 10.3390/a13100249

CrossRef Full Text | Google Scholar

Attaway, A. H., Scheraga, R. G., Bhimraj, A., Biehl, M., Hatipoğlu, U. (2021). Severe Covid-19 Pneumonia: Pathogenesis and Clinical Management. BMJ 372, n436. doi: 10.1136/bmj.n436

PubMed Abstract | CrossRef Full Text | Google Scholar

Bartlett, J. G., Mundy, L. M. (1995). Community-Acquired Pneumonia. N. Engl. J. Med. 333 (24), 1618–1624. doi: 10.1056/NEJM199512143332408

PubMed Abstract | CrossRef Full Text | Google Scholar

Bastug, A., Bodur, H., Erdogan, S., Gokcinar, D., Kazancioglu, S., Kosovali, B. D., et al. (2020). Clinical and Laboratory Features of COVID-19: Predictors of Severe Prognosis. Int. Immunopharmacol. 88, 106950. doi: 10.1016/j.intimp.2020.106950

PubMed Abstract | CrossRef Full Text | Google Scholar

Bermejo-Martin, J. F., Almansa, R., Menéndez, R., Mendez, R., Kelvin, D. J., Torres, A. (2020). Lymphopenic Community Acquired Pneumonia as Signature of Severe COVID-19 Infection. J. Infect 80 (5), e23. doi: 10.1016/j.jinf.2020.02.029

PubMed Abstract | CrossRef Full Text | Google Scholar

Biau, G., Scornet, E. (2016). A Random Forest Guided Tour. Test 25 (2), 197–227. doi: 10.1007/s11749-016-0481-7

CrossRef Full Text | Google Scholar

Bouyer, J. (1991). Logistics Regression in Epidemiology. II. Rev. d’epidemiologie Sante Publique 39 (2), 183–196.

Google Scholar

Cataudella, E., Giraffa, C. M., Di Marca, S., Pulvirenti, A., Alaimo, S., Pisano, M., et al. (2017). Neutrophil-To-Lymphocyte Ratio: An Emerging Marker Predicting Prognosis in Elderly Adults With Community-Acquired Pneumonia. J. Am. Geriatrics Soc. 65 (8), 1796–1801. doi: 10.1111/jgs.14894

CrossRef Full Text | Google Scholar

Chen, T., He, T., Benesty, M., Khotilovich, V., Tang, Y., Cho, H., et al. (2015). Xgboost: Extreme Gradient Boosting. R Package Version 0.4-2 1 (4), 1–4. doi: 10.1145/2939672.2939785

CrossRef Full Text | Google Scholar

Chen, G., Wu, D., Guo, W., Cao, Y., Huang, D., Wang, H., et al. (2020). Clinical and Immunological Features of Severe and Moderate Coronavirus Disease 2019. J. Clin. Invest. 130 (5), 2620–2629. doi: 10.1172/JCI137244

PubMed Abstract | CrossRef Full Text | Google Scholar

Choi, E., Bahadori, M. T., Kulas, J. A., Schuetz, A., Stewart, W. F., Sun, J. (2016). Retain: An Interpretable Predictive Model for Healthcare Using Reverse Time Attention Mechanism. arXiv preprint arXiv:1608.05745. 29, 3504–3512. doi: 10.48550/arXiv.1608.05745

CrossRef Full Text | Google Scholar

Cui, Y., Jia, M., Lin, T.-Y., Song, Y., Belongie, S. (2019). “Class-Balanced Loss Based on Effective Number of Samples”, In Proceedings of the IEEE/CVF Conference on Computer Vision and Pattern Recognition, (Long Beach, CA: IEEE) 9268–9277.

Google Scholar

de Jager, C. P., Wever, P. C., Gemen, E. F., Kusters, R., van Gageldonk-Lafeber, A. B., van der Poll, T., et al. (2012). The Neutrophil-Lymphocyte Count Ratio in Patients With Community-Acquired Pneumonia. PloS One 7 (10), e46561. doi: 10.1371/journal.pone.0046561

PubMed Abstract | CrossRef Full Text | Google Scholar

Falcone, M., Russo, A., Shindo, Y., Farcomeni, A., Pieralli, F., Cangemi, R., et al. (2019). A Hypothesis-Generating Study of the Combination of Aspirin Plus Macrolides in Patients With Severe Community-Acquired Pneumonia. Antimicrob. Agents Chemother 63 (2), e01556–e01518. doi: 10.1128/AAC.01556-18

PubMed Abstract | CrossRef Full Text | Google Scholar

Fauci, A. S., Lane, H. C., Redfield, R. R. (2020). Covid-19—Navigating the Uncharted”. Mass Med. Soc. 382, 1268–1269. doi: 10.1056/NEJMe2002387

CrossRef Full Text | Google Scholar

Gasparyan, A. Y., Watson, T., Lip, G. Y. (2008). The Role of Aspirin in Cardiovascular Prevention: Implications of Aspirin Resistance. J. Am. Coll. Cardiol. 51 (19), 1829–1843. doi: 10.1016/j.jacc.2007.11.080

PubMed Abstract | CrossRef Full Text | Google Scholar

González-Nóvoa, J. A., Busto, L., Rodríguez-Andina, J. J., Fariña, J., Segura, M., Gómez, V., et al. (2021). Using Explainable Machine Learning to Improve Intensive Care Unit Alarm Systems. Sensors 21 (21), 7125. doi: 10.3390/s21217125

PubMed Abstract | CrossRef Full Text | Google Scholar

Goodwin, T. R., Demner-Fushman, D. (2019). “Deep Learning From Incomplete Data: Detecting Imminent Risk of Hospital-Acquired Pneumonia in ICU Patients,” in AMIA Annual Symposium Proceedings: American Medical Informatics Association. (Washington: Springer) 467.

Google Scholar

Hamilton, F., Arnold, D., Henley, W., Payne, R. A. (2021). Aspirin Reduces Cardiovascular Events in Patients With Pneumonia: A Prior Event Rate Ratio Analysis in a Large Primary Care Database. Eur. Respir. J. 57 (2), 2002795. doi: 10.1183/13993003.02795-2020

PubMed Abstract | CrossRef Full Text | Google Scholar

Harutyunyan, H., Khachatrian, H., Kale, D. C., Ver Steeg, G., Galstyan, A. (2019). Multitask Learning and Benchmarking With Clinical Time Series Data. Sci. Data 6 (1), 1–18. doi: 10.1038/s41597-019-0103-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Higuera-De-La-Tijera, F., Servín-Caamaño, A., Reyes-Herrera, D., Flores-López, A., Robiou-Vivero, E. J., Martínez-Rivera, F., et al. (2021). The Age-AST-D Dimer (AAD) Regression Model Predicts Severe COVID-19 Disease. Dis. Markers 2021, 6658270. doi: 10.1155/2021/6658270

PubMed Abstract | CrossRef Full Text | Google Scholar

Hussain, L., Nguyen, T., Li, H., Abbasi, A. A., Lone, K. J., Zhao, Z., et al. (2020). Machine-Learning Classification of Texture Features of Portable Chest X-Ray Accurately Classifies COVID-19 Lung Infection. BioMed. Eng. Online 19 (1), 88. doi: 10.1186/s12938-020-00831-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Janson, C., Stratelis, G., Miller-Larsson, A., Harrison, T. W., Larsson, K. (2017). Scientific Rationale for the Possible Inhaled Corticosteroid Intraclass Difference in the Risk of Pneumonia in COPD. Int. J. Chronic Obstruct. Pulmon. Dis. 12, 3055. doi: 10.2147/COPD.S143656

CrossRef Full Text | Google Scholar

Johnson, A. E., Pollard, T. J., Shen, L., Li-Wei, H. L., Feng, M., Ghassemi, M., et al. (2016). MIMIC-III, a Freely Accessible Critical Care Database. Sci. Data 3 (1), 1–9. doi: 10.1038/sdata.2016.35

CrossRef Full Text | Google Scholar

Jones, B. E., Jones, J., Bewick, T., Lim, W. S., Aronsky, D., Brown, S. M., et al. (2011). CURB-65 Pneumonia Severity Assessment Adapted for Electronic Decision Support. Chest 140 (1), 156–163. doi: 10.1378/chest.10-1296

PubMed Abstract | CrossRef Full Text | Google Scholar

Kaji, D. A., Zech, J. R., Kim, J. S., Cho, S. K., Dangayach, N. S., Costa, A. B., et al. (2019). An Attention Based Deep Learning Model of Clinical Events in the Intensive Care Unit. PloS One 14 (2), e0211057. doi: 10.1371/journal.pone.0211057

PubMed Abstract | CrossRef Full Text | Google Scholar

Kang, J., Chen, T., Luo, H., Luo, Y., Du, G., Jiming-Yang, M. (2021). Machine Learning Predictive Model for Severe COVID-19. Infect Genet. Evol. 90, 104737. doi: 10.1016/j.meegid.2021.104737

PubMed Abstract | CrossRef Full Text | Google Scholar

Kim, D., Lee, J.-Y., Yang, J.-S., Kim, J. W., Kim, V. N., Chang, H. (2020). The Architecture of SARS-CoV-2 Transcriptome. Cell 181 (4), 914–921. e910. doi: 10.1016/j.cell.2020.04.011

PubMed Abstract | CrossRef Full Text | Google Scholar

Lefrançais, E., Mallavia, B., Zhuo, H., Calfee, C. S., Looney, M. R. (2018). Maladaptive Role of Neutrophil Extracellular Traps in Pathogen-Induced Lung Injury. JCI Insight 3 (3), e98178. doi: 10.1172/jci.insight.98178

CrossRef Full Text | Google Scholar

Liapikou, A., Ferrer, M., Polverino, E., Balasso, V., Esperatti, M., Piñer, R., et al. (2009). Severe Community-Acquired Pneumonia: Validation of the Infectious Diseases Society of America/American Thoracic Society Guidelines to Predict an Intensive Care Unit Admission. Clin. Infect. Dis. 48 (4), 377–385. doi: 10.1086/596307

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, W. T., Ma, J., Shende, N., Castaneda, G., Chakladar, J., Tsai, J. C., et al. (2020). Using Machine Learning of Clinical Data to Diagnose COVID-19: A Systematic Review and Meta-Analysis. BMC Med. Inform. Decis. Mak. 20 (1), 247. doi: 10.1186/s12911-020-01266-z

PubMed Abstract | CrossRef Full Text | Google Scholar

Lippi, G., Plebani, M. (2020). Laboratory Abnormalities in Patients With COVID-2019 Infection. Clin. Chem. Lab. Med. (CCLM) 58 (7), 1131–1134. doi: 10.1515/cclm-2020-0198

CrossRef Full Text | Google Scholar

Liu, Y., Du, X., Chen, J., Jin, Y., Peng, L., Wang, H. H., et al. (2020). Neutrophil-To-Lymphocyte Ratio as an Independent Risk Factor for Mortality in Hospitalized Patients With COVID-19. J. Infect 81 (1), e6–e12. doi: 10.1016/j.jinf.2020.04.002

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, S., Su, X., Pan, P., Zhang, L., Hu, Y., Tan, H., et al. (2016). Neutrophil Extracellular Traps Are Indirectly Triggered by Lipopolysaccharide and Contribute to Acute Lung Injury. Sci. Rep. 6 (1), 1–8. doi: 10.1038/srep37252

PubMed Abstract | CrossRef Full Text | Google Scholar

Lutfiyya, M. N., Henley, E., Chang, L. F., Reyburn, S. W. (2006). Diagnosis and Treatment of Community-Acquired Pneumonia. Am. Family physician 73 (3), 442–450. doi: 10.1016/j.annepidem.2003.08.001

CrossRef Full Text | Google Scholar

Narasaraju, T., Yang, E., Samy, R. P., Ng, H. H., Poh, W. P., Liew, A.-A., et al. (2011). Excessive Neutrophils and Neutrophil Extracellular Traps Contribute to Acute Lung Injury of Influenza Pneumonitis. Am. J. Pathol. 179 (1), 199–210. doi: 10.1016/j.ajpath.2011.03.013

PubMed Abstract | CrossRef Full Text | Google Scholar

Ng, K., Steinhubl, S. R., deFilippi, C., Dey, S., Stewart, W. F. (2016). Early Detection of Heart Failure Using Electronic Health Records: Practical Implications for Time Before Diagnosis, Data Diversity, Data Quantity, and Data Density. Circ: Cardiovasc. Qual. Outcomes 9 (6), 649–658. doi: 10.1161/CIRCOUTCOMES.116.002797

PubMed Abstract | CrossRef Full Text | Google Scholar

Nguyen, Y., Corre, F., Honsel, V., Curac, S., Zarrouk, V., Fantin, B., et al. (2020). Applicability of the CURB-65 Pneumonia Severity Score for Outpatient Treatment of COVID-19. J. Infect. 81 (3), e96–e98. doi: 10.1016/j.jinf.2020.05.049

PubMed Abstract | CrossRef Full Text | Google Scholar

Niazkar, H. R., Niazkar, M. (2020). Application of Artificial Neural Networks to Predict the COVID-19 Outbreak. Global Health Res. Policy 5 (1), 1–11. doi: 10.1186/s41256-020-00175-y

CrossRef Full Text | Google Scholar

Oblokulov, A. R., Husenova, Z. Z., Ergashev, M. M. (2021). Procalcitonin as an Indicator of Antibacterial Therapy in Covid-19. Ann. Romanian Soc. Cell Biol. 25, 5220–5224.

Google Scholar

Ogale, S. S., Lee, T. A., Au, D. H., Boudreau, D. M., Sullivan, S. D. (2010). Cardiovascular Events Associated With Ipratropium Bromide in COPD. Chest 137 (1), 13–19. doi: 10.1378/chest.08-2367

PubMed Abstract | CrossRef Full Text | Google Scholar

Ou, M., Zhu, J., Ji, P., Li, H., Zhong, Z., Li, B., et al. (2020). Risk Factors of Severe Cases With COVID-19: A Meta-Analysis. Epidemiol. Infect. 148, e175–e175. doi: 10.1017/S095026882000179X

PubMed Abstract | CrossRef Full Text | Google Scholar

Podder, P., Bharati, S., Mondal, M. R. H., Kose, U. (2021). Application of Machine Learning for the Diagnosis of COVID-19. Data Science for COVID-19: Academic Press. (Bandung: IEEE) 175–194.

Google Scholar

Podder, P., Mondal, M. R. H. (2020). “Machine Learning to Predict COVID-19 and ICU Requirement,” in 2020 11th International Conference on Electrical and Computer Engineering (ICECE). (Bandung, IEEE).

Google Scholar

Pollard, T., Johnson, A., Raffa, J., Celi, L. A., Badawi, O., Mark, R. (2019). ICU Collaborative Research Database (Version 2.0). PhysioNet. Published online. doi: 10.13026/C2WM1R

CrossRef Full Text | Google Scholar

Satici, C., Demirkol, M. A., Altunok, E. S., Gursoy, B., Alkan, M., Kamat, S., et al. (2020). Performance of Pneumonia Severity Index and CURB-65 in Predicting 30-Day Mortality in Patients With COVID-19. Int. J. Infect. Dis. 98, 84–89. doi: 10.1016/j.ijid.2020.06.038

PubMed Abstract | CrossRef Full Text | Google Scholar

Stark, J. L. (1980). BUN/creatinine: Your Keys to Kidney Function. Nursing2021 10 (5), 33–38. doi: 10.1097/00152193-198005000-00007

CrossRef Full Text | Google Scholar

Tan, L., Wang, Q., Zhang, D., Ding, J., Huang, Q., Tang, Y.-Q., et al. (2020). Lymphopenia Predicts Disease Severity of COVID-19: A Descriptive and Predictive Study. Signal Transduction Targeted Ther. 5 (1), 1–3. doi: 10.1038/s41392-020-0148-4

CrossRef Full Text | Google Scholar

Tatum, D., Taghavi, S., Houghton, A., Stover, J., Toraih, E., Duchesne, J. (2020). Neutrophil-To-Lymphocyte Ratio and Outcomes in Louisiana Covid-19 Patients. Shock (Augusta Ga.). 54, 652–658. doi: 10.1097/SHK.0000000000001585

CrossRef Full Text | Google Scholar

Taylor, J., Kotch, A., Rice, K., Ghafouri, M., Kurland, C. L., Fagan, N. M., et al. (2001). Ipratropium Bromide Hydrofluoroalkane Inhalation Aerosol is Safe and Effective in Patients With COPD. Chest 120 (4), 1253–1261. doi: 10.1378/chest.120.4.1253

PubMed Abstract | CrossRef Full Text | Google Scholar

Velavan, T. P., Meyer, C. G. (2020). The COVID-19 Epidemic. Trop. Med. Int. Health 25 (3), 278. doi: 10.1111/tmi.13383

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, D., Hu, B., Hu, C., Zhu, F., Liu, X., Zhang, J., et al. (2020). Clinical Characteristics of 138 Hospitalized Patients With 2019 Novel Coronavirus–Infected Pneumonia in Wuhan, China. Jama 323 (11), 1061–1069. doi: 10.1001/jama.2020.1585

PubMed Abstract | CrossRef Full Text | Google Scholar

Ward, M. P., Li, X., Tian, K. (2020). Novel Coronavirus 2019, an Emerging Public Health Emergency. Transboundary emerging Dis. 67 (2), 469. doi: 10.1111/tbed.13509

CrossRef Full Text | Google Scholar

Watkins, J. (2020). Preventing a Covid-19 Pandemic. BMJ (Clinical Research Ed.) 368, m810. doi: 10.1136/bmj.m810

PubMed Abstract | CrossRef Full Text | Google Scholar

Wei, P. F. (2020). Diagnosis and Treatment Protocol for Novel Coronavirus Pneumonia (Trial Version 7). Chinese Med. J. 133 (9), 1087–1095. doi: 10.1097/CM9.0000000000000819

CrossRef Full Text | Google Scholar

Wunderink, R. G., Niederman, M. S., Kollef, M. H., Shorr, A. F., Kunkel, M. J., Baruch, A., et al. (2012). Linezolid in Methicillin-Resistant Staphylococcus Aureus Nosocomial Pneumonia: A Randomized, Controlled Study. Clin. Infect. Dis. 54 (5), 621–629. doi: 10.1093/cid/cir895

PubMed Abstract | CrossRef Full Text | Google Scholar

Yamasaki, Y., Ooka, S., Tsuchida, T., Nakamura, Y., Hagiwara, Y., Naitou, Y., et al. (2020). The Peripheral Lymphocyte Count as a Predictor of Severe COVID-19 and the Effect of Treatment With Ciclesonide. Virus Res. 290, 198089. doi: 10.1016/j.virusres.2020.198089

PubMed Abstract | CrossRef Full Text | Google Scholar

Yan, L., Zhang, H.-T., Goncalves, J., Xiao, Y., Wang, M., Guo, Y., et al. (2020). An Interpretable Mortality Prediction Model for COVID-19 Patients. Nat. Mach. Intell. 2 (5), 283–288. doi: 10.1038/s42256-020-0180-7

CrossRef Full Text | Google Scholar

Zhang, Y., He, L., Chen, H., Lu, S., Xiong, Y., Liu, J., et al. (2020). Manifestations of Blood Coagulation and its Relation to Clinical Outcomes in Severe COVID-19 Patients: Retrospective Analysis. Int. J. Lab. Hematol. 42 (6), 766–772. doi: 10.1111/ijlh.13273

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, S., Zhang, K., Yu, Y., Tian, B., Cui, W., Zhang, G. (2019). A New Prediction Model for Assessing the Clinical Outcomes of ICU Patients With Community-Acquired Pneumonia: A Decision Tree Analysis. Ann. Med. 51 (1), 41–50. doi: 10.1080/07853890.2018.1518580

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, Z., Zhu, C., Mo, L., Hong, Y. (2018). Effectiveness of Sodium Bicarbonate Infusion on Mortality in Septic Patients With Metabolic Acidosis. Intensive Care Med. 44 (11), 1888–1895. doi: 10.1007/s00134-018-5379-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhao, D., Yao, F., Wang, L., Zheng, L., Gao, Y., Ye, J., et al. (2020). A Comparative Study on the Clinical Features of Coronavirus 2019 (COVID-19) Pneumonia With Other Pneumonias. Clin. Infect. Dis. 71 (15), 756–761. doi: 10.1093/cid/ciaa247

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhou, H., Cheng, C., Lipton, Z. C., Chen, G. H., Weiss, J. C. (2020). Predicting Mortality Risk in Viral and Unspecified Pneumonia to Assist Clinicians With COVID-19 ECMO Planning. arXiv preprint arXiv:2006.01898. doi: 10.48550/arXiv.2006.01898

CrossRef Full Text | Google Scholar

Zimlichman, E., Henderson, D., Tamir, O., Franz, C., Song, P., Yamin, C. K., et al. (2013). Health Care–Associated Infections: A Meta-Analysis of Costs and Financial Impact on the US Health Care System. JAMA Internal Med. 173 (22), 2039–2046. doi: 10.1001/jamainternmed.2013.9763

CrossRef Full Text | Google Scholar

Keywords: the coronavirus disease 2019, pneumonia, statistical analysis, machine learning, attention mechanism, clinical indicators

Citation: Zhao Y, Zhang R, Zhong Y, Wang J, Weng Z, Luo H and Chen C (2022) Statistical Analysis and Machine Learning Prediction of Disease Outcomes for COVID-19 and Pneumonia Patients. Front. Cell. Infect. Microbiol. 12:838749. doi: 10.3389/fcimb.2022.838749

Received: 18 December 2021; Accepted: 07 March 2022;
Published: 19 April 2022.

Edited by:

Christoph M. Friedrich, Dortmund University of Applied Sciences and Arts, Germany

Reviewed by:

Prajoy Podder, Bangladesh University of Engineering and Technology, Bangladesh
Sina Ardabili. University of Mohaghegh Ardabili, Iran

Copyright © 2022 Zhao, Zhang, Zhong, Wang, Weng, Luo and Chen. 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: Zuquan Weng, d2VuZ3pxQGZ6dS5lZHUuY24=; Heng Luo, aGVuZ2x1bzg4QGdtYWlsLmNvbQ==; Cunrong Chen, Y2hjcjc4OUAxMzkuY29t

These authors have contributed equally to this work

Disclaimer: 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.