Advancing Bloodstream Infection Prediction Using Explanable Artificial Intelligence Framework

1. INTRODUCTION

Bloodstream infections (BSIs) are a significant public health challenge, often leading to severe clinical conditions such as sepsis and septic shock, particularly when unrecognized or untreated. The rapid progression of these infections, coupled with their association with high morbidity, mortality, and healthcare costs, renders BSIs a critical challenge in clinical care[1, 2]. The accuracy of available clinical decision tools for BSI and sepsis is in need of improvement. Currently, most of them are mostly based on changes in vital signs and abnormal blood test results[3, 4]. Improved prediction can lead to more efficient allocation of resources and reduced health costs. A refined initial stratification would aid in allocating resources to patients with a high risk of a BSI and reduce needless testing of patients with low risk. BC may yield not only relevant bacteria causing disease but also growth of contaminating microbes. Ambigous culture results may be difficult to interpret and lead to clinical uncertainty, often resulting in longer hospital stay and unnecessary administration of antibiotics [5, 6]. Therefore, reduced collection of BC in patients with a low risk of BSI can lead to a higher positive predictive value of blood cultures as well as reduced use of antibiotics. Prompt and early identification of high-risk versus low-risk patients is thus imperative for adequate and effective initial handling of patients with suspected BSI, optimized allocation of healthcare resources and reduction of associated costs.

The advent of Artificial Intelligence (AI) allows innovative methodologies for BSI diagnostics, showcasing the potential to enhance or even surpass human expertise in diagnostic endeavors [7]. Despite its demonstrated efficacy, the integration of AI into clinical workflows remains limited [8, 9]. Facilitating this integration may involve leveraging AI models to augment guidelines-based Clinical Decision Support Systems (CDSS) rather than striving for the development of fully autonomous AI-based CDSS [10, 11]. Furthermore, adopting a stance grounded in Predictive, Preventive, and Personalized Medicine (PPPM) principles could refine the utilization of AI, emphasizing the analysis of historical rather than real-time data variables [12, 13]. The primary aim of this study was to develop and validate an eXplainable Artificial Intelligence (XAI) framework that utilizes historical electronic health records (EHRs) to enhance the prediction of BSIs, thus facilitating early and precise intervention. By integrating a rich dataset with advanced machine learning techniques, we sought to overcome the limitations of current real-time diagnostic methods. Our findings confirm that the XAI framework not only meets but exceeds traditional diagnostic performance, demonstrating enhanced predictive accuracy and operational efficiency in BSI management. This achievement marks a significant step forward in the application of AI in medical diagnostics, potentially transforming patient outcomes in healthcare settings.

2. LITERATURE REVIEW

BSI are a critical precursor to sepsis, a severe and potentially life-threatening condition. Early detection and management of BSI can significantly mitigate the risk of progressing to sepsis. In this context, leveraging AI models to predict and manage BSI presents a promising avenue. Our literature review sought to capture the breadth and heterogeneity of recent advances in machine learning ML-based models for BSI prediction. We compiled and analyzed thirty studies published in the last five years, focusing on various healthcare settings and patient demographics [14–44]. These studies predominantly focused on inpatient settings, constituting 56% (n=17) of the research, followed by emergency department (ED) settings at 23% (n=7), and ICU settings at 20% (n=6). Within the inpatient group, the studies varied, with nine examining general populations [14–22], two targeting inpatients with central venous catheters (CVC) [23, 24], and others focusing on specific patient categories such as, hemodialysis (HD) patients [26], cancer patients [27], maternity patients [28], patients with low procalcitonin levels (PCT ≤2.0 ng/ml) [29], and HIV patients [30]. Bacteremia was the primary condition under study in 24 articles, including varied focuses such as fungemia [14] and Candidemia [27, 43], while three studies aimed at predicting central line-associated bloodstream infections (CLABSIs) [33, 34, 41] and one on hospital-acquired BSI (HA-BSI) [22]. All articles reported high prediction performance (AUROC > 0.7) except for one article [20] which reported poor performance metrics. Most studies were single centered, three articles sourced their data from two hospitals [14, 16, 38], and two articles used data from multiple centres [35, 43]. The key characteristics such as, target condition, number of patients or samples, data source, prevalence, ML models, and top predictors for each study grouped by settings are shown in Table 1 on Page No. 28.

