Machine Learning Predicts Unplanned Care Escalations for Post-Anesthesia Care Unit Patients during the Perioperative Period: A Single-Center Retrospective Study

Background: Despite low mortality for elective procedures in the United States and developed countries, some patients have unexpected care escalations (UCE) following post-anesthesia care unit (PACU) discharge. Studies indicate patient risk factors for UCE, but determining which factors are most important is unclear. Machine learning (ML) can predict clinical events. We hypothesized that ML could predict patient UCE after PACU discharge in surgical patients and identify specific risk factors. Methods: We conducted a single center, retrospective analysis of all patients undergoing non-cardiac surgery (elective and emergent). We collected data from pre-operative visits, intra-operative records, PACU admissions, and the rate of UCE. We trained a ML model with this data and tested the model on an independent data set to determine its efficacy. Finally, we evaluated the individual patient and clinical factors most likely to predict UCE risk. Results: Our study revealed that ML could predict UCE risk which was approximately 5% in both the training and testing groups. We were able to identify patient risk factors such as patient vital signs, emergent procedure, ASA Status, and non-surgical anesthesia time as significant variable. We plotted Shapley values for significant variables for each patient to help determine which of these variables had the greatest effect on UCE risk. Of note, the UCE risk factors identified frequently by ML were in alignment with anesthesiologist clinical practice and the current literature. Conclusions: We used ML to analyze data from a single-center, retrospective cohort of non-cardiac surgical patients, some of whom had an UCE. ML assigned risk prediction for patients to have UCE and determined perioperative factors associated with increased risk. We advocate to use ML to augment anesthesiologist clinical decision-making, help decide proper disposition from the PACU, and ensure the safest possible care of our patients. Supplementary Information The online version contains supplementary material available at 10.1007/s10916-024-02085-9.


Introduction
Since the advent of general anesthesia, the risk of mortality for elective procedures in the United States and other developed countries has never been lower [1,2].Despite improvements in perioperative medicine, a fraction of patients experience unexpected care escalations (UCE) following post-anesthesia care unit (PACU) discharge [3,4].UCE can be defined as the unforeseen need for higher acuity care, i.e. unanticipated transfer from a standard inpatient unit to an intermediate level unit, from a standard or intermediate level unit to an intensive care unit (ICU), or intervention by a Medical Emergency Team (MET).These UCE increase morbidity for patients, increase medical cost and length of hospital stay, and slow rehabilitation and recovery for patients undergoing elective surgery [3,5].
There are a variety of clinical risk factors that may increase a patient's risk of UCE.For example, having emergency surgery increases risk for ICU admission during the perioperative period [6,7].Furthermore, patients with risk factors for respiratory insufficiency after surgery, such as obesity, a diagnosis of asthma or chronic obstructive pulmonary disease (COPD), or insufficient paralytic reversal are more likely to have post-operative pulmonary complications [8,9].Finally, pre-operative co-morbidities such as poor nutrition, end-stage renal disease (ESRD), and chronic anemia are risk factors for poorer outcomes after surgery and anesthesia [9,10].While these studies indicate that there are important risk factors that may predict UCE after surgery and anesthesia, it may be difficult to determine how a combination of multiple factors conspire to predict occurrence of an UCE.
Machine learning (ML) is gaining traction in medicine and has been used to perform some tasks as well or better than well-trained physicians.Two of these contexts include interpretation of radiologic imaging [11,12] and determination of cancer from pathology samples [13,14].Recently, ML has turned its attention to "big data" and the ability to predict clinical events before they occur [15,16].For example, studies indicate that ML has the ability to predict hypotension in the PACU [17], to predict the occurrence of ventilator-associated pneumonia after traumatic brain injury [18], and to detect heterogeneity of treatment effect in acute respiratory distress syndrome [19].
Therefore, we conducted a single center, retrospective study hypothesizing that a ML algorithm could detect patients at risk for UCE and provide risk factors that were important to UCE risk.We collected data from pre-operative visits, intra-operative records, and PACU admissions and determined which patients had an UCE.We trained a supervised ML model with this data and, after training, tested the ML model on a separate data set to determine efficacy and patient clinical factors that may be important in determining patient UCE risk.

