Racial Disparities in Pulse Oximetry, in COVID-19 and ICU Settings.

OBJECTIVES BACKGROUND: This study aimed to assess the impact of race on pulse oximetry reliability, taking into account Spo2 ranges, COVID-19 diagnosis, and ICU admission. DESIGN: Retrospective cohort study covering admissions from January 2020 to April 2024. SETTING: National COVID Cohort Collaborative (N3C) database, consisting of electronic health records from 80 U.S. institutions. PATIENTS/SUBJECTS: Patients were selected from the N3C database based on the availability of data on self-identified race and both pulse oximetry estimated Spo2 and Sao2. Subgroups included patients in ICU and non-ICU settings, with or without a diagnosis of COVID-19 disease. INTERVENTIONS: None. MEASUREMENTS AND MAIN RESULTS: The agreement between Spo2 and Sao2 was assessed across racial groups (American Indian or Alaska Native, Asian, Black, Hispanic or Latino, Pacific Islander, and White). Each patient's initial Sao2 measurement was matched with the closest Spo2 values recorded within the preceding 10-minute time frame. The risk of hidden hypoxemia (Spo2 ≥ 88% but Sao2 < 88%) was determined for various Spo2 ranges, races, and clinical scenarios. We used a generalized logistic mixed-effects model to evaluate the impact of relevant variables, such as COVID-19, ICU admission, age, sex, race, and Spo2, on the risk of hidden hypoxemia, while accounting for the random effects within each hospital. A total of 80,541 patients were included, consisting of 596 American Indian or Alaska Native, 2,729 Asian, 11,889 Black, 13,154 Hispanic or Latino, 221 Pacific Islander, and 51,952 White individuals. Discrepancies between Spo2 and Sao2 were observed across all racial groups, with the most pronounced bias in Black patients. Hidden hypoxemia rates were higher in Black patients across all Spo2 subgroups, for all clinical scenarios. The odds of hidden hypoxemia were higher for Black and Hispanic or Latino patients and for those with COVID-19 disease. CONCLUSIONS: Race significantly impacts pulse oximetry reliability. Not only Black and Hispanic or Latino patients were at higher risk for hidden hypoxemia, but also those admitted with a COVID-19 diagnosis. Future in-depth explorations into the underlying causes and potential solutions are needed.

Augmented intelligence facilitates concept mapping across different electronic health records.

INTRODUCTION: With the advent of artificial intelligence, the secondary use of routinely collected medical data from electronic healthcare records (EHR) has become increasingly popular. However, different EHR systems typically use different names for the same medical concepts. This obviously hampers scalable model development and subsequent clinical implementation for decision support. Therefore, converting original parameter names to a so-called ontology, a standardized set of predefined concepts, is necessary but time-consuming and labor-intensive. We therefore propose an augmented intelligence approach to facilitate ontology alignment by predicting correct concepts based on parameter names from raw electronic health record data exports. METHODS: We used the manually mapped parameter names from the multicenter "Dutch ICU data warehouse against COVID-19" sourced from three types of EHR systems to train machine learning models for concept mapping. Data from 29 intensive care units on 38,824 parameters mapped to 1,679 relevant and unique concepts and 38,069 parameters labeled as irrelevant were used for model development and validation. We used the Natural Language Toolkit (NLTK) to preprocess the parameter names based on WordNet cognitive synonyms transformed by term-frequency inverse document frequency (TF-IDF), yielding numeric features. We then trained linear classifiers using stochastic gradient descent for multi-class prediction. Finally, we fine-tuned these predictions using information on distributions of the data associated with each parameter name through similarity score and skewness comparisons. RESULTS: The initial model, trained using data from one hospital organization for each of three EHR systems, scored an overall top 1 precision of 0.744, recall of 0.771, and F1-score of 0.737 on a total of 58,804 parameters. Leave-one-hospital-out analysis returned an average top 1 recall of 0.680 for relevant parameters, which increased to 0.905 for the top 5 predictions. When reducing the training dataset to only include relevant parameters, top 1 recall was 0.811 and top 5 recall was 0.914 for relevant parameters. Performance improvement based on similarity score or skewness comparisons affected at most 5.23% of numeric parameters. CONCLUSION: Augmented intelligence is a promising method to improve concept mapping of parameter names from raw electronic health record data exports. We propose a robust method for mapping data across various domains, facilitating the integration of diverse data sources. However, recall is not perfect, and therefore manual validation of mapping remains essential.