Among studies with study design of inpatient settings, Bhavani et al. (2020) used EHRs from two academic tertiary medical centers between 2007 and 2018. Data types included Demographic data, International Classification of Diseases billing codes, clinician orders, blood culture results, vital signs, nursing assessments, and laboratory values [14]. The LR and GBM models were developed. GBM models demonstrated superior performance in predicting bacteremia and fungemia with area under the receiver operating characteristic curve (AUC) scores of 0.78 for bacteremia and 0.88 for fungemia. Lee et al. (2019) investigated the early detection of bacteraemia using an artificial neural network (ANN) model, specifically a multi-layer perceptron (MLP). The study utilized data from 13,402 patients at Gangnam Severance Hospital in South Korea, identifying 1,260 episodes of bacteraemia from blood cultures [15]. Data included 20 clinical variables such as vital signs and various laboratory data. The study highlighted the effectiveness of MLP models, showing remarkable sensitivity in identifying bacteraemia episodes based on a well-curated set of clinical variables. Lee et al. (2022) expanded on earlier research, utilizing an extensive dataset, also from Gangnam Severance Hospital, consisting of a larger patient sample size and covering more recent patient admissions from 2013 to 2018 [16]. This study employed enhanced ANN models alongside other ML techniques like Random Forest (RF) and Support Vector Machines (SVM) to predict bacteraemia from clinical and laboratory data. The study by Mahmoud et al. (2021) developed a predictive model for bacteremia using data from 7,157 adult patients admitted to King Abdulaziz Medical City in Riyadh. This retrospective cohort study utilized EHRs from July 2017 to July 2019 to analyze 36,405 blood culture tests [20]. The dataset included demographic details, clinical variables such as vital signs (e.g., temperature, heart rate, blood pressure), and laboratory test results (e.g., white blood cell count, platelet count, creatinine level, lactic acid level, C-reactive protein (CRP), and procalcitonin levels). Several machine learning models were employed to determine the best predictor of positive blood cultures, including Neural Networks (NN), Random Forest (RF), Logistic Regression (LR), Decision Trees (DT), Naive Bayes (NB), and Support Vector Machines (SVM) with a Radial Basis Function (RBF) kernel. The study by Garnice et al. (2021) utilized a dataset from the Hospital Universitario de Fuenlabrada in Madrid, comprising 4,357 patients with 117 features per patient [21]. The features included patient demographics, medical history, clinical analysis, comorbidities, and the results of blood cultures, differentiated into cases of bacteraemia and no bacteraemia. Three supervised machine learning classifiers were implemented: Support Vector Machine (SVM), Random Forest (RF), and K-Nearest Neighbours (KNN). Each classifier was employed in two scenarios: using only the features available at the time of blood extraction, and a second scenario that included additional features revealed during the blood culture process. These studies highlight the potential of ML to enhance diagnostic accuracy and patient outcomes in detecting critical infections.

Among studies in the ED settings, Choi et al. (2022) developed ML models to predict bacteremia in emergency department (ED) settings at Seoul National University Hospital [31]. The study used data from 24,768 adult patients collected between 2016 and 2018. The models utilized demographic information, chief complaints, vital signs, and laboratory data collected during ED triage and disposition. Two primary models were developed: the Triage XGB model and the Disposition XGB model. In a subsequent study Choi et al. (2023) aimed to refine the predictive accuracy of ED triage-based bacteremia identification using an advanced ensemble of ML techniques. The study analyzed data from over 30,000 ED visits, employing a wide array of clinical inputs including detailed symptom descriptions, vital signs, and initial lab results [32]. The developed model incorporated a Gradient Boosting Machine (GBM) framework that effectively integrated the diverse dataset to predict bacteremia risk. The study by Schinkel et al. (2022) harnessed data from EHRs of 44,123 unique emergency department visits across four hospitals: Amsterdam UMC, Zaans Medical Center, and Beth Israel Deaconess Medical Center, covering the period from 2011 to 2021 [35]. The model employed was a hybrid of LR and (eXtreme Gradient Boosting) XGBoost, with the latter outperforming in predictive accuracy. The data utilized included patient demographics, vital sign measurements, and common laboratory results such as temperature, creatinine, and C-reactive protein, which were identified as the most influential predictors. This predictive model was integrated into the VUMC’s EHRs system for real-time prospective evaluation, affirming its practical utility by potentially reducing unnecessary blood culture analyses by at least 30%. The study by Boerman et al. (2022) utilized a single-center, retrospective observational design, the study encompassed data from 51,399 ED visits at the Amsterdam UMC from September 2018 to June 2020 [33]. Data included demographics, vital signs, laboratory and radiology results, and medications administered during ED visits. The study employed two predictive models: a LR model and a gradient boosted tree (GBT) model, both of which demonstrated good predictive performance with an area under the receiver operating characteristic curve (AUROC) of approximately 0.77 and 0.78 respectively. Notably, the gradient boosted tree model was optimized to predict 69% of blood cultures as negative, with a negative predictive value exceeding 94%, indicating its utility in potentially reducing unnecessary blood cultures and associated healthcare costs. The models harnessed a comprehensive array of features, including commonly available clinical data such as C-reactive protein levels and white cell counts, to predict the likelihood of bacteraemia. These studies illustrate how integrating ML models into ED workflows can improve the speed and accuracy of BSI detection, potentially reducing unnecessary interventions and optimizing resource allocation. However, all the model focused on recent data pertaining to the patients and none of the study utilized the predictors from the complete medical history of their patients, apart from demographic and information on co-morbidities. In our previous work through innovative feature engineering from historical medical records and employing an array of machine learning classifiers, we showcased the efficacy of eXtreme Gradient Boosting (XGBoost) model, in predicting 30-day mortality using EHR trajectory features [56].