Methods
This study was conducted in accordance with the 1964 Declaration of Helsinki and the protocol was approved by The University of Alabama at Birmingham Institutional Review Board (Protocol Number-300007621).Consent was waived by the IRB.This manuscript adheres to the STROBE guidelines.

Study Design and Clinical Parameters
We conducted a single center, retrospective analysis of all patients undergoing non-cardiac surgery at the University of Alabama at Birmingham from January 2016 to December 2019.All surgical patients were included if they had noncardiac surgery (elective or emergent), were admitted to the PACU after surgery, and were admitted to an inpatient floor or intermediate care after PACU discharge.Patients were excluded if they had cardiac surgery, had outpatient surgery (discharged home after PACU stay), were directly admitted to an ICU, were < 18 years-old, or we could not determine whether they had an escalation event.We collected clinical parameters and data from patients in the form of labs collected before surgery, patient characteristics, intraoperative data, and data collected during their PACU admission from Cerner PowerChart and CompuRecord electronic medical records.These features were chosen because they represent pre-operative (hemoglobin, bicarbonate, etc.), intra-operative (blood loss, transfusion, etc.), or post-operative (vital signs, urine output (UOP), etc.) measures of end-organ dysfunction and/or morbidity.These data are completely listed in Table 1 and Supplemental Table 1.We defined an UCE as an increase in care from the inpatient floor to intermediate care or Surgical ICU (SICU), an increase in care from intermediate care to the SICU, or a MET response within 3 midnights of PACU discharge.We made this definition based on Major Adverse Cardiac Events (MACE) criteria that myocardial infarction risk (a major perioperative clinical outcome) decreases significantly after 3 midnights [20][21][22].Therefore, after this timeframe, UCE is less likely to be secondary to perioperative events.For some patients, there was a MET response without escalation or delayed escalation.If this was the case, the entire spectrum of such events counted as a single UCE event.

Machine Learning Algorithm
A "credit scorecard" modeling system [23] was applied as a set of pre-and post-processing steps layered on top of regularized logistic regression.This modeling methodology has exhibited some recent success in predicting COVID-19 risk from electronic health record data [24][25][26].Variables employed are listed in Table 1 and Supplemental Table 1.80% of patients were chosen randomly as a training cohort for ML algorithm "learning" how to determine factors that attributed UCE risk to a patient.Our supervised ML model was trained with examples where the correct answer (patient had an UCE or not) is known and used in the training process.After training, the remaining 20% of patients were used as a test cohort to determine how well the ML algorithm had learned.Continuous variables were binned into discrete categories and extant categorical variables were grouped using weight-of-evidence (WoE) binning via the Python package OptBinning [27,28].WoE transformations are appropriate for electronic health record data, which often contains missing and (potentially erroneous) extreme values -situations for which WoE is particularly well suited [29].Post-regression, WoE values and model coefficients are used to calculate scorecard points [23].While other machine learning models can handle continuous variables, the "credit scorecard" model we used requires binned variables to assign specific score points per bin.By binning the continuous variables, we provide a model that is not only predictive but also intuitively understandable and usable by healthcare providers in a real-world clinical setting.Additionally, WoE binning process of continuous variables also has advantages in terms of robustness and stability.Furthermore, it has the "nice" theoretical property of ensuring a monotonic relationship of predictors and outcomes.Additionally, it can mitigate the influence of outliers and less frequent or missing values which might skew the results in a model that uses continuous values without binning [26,29]. is a physical status classification system to score patient fitness before surgery.Of note, in the p−value column, "C" denotes use of Chi−squared test and "F" denotes use of Fisher−exact test right on the plot indicates an increase in risk probability, where some lines move left across the variables to indicate a decrease in probability of UCE risk.The variables are ordered by magnitude of importance.
Continuous variables were compared using a two-sample t-test and categorical variables (including Boolean variables) were compared with Chi-squared tests or Fisher's exact test, as applicable.Correction for multiple comparisons was performed using Hommel's method [35], as implemented in the Python package "TableOne" [36].A p-value of < 0.05 was considered significantly different.