Determining and assessing characteristics of data element names impacting the performance of annotation using Usagi.

INTRODUCTION: Hospitals generate large amounts of data and this data is generally modeled and labeled in a proprietary way, hampering its exchange and integration. Manually annotating data element names to internationally standardized data element identifiers is a time-consuming effort. Tools can support performing this task automatically. This study aimed to determine what factors influence the quality of automatic annotations. METHODS: Data element names were used from the Dutch COVID-19 ICU Data Warehouse containing data on intensive care patients with COVID-19 from 25 hospitals in the Netherlands. In this data warehouse, the data had been merged using a proprietary terminology system while also storing the original hospital labels (synonymous names). Usagi, an OHDSI annotation tool, was used to perform the annotation for the data. A gold standard was used to determine if Usagi made correct annotations. Logistic regression was used to determine if the number of characters, number of words, match score (Usagi's certainty) and hospital label origin influenced Usagi's performance to annotate correctly. RESULTS: Usagi automatically annotated 30.5% of the data element names correctly and 5.5% of the synonymous names. The match score is the best predictor for Usagi finding the correct annotation. It was determined that the AUC of data element names was 0.651 and 0.752 for the synonymous names respectively. The AUC for the individual hospital label origins varied between 0.460 to 0.905. DISCUSSION: The results show that Usagi performed better to annotate the data element names than the synonymous names. The hospital origin in the synonymous names dataset was associated with the amount of correctly annotated concepts. Hospitals that performed better had shorter synonymous names and fewer words. Using shorter data element names or synonymous names should be considered to optimize the automatic annotating process. Overall, the performance of Usagi is too poor to completely rely on for automatic annotation.

Age-related changes in plasma biomarkers and their association with mortality in COVID-19.

BACKGROUND: Coronavirus disease 2019 (COVID-19)-induced mortality occurs predominantly in older patients. Several immunomodulating therapies seem less beneficial in these patients. The biological substrate behind these observations is unknown. The aim of this study was to obtain insight into the association between ageing, the host response and mortality in patients with COVID-19. METHODS: We determined 43 biomarkers reflective of alterations in four pathophysiological domains: endothelial cell and coagulation activation, inflammation and organ damage, and cytokine and chemokine release. We used mediation analysis to associate ageing-driven alterations in the host response with 30-day mortality. Biomarkers associated with both ageing and mortality were validated in an intensive care unit and external cohort. RESULTS: 464 general ward patients with COVID-19 were stratified according to age decades. Increasing age was an independent risk factor for 30-day mortality. Ageing was associated with alterations in each of the host response domains, characterised by greater activation of the endothelium and coagulation system and stronger elevation of inflammation and organ damage markers, which was independent of an increase in age-related comorbidities. Soluble tumour necrosis factor receptor 1, soluble triggering receptor expressed on myeloid cells 1 and soluble thrombomodulin showed the strongest correlation with ageing and explained part of the ageing-driven increase in 30-day mortality (proportion mediated: 13.0%, 12.9% and 12.6%, respectively). CONCLUSIONS: Ageing is associated with a strong and broad modification of the host response to COVID-19, and specific immune changes likely contribute to increased mortality in older patients. These results may provide insight into potential age-specific immunomodulatory targets in COVID-19.

Evolution of Clinical Phenotypes of COVID-19 Patients During Intensive Care Treatment: An Unsupervised Machine Learning Analysis.

BACKGROUND: Identification of clinical phenotypes in critically ill COVID-19 patients could improve understanding of the disease heterogeneity and enable prognostic and predictive enrichment. However, previous attempts did not take into account temporal dynamics with high granularity. By including the dimension of time, we aim to gain further insights into the heterogeneity of COVID-19. METHODS: We used granular data from 3202 adult COVID patients in the Dutch Data Warehouse that were admitted to one of 25 Dutch ICUs between February 2020 and March 2021. Parameters including demographics, clinical observations, medications, laboratory values, vital signs, and data from life support devices were selected. Twenty-one datasets were created that each covered 24 h of ICU data for each day of ICU treatment. Clinical phenotypes in each dataset were identified by performing cluster analyses. Both evolution of the clinical phenotypes over time and patient allocation to these clusters over time were tracked. RESULTS: The final patient cohort consisted of 2438 COVID-19 patients with a ICU mortality outcome. Forty-one parameters were chosen for cluster analysis. On admission, both a mild and a severe clinical phenotype were found. After day 4, the severe phenotype split into an intermediate and a severe phenotype for 11 consecutive days. Heterogeneity between phenotypes appears to be driven by inflammation and dead space ventilation. During the 21-day period, only 8.2% and 4.6% of patients in the initial mild and severe clusters remained assigned to the same phenotype respectively. The clinical phenotype half-life was between 5 and 6 days for the mild and severe phenotypes, and about 3 days for the medium severe phenotype. CONCLUSIONS: Patients typically do not remain in the same cluster throughout intensive care treatment. This may have important implications for prognostic or predictive enrichment. Prominent dissimilarities between clinical phenotypes are predominantly driven by inflammation and dead space ventilation.

Predicting responders to prone positioning in mechanically ventilated patients with COVID-19 using machine learning.

BACKGROUND: For mechanically ventilated critically ill COVID-19 patients, prone positioning has quickly become an important treatment strategy, however, prone positioning is labor intensive and comes with potential adverse effects. Therefore, identifying which critically ill intubated COVID-19 patients will benefit may help allocate labor resources. METHODS: From the multi-center Dutch Data Warehouse of COVID-19 ICU patients from 25 hospitals, we selected all 3619 episodes of prone positioning in 1142 invasively mechanically ventilated patients. We excluded episodes longer than 24 h. Berlin ARDS criteria were not formally documented. We used supervised machine learning algorithms Logistic Regression, Random Forest, Naive Bayes, K-Nearest Neighbors, Support Vector Machine and Extreme Gradient Boosting on readily available and clinically relevant features to predict success of prone positioning after 4 h (window of 1 to 7 h) based on various possible outcomes. These outcomes were defined as improvements of at least 10% in PaO2/FiO2 ratio, ventilatory ratio, respiratory system compliance, or mechanical power. Separate models were created for each of these outcomes. Re-supination within 4 h after pronation was labeled as failure. We also developed models using a 20 mmHg improvement cut-off for PaO2/FiO2 ratio and using a combined outcome parameter. For all models, we evaluated feature importance expressed as contribution to predictive performance based on their relative ranking. RESULTS: The median duration of prone episodes was 17 h (11-20, median and IQR, N = 2632). Despite extensive modeling using a plethora of machine learning techniques and a large number of potentially clinically relevant features, discrimination between responders and non-responders remained poor with an area under the receiver operator characteristic curve of 0.62 for PaO2/FiO2 ratio using Logistic Regression, Random Forest and XGBoost. Feature importance was inconsistent between models for different outcomes. Notably, not even being a previous responder to prone positioning, or PEEP-levels before prone positioning, provided any meaningful contribution to predicting a successful next proning episode. CONCLUSIONS: In mechanically ventilated COVID-19 patients, predicting the success of prone positioning using clinically relevant and readily available parameters from electronic health records is currently not feasible. Given the current evidence base, a liberal approach to proning in all patients with severe COVID-19 ARDS is therefore justified and in particular regardless of previous results of proning.

Assess and validate predictive performance of models for in-hospital mortality in COVID-19 patients: A retrospective cohort study in the Netherlands comparing the value of registry data with high-granular electronic health records.

PURPOSE: To assess, validate and compare the predictive performance of models for in-hospital mortality of COVID-19 patients admitted to the intensive care unit (ICU) over two different waves of infections. Our models were built with high-granular Electronic Health Records (EHR) data versus less-granular registry data. METHODS: Observational study of all COVID-19 patients admitted to 19 Dutch ICUs participating in both the national quality registry National Intensive Care Evaluation (NICE) and the EHR-based Dutch Data Warehouse (hereafter EHR). Multiple models were developed on data from the first 24 h of ICU admissions from February to June 2020 (first COVID-19 wave) and validated on prospective patients admitted to the same ICUs between July and December 2020 (second COVID-19 wave). We assessed model discrimination, calibration, and the degree of relatedness between development and validation population. Coefficients were used to identify relevant risk factors. RESULTS: A total of 1533 patients from the EHR and 1563 from the registry were included. With high granular EHR data, the average AUROC was 0.69 (standard deviation of 0.05) for the internal validation, and the AUROC was 0.75 for the temporal validation. The registry model achieved an average AUROC of 0.76 (standard deviation of 0.05) in the internal validation and 0.77 in the temporal validation. In the EHR data, age, and respiratory-system related variables were the most important risk factors identified. In the NICE registry data, age and chronic respiratory insufficiency were the most important risk factors. CONCLUSION: In our study, prognostic models built on less-granular but readily-available registry data had similar performance to models built on high-granular EHR data and showed similar transportability to a prospective COVID-19 population. Future research is needed to verify whether this finding can be confirmed for upcoming waves.

Right dose, right now: bedside, real-time, data-driven, and personalised antibiotic dosing in critically ill patients with sepsis or septic shock-a two-centre randomised clinical trial.

BACKGROUND: Adequate antibiotic dosing may improve outcomes in critically ill patients but is challenging due to altered and variable pharmacokinetics. To address this challenge, AutoKinetics was developed, a decision support system for bedside, real-time, data-driven and personalised antibiotic dosing. This study evaluates the feasibility, safety and efficacy of its clinical implementation. METHODS: In this two-centre randomised clinical trial, critically ill patients with sepsis or septic shock were randomised to AutoKinetics dosing or standard dosing for four antibiotics: vancomycin, ciprofloxacin, meropenem, and ceftriaxone. Adult patients with a confirmed or suspected infection and either lactate > 2 mmol/L or vasopressor requirement were eligible for inclusion. The primary outcome was pharmacokinetic target attainment in the first 24 h after randomisation. Clinical endpoints included mortality, ICU length of stay and incidence of acute kidney injury. RESULTS: After inclusion of 252 patients, the study was stopped early due to the COVID-19 pandemic. In the ciprofloxacin intervention group, the primary outcome was obtained in 69% compared to 3% in the control group (OR 62.5, CI 11.4-1173.78, p < 0.001). Furthermore, target attainment was faster (26 h, CI 18-42 h, p < 0.001) and better (65% increase, CI 49-84%, p < 0.001). For the other antibiotics, AutoKinetics dosing did not improve target attainment. Clinical endpoints were not significantly different. Importantly, higher dosing did not lead to increased mortality or renal failure. CONCLUSIONS: In critically ill patients, personalised dosing was feasible, safe and significantly improved target attainment for ciprofloxacin. TRIAL REGISTRATION: The trial was prospectively registered at Netherlands Trial Register (NTR), NL6501/NTR6689 on 25 August 2017 and at the European Clinical Trials Database (EudraCT), 2017-002478-37 on 6 November 2017.

INCIDENCE, RISK FACTORS, AND OUTCOME OF SUSPECTED CENTRAL VENOUS CATHETER-RELATED INFECTIONS IN CRITICALLY ILL COVID-19 PATIENTS: A MULTICENTER RETROSPECTIVE COHORT STUDY.

ABSTRACT: Background: Aims of this study were to investigate the prevalence and incidence of catheter-related infection, identify risk factors, and determine the relation of catheter-related infection with mortality in critically ill COVID-19 patients. Methods: This was a retrospective cohort study of central venous catheters (CVCs) in critically ill COVID-19 patients. Eligible CVC insertions required an indwelling time of at least 48 hours and were identified using a full-admission electronic health record database. Risk factors were identified using logistic regression. Differences in survival rates at day 28 of follow-up were assessed using a log-rank test and proportional hazard model. Results: In 538 patients, a total of 914 CVCs were included. Prevalence and incidence of suspected catheter-related infection were 7.9% and 9.4 infections per 1,000 catheter indwelling days, respectively. Prone ventilation for more than 5 days was associated with increased risk of suspected catheter-related infection; odds ratio, 5.05 (95% confidence interval 2.12-11.0). Risk of death was significantly higher in patients with suspected catheter-related infection (hazard ratio, 1.78; 95% confidence interval, 1.25-2.53). Conclusions: This study shows that in critically ill patients with COVID-19, prevalence and incidence of suspected catheter-related infection are high, prone ventilation is a risk factor, and mortality is higher in case of catheter-related infection.

Rapid Evaluation of Coronavirus Illness Severity (RECOILS) in intensive care: Development and validation of a prognostic tool for in-hospital mortality.

BACKGROUND: The prediction of in-hospital mortality for ICU patients with COVID-19 is fundamental to treatment and resource allocation. The main purpose was to develop an easily implemented score for such prediction. METHODS: This was an observational, multicenter, development, and validation study on a national critical care dataset of COVID-19 patients. A systematic literature review was performed to determine variables possibly important for COVID-19 mortality prediction. Using a logistic multivariable model with a LASSO penalty, we developed the Rapid Evaluation of Coronavirus Illness Severity (RECOILS) score and compared its performance against published scores. RESULTS: Our development (validation) cohort consisted of 1480 (937) adult patients from 14 (11) Dutch ICUs admitted between March 2020 and April 2021. Median age was 65 (65) years, 31% (26%) died in hospital, 74% (72%) were males, average length of ICU stay was 7.83 (10.25) days and average length of hospital stay was 15.90 (19.92) days. Age, platelets, PaO2/FiO2 ratio, pH, blood urea nitrogen, temperature, PaCO2, Glasgow Coma Scale (GCS) score measured within +/-24 h of ICU admission were used to develop the score. The AUROC of RECOILS score was 0.75 (CI 0.71-0.78) which was higher than that of any previously reported predictive scores (0.68 [CI 0.64-0.71], 0.61 [CI 0.58-0.66], 0.67 [CI 0.63-0.70], 0.70 [CI 0.67-0.74] for ISARIC 4C Mortality Score, SOFA, SAPS-III, and age, respectively). CONCLUSIONS: Using a large dataset from multiple Dutch ICUs, we developed a predictive score for mortality of COVID-19 patients admitted to ICU, which outperformed other predictive scores reported so far.

Predictors for extubation failure in COVID-19 patients using a machine learning approach.

INTRODUCTION: Determining the optimal timing for extubation can be challenging in the intensive care. In this study, we aim to identify predictors for extubation failure in critically ill patients with COVID-19. METHODS: We used highly granular data from 3464 adult critically ill COVID patients in the multicenter Dutch Data Warehouse, including demographics, clinical observations, medications, fluid balance, laboratory values, vital signs, and data from life support devices. All intubated patients with at least one extubation attempt were eligible for analysis. Transferred patients, patients admitted for less than 24 h, and patients still admitted at the time of data extraction were excluded. Potential predictors were selected by a team of intensive care physicians. The primary and secondary outcomes were extubation without reintubation or death within the next 7 days and within 48 h, respectively. We trained and validated multiple machine learning algorithms using fivefold nested cross-validation. Predictor importance was estimated using Shapley additive explanations, while cutoff values for the relative probability of failed extubation were estimated through partial dependence plots. RESULTS: A total of 883 patients were included in the model derivation. The reintubation rate was 13.4% within 48 h and 18.9% at day 7, with a mortality rate of 0.6% and 1.0% respectively. The grandient-boost model performed best (area under the curve of 0.70) and was used to calculate predictor importance. Ventilatory characteristics and settings were the most important predictors. More specifically, a controlled mode duration longer than 4 days, a last fraction of inspired oxygen higher than 35%, a mean tidal volume per kg ideal body weight above 8 ml/kg in the day before extubation, and a shorter duration in assisted mode (< 2 days) compared to their median values. Additionally, a higher C-reactive protein and leukocyte count, a lower thrombocyte count, a lower Glasgow coma scale and a lower body mass index compared to their medians were associated with extubation failure. CONCLUSION: The most important predictors for extubation failure in critically ill COVID-19 patients include ventilatory settings, inflammatory parameters, neurological status, and body mass index. These predictors should therefore be routinely captured in electronic health records.

Some Patients Are More Equal Than Others: Variation in Ventilator Settings for Coronavirus Disease 2019 Acute Respiratory Distress Syndrome.

OBJECTIVES: As coronavirus disease 2019 is a novel disease, treatment strategies continue to be debated. This provides the intensive care community with a unique opportunity as the population of coronavirus disease 2019 patients requiring invasive mechanical ventilation is relatively homogeneous compared with other ICU populations. We hypothesize that the novelty of coronavirus disease 2019 and the uncertainty over its similarity with noncoronavirus disease 2019 acute respiratory distress syndrome resulted in substantial practice variation between hospitals during the first and second waves of coronavirus disease 2019 patients. DESIGN: Multicenter retrospective cohort study. SETTING: Twenty-five hospitals in the Netherlands from February 2020 to July 2020, and 14 hospitals from August 2020 to December 2020. PATIENTS: One thousand two hundred ninety-four critically ill intubated adult ICU patients with coronavirus disease 2019 were selected from the Dutch Data Warehouse. Patients intubated for less than 24 hours, transferred patients, and patients still admitted at the time of data extraction were excluded. MEASUREMENTS AND MAIN RESULTS: We aimed to estimate between-ICU practice variation in selected ventilation parameters (positive end-expiratory pressure, Fio2, set respiratory rate, tidal volume, minute volume, and percentage of time spent in a prone position) on days 1, 2, 3, and 7 of intubation, adjusted for patient characteristics as well as severity of illness based on Pao2/Fio2 ratio, pH, ventilatory ratio, and dynamic respiratory system compliance during controlled ventilation. Using multilevel linear mixed-effects modeling, we found significant (p ≤ 0.001) variation between ICUs in all ventilation parameters on days 1, 2, 3, and 7 of intubation for both waves. CONCLUSIONS: This is the first study to clearly demonstrate significant practice variation between ICUs related to mechanical ventilation parameters that are under direct control by intensivists. Their effect on clinical outcomes for both coronavirus disease 2019 and other critically ill mechanically ventilated patients could have widespread implications for the practice of intensive care medicine and should be investigated further by causal inference models and clinical trials.

The Dutch Data Warehouse, a multicenter and full-admission electronic health records database for critically ill COVID-19 patients.

BACKGROUND: The Coronavirus disease 2019 (COVID-19) pandemic has underlined the urgent need for reliable, multicenter, and full-admission intensive care data to advance our understanding of the course of the disease and investigate potential treatment strategies. In this study, we present the Dutch Data Warehouse (DDW), the first multicenter electronic health record (EHR) database with full-admission data from critically ill COVID-19 patients. METHODS: A nation-wide data sharing collaboration was launched at the beginning of the pandemic in March 2020. All hospitals in the Netherlands were asked to participate and share pseudonymized EHR data from adult critically ill COVID-19 patients. Data included patient demographics, clinical observations, administered medication, laboratory determinations, and data from vital sign monitors and life support devices. Data sharing agreements were signed with participating hospitals before any data transfers took place. Data were extracted from the local EHRs with prespecified queries and combined into a staging dataset through an extract-transform-load (ETL) pipeline. In the consecutive processing pipeline, data were mapped to a common concept vocabulary and enriched with derived concepts. Data validation was a continuous process throughout the project. All participating hospitals have access to the DDW. Within legal and ethical boundaries, data are available to clinicians and researchers. RESULTS: Out of the 81 intensive care units in the Netherlands, 66 participated in the collaboration, 47 have signed the data sharing agreement, and 35 have shared their data. Data from 25 hospitals have passed through the ETL and processing pipeline. Currently, 3464 patients are included in the DDW, both from wave 1 and wave 2 in the Netherlands. More than 200 million clinical data points are available. Overall ICU mortality was 24.4%. Respiratory and hemodynamic parameters were most frequently measured throughout a patient's stay. For each patient, all administered medication and their daily fluid balance were available. Missing data are reported for each descriptive. CONCLUSIONS: In this study, we show that EHR data from critically ill COVID-19 patients may be lawfully collected and can be combined into a data warehouse. These initiatives are indispensable to advance medical data science in the field of intensive care medicine.

Predicting mortality of individual patients with COVID-19: a multicentre Dutch cohort.

OBJECTIVE: Develop and validate models that predict mortality of patients diagnosed with COVID-19 admitted to the hospital. DESIGN: Retrospective cohort study. SETTING: A multicentre cohort across 10 Dutch hospitals including patients from 27 February to 8 June 2020. PARTICIPANTS: SARS-CoV-2 positive patients (age ≥18) admitted to the hospital. MAIN OUTCOME MEASURES: 21-day all-cause mortality evaluated by the area under the receiver operator curve (AUC), sensitivity, specificity, positive predictive value and negative predictive value. The predictive value of age was explored by comparison with age-based rules used in practice and by excluding age from the analysis. RESULTS: 2273 patients were included, of whom 516 had died or discharged to palliative care within 21 days after admission. Five feature sets, including premorbid, clinical presentation and laboratory and radiology values, were derived from 80 features. Additionally, an Analysis of Variance (ANOVA)-based data-driven feature selection selected the 10 features with the highest F values: age, number of home medications, urea nitrogen, lactate dehydrogenase, albumin, oxygen saturation (%), oxygen saturation is measured on room air, oxygen saturation is measured on oxygen therapy, blood gas pH and history of chronic cardiac disease. A linear logistic regression and non-linear tree-based gradient boosting algorithm fitted the data with an AUC of 0.81 (95% CI 0.77 to 0.85) and 0.82 (0.79 to 0.85), respectively, using the 10 selected features. Both models outperformed age-based decision rules used in practice (AUC of 0.69, 0.65 to 0.74 for age >70). Furthermore, performance remained stable when excluding age as predictor (AUC of 0.78, 0.75 to 0.81). CONCLUSION: Both models showed good performance and had better test characteristics than age-based decision rules, using 10 admission features readily available in Dutch hospitals. The models hold promise to aid decision-making during a hospital bed shortage.

Risk factors for adverse outcomes during mechanical ventilation of 1152 COVID-19 patients: a multicenter machine learning study with highly granular data from the Dutch Data Warehouse.

BACKGROUND: The identification of risk factors for adverse outcomes and prolonged intensive care unit (ICU) stay in COVID-19 patients is essential for prognostication, determining treatment intensity, and resource allocation. Previous studies have determined risk factors on admission only, and included a limited number of predictors. Therefore, using data from the highly granular and multicenter Dutch Data Warehouse, we developed machine learning models to identify risk factors for ICU mortality, ventilator-free days and ICU-free days during the course of invasive mechanical ventilation (IMV) in COVID-19 patients. METHODS: The DDW is a growing electronic health record database of critically ill COVID-19 patients in the Netherlands. All adult ICU patients on IMV were eligible for inclusion. Transfers, patients admitted for less than 24 h, and patients still admitted at time of data extraction were excluded. Predictors were selected based on the literature, and included medication dosage and fluid balance. Multiple algorithms were trained and validated on up to three sets of observations per patient on day 1, 7, and 14 using fivefold nested cross-validation, keeping observations from an individual patient in the same split. RESULTS: A total of 1152 patients were included in the model. XGBoost models performed best for all outcomes and were used to calculate predictor importance. Using Shapley additive explanations (SHAP), age was the most important demographic risk factor for the outcomes upon start of IMV and throughout its course. The relative probability of death across age values is visualized in Partial Dependence Plots (PDPs), with an increase starting at 54 years. Besides age, acidaemia, low P/F-ratios and high driving pressures demonstrated a higher probability of death. The PDP for driving pressure showed a relative probability increase starting at 12 cmH2O. CONCLUSION: Age is the most important demographic risk factor of ICU mortality, ICU-free days and ventilator-free days throughout the course of invasive mechanical ventilation in critically ill COVID-19 patients. pH, P/F ratio, and driving pressure should be monitored closely over the course of mechanical ventilation as risk factors predictive of these outcomes.

[Clinical course of COVID-19 in the Netherlands: an overview of 2607 patients in hospital during the first wave]. / Klinisch beloop van covid-19 in Nederland.

OBJECTIVE: To systematically collect clinical data from patients with a proven COVID-19 infection in the Netherlands. DESIGN: Data from 2579 patients with COVID-19 admitted to 10 Dutch centers in the period February to July 2020 are described. The clinical data are based on the WHO COVID case record form (CRF) and supplemented with patient characteristics of which recently an association disease severity has been reported. METHODS: Survival analyses were performed as primary statistical analysis. These Kaplan-Meier curves for time to (early) death (3 weeks) have been determined for pre-morbid patient characteristics and clinical, radiological and laboratory data at hospital admission. RESULTS: Total in-hospital mortality after 3 weeks was 22.2% (95% CI: 20.7% - 23.9%), hospital mortality within 21 days was significantly higher for elderly patients (> 70 years; 35, 0% (95% CI: 32.4% - 37.8%) and patients who died during the 21 days and were admitted to the intensive care (36.5% (95% CI: 32.1% - 41.3%)). Apart from that, in this Dutch population we also see a risk of early death in patients with co-morbidities (such as chronic neurological, nephrological and cardiac disorders and hypertension), and in patients with more home medication and / or with increased urea and creatinine levels. CONCLUSION: Early death due to a COVID-19 infection in the Netherlands appears to be associated with demographic variables (e.g. age), comorbidity (e.g. cardiovascular disease) but also disease char-acteristics at admission.

Transatlantic transferability of a new reinforcement learning model for optimizing haemodynamic treatment for critically ill patients with sepsis.

INTRODUCTION: In recent years, reinforcement learning (RL) has gained traction in the healthcare domain. In particular, RL methods have been explored for haemodynamic optimization of septic patients in the Intensive Care Unit. Most hospitals however, lack the data and expertise for model development, necessitating transfer of models developed using external datasets. This approach assumes model generalizability across different patient populations, the validity of which has not previously been tested. In addition, there is limited knowledge on safety and reliability. These challenges need to be addressed to further facilitate implementation of RL models in clinical practice. METHOD: We developed and validated a new reinforcement learning model for hemodynamic optimization in sepsis on the MIMIC intensive care database from the USA using a dueling double deep Q network. We then transferred this model to the European AmsterdamUMCdb intensive care database. T-Distributed Stochastic Neighbor Embedding and Sequential Organ Failure Assessment scores were used to explore the differences between the patient populations. We apply off-policy policy evaluation methods to quantify model performance. In addition, we introduce and apply a novel deep policy inspection to analyse how the optimal policy relates to the different phases of sepsis and sepsis treatment to provide interpretable insight in order to assess model safety and reliability. RESULTS: The off-policy evaluation revealed that the optimal policy outperformed the physician policy on both datasets despite marked differences between the two patient populations and physician's policies. Our novel deep policy inspection method showed insightful results and unveiled that the model could initiate therapy adequately and adjust therapy intensity to illness severity and disease progression which indicated safe and reliable model behaviour. Compared to current physician behavior, the developed policy prefers a more liberal use of vasopressors with a more restrained use of fluid therapy in line with previous work. CONCLUSION: We created a reinforcement learning model for optimal bedside hemodynamic management and demonstrated model transferability between populations from the USA and Europe for the first time. We proposed new methods for deep policy inspection integrating expert domain knowledge. This is expected to facilitate progression to bedside clinical decision support for the treatment of critically ill patients.

Why we should sample sparsely and aim for a higher target: Lessons from model-based therapeutic drug monitoring of vancomycin in intensive care patients.

AIMS: To explore the optimal data sampling scheme and the pharmacokinetic (PK) target exposure on which dose computation is based in the model-based therapeutic drug monitoring (TDM) practice of vancomycin in intensive care (ICU) patients. METHODS: We simulated concentration data for 1 day following four sampling schemes, Cmin , Cmax + Cmin , Cmax + Cmid-interval + Cmin , and rich sampling where a sample was drawn every hour within a dose interval. The datasets were used for Bayesian estimation to obtain PK parameters, which were used to compute the doses for the next day based on five PK target exposures: AUC24 = 400, 500, and 600 mg·h/L and Cmin = 15 and 20 mg/L. We then simulated data for the next day, adopting the computed doses, and repeated the above procedure for 7 days. Thereafter, we calculated the percentage error and the normalized root mean square error (NRMSE) of estimated against "true" PK parameters, and the percentage of optimal treatment (POT), defined as the percentage of patients who met 400 ≤ AUC24 ≤ 600 mg·h/L and Cmin ≤ 20 mg/L. RESULTS: PK parameters were unbiasedly estimated in all investigated scenarios and the 6-day average NRMSE were 32.5%/38.5% (CL/V, where CL is clearance and V is volume of distribution) in the trough sampling scheme and 27.3%/26.5% (CL/V) in the rich sampling scheme. Regarding POT, the sampling scheme had marginal influence, while target exposure showed clear impacts that the maximum POT of 71.5% was reached when doses were computed based on AUC24 = 500 mg·h/L. CONCLUSIONS: For model-based TDM of vancomycin in ICU patients, sampling more frequently than taking only trough samples adds no value and dosing based on AUC24 = 500 mg·h/L lead to the best POT.

Optimizing Predictive Performance of Bayesian Forecasting for Vancomycin Concentration in Intensive Care Patients.

PURPOSE: Bayesian forecasting is crucial for model-based dose optimization based on therapeutic drug monitoring (TDM) data of vancomycin in intensive care (ICU) patients. We aimed to evaluate the performance of Bayesian forecasting using maximum a posteriori (MAP) estimation for model-based TDM. METHODS: We used a vancomycin TDM data set (n = 408 patients). We compared standard MAP-based Bayesian forecasting with two alternative approaches: (i) adaptive MAP which handles data over multiple iterations, and (ii) weighted MAP which weights the likelihood contribution of data. We evaluated the percentage error (PE) for seven scenarios including historical TDM data from the preceding day up to seven days. RESULTS: The mean of median PEs of all scenarios for the standard MAP, adaptive MAP and weighted MAP method were - 7.7%, -4.5% and - 6.7%. The adaptive MAP also showed the narrowest inter-quartile range of PE. In addition, regardless of MAP method, including historical TDM data further in the past will increase prediction errors. CONCLUSIONS: The proposed adaptive MAP method outperforms standard MAP in predictive performance and may be considered for improvement of model-based dose optimization. The inclusion of historical data beyond either one day (standard MAP and weighted MAP) or two days (adaptive MAP) reduces predictive performance.