2. MATERIALS AND METHODS

2.2. XBSI Framework

The XBSI framework aggregated various data types from the raw EHRs including demographics, laboratory tests including microbiology tests, discharge summaries, and ICU admissions, as depicted in Figure 1. This dataset underwent preprocessing, event log creation for sequential ML models, and feature engineering for static ML models before being transformed and scaled to facilitate the respective model development pipelines. Following subsections details of the various steps in the XBSI framework.

• Download figure
• Open in new tab

Figure 1. Schematic Overview of the XBSI framework.

The diagram illustrates the XBSI framework applied in the study, starting with the extraction of electronic health records (EHRs), which include demographics, laboratory tests, discharge summaries, microbiology tests, and ICU stays. The workflow bifurcates into two parallel processes: Sequence Creation and Feature Engineering, both incorporating Anomaly Scores derived from the data. Subsequent steps include Data Transformation and Scaling, with the data split into training (at 65% and 80%, respectively), validation (at 15%), and testing (at 20%) subsets. The bottom layer of the workflow depicts the range of ML models, the performance metrics used to assess them, and the model explanations generated to study the feature importances.

2.2.1. Data Preprocessing and Transformation

Utilizing Python libraries such as Pandas and NumPy, the CSV files comprising raw EHRs, were loaded into dataframes for further processing. The discharge summaries required several data cleaning steps to ensure quality and relevance of the data: Relevant patient information such as identifiers, admission and discharge times, and diagnostic codes were retained. Instances of missing identifiers were addressed by replacing empty strings with NaN values and subsequently removing these records. Data were organized by patient identifier and admission/discharge times to maintain coherent episode tracking. Non-standard characters within diagnostic codes, such as semicolons and commas, were standardized to spaces, and any duplicates were removed. The timestamps were converted into datetime format, facilitating the calculation of the length of stay in hours for each episode. The duration of each ICU stay was calculated in hours, along with the total count of each type of hospital admissions per patient. The request dates in both laboratory test results and microbiology test results were standardized to datetime objects and used to create eventlogs of tests per patient. Specialized functions were created to clean the laboratory an microbiology test table entries. This function performed tasks to remove any non-numeric characters, which could represent encoding errors or artifacts from data entry. It standardized decimal point characters by replacing commas with periods, which is necessary for consistent numerical representation across different regions that may use varying formats for decimal points. The microbiology test table was filtered to identify suspected BSI episodes, with an emphasis on blood culture tests. The results column was processed to standardize and clean the values, categorizing them as ‘positive’, ‘negative’, or ‘contaminant’ based on the results column. For each patient ID, four event logs were created, from discharge summaries, ICU admissions, laboratory tests and microbiology tests, following subsection describes the event logs.

2.2.2. Event log description

Discharge Summaries Event Log: Captures patient discharge information, including admission and discharge times, diagnostic and procedural codes, urgency and care level code.

ICU Admissions Event Log: Records details of each ICU stay, including the duration in hours and the total count of ICU admissions per patient and total length of ICU stays per patient.

Laboratory Tests Event Log: Includes results of various laboratory tests standardized and organized chronologically for each patient.