Results
We conducted a retrospective analysis of all surgical inpatients at the University of Alabama at Birmingham from January 2016 and December 2019.Patients were included in the study and had UCE determined as described in the Methods.We screened 138,929 patients during the study timeframe and 79,998 were excluded per the Methods (Fig. 1).A total of 58,931 patients were included in the study and 47,144 (80%) were randomly allocated to the training group and the remainder, 11,787, were allocated to the testing group.The training and testing groups had an UCE incidence rate of 5.12% and 5.37% (p > 0.99), respectively, indicating no significant differences.Additionally, the reason(s) for UCE (MET, ICU admission, and intermediate care admission) occurred at similar rates between the training and testing groups.
Table 1 shows patient demographics and clinical characteristics for the overall cohort (n = 58,931), training group (n = 47,144), and testing group (n = 11,787).For the overall cohort, 46.4% were male, mean age was 51.7 years, and mean body-mass index was 30.4.Roughly 50% of the cohort was Caucasian and 32.8% were African-American.At the time of surgery, 12.3% had chronic kidney disease and 3.8% had ESRD; 18.3% developed acute kidney injury (AKI) during the preoperative period.Before surgery, 3.3% of patients had COPD, 4.5% had asthma, and 12.5% had diabetes listed on their pre-operative problem list.The mean score for MACE was 0.9, the most prevalent ASA Status score was 3 (69.8%),and 9.3% of cases were emergencies.UCE rate was 5.2%.There were no differences in demographics or clinical characteristics between the training or testing groups.A more complete list of variables collected for the overall cohort, training group, and testing group is available (Supplemental Table 1).In both Table 1 and Supplemental Table 1, we include a column that lists the number of events that are missing from a discrete category.For example, there are no data missing for sex or age, but there are 1,353 events missing for the overall cohort for race.
The resulting model is straightforward to interpret and easy to implement.
Elastic-net regularized logistic regression served as both a variable selection method and hedge for colinear variables [30].The optimal regularization strength and mixing parameter for elastic-net regression was chosen using 10-fold cross-validation.The final scorecard model range of values was enforced upon the model by solving a mixedinteger programming problem constrained by the WoE bins and results of the regularized logistic regression model.This final step enforces whole-number point values, avoiding decimals in the final scorecard -adding to the ease of interpretation.
The final total score needs to accurately reflect the patient's risk of UCE, but for interpretability and ease of calculation, we wanted this to be a number between 0 and 100.This goal of getting a score that both reflects the risk accurately and stays within a specific range is called the objective function in a mixed integer programming nomenclature [31].
The mixed integer programming model adjusts the points assigned to each patient factor (called the decision variables) to optimize the scoring within the set constraints.The constraints themselves, such as the total score needing to be between 0 and 100, are the parameters.The key outcome of this process is a set of scorecard points for each factor, which when added together for a patient, gives a total score that predicts the risk of UCE within the desired range of 0 to 100.
Using mixed integer programming here is akin to finely tuning a balance scale.The algorithm adjusts the weights (the points for each factor) to achieve a balanced and accurate reading (the total score).This must be done so that the final scores not only remain within the desired scoring range of 0 to 100 but also as closely as possible reflect the true risk of UCE based on the underlying logistic regression model.
For those wishing to reproduce our work, a copy of the code and parameters used to produce our scorecard model is available at https://github.com/UABPeriopAI/post-pacu-uce.