Microbiology Tests Event Log: Consists of microbiology test results, grouped by collection sample type, and categorizes them based on outcomes such as ‘positive’, ‘negative’, or ‘contaminant’. Groups of microbiology tests categorized by collection sample type are given in Supplementary List 1.

2.2.3. Sequence creation

The sequence creation process was implemented using the create_sequences function, designed to compile a comprehensive view of a patient’s medical history over their entire recorded history. The create_sequences function systematically constructs a timeline of medical events for each patient, for each BC test. The medical events were merged from the four event logs. The filtered event logs are merged into a single dataframe, ensuring no information is lost. This step involves an outer join on patient ID and date, maintaining all records from each event log.

2.2.4. Feature engineering

This approach involved creating a dataset with attributes derived from ICD codes, procedure codes, laboratory, and microbiology test results. Laboratory test results were organized using pivot tables, ensuring a structured format for analysis. Tests such as ‘bilirubin’ (total, conjugated, and unconjugated), ‘C-Reactive Protein’ (CRP), and ‘lactate’ (various measurements) were included, alongside white blood cell count (leukocytes), platelet count (thrombocytes), and blood gas measurements (pH, PO2). Similarly, microbiology test results were consolidated to reflect various sample types such as blood, urine, and other fluids, employing a dictionary mapping to streamline similar types. The resulting pivot table included columns for diverse samples, ranging from ‘blood culture test’ to ‘urine’, ‘feces’, and ‘nasal swabs’. A function was created to determine the aggregated result of blood culture tests for each patient group, considering the possibility of concurrent positive and contaminant results. This aggregation provided a comprehensive view of the infection status per patient. Moreover, the history of prior positive, negative, and contaminant results was calculated and added to the dataset, offering a valuable perspective on the patient‘s previous encounters with BSIs. Comorbidities were extracted and processed to identify unique diseases from patient records. New columns were created for each disease, and the counts were updated based on patient history. Finally, aggregate columns were added for both ICD and procedure codes to calculate cumulative sums. This method enabled the capturing of the cumulative history of medical conditions and procedures for each patient. All the NaN values representing absence of a condition or measurement were filled with zero and “0” is not interpreted as a value.

2.2.5 Model Development

This section outlines the development and evaluation of ML models. The process bifurcates into the sequence creation and feature engineering pipelines for the temporal and static ML pipelines respectively. An Isolation Forest model was implemented to detect anomalies within the data. This model was trained exclusively on normal data (negative BC episodes) and then used to compute anomaly scores for the training, validation, and test sets. These scores were normalized and appended to the original dataset to serve as additional features, enhancing the model’s ability to distinguish between normal and anomalous patterns. Given the sequential and tabular nature of the input data for the respective ML pipelines, the necessary transformation steps were employed to prepare it for the corresponding learning algorithms. The data was reshaped and scaled using the StandardScaler from sklearn.preprocessing. This normalization step is crucial for models that are sensitive to the scale of input features. After scaling, the data was reshaped back to its original form, ensuring compatibility with the machine learning models used. Temporal ML models included Long Short-Term Memory (LSTM), Gated Recurrent Unit (GRU), Convolutional Neural Network-Long Short-Term Memory (CNN-LSTM), Convolutional Neural Network-Gated Recurrent Unit (CNN-GRU), and Transformer models. These models were trained on sequences of medical events to capture temporal dependencies within the data. On the other hand, the static ML models, which treat data points as independent and identically distributed, comprised of Light Gradient Boosting Machine (LightGBM), eXtreme Gradient Boosting (XGBoost), Categorical Boosting (CatBoost), Artificial Neural Network (ANN), Random Forest (RF), and Logistic Regression (LR). Detailed information of each ML model is given in Supplementary Methods, Section 1.1., 1.2.

2.4.6 Model Validation

For cross validation, we chose a 5-fold cross-validation scheme. This technique involves dividing the entire dataset into ‘k’ equally sized subsets or folds. The model is then trained on ‘k-1’ folds and tested on the remaining fold. This process is repeated ‘k’ times, each time with a different fold as the validation set. In addition to cross-validation, we also validated our model on an independent test set, which was separated from the dataset at the outset and not used during the training phase. The test set comprised of the most recent 20% of the BC episodes while the earliest 80% of BC episodes were assigned as training cohort. Further, the training data was divided into a smaller training set and a validation set, with 15% of the training data allocated for validation. This additional split allowed for the tuning of hyperparameters and the assessment of the model’s performance during training. Each model was evaluated using a variety of metrics to evaluate model performance, including accuracy, precision, recall, F1-score, specificity, the Area Under Precision-Recall (PR) Curve (AUPRC) and the Area Under the Receiver Operating Characteristic (AUROC). Detailed information on each performance metric is given in Supplementary Methods, Section 1.3.

2.4.7 Model Explanations

For interpreting the outputs of our ML models, we employ both global and local explanation methods to provide comprehensive insights into feature importances which guide clinical decision-making. To evaluate the overall impact of features across all predictions, we generate SHAP (SHapley Additive exPlanations) summary plots [55]. These plots aggregate SHAP values to illustrate the average influence of each feature on the model output, ranked by significance. This method, rooted in game theory, decomposes a prediction into the contribution of each feature, providing a transparent view into the predictive process. SHAP’s suitability for tree-based models was enhanced by its efficient computation of exact SHAP values using the TreeExplainer algorithm [57], which significantly reduced computational complexity by exploiting the structural properties of decision trees [58]. By integrating SHAP summary plots, we can discern which features are most instrumental globally in predicting bloodstream infections, facilitating an understanding of model behavior over the entire dataset On an individual prediction level, we utilize waterfall plots and force plots to detail the contribution of each feature to specific predictive outcomes. Waterfall plots provide a step-by-step breakdown of how each feature’s value contributes to the final prediction, starting from the base value (the average model output across all data points) and adding the effect of each feature sequentially. This visualization helps in understanding the decision-making process for individual predictions, which is crucial for clinical validation and personalized patient insights. Force plots, another type of local interpretability tool, display how each feature’s value pushes the model’s prediction higher or lower, which is particularly useful for individual patient assessments. These plots highlight the positive or negative contribution of each feature towards the final prediction, allowing healthcare providers to grasp the underlying reasons for a model’s decision on a case-by-case basis.

3. RESULTS

3.1. Patient Characteristics

There were a total of 72,495 BC episodes in the dataset. Following the exclusion of pediatrics and outpatient BC episodes 65,975 adult inpatient BC episodes were included in the analysis. Of the BC episodes, 5,288 (8%) were classified as positives. Please see the flow chart provided in Figure 2. The differences in the mean values for all feature across the two classes along with the T-statistic and p-value is given in Table (Supplementary Results, Supplementary Table 2). The top 25 most significant features and the comparison between their mean values across the two classes, along with the T-statistic and p-value are given as Table 1. There were significant differences between the positive and negative blood culture groups regarding the occurrence of ICD codes starting with A, B, and N in the medical records. The ICD Chapter I (A00-B99) describes certain infectious diseases and parasitic diseases and ICD Chapter XIV (N00-N99) concerns kidney and urinary disorders. From laboratory markers, bilirubin, creatinine, and CRP showed the most significant differences between the groups. The positive BC group also had higher mean age, higher total length of stay (LOS) till the time of BC, shorter time to the last, most recent episode, and a higher number of previous positive BC test results. The correlation matrix heatmap depicts a visual representation of the strength and directionality of correlations between various clinical features within the EPRs, highlighting the complex interdependencies relevant to BSI prediction, and is presented in Supplementary Figure 5. The top three most correlated coefficients among significant features is given in Supplementary Table 3. The list of all the features, their description, and the mean value across the dataset is given in Supplementary Table 4.

• Download figure
• Open in new tab

Figure 2. Flow chart

The flowchart depicts the categorization of the BC episodes.

3.3. Global Feature Importances

The SHAP summary plots in Figure 3. reveal the feature importances derived for the XGBoost model. The bar plot indicates that the count of ICD codes starting with letters A, B, and N, bilirubin (BILIRUBIN TOTAL), creatinine (KREATININ), leukocytes (LEUKOCYTTER), thrombocytes (TROMBOCYTTER), age, number of previous positive blood culture results, and the LOS of current or recent hospital episode were the top ten most important features for predicting positive blood culture. The beeswarm plot further elucidates the directional impact of each feature’s value on the model’s predictions. Higher values of bilirubin, creatinine, and leukocytes (as indicated by red dots on the right side of the x-axis) lead to a greater increase in the model’s prediction values, potentially signifying a higher risk whereas a lower value for thrombocytes signifies a higher risk of bloodstream infection. Other significant factors towards a positive prediction were higher cumulative hospital length of stay (total_los), lower time to the recent episodes, and higher CRP.