Statistical Analysis
Shapley Additive exPlainations (SHAP) is a framework built around Shapley values to explain ML models [32][33][34].By iteratively removing each variable from the model and looking at the effects on all other variables, Shapley values can quantify the influence of each variable.Decision plots are examples of how Shapely values influence the estimated probability.Each line in the plot represents a patient and the color gradient goes from red to blue as highest to lowest estimated probability of UCE.The line moving further second bin (Supplemental Fig. 1).These data indicate that the model functions as expected within a dataset collected from the same institution.
After model training, we examined how the model predicted risk for each patient characteristic or clinical parameter that was provided.For each datum, the model created bins for continuous variables and assigned risk with points for each bin (Table 2).For example, non-surgical anesthesia time is a risk factor for UCE and the longer this time is the higher the risk the model assigns.Additionally, the model applies appropriate risk to heart rates that are too low or too high in the last hour of PACU care.Finally, low SpO 2 values are accorded more risk by the model.
After acquiring available patient data, we trained a supervised ML model to detect differences in patients who had an UCE within 3 midnights of PACU discharge with those who did not.We used the Credit Card modeling system as it incorporates missing data and includes it as a risk factor.Using our training data, the model was able to predict UCE risk with an AUC of 0.757 (Fig. 2A).After the model had learned using 80% of patients in the study, it was applied to a test group of patients and achieved an AUC of 0.752 (Fig. 2B).These are both acceptable AUC scores with ML models [37].Across the training and test groups, the model found similar percentages of patients within the risk bins it created, and within each of these bins, the patient's UCE risk is not significantly different, with the exception of the We next sought to understand the variables that had the greatest effect on patient UCE risk.We found the 10 most influential variables for patients whose estimated probability of UCE was greater than 20% (Fig. 3).For each   3 Shapley values that most affect patient UCE risk.On the left, the 10 most influential variables for patients whose estimated probability of UCE was greater than 20% are listed.Each patient is visualized as a path moving from lower risk (left-hand side) to higher risk (righthand side) as variables of increasing importance (moving from bottom to top) are evaluated for that patient.The greater a patient moves to the right, the more risk that variable provided for that patient.If they move to the left, the variable caused a decrease in risk of UCE for that patient.Of note, any value labeled as "Maximum PACU" (e.g.HR, RR, Aldrete, SBP) is measured in the last hour of the PACU stay find these patients and identify specific risk factors associated with UCE.Our study obtained the following significant results: (1) the rate of UCE was ~ 5%; (2) ML can be used to find patients at risk for UCE; (3) there are specific risk factors associated with UCE; and (4) after training, the algorithm identifies risk factors indicative of increased UCE risk similar to that known by anesthesiologists and the literature without prompting.
The UCE rate was approximately 5% in the training and testing groups.This is an increased incidence rate compared to other studies that have examined unplanned ICU admission [3].It is generally accepted that the risk of post-operative MACE decreases three days after surgery [20][21][22].We therefore used "three midnights" as our time-frame for the assessment of post-operative UCE.Furthermore, we define UCE as an incident of increased patient care level or MET response; including from the inpatient floor to intermediate care or an unplanned ICU admission.We included MET response within our definition as this typically indicates a point of clinical decline requiring medical intervention.These broader definitions are a likely reason we demonstrate an increased UCE incidence.
Numerous ML algorithms detect risk within specific patient populations or problems better than physicians.A recent example includes determining proper endotracheal tube depth [39].Our algorithm can analyze individual patient data prior to PACU discharge and assign UCE risk levels within the following three midnights to individual patients.A primary motivation for our a priori approach (scorecard modeling) was to enhance the interpretability and clinical applicability of our model.While other machine learning models can handle continuous variables, the "credit scorecard" model we used requires binned variables to assign specific score points per bin.The scorecard model was selected because it can provide a straightforward, easyto-understand risk score for clinicians.The Aldrete Score as determined by SHAP values were non-surgical anesthesia time, AKI, maximum heart rate in the last hour before PACU discharge, ASA status, and maximum respiratory rate in the last hour before PACU discharge.Complete analysis of all variables and their influence in described in Supplemental Fig. 2. To further understand which specific bins within each variable provided the highest patient UCE risk, we examined any bin variable that provided 5 or more risk points from the ML algorithm (Table 3).For reference, the maximum score of the scorecard is 100, while the minimum is 0. For a visualization of the association of a scorecard score from this model with actual risk, see Supplemental Fig. 2.There are thirteen bins that provide 5 or more risk points: non-surgical anesthesia time, SpO 2 , systolic blood pressure, heart rate, respiratory rate, emergency case, ASA status > 3, and Aldrete score.Interestingly, non-surgical anesthesia time was the most important variable identified by SHAP values.Any duration longer than 53.5 min scored 5 or more points with longer durations equating to larger point totals.A non-surgical anesthesia time duration longer than 88.5 min was the largest single contributor to the ML model with 10 points.Finally, and intriguingly, a missing value that should have been measured in the last hour before PACU discharge also plays an important role in the risk probability for an UCE.For example, a missing SpO 2 measurement, systolic blood pressure measurement, or Aldrete score in the last hour before PACU discharge all caused the ML model to score 5 points to patient UCE risk.