• Download figure
• Open in new tab

Figure 3. SHAP Summary Plots for XGBoost Model Feature Importance.

Bar plot on the left illustrating the global feature importance ranked by the sum of SHAP values across all samples. On the right is the Beeswarm plot detailing the individual SHAP values for each feature and their impact on the model’s output.

3.4. Local Feature Importances

The figure 4. on Page No. 32. Details the waterfall and force plots for the first three prediction tasks. In the first prediction task, the feature ICD_A (representing infectious and parasitic diseases) has a significantly high positive SHAP value, indicating that an increase in counts of ICD-10 Chapter I codes strongly sways the model towards predicting a bloodstream infection. The ‘anomaly_score’ feature contributes negatively, suggesting that anomalies detected by the model are not predictive of the outcome. Similarly, creatinine level and ICD_N (genitourinary diseases) show negative SHAP values, implying that higher levels of creatinine and the presence of genitourinary diseases are linked to a lower probability of bloodstream infection in this instance. For the second prediction, counts of prior positive BC results has the most substantial positive impact, which aligns with clinical reasoning that past positive tests could indicate a higher risk for current infection. Additionally, ICD_B (infectious and parasitic diseases) and ICD_R has a positive effect, and total bilirubin has a positive contribution, suggesting high risk of BSI. In the third prediction, similar to the second prediction, counts of prior positive BC results again shows a positive influence. Additionally, ‘ICD_C_aggregate’ (aggregate count of cancer-related codes) exerts a negative influence on the prediction outcome, whereas ‘procedure_W_aggregate’ (count of procedures on female reproductive organs in the history) is associated with an increase in the risk of BSI.

• Download figure
• Open in new tab

Figure 4.

Waterfall and Force plots for the first three predictions from the test set.

4. DISCUSSIONS

In this study we presented the XBSI framework for early prediction of BSI in the hospitals. For the sequential models, the AUC ranged from 0.6961 to 0.7182, indicating a modest discriminatory capacity, with CNN-LSTM marginally leading in performance. Notably, the Transformer model demonstrated the highest sensitivity 0.4688. In contrast, static models depicted better AUC values ranging from 0.7751 to 0.8410. Our findings underscore the superior performance of static ML models, particularly those employing tree-based algorithms such as CatBoost, which achieved a notable AUC of 0.8167 [52]. The statistical analysis and global and local SHAP value interpretations reveal that the count of diagnostic and procedural ICD codes and routine laboratory significantly contribute to the model’s predictions. These findings align with existing literature, underscoring the clinical relevance of laboratory tests and medical history for early BSI prediction [14, 15, 16, 20, 21, 51]. The demonstrated efficacy of our XBSI framework highlights the predictive capability of combining medical history with routine laboratory tests for early BSI detection. The modular design of our framework reinforces robustness to all types of medical data and emphasize the utility of integrating diagnostic and procedural ICD codes into predictive models. This approach allows for a nuanced understanding of patient profiles, which is pivotal for implementing predictive healthcare solutions effectively.