Discussion
We hypothesized that ML could predict patients that may experience UCE and identify factors that increased UCE risk.We compared patients who experienced an UCE to those that did not to determine whether the model could made.Furthermore, testing this variable prospectively will improve the algorithm and, perhaps, find other risk factors for patient UCE that can be intensely studied.Both of the above findings are sources of future investigation.
Data from the model indicate some demographic and/ or clinical factors that increase UCE risk are known.This is paramount and beneficial to this study as the ML model applied scores, without assistance, to these variables and agree with changes in patient factors that would similarly concern anesthesiologists regarding patient clinical status.Knowing that the ML model views clinical patient factors similar to that found in clinical practice, should improve trust between anesthesiologists and AI and should be viewed as augmentation of our discipline's clinical skillset [47][48][49].For example, many of the diagnoses-emergency, increased heart rate, respiratory rate, AKI, etc.-are known to cause perioperative mortality and morbidity for patients [6][7][8][9][10].However, most studies examine these factors as individual entities rather than in tandem.It would be difficult for anesthesiologists to evaluate the number of variables that the ML model evaluates for every patient discharged from the PACU at large medical centers.However, a ML algorithm can do this for every patient in the PACU quickly and alert anesthesiologists to patients that meet predetermined UCE risk levels and indicate factors causing the ML algorithm to attribute risk.Identification of patient UCE risk will augment anesthesiologist clinical practice as they can examine select patient charts for specific risk elements and make informed decisions about the most appropriate patient care levels (inpatient floor, intermediate care, or ICU).Furthermore, this will improve communication with surgical teams to alert them of patient UCE risk factors and concern for closer monitoring.In our study where ~ 5% of patients had an UCE, this became an average of ~ 600 patients/year or 2 patients/day that would require evaluation at our specific institution.The potential benefit to the patient should outweigh the cost to anesthesiologists, streamline analysis of multiple variables, and accentuate our perioperative practice.
Weaknesses of the study include inclusion of only a single center, no prospective testing of the model, and only non-cardiac surgery patients.While this is a single center study, it employs nearly 60,000 patient records and, as can be seen in Table 1, we have rich diversity among age, race, gender, and co-morbidities of patients, improving the generalizability of our results.Despite no prospective testing of the model, our model did test well after training and was able to find risk factors that are similar to that known by anesthesiologists and the literature without prompting.This is a positive indicator for the algorithm as it is not frequently including obsequious risk factors.Cardiac surgery patients were excluded from our study because nearly 100% of them was developed to help determine when a patient is ready to be discharged from the PACU [40,41].While highly useful to determine readiness for discharge and recovery from anesthesia, the score is not designed to determine illness severity or proper unit disposition.A model that can analyze larger pools of data quickly and attribute UCE risk to surgical patients would detect patients at risk of clinical decompensation and augment anesthesiologist clinical decision making regarding appropriate disposition.Furthermore, these results could lead to a new scoring system that is better able to decide appropriate disposition for post-surgical PACU patients.Other published algorithms also consider patient escalation and predict when patients will do poorly [42][43][44].These studies consider UCE throughout the hospital stay or consider patient events 1 year post-discharge, but are not specific to the surgical patient or the perioperative period.Future studies should focus on comparison with other warning system scores [45,46] for greater generalizability to non-surgical hospital patients.
Our ML algorithm was able to identify specific risk factors for patient UCE (Fig. 3; Table 2).Some of these risks are already known-AKI, maximal respiratory rate, etc.However, two new variables were non-surgical anesthesia time and "missing" data.First, increased non-surgical anesthesia time incorporates both a timeframe between anesthesia start time and surgical procedure start time and another timeframe between the surgical procedure finish time and anesthesia finish time.Importantly, different events occur during each timeframe.Interestingly, increased non-surgical anesthesia time is likely due to patient characteristics, procedural technicalities, and process variables.For example, before the surgical procedure, there may be placement of an arterial line or central venous line, an unanticipated difficult airway, or difficulty in patient positioning.After the procedure, the patient could have a prolonged awakening time after anesthesia and/or co-morbidities, such as age or cardiac/respiratory diagnoses, that prolong extubation.It is currently unclear whether the patient risk factors lead to needed procedures or the procedures and other processes themselves are responsible for increased non-surgical anesthesia time.Second, "missing" data in the last hour before PACU discharge increased UCE risk.The clinical concern is that patient data is not available for patient evaluation before PACU discharge.Whether missing data is normal or abnormal and how this data would be interpreted by the ML model to attribute UCE risk is unclear.Unfortunately, patient decisions were made without this data in the past.However, in future prospective studies, we hope that because our algorithm indicates a higher UCE risk score when critical patient data is missing, a physician or PACU nurse will be alerted, the data can be added to the algorithm, and more accurate evaluations of patient UCE risk can be Standards of Reporting This manuscript adheres to the STROBE guidelines.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made.The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material.If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder.To view a copy of this licence, visit http://creativecommons. org/licenses/by/4.0/.go immediately to an ICU after surgery.Determining risk of re-entry, post-operative bleeding, etc. in cardiac surgery patients will require a different analysis and algorithm.These three limitations are the goals of future studies.
In conclusion, this ML algorithm can be used to predict patients that are at increased UCE risk.While it cannot absolutely tell anesthesiologists where to place patients within the hospital, it does provide us with increased data to predict patients that will have a problem within the next three midnights.This data will be useful as it allows us to closely evaluate these patients prior to PACU discharge, speak with their surgical teams, and provide interventions that may reduce morbidity for these patients during their hospital stay.Future studies are needed to investigate: (1) prospective use of the model; (2) testing against an external cohort; and (3) the possibility of a new scoring system to better evaluate PACU discharge and appropriate care level.These findings are exciting as they use ML to augment anesthesiologist decision-making, further our role in precision medicine and the perioperative period, and provide the safest possible care of our patients.

Fig. 2 Fig. 1
Fig. 2 ML algorithm learns and can predict UCE in post-PACU discharge patients within three midnights.(A) Testing and (B) training receiver-operator curves are displayed.For both figures, what a random ML model would produce is shown with a dashed line.The AUC for training (A) is 0.757 and testing (B) is 0.752, both of which are acceptable values for ML Fig. 1 Study Design Flow Chart.The total number of patients assessed, excluded, included, and allocated to the different arms of the study are shown.Furthermore, the number of patients within in each arm that had UCE are indicated

Fig.
Fig.3 Shapley values that most affect patient UCE risk.On the left, the 10 most influential variables for patients whose estimated probability of UCE was greater than 20% are listed.Each patient is visualized as a path moving from lower risk (left-hand side) to higher risk (righthand side) as variables of increasing importance (moving from bottom

Table 1
Body Mass Index; CKD: Chronic Kidney Disease; AKI: Acute Kidney Injury; ESRD: End Stage Renal Disease; COPD: Chronic Obstructive Pulmonary Disease; MACE: Major Adverse Cardiac Events; ASA Status

Table 2
Breakdown of example bins with percentage of events within each bin and points scored by machine learning algorithm for being within HR: Heart Rate; PACU: Post−Anesthesia Care Unit

Table 3
Individual variables that provide the highest risk points in the ML model These values are measured in the final hour of the PACU stay.HR: Heart Rate; PACU: Post−Anesthesia Care Unit; ASA: American Society of Anesthesiologists; SBP: Systolic Blood Pressure; RR: Respiratory Rate