Comparing the recent studies on predicting BSI in hospital settings, Bhavani et al. (2020) utilized LR and GBM models to predict bacteremia and fungemia from EHRs, achieving AUCs of 0.73 and 0.88 respectively. The key predictors identified in the study were time from admission to BC, temperature, age, heart rate, prior bacteremia/fungemia, white blood cells (WBC), blood urea nitrogen (BUN), glucose, diastolic blood pressure (DBP), and systolic blood pressure (SBP) [14]. Lee et al. (2019) compared multiple ML algorithms, including MLP, SVM, and RF, achieving their best AUC of 0.732 with the RF model [15]. Our study surpasses this benchmark, suggesting that the sequential analysis of patient data may capture the temporal dynamics of BSI development more effectively than models relying solely on static point-in-time data. The key predictors reported by the study included, alkaline phosphate (ALP), platelet, maximum body temperature, SBP, WBC, CRP, ICU stay, hospital day to blood culture, age, heart rate, prothrombin time, and albumin. In a subsequent study, Lee et al. (2022) further explored MLP, RF, and XGBoost among others, focusing on their application over a long-term dataset. Their models achieved AUCs of 0.762 for MLP and 0.758 for RF in the 12-hour data group, which are slightly lower than our CatBoost model’s performance. the key predictors in this study were monocyte, platelet, hospital stay, neutrophil, total bilirubin, BUN, albumin, ALP, WBC, CRP, creatinine, pulse rate, and chloride [16]. The study by Mahmoud et al. used data from a tertiary care center comprising patient demographics, length of stay before blood culture collection, presence of central line, vital signs, laboratory results, and SIRS and qSOFA scores [20]. They employed various models including NN and LR, with their best models achieving modest performance metrics (highest specificity at 89% but with low sensitivity). This study extensively used vital signs and other real-time clinical parameters such as temperature and heart rate as predictors. Despite this, the highest sensitivity achieved was only 31% with Logistic Regression, and even though some of their models achieved high specificity, they struggled with very low sensitivity, limiting their practical utility in clinical settings. Our model’s ability to outperform these studies without the need for immediate clinical data or vital signs not only showcases the robustness and efficiency of our approach but also suggests a potential for earlier and simpler implementation in clinical workflows. This is particularly advantageous in healthcare settings where immediate comprehensive data collection is challenging, offering a powerful tool for early BSI detection that is less dependent on the specific timing of clinical data acquisition. Similar to our study, the study by Garnica et al. (2021) used SVM, RF, and K-nearest neighbors (KNN), with a combination of RF and SVM yielding the most robust performance metrics [21]. The number of days in ICU before BC extraction, presence of catheters, age, chronic respiratory disease, fever, and CRP were reported as key predictors. Their better perfomance in comparison can also be attributed to the fact that they used a BSI dataset with prevalence rate of 51.3% which is highest reported among all the studies reporting ML based BSI prediction models.

Our study’s findings also contribute to the evolving landscape of ML applications for BSI prediction in ED settings. Schinkel et al. employed LR and XGBoost models using vital signs, laboratory results, and demographics [35]. Their XGBoost model achieved an AUROC of 0.81 (95% CI 0.78-0.83), slightly higher compared to our XGBoost model AUROC of 0.7951 (95% CI 0.7833-0.8069) and an AUPRC of 0.34 (95% CI 0.29-0.38) slightly lower than our XGBoost model AUPRC of 0.3942. The key predictors brought out by this study were temperature, creatinine, CRP, lymphocytes, DBP, bilirubin, thrombocytes, neutrophils, APL, heart rate, SBP, leukocytes, glucose, age, potassium, BUN, sodium, monocytes. Their model effectively reduced unnecessary blood cultures by approximately 30% during real-time prospective evaluation, which aligns with our goals of enhancing diagnostic efficiency and reducing healthcare costs. Our study builds on this foundation by implementing a similar ML approach but extends its application by incorporating a wider array of clinical variables derived from historical EHRs and employing a novel algorithmic configuration that may provide improved predictive performance. Similar to our study, Boerman et al. focused on the ED setting of a large teaching hospital, developing predictive models specifically for BSI outcomes based on data available at the end of ED visits, such as demographics, vital signs, administered medications, and laboratory results. They reported an AUC of 0.77 for their gradient-boosted tree (GBT) model and 0.78 for their LR model, indicating good performance in predicting bacteremia in ED [33]. The key predictors in their models were bilirubin, urea, lymphocytes, pulse rate, CRP, neutrophil, age, Temperature, DBP, potassium, glucose, thrombocytes, creatinine, ALP, SBP, and organ damage. The study highlighted the ability of their models to significantly reduce unnecessary blood cultures by predicting negative outcomes with a high degree of accuracy, reflected by a negative predictive value of over 94%. In parallel, our model not only confirms these findings but also demonstrates improved predictive accuracy and efficiency in differentiating between positive and negative BSI outcomes, which could further optimize the use of resources in ED settings. In another study, Choi et al. (2022) demonstrated the effectiveness of XGBoost model, to predict bacteremia at different stages of patient care in the ED, achieving an AUROC up to 0.853 [31]. The key predictors in the study were chief complaint, age, temperature, heart rate, and DBP at triage stage, and neutrophils, platelets, CRP, chief complaint, and creatinine at disposition stage. Their phased approach using predictions at both triage and disposition stages aligns with our methodology of employing dynamic modeling to adapt predictions based on real-time data updates. However, our model differentiates itself by utilizing a more complex array of inputs and a refined computational technique that enhances predictive generalizability. By extending the predictive capabilities beyond structured data, our study underscores the importance of comprehensive data utilization in developing models that can more effectively guide early clinical decisions in real-time settings. Our model’s ability to accommodate a broader array of clinical variables may explain any improvements in prediction robustness compared to the framework used by Choi et. al. (2023) [32]. The key predictors in their study were age, vital signs, history of chills, ambulance use. The collective insights from these comparisons suggest that while our model shares common ground with existing approaches, it also explores additional layers of complexity such as patient recent hospital interactions, comorbidities, and previous history of infections, which may influence the generalizability and effectiveness of the model across diverse healthcare environments, capturing more subtle nuances of BSI risk factors and could be pivotal in reducing BSI misdiagnosis.

Our framework’s modular design facilitates the inclusion of additional data types as they become available, enhancing its adaptability across various clinical settings. This feature is particularly valuable in current healthcare environments where flexibility and comprehensive data utilization are crucial for advancing diagnostic accuracies. Moreover, the intuitive nature of the XBSI framework ensures that it can be seamlessly integrated into existing clinical workflows, making it a practical tool for clinicians seeking to leverage AI for improved patient outcomes. The principal clinical value of our approach lies in the ability to identify patients at low risk of a positive BC, at the time of suspicion of BSI without the need for waiting to capture vital signs at the moment or within a specified time window, which could increase patient risk and stay [53, 54]. Integrating our proposed framework in hospitals as a pre-emptive BSI prediction tool can reduce BC ordering and its resulting costs and harms [49]. As it is reported in the literature that the use of data that is not routinely captured in clinical practice is the main reason why none of the prediction models have been implemented in the clinical practice yet [50]. Moreover, in the context of ICU patient monitoring, the application of AI to real-time data may seem redundant. Patients in critical care are already under intensive surveillance, and the utilization of AI for immediate alert systems could potentially clutter the workflow rather than enhance it. Instead, a more strategic use of AI lies in its ability to predict a patient’s worsening condition well before critical thresholds are reached. By analyzing historical EHRs, AI can identify subtle patterns and indicators of decline that precede acute episodes, thereby enabling preemptive medical interventions. Moreover, as the complexity of predicting health events increases the earlier the prediction occurs, achieving perfect accuracy, as often highlighted through metrics like AUC in ICU settings [38–43], may not be as critical as maintaining reasonable predictability at the initial stages of patient contact. This approach could shift the focus from crisis response to proactive patient care management, optimizing outcomes through early and targeted intervention. The only limitation of this study is its reliance on data from a single center, which may not represent the diverse patient demographics. Despite these constraints, our model’s ability to integrate a broader spectrum of data and apply a real-time learning framework supports its potential utility in clinical settings, promising reductions in unnecessary interventions and improvements in patient management. Furthermore, our use of historical patient data and variety of ML techniques allowed for a more flexible and scalable prediction framework. This aspect of our model enhances its applicability and provides a significant contribution towards advancing BSI prediction using AI.

5. CONCLUSIONS

This study presents a simple and intuitive XAI framework which comprehensively captures the complete medical history of a patient to accurately predict the risk of BSIs. The XBSI framework is the first of its kind to significantly enhances the predictive analysis by integrating information stored as diagnostic and procedural ICD-10 codes. This novel approach diverges significantly from traditional real-time monitoring systems, emphasizing incorporating comprehensive historical patient data. The key predictors include the count of ICD-10 codes for infectious diseases, and kidney and urinary disorders, total bilirubin, creatinine, leukocytes, thrombocytes, age, prior positive blood culture tests, and the total length of hospital stay till the blood culture test. To further enhance the generalizability of our findings and adoption of our framework, there is a need to validate our framework on administrative datasets of hospitals outside Norway and include more diverse data sources, such as biomarkers, genotypes, and phenotypes in future works.

Author’s contributions

RB conceptualized and designed the study, with inputs from ØN, JEA, BE and JKD. ØN and LGT secured the data and funding for this project. RB performed data preprocessing, model development and interpretation. RB wrote the initial draft of the paper, to which the remaining authors provided comments and feedback. All authors reviewed and approved the final manuscript.

Availability of data and materials

The jupyter notebooks implementing the data preprocessing, sequence creation, feature engineering and model development pipeline are available through the corresponding author. The processed and transformed final datasets, and the list of derived features are available through the corresponding author on a reasonable request.

Ethics

The use of the EPRs data in this project has been approved by the Regional Commitees for Medical and Health Research Ethics (REK) in Central Norway by REK no. 2020/26184.

Conflict of interests

None declared.

Funding

Financial support for this study was provided by the Computational Sepsis Mining and Modelling project through the Norwegian University of Science and Technology Health Strategic Area.