This paper is in the following e-collection/theme issue:

Reviews in Medical Informatics (93) Decision Support for Health Professionals (1210) Machine Learning (1501)

Published on in Vol 10, No 5 (2022): May

Preprints (earlier versions) of this paper are available at https://preprints.jmir.org/preprint/35293, first published .
Comparison of Severity of Illness Scores and Artificial Intelligence Models That Are Predictive of Intensive Care Unit Mortality: Meta-analysis and Review of the Literature

Comparison of Severity of Illness Scores and Artificial Intelligence Models That Are Predictive of Intensive Care Unit Mortality: Meta-analysis and Review of the Literature

Comparison of Severity of Illness Scores and Artificial Intelligence Models That Are Predictive of Intensive Care Unit Mortality: Meta-analysis and Review of the Literature

Authors of this article:

Cristina Barboi1 Author Orcid Image ;   Andreas Tzavelis2, 3 Author Orcid Image ;   Lutfiyya NaQiyba Muhammad4 Author Orcid Image
Article Authors Cited by (14) Tweetations (7) Metrics

    Review

    1Indiana University Purdue University, Regenstrief Institue, Indianapolis, IN, United States

    2Medical Scientist Training Program, Feinberg School of Medicine, Chicago, IL, United States

    3Department of Biomedical Engineering, Northwestern University, Chicago, IL, United States

    4Department of Preventive Medicine and Biostatistics, Northwestern University, Evanston, IL, United States

    Corresponding Author:

    Cristina Barboi, MSc, MD

    Indiana University Purdue University

    Regenstrief Institue

    1101 W 10th St,

    Indianapolis, IN, 46202

    United States

    Phone: 1 2628538872

    Email: cbarboi@iu.edu


    Background: Severity of illness scores—Acute Physiology and Chronic Health Evaluation, Simplified Acute Physiology Score, and Sequential Organ Failure Assessment—are current risk stratification and mortality prediction tools used in intensive care units (ICUs) worldwide. Developers of artificial intelligence or machine learning (ML) models predictive of ICU mortality use the severity of illness scores as a reference point when reporting the performance of these computational constructs.

    Objective: This study aimed to perform a literature review and meta-analysis of articles that compared binary classification ML models with the severity of illness scores that predict ICU mortality and determine which models have superior performance. This review intends to provide actionable guidance to clinicians on the performance and validity of ML models in supporting clinical decision-making compared with the severity of illness score models.

    Methods: Between December 15 and 18, 2020, we conducted a systematic search of PubMed, Scopus, Embase, and IEEE databases and reviewed studies published between 2000 and 2020 that compared the performance of binary ML models predictive of ICU mortality with the performance of severity of illness score models on the same data sets. We assessed the studies' characteristics, synthesized the results, meta-analyzed the discriminative performance of the ML and severity of illness score models, and performed tests of heterogeneity within and among studies.

    Results: We screened 461 abstracts, of which we assessed the full text of 66 (14.3%) articles. We included in the review 20 (4.3%) studies that developed 47 ML models based on 7 types of algorithms and compared them with 3 types of the severity of illness score models. Of the 20 studies, 4 (20%) were found to have a low risk of bias and applicability in model development, 7 (35%) performed external validation, 9 (45%) reported on calibration, 12 (60%) reported on classification measures, and 4 (20%) addressed explainability. The discriminative performance of the ML-based models, which was reported as AUROC, ranged between 0.728 and 0.99 and between 0.58 and 0.86 for the severity of illness score–based models. We noted substantial heterogeneity among the reported models and considerable variation among the AUROC estimates for both ML and severity of illness score model types.

    Conclusions: ML-based models can accurately predict ICU mortality as an alternative to traditional scoring models. Although the range of performance of the ML models is superior to that of the severity of illness score models, the results cannot be generalized due to the high degree of heterogeneity. When presented with the option of choosing between severity of illness score or ML models for decision support, clinicians should select models that have been externally validated, tested in the practice environment, and updated to the patient population and practice environment.

    Trial Registration: PROSPERO CRD42021203871; https://tinyurl.com/28v2nch8

    JMIR Med Inform 2022;10(5):e35293

    doi:10.2196/35293

    Keywords

    artificial intelligence; machine learning; intensive care unit mortality; severity of illness models 


    Background

    In the United States, intensive care unit (ICU) care costs account for 1% of the US gross domestic product, underscoring the need to optimize its use to attenuate the continued increase in health care expenditures [1]. Models that characterize the severity of illnesses of patients who are critically ill by predicting complications and ICU mortality risk can guide organizational resource management and planning, implementation and support of critical clinical protocols, and benchmarking and are proxies for resource allocation and clinical performance [2]. Although the medical community values the information provided by such models, they are not consistently used in practice because of their complexity, marginal predictive capacity, and limited internal or external validation [2-5].

    Severity of illness score models require periodic updates and customizations to reflect changes in medical care and regional case pathology [6]. Scoring models are prone to high interrater variability, are less accurate for patients with increased severity of illness score or specific clinical subgroups, are not designed for repeated applications, and cannot represent patients’ status trends [7]. The Acute Physiology and Chronic Health Evaluation (APACHE)-II (APACHE-II) and Simplified Acute Physiology Score (SAPS), developed in the 80s, are still in use [8]. The underlying algorithms for APACHE-IV are in the public domain and are available at no cost; however, their use is time intensive and is facilitated by software that requires payments for licensing implementation and maintenance [9]. Compared with SAPS-III, which uses data exclusively obtained within the first hour of ICU admission [10], APACHE-IV uses data from the first day (24 hours) [11]. Although the Sequential Organ Failure Assessment (SOFA) is an organ dysfunction score that detects differences in the severity of illness and is not designed to predict mortality, it is currently used to estimate mortality risk based on the mean, highest, and time changes accrued in the score during the ICU stay [11].

    The availability of machine-readable data from electronic health records enables the analysis of large volumes of medical data using machine learning (ML) methods. ML algorithms enable the exploration of high-dimensional data and the extraction of features to develop models that solve classification or regression problems. These algorithms can fit linear and nonlinear associations and interactions between predictive variables and relate all or some of the predictive variables to an outcome. The increased flexibility of ML models comes with the risk of overfitting training data; therefore, model testing on external data is essential to ensure adequate performance on previously unseen data. In model development, the balance between the model’s accuracy and generalizability, or bias and variance, is achieved through model training on a training set and hyperparameter optimization on a tuning set. Once a few models have been trained, they can be internally validated on a split-sample data set or cross-validated; the candidate model chosen is then validated on an unseen test data set to calculate its performance metrics and out of sample error [12]. The choice of algorithm is critical for providing a balance between interpretability, accuracy, and susceptibility to bias and variance [13]. Compared with the severity of illness scores, ML models can incorporate large numbers of covariates and temporal data, nonlinear predictors, trends in measured variables, and complex interactions between variables [14]. Numerous ML algorithms have been integrated into ICU predictive models, such as artificial neural networks (NNs), deep reinforcement learning, support vector machines (SVMs), random forest models, genetic algorithms, clinical trajectory models, gradient boosting models, k-nearest neighbor, naive Bayes, and the Ensemble approach [15]. Despite the rapidly growing interest in using ML methods to support clinical care, modeling processes and data sources have been inadequately described [16,17]. Consequently, the ability to validate and generalize the current literature’s results is questionable.

    Objectives

    This study aims to systematically review and meta-analyze studies that compare binary classification ML models with the severity of illness scores for predicting ICU mortality and determine which models have superior performance. This review intends to provide actionable guidance to clinicians on the prognostic value of ML models compared with the severity of illness scores in supporting clinical decision-making, as well as on their performance, in the context of the current guidelines [18] and recommendations for reporting ML analysis in clinical research [19] (Table 1).

    Table 1. Recommended structure for reporting MLa models.
    Research question and ML justificationData sources and preprocessing (feature selection)Model training and validation
    Clinical questionPopulationHardware, software, and packages used
    Intended use of the resultSample record and measurement characteristicsEvaluation (calibration and discrimination)
    Defined problem typeData collection and qualityConfiguration (parameters and hyperparameters)
    Available dataData structure and typesModel optimization and generalization (hyperparameter tuning and parameter limits)
    Defined ML method and rationaleDifferences between evaluation and validation setsValidation method and data split and cross-validation
    Defined evaluation measures, training protocols, and validationData preprocessing (data aggregation, missing data, transformation, and label source)Validation method performance metrics on an external data set
    N/AbInput configurationReproducibility, code reuse, and explainability

    aML: machine learning.

    bN/A: not applicable.


    We conducted a systematic review of the relevant literature. The research methods and reporting followed the PRISMA (Preferred Reporting Items for Systematic Reviews and Meta-Analyses) 2020 statement and guide to review and meta-analysis of prediction models [20,21].

    Information Sources and Search Strategy

    Between December 15 and 18, 2020, we performed a comprehensive search in the bibliographic databases PubMed, Scopus, Embase, and IEEE of the literature published between December 2000 and December 15, 2020. These databases were available free of charge from the university library. We selected PubMed for its significance in biomedical electronic research; Scopus for its wide journal range, keyword search, and citation analysis; Embase because of its European union literature coverage; and IEEE Xplore for its access to engineering and computer science literature.

    The search terms included control terms (Medical Subject Headings and Emtree) and free-text terms. The filters applied during the search of all 4 databases were Humans and Age:Adult. A search of the PubMed database using the terms (AI artificial intelligence) OR (machine learning) AND (intensive care unit) AND (mortality) identified 125 articles. The Scopus database was searched using the terms KEY (machine learning) OR KEY (artificial-intelligence) AND KEY (intensive care unit) AND KEY (mortality) revealed 182 articles. The Embase database queries using the terms (AI Artificial Intelligence) OR (machine learning) AND (intensive care unit) AND (mortality) resulted in 103 articles. The IEEE database search using the terms (machine learning) OR (artificial intelligence) AND (intensive care unit) AND (mortality) produced 51 citations.

    A total of 2 authors (CB and AT) screened titles and abstracts and recorded the reasons for exclusion. The same authors (CB and AT) independently reviewed the previously selected full-text articles to determine their eligibility for quantitative and qualitative assessments. Both authors revisited the discrepancies to guarantee database accuracy and checked the references of the identified articles for additional papers. A third researcher (LNM) was available to resolve any disagreements.

    Eligibility Criteria and Study Selection

    We included studies that compared the predictive performance of newly developed ML classification models predictive of ICU mortality with the severity of illness score models on the same data sets in the adult population. To be included in the review, the studies had to provide information on the patient cohort, model development and validation, and performance metrics. Both prospective and retrospective studies were eligible for inclusion.

    Data Collection Process

    Data extraction was performed by CB, reviewed by AT, and guided by the CHARMS (Critical Appraisal and Data Extraction for Systematic Reviews of Prediction Modeling Studies) checklist [22] specifically designed for systematic reviews of prognostic prediction models. The methodological qualities of the included studies were appraised with guidance from the Prediction model Risk of Bias (ROB) Assessment Tool (PROBAST) [23]. The reported features of the ML models are shown in Table 2.

    Table 2. CHARMS (Checklist for Critical Appraisal and Data Extraction for Systematic Reviews of Prediction Modeling Studies) checklist.
    AuthorData source (description)Outcome mortalityData preparationModel trainingPredictive performanceGeneralizability



    		AaBbCcDdEeFfGgHhIiJjKkLlMmNn
    Pirracchio et al [1]MIMICo 2Hospital


    		✓✓✓✓

    		✓
    		✓
    		✓
    Nielsen et al [24]Danish ICUpHospital 30/90 days✓✓✓✓✓
    		✓✓
    		✓

    		✓
    Nimgaonkar et al [25]ICU IndiaHospital


    		✓✓✓✓






    Xia et al [26]MIMIC 328 days/hospital✓✓
    		✓✓
    		✓✓



    		✓✓
    Purushotham et al [27]MIMIC 3Hospital, 2 days, 3 days, 30 days, 1 year✓✓
    		✓✓
    		✓✓
    		✓

    		✓✓
    Nanayakkara et al [28]ANZICSqHospital✓✓✓✓✓✓✓✓

    		✓✓✓✓
    Meyer et al [29]GermanyHospital
    		✓
    		✓✓
    		✓✓
    		✓
    		✓
    		✓
    Meiring et al [7]CCHICr United KingdomHospital✓✓✓✓✓
    		✓




    		✓✓
    Lin et al [30]MIMIC 3Hospital✓✓
    		✓✓✓✓✓





    Krishnan et al [31]MIMIC 3ICU✓✓
    		✓✓
    		✓✓


    		✓

    Kang et al [32]KoreaHospital✓
    		✓✓✓✓✓✓

    		✓


    Johnson et al [33]United KingdomICU and hospital✓✓
    		✓✓✓✓
    		✓✓
    		✓

    Holmgren et al [34]SwedenHospital and 30 days

    		✓✓✓✓✓✓



    		✓✓
    Garcia-Gallo et al [35]MIMIC 3Hospital and 1 year✓✓✓✓✓✓✓✓



    		✓✓
    El-Rashidy et al [36]MIMIC 3ICU and hospital✓✓
    		✓✓
    		✓✓
    		✓

    		✓✓
    Silva et al [37]EURICUSs 2ICU✓✓✓✓✓
    		✓✓✓
    		✓✓

    Caicedo-Torres et al [38]MIMIC 3ICU✓✓
    		✓✓
    		✓✓

    		✓


    Deshmukh et al [39]eICU-CRDtICU✓✓✓✓✓
    		✓✓

    		✓
    		✓
    Ryan et al [40]MIMIC 2ICU and hospital✓✓✓✓✓
    		✓✓
    		✓

    		✓✓
    Mayaud et al [41]MIMIC 2Hospital✓✓
    		✓✓✓✓
    		✓




    aData normalization/outlier addressed.

    bMissing data addressed.

    cHyperparameter optimization addressed.

    dOverfitting/shrinkage and cross-validation addressed.

    ePredictor selection, full model versus backward elimination.

    fCalibration assessed (Brier, Hosmer-Lemeshow, and calibration plot).

    gDiscrimination/reclassification performed (net reclassification improvement/integrated discrimination improvement).

    hClassification reported.

    iRecalibration performed.

    jExternally validated.

    kExplainability addressed/decision curve analysis.

    lClinical applicability addressed.

    mPrediction span defined.

    nIntended moment of use reported.

    oMIMIC: Medical Information Mart for Intensive Care.

    pICU: intensive care unit.

    qANZICS: Australia New Zealand Intensive Care Unit Society.

    rCCHIC: Critical Care Health Informatics Collaborative.

    sEURICUS: European ICU studies.

    teICU CRD: Electronic ICU Collaborative Research Database.

    Assessment of the ROB and Quality of Reviewed Studies

    The reviewers used the PROBAST tool to assess the methodological quality of each study for ROB and concerns regarding applicability in 4 domains: study participants, predictors, outcome, and analysis [23]. The reviewers evaluated the applicability of the selected studies by assessing the extent to which the studied outcomes matched the goals of the review in the 4 domains. We evaluated the ROB by assessing the primary study design and conduct, predictor selection process, outcome definition, and performance analysis. The ROB in the reporting models’ performance was appraised by exploring the reported measures of calibration (model’s predicted risk of mortality vs the observed risk), discrimination (model’s ability to discriminate between patients who are alive or expired), classification (sensitivity and specificity), and reclassification (net reclassification index). The performance of the models on internal data sets not used for model development—internal validation—and on data sets originating from an external patient population–external validation—were weighted in the ROB assignment. The ROB and applicability were assigned as low risk, high risk, or unclear risk according to PROBAST recommendations [42].

    Meta-analysis and Performance Metrics

    The C statistic–area under the receiver operating curve (AUROC) is the most commonly reported estimate of discriminative performance for binary outcomes [43-46] and the pragmatic performance measure of ML and severity of illness score models previously used in the medical literature to compare models based on different computational methods [21,45-47]. It is generally interpreted as follows: an AUROC of 0.5 suggests no discrimination, 0.7 to 0.8 is considered acceptable performance, 0.8 to 0.9 is considered excellent performance, and >0.9 is considered outstanding performance [48]. We included the performance of models developed using similar algorithms in forest plots and performed heterogeneity diagnostics and investigations without calculating a pooled estimate [49]. The results were pooled only for studies that followed a consistent methodology that included the external validation or benchmarking of the models. Random-effects meta-analyses computed the pooled AUROC for the following subgroups of ML algorithms—NNs and Ensemble—and the following subgroups of scoring models—SAPS II, APACHE II and SOFA. The AUROC for each model type was weighted using the inverse of its variance. Pooled AUROC estimates for each model were meta-analyzed along with 95% CIs of the estimates and were reported in forest plots together with the associated heterogeneity statistics (I2, τ2, and Cochran Q). Cochran Q statistic (also known as the chi-square statistic) determines the within-study variation, τ2 determines the between-study variability, and I2 represents the percentage of variability from the AUROC estimate not caused by sampling error [36]. The Cochran Q P value is denoted as P. Meta-analyses were conducted in R (version 3.6.1) [37] (see Multimedia Appendix 1 for scripts).


    Selection Process

    Of the 461 screened abstracts, we excluded 372 (80.7%) because of relevance (models not developed to predict ICU mortality), 9 (2%) duplicates, 6 (1.3%) reviews, and 8 (1.7%) conference proceedings (not intended for clinical application). We assessed the full text of 66 articles; the most common performance method reported to allow comparison between all models and a meta-analysis was the C statistic–AUROC. Of the 66 articles, we excluded 12 (18%) articles because of limited information on model development, 22 (33%) articles because of a lack of comparison with clinical scoring models, and 12 (18%) articles as the AUROC was not reported. The search strategy and selection process are illustrated in Figure 1.

    ‎
    Figure 1. Search strategy and selection process. AUROC: area under the receiver operating curve; ICU: intensive care unit.
    View this figure

    Assessment of the Prediction Model Development

    The 20 studies reported 47 ML models that were developed based on 7 types of algorithms and compared them with 3 severity of illness score models. All ML models were developed through a retrospective analysis of the ICU data sets. Of the 20 studies, 10 (50%) used data from the publicly available Medical Information Mart for Intensive Care database (Beth Israel Deaconess Medical Center in the United States) at different stages of expansion. Of the 20 studies, 10 (50%) used national health care databases (Danish, Australia-New Zealand, United Kingdom, and Sweden) or ICU-linked databases (Korea, India, and the United Kingdom). One of the studies included data from >80 ICUs belonging to >40 hospitals [33], and one of the studies’ ICU-linked database collected data from 9 European countries [37]. The cohorts generating the data sets used for model development and internal testing ranged from 1571 to 217,289 patients, with a median of 15,789 patients. Of the 20 studies, 10 (50%) used data from patients admitted to general ICUs, whereas 10 (50%) studies used data from patients who were critically ill with specific pathologies: gastrointestinal bleeds [39], COVID-19 and pneumonia–associated respiratory failure [40], postcardiac arrest [28], postcardiac surgery [29,36], acute renal insufficiency [30,32], sepsis [35,41], or neurological pathology [25]. The lower age thresholds for study inclusion ranges were 12 years [25], 15 years [26,27], 16 years [33,35,38], 18 years [24,29,40], and 19 years [30]. Within the studied cohorts, mortality ranged from 0.08 to 0.5 [29,32,36].

    The processes and tools used for the selection of predicting variables were described in 65% (13/20) of studies and included the least absolute shrinkage and selection operator, stochastic gradient boosting [33,35], genetic algorithms, and particle swarm optimization [33]. Approximately 15% (3/20) of studies [25,26,35] reported multiple models developed on variable predictor sets, which were subsequently tested for the best performance, validation, and calibration. The number of predictive variables used in the final models varied between 1 and 80, with a median of 21. The most common predicting variables are shown in Figure 2 and are grouped by the frequency of occurrence in the studies.

    ‎
    Figure 2. Frequency and type of ML model input variables (x-axis: number of studies using the input variables; y-axis: input variable). ASA: American Society of Anesthesiology; COPD: chronic obstructive pulmonary disease; CVA: cerebral vascular accident; FIO2: fraction of inspired oxygen; ICU: intensive care unit; LFT: liver function test; ML: machine learning; RBC: red blood cell; SpO2: oxygen saturation; PaO2: arterial oxygen pressure; PaCO2: arterial CO2 pressure.
    View this figure

    All studies developed models on 24-hour data; furthermore, ML models were developed on the first hour of ICU data [34]; the first 48-hour data [27,38,41]; 3-day data [40]; 5-day data [7]; 10-day data [26]; or on patients’ prior medical history collected from 1 month, 3 months, 6 months, 1 year, 2.5 years, 5 years, 7.5 years, 10 years, and 23 years [24]. The frequency of data collection ranged from every 30 minutes [29], 1 hour [1,25,27,37], 3 hours, 6 hours, 12 hours, 15 hours [38], and 24 hours [7,36] to every 27 hours, 51 hours, and 75 hours [40].

    Researchers handled missing data and continuous and fixed variables differently. A total of 6 model developers provided no information on missing data [1,25,29,32,34,39], and 1 [27] addressed data cleaning. Researchers [24,30,33,35-37,40,41] removed the records with missing values ranging from 1 missing value per admission to 30%, 50%, and 60% of missing data. One of the studies [29] included only variables documented for at least 50% of the patients and imputed the missing values with the last measured value for the feature. Missing values (up to 60%) were forward-filled; backward-filled; or replaced with means (continuous variables) or modes (categorical variables), normal values, averages [24,28,36,38,40], predictive mean matching [7], or linear interpolation imputation method [26]. The data were normalized using the minimum-maximum normalization technique. The time prediction of hospital mortality was undefined in 45% (9/20) of studies and varied from 2 or 3 days to 28 days, 30 days [26], 90 days [24], and up to 1 year [24] in the others.

    There was a wide range in the prevalence of mortality among studies (0.08-0.56), creating a class imbalance in the data sets. In studies with low investigated outcome mortality, few researchers addressed the problem of class imbalance (survivors vs nonsurvivors) through balanced training [24,37], random resampling [29], undersampling [36], or class penalty and reweighting schemes [38]. A breakdown of the model characteristics is presented in Table 3.

    Table 3. Information on the MLa prediction model development, validation, and performance, and on the severity of illness score performance.
    AuthorML model type (AUROCb test)Data training/test (split %)FeaturesK-fold/validationExternal validation data setML AUROC externalSeverity of illness score model type (AUROC)
    Pirracchio et al, 2015 [1]
    • Ensemble SICULAc (0.85)
      • ‎
    		24,508175-fold cross-validation2000.94
    • SAPSd-II (0.78)
      • ‎
    • APACHEe-II (0.83)
      • ‎
    • SOFf (0.71)
      • ‎
    Nielsen et al, 2019 [24]
    • NNg (0.792)
      • ‎
    		10,368 (80/20)445-fold cross-validation15280.773
    • SAPS-II (0.74)
      • ‎
    • APACHE-II (0.72)
      • ‎
    Nimgaonkar et al, 2004 [25]
    • NN (0.88)
      • ‎
    		2962 (70/30)15N/AhN/AN/A
    • APACHE-II (0.77)
      • ‎
    Xia et al, 2019 [26]
    • Ensemble-LSTMi (0.85)
      • ‎
    • LSTM (0.83)
      • ‎
    • DTj (0.82)
      • ‎
    		18,415 (90/10)50Bootstrap and RSMkN/AN/A
    • SAPS-II (0.77)
      • ‎
    • SOFA (0.73)
      • ‎
    • APACHE-II (0.74)
      • ‎
    Purushotham et al, 2018 [27]
    • NN (0.87)
      • ‎
    • Ensemble (0.84)
      • ‎
    		35,62717/22/ 1365-fold cross-validationExternal benchmarkN/A
    • SAPS-II (0.80)
      • ‎
    • SOFA (0.73)
      • ‎
    Nanayakkara et al, 2018 [28]
    • DT (0.86)
      • ‎
    • SVMl (0.86)
      • ‎
    • NN (0.85)
      • ‎
    • Ensemble (0.87)
      • ‎
    • GBMm (0.87)
      • ‎
    		39,560 (90/10)295-fold cross-validationN/AN/A
    • APACHE-III (0.8)
      • ‎
    Meyer et al, 2018 [29]
    • NN (0.95)
      • ‎
    		5898 (90/10)5210-fold cross-validation59890.81
    • SAPS-II (0.71)
      • ‎
    Meiring et al, 2018 [7]
    • DT (0.85)
      • ‎
    • NN (0.86)
      • ‎
    • SVM (0.86)
      • ‎
    		80/202521,911
    LOOn    		N/AN/A
    • APACHE-II (0.83)
      • ‎
    Lin et al, 2019 [30]
    • DT (0.86)
      • ‎
    • NN (0.83)
      • ‎
    • SVM (0.86)
      • ‎
    		19,044155-fold cross-validationN/AN/A
    • SAPS-II (0.79)
      • ‎
    Krishnan et al, 2018 [31]
    • NN-ELMo (0.99)
      • ‎
    		10,155 (75/25)110-fold cross-validationN/AN/A
    • SAPS (0.80)
      • ‎
    • SOFA (0.73)
      • ‎
    • APSp-III (0.79)
      • ‎
    Kang et al, 2020 [32]
    • SVM (0.77)
      • ‎
    • DT (0.78)
      • ‎
    • NN (0.776)
      • ‎
    • k-NNq (0.76)
      • ‎
    		1571 (70/30)3310-fold cross-validationN/AN/A
    • SOFA (0.66)
      • ‎
    • APACHE-II (0.59)
      • ‎
    Johnson et al, 2013 [33]
    • LRr univariate (0.902)
      • ‎
    • LR multivariate (0.876)
      • ‎
    		39,070 (80/20)1010-fold cross-validation23,6180.837 (univariate); 0.868 (multivariate)
    • APS-III (0.86)
      • ‎
    Holmgren et al, 2019 [34]
    • NN (0.89)
      • ‎
    		217,289 (80/20)85-fold cross-validationN/AN/A
    • SAPS-III (0.85)
      • ‎
    Garcia-Gallo et al, 2020 [35]
    • SGB-LASSOs (0.803)
      • ‎
    		5650 (70/30)18 140 3710-fold cross-validationN/AN/A
    • SOFA (0.58)
      • ‎
    • SAPS (0.70)
      • ‎
    El-Rashidy et al, 2020 [36]
    • Ensemble (0.93)
      • ‎
    		10,664 (75/25)8010-fold cross-validationExternal benchmarkN/A
    • APACHE-II (0.73)
      • ‎
    • SAPS-II (0.81)
      • ‎
    • SOFA-II (0.78)
      • ‎
    Silva et al, 2006 [37]
    • NN (0.85)
      • ‎
    		13,164 (66/33)12Hold outN/AN/A
    • SAPS- II (0.8)
      • ‎
    Caicedo-Torres et al, 2019 [38]
    • NN (0.87)
      • ‎
    		22,413225-fold cross-validationN/AN/A
    • SAPS-II (0.73)
      • ‎
    Deshmukh et al, 2020 [39]
    • XGBt (0.85)
      • ‎
    		5691 (80/20)345-fold cross-validationN/AN/A
    • APACHE-IV (0.8)
      • ‎
    Ryan et al, 2020 [40]
    • DT (0.86)
      • ‎
    		35,061 (80/20)125-fold cross-validation1140.91
    • qSOFAu (0.76)
      • ‎
    Mayaud et al, 2013 [41]
    • GAv+LR (0.82)
      • ‎
    		2113 (70/30)25BBCCVwN/AN/A
    • APACHE-III (0.68)
      • ‎

    aML: machine learning.

    bAUROC: area under the receiver operating curve.

    cSICULA: Super ICU Learner Algorithm.

    dSAPS: Simplified Acute Physiology Score.

    eAPACHE: Acute Physiology and Chronic Health Evaluation.

    fSOFA: Sequential Organ Failure Assessment.

    gNN: neural network.

    hN/A: not applicable.

    iLSTM: long short-term memory.

    jDT: decision tree.

    kRSM: random subspace method.

    lSVM: support vector machine.

    mGBM: gradient boosting machine.

    nLOO: leave one out.

    oELM: extreme learning machine.

    pAPS: Acute Physiology Score.

    qk-NN: k-nearest neighbor.

    rLR: logistic regression.

    sSGB-LASSO: stochastic gradient boosting least absolute shrinkage and selection operator.

    tXGB: extreme gradient boosting.

    uqSOFA: Quick Sequential Organ Failure Assessment.

    vGA: genetic algorithm.

    wBBCV: bootstrap bias–corrected cross-validation.

    Overview of ML Algorithms and Model Validation

    The reviewers recorded the ML model types based on the final trained model structure rather than on the algorithm used for fitting the model (Table 3). The reviewers noted a diversity of strategies in model fitting, although the implemented models defined the operating functions and transformations. Of the 20 studies, NNs were applied in 13 (65%) [7,24-32,34,37,38], decision trees in 8 (40%) [7,26,28,30,32,35,39,40], SVM in 4 (20%) [7,28,30,32], and Ensemble of algorithms in 4 (20%) [1,27,28,36]. The types of algorithms used in the same study varied between 1 and 5. All studies provided information on data training and internal testing (see Table 2 for k-fold validation and data splitting). Of the 20 studies, 5 (25%) [1,24,29,33,40] performed validation on external data sets ranging from 114 to 23,618 patients, and 2 (10%) studies [27,36] benchmarked the ML model performance against existing ML mortality prediction models; 14 (70%) studies reported CIs for the measure of discrimination AUROC, 9 (45%) studies reported on calibration (Hosmer-Lemeshow, calibration curve, or Brier score), and 12 (60%) studies reported on classification measures (Table 4). Approximately 10% (2/20) of studies were available for use in clinical practice [1,33]; the models’ decisions were explained with local interpretable model-agnostic explanations [28] or the Shapley additive explanations method (SHAP) [39]. The AUROC of the ML models ranged from 0.728 to 0.99 for predicting mortality.

    Table 4. Reported performance measures of the MLa models.
    Author and ML modelClassification measurementsCalibration measurementsOther

    		SpecificityPPVb/precisionRecall/sensitivityF1 scoreAccuracyHLc scoreBrier scoreCalibration curve
    Pirrachio et al[1]

    		Ensemble SLd-1N/AeN/AN/AN/AN/AN/A0.079Uf=0.0007 (calibration plot)DSg=0.21

    		Ensemble SL-2N/AN/AN/AN/AN/AN/A0.079U=0.006 (calibration plot)DS=0.26
    Nielsen et al[24]

    		NNhN/A0.388N/AN/AN/AN/AN/AN/AMathews correlation coefficient
    Purushotham et al[27]

    		NNN/AN/AN/AN/AN/AN/AN/AN/A0.491 (AUPRCi)

    		EnsembleN/AN/AN/AN/AN/AN/AN/AN/A0.435 (AUPRC)
    Nimgaonkar etal [25]

    		NN-15 featuresN/AN/AN/AN/AN/A27.7N/ACalibration plotN/A

    		NN-22 featuresN/AN/AN/AN/AN/A22.4N/ACalibration plotN/A
    Xia et al[26]

    		Ensemble-LSTMj0.75030.2940.77580.42620.7533N/AN/AN/AN/A

    		LSTM0.77460.3050.73840.43170.7703N/AN/AN/AN/A

    		RFk0.78070.3060.711970.42900.7734N/AN/AN/AN/A
    Nanayakkara et al[28]

    		RF0.790.750.76N/A0.78N/A0.156Calibration plot0.47 (log loss)

    		SVCl0.810.770.75N/A0.78N/A0.153Calibration plot0.47 (log loss)

    		GBMm0.780.750.8N/A0.79N/A0.147Calibration plot0.45 (log loss)

    		NN0.720.710.82N/A0.77N/A0.158Calibration plot0.48 (log loss)

    		Ensemble0.810.770.77N/A0.79N/A0.148Calibration plot0.45 (log loss)
    Meyer et al[29]

    		RNNn0.910.90.850.880.88N/AN/AN/AN/A
    Meiring et al[7]

    		DTo, NN, SVMpN/AN/AN/AN/AN/AN/AN/AN/AN/A
    Lin et al[30]

    		RFN/AN/AN/A0.4590.728N/A0.085Calibration plotN/A

    		NNN/AN/AN/A0.4060.666N/A0.091Calibration plotN/A

    		SVMN/AN/AN/A0.4600.729N/A0.086Calibration plotN/A
    Krishnan et al[31]

    		ANN-ELMqN/AN/A0.980.980.98N/AN/AN/AMathews correlation coefficient
    Kang et al[32]

    		k-NNrN/AN/AN/A0.7450.673N/AN/ACalibration plotN/A

    		SVMN/AN/AN/A0.7520.696N/AN/ACalibration plotN/A

    		RFN/AN/AN/A0.7620.69N/AN/ACalibration plotN/A

    		XGBsN/AN/AN/A0.7630.711N/AN/ACalibration plotN/A

    		NNN/AN/AN/A0.749
    		N/AN/ACalibration plotN/A
    Johnson et al[33]

    		LRt univariateN/AN/AN/AN/AN/A220.051N/AN/A

    		LR multivariateN/AN/AN/AN/AN/A19.60.048N/AN/A
    Holmgren et al[34]

    		NNN/AN/AN/AN/AN/AN/A0.106Calibration plotN/A
    Garcia-Gallo et al[35]

    		SGBuN/AN/AN/AN/A0.7250.0916N/ACalibration plotN/A

    		SGB-LASSOvN/AN/AN/AN/A0.7120.0916N/ACalibration plotN/A
    El-Rashidy et al[36]

    		Ensemble0.94N/A0.9110.9370.944N/AN/AN/AN/A
    Silva et al[37]

    		NN0.79N/A0.78N/A0.7921N/AN/AN/AN/A
    Caicedo-Torres et al[38]

    		NN0.827N/A0.75N/AN/AN/AN/AN/AN/A
    Deshmukh et al[39]

    		XGB0.27N/A1N/AN/AN/AN/AN/AN/A
    Ryan et al[40]

    		XGB0.75N/A0.8010.3780.75N/AN/AN/AN/A
    Mayaud et al[41]

    		GAw+LRN/AN/AN/AN/AN/A10.43N/ACalibration plotN/A

    aML: machine learning.

    bPPV: positive predictive value.

    cHL: Hosmer-Lemeshow.

    dSL: super learner.

    eN/A: not available.

    fU statistics.

    gDS: discrimination slope.

    hNN: neural network.

    iAUPRC: area under the precison-recall curve.

    jLSTM: long short-term memory.

    kRF: random forest.

    lSVC: support vector classifier.

    mGBM: gradient boosting machine.

    nRNN: recurrent neural network.

    oDT: decision tree.

    pSVM: support vector machine.

    qANN-ELM: artificial neural network extreme learning machine.

    rk-NN: k-nearest neighbor.

    sXGB: extreme gradient boosting.

    tLR: logistic regression.

    uSGB: stochastic gradient boosting.

    vLASSO: least absolute shrinkage and selection operator.

    wGA: genetic algorithm.

    The performance of the ML models was compared with that of the following severity of illness scoring models: APACHE-II (6/20, 30%), APACHE-III (2/20, 10%), APACHE-IV (1/20, 5%), SAPS-II (11/20, 55%), SAPS-III (1/20, 5%), SOFA (9/20, 45%), and Acute Physiology Score-3 (2/20, 10%; Table 3). The severity of illness scores’ discrimination reported as AUROC was associated with a CI in 65% (13/20) of studies. Calibration of the severity of illness score models was reported in 30% (6/20) of studies. Approximately 60% (12/20) of studies reported binary classification results. The severity of illness scores used for comparison and associated AUROCs were 0.70 to 0.803 for SAPS, 0.588 to 0.782 for SOFA, and 0.593 to 0.86 for APACHE (Table 3).

    Analysis of ROB and Applicability

    The results of the analysis of the ROB in the selection of the study population, predictors, outcome definition, and performance reporting are presented in Table 5. The results of the assessment of the developed ML models’ applicability regarding the study participants and setting, the predictors used in the ML models’ development and their timing, the outcome definition and prediction by the models, and the analysis that reports the models’ performance are also presented in Table 5. Of the 47 models, 4 (9%) models [1,17,23,29] were identified as having a low risk, and 3 (6%) models were rated as having an uncertain ROB and applicability model development [24,29,40]. The main reason for the high ROB in the overall judgment of the study was the lack of external validation, which was identified in 28% (13/47) of the models.

    Table 5. Assessment for ROBa and applicability for prognostic models with the Prediction model ROB Assessment Tool checklist.
    AuthorsROB and applicability

    		ParticipantsPredictorsOutcomeAnalysisOverall judgment

    		ROBApplicabilityROBApplicabilityROBApplicabilityROBROBApplicability
    Pirracchio et al [1]LowbLowLowLowLowLowLowLowLow
    Nielsen et al [24]LowLowUnclearcLowUnclearLowLowUnclearLow
    Nimgaonkar et al [25]LowUnclearLowLowLowLowHighdHighUnclear
    Xia et al [26]LowLowLowLowUnclearLowHighHighLow
    Purushotham et al [27]LowLowLowLowLowLowLowLowLow
    Nanayakkara et al [28]LowUnclearLowLowLowLowHighHighUnclear
    Meyer et al [29]LowUnclearLowLowLowLowLowLowUnclear
    Meiring et al [7]LowLowLowLowLowLowHighHighLow
    Lin et al [30]LowUnclearLowLowLowLowHighHighUnclear
    Krishnan et al [31]LowLowLowLowLowLowHighHighLow
    Kang et al [32]LowUnclearLowLowLowLowHighHighUnclear
    Johnson et al [33]LowLowLowLowLowLowLowLowLow
    Holmgren et al [34]LowLowLowLowUnclearLowHighHighLow
    Garcia-Gallo et al [35]LowUnclearLowLowLowLowHighHighUnclear
    El-Rashidy et al [36]LowLowLowLowLowLowLowLowLow
    Silva et al [37]LowLowLowLowLowLowHighHighLow
    Caicedo-Torres et al [38]LowLowLowLowLowLowHighHighLow
    Deshmukh et al [39]LowUnclearLowLowLowLowHighHighUnclear
    Ryan et al [40]LowLowUnclearLowLowLowLowUnclearLow
    Mayaud et al [41]LowUnclearLowUnclearLowLowHighHighUnclear

    aROB: risk of bias.

    bLow risk: no relevant shortcomings in ROB assessment.

    cUnclear risk: unclear ROB in at least one domain and all other domains at low ROB.

    dHigh risk: relevant shortcomings in the ROB assessment, at least one domain with high ROB, or model developed without external validation.

    Meta-analysis

    Forest plots for the NN, Ensemble, SOFA, SAPS II, and APACHE-II models and the associated heterogeneity tests are shown in Figures 3-7. The forest plots and tests of heterogeneity for SVM, NN, DT, and Ensemble models that were not externally validated can be seen in Multimedia Appendix 2. The AUROC for each model type was weighted using the inverse of its variance. Most of the 95% CIs of AUROC estimates from various studies did not overlap within the forest plot; considerable variation among AUROC estimates for both ML and severity of illness score model types was noted. Regrading tests of heterogeneity, I2 varied between 99% and 100%, τ2 ranged from 0.0003 to 0.0034, and P was consistently <.01. In Figures 3-7 and Multimedia Appendix 2, the gray boxes represent the weight estimates of the AUROC value from each study. The horizontal line through each gray box illustrates the 95% CI of the AUROC value from that study. Black horizontal lines through a gray box indicate that the CI limits exceeded the length of the gray box. White horizontal lines represent the CI limits that are within the length of the gray box. I2, τ2, and Cochran Q P value (denoted as P) are heterogeneity tests.

    Random-effects meta-analysis results of the computed pooled AUROC of the ML subgroup models that were externally validated or benchmarked NNs and Ensemble are shown in Figure 3 and Figure 4.

    ‎
    Figure 3. Meta-analysis results: pooled AUROC for externally validated Ensemble models. Gray boxes represent the fixed weight estimates of the AUROC value from each study. Larger gray boxes represent larger fixed weight estimates of the AUROC values. The horizontal line through each gray box illustrates the 95% CI of the AUROC value from that study. Black horizontal lines through a gray box indicate that the CI limits exceed the length of the gray box. White horizontal lines represent CI limits that are within the length of the gray box. The vertical dashed lines in the forest plot are the estimated random pooled effect of the AUROC value from the random-effects meta-analysis. The gray diamonds illustrate the 95% CI for the random pooled effects. Tests of heterogeneity included I2, τ2, and Cochran Q P value (denoted as P). AUROC: area under the receiver operating curve;.
    View this figure
    ‎
    Figure 4. Meta-analysis results: pooled AUROC for externally validated NN models. Gray boxes represent the fixed weight estimates of the AUROC value from each study. Larger gray boxes represent larger fixed weight estimates of the AUROC values. The horizontal line through each gray box illustrates the 95% CI of the AUROC value from that study. Black horizontal lines through a gray box indicate that the CI limits exceed the length of the gray box. White horizontal lines represent CI limits that are within the length of the gray box. The vertical dashed lines in the forest plot are the estimated random pooled effect of the AUROC value from the random-effects meta-analysis. The gray diamonds illustrate the 95% CI for the random pooled effects. Tests of heterogeneity included I2, τ2, and Cochran Q P value (denoted as P). AUROC: area under the receiver operating curve; NN: neural network.
    View this figure
    ‎
    Figure 5. Meta-analysis results: pooled AUROC for SAPS-II. Gray boxes represent the fixed weight estimates of the AUROC value from each study. Larger gray boxes represent larger fixed weight estimates of the AUROC values. The horizontal line through each gray box illustrates the 95% CI of the AUROC value from that study. Black horizontal lines through a gray box indicate that the CI limits exceed the length of the gray box. White horizontal lines represent CI limits that are within the length of the gray box. The vertical dashed lines in the forest plot are the estimated random pooled effect of the AUROC value from the random-effects meta-analysis. The gray diamonds illustrate the 95% CI for the random pooled effects. Tests of heterogeneity included I2, τ2, and Cochran Q P value (denoted as P). AUROC: area under the receiver operating curve; SAPS-II: Simplified Acute Physiology Score II.
    View this figure
    ‎
    Figure 6. Meta-analysis results: pooled AUROC for SOFA. Gray boxes represent the fixed weight estimates of the AUROC value from each study. Larger gray boxes represent larger fixed weight estimates of the AUROC values. The horizontal line through each gray box illustrates the 95% CI of the AUROC value from that study. Black horizontal lines through a gray box indicate that the CI limits exceed the length of the gray box. White horizontal lines represent CI limits that are within the length of the gray box. The vertical dashed lines in the forest plot are the estimated random pooled effect of the AUROC value from the random-effects meta-analysis. The gray diamonds illustrate the 95% CI for the random pooled effects. Tests of heterogeneity included I2, τ2, and Cochran Q P value (denoted as P). AUROC: area under the receiver operating curve; SOFA: Sequential Organ Failure Assessment.
    View this figure
    ‎
    Figure 7. Meta-analysis results: pooled AUROC for APACHE-II. Gray boxes represent the fixed weight estimates of the AUROC value from each study. Larger gray boxes represent larger fixed weight estimates of the AUROC values. The horizontal line through each gray box illustrates the 95% CI of the AUROC value from that study. Black horizontal lines through a gray box indicate that the CI limits exceed the length of the gray box. White horizontal lines represent CI limits that are within the length of the gray box. The vertical dashed lines in the forest plot are the estimated random pooled effect of the AUROC value from the random-effects meta-analysis. The gray diamonds illustrate the 95% CI for the random pooled effects. Tests of heterogeneity included I2, τ2, and Cochran Q P value (denoted as P). APACHE-II: Acute Physiology and Chronic Health Evaluation-II; AUROC: area under the receiver operating curve;.
    View this figure

    The results of heterogeneity for the NN models were as follows: τ2= 0.0043 (95% CI 0.0014-0.2100), I2=99.9% (95% CI 99.8%-99.9%), P<.01. The results of heterogeneity for the Ensemble models were as follows:

    τ2=0.0015 (95% CI 0.0005-0.0223), I2=99.7% (95% CI 99.6%-99.8%), P<.01. The results were synthesized, and the models are presented in Figure 3 and Figure 4. The results of heterogeneity for the APACHE-2 models were as follows: τ2=0.0046 (95% CI 0.0011-0.1681), I2=99.7% (95% CI 99.6%-99.8%), P<.01. The results of heterogeneity for the SAPS-II models were as follows: τ2=0.0012 (95% CI 0.0005-0.0133), I2=99.2% (95% CI 98.9%-99.4%), P<.01. The results of heterogeneity for the SOFA models were as follows: τ2=0.0009 (95% CI 0.0003-0.0461), I2=99.1% (95% CI 98.5%-99.4%), P<.01 (Figures 5-7).


    Principal Findings

    This is the first study to critically appraise the literature comparing the ML and severity of illness score models to predict ICU mortality. In the reviewed articles, the AUROC of the ML models demonstrated very good discrimination. The range of the ML model AUROC was superior to that of the severity of illness score AUROC. The meta-analysis demonstrated a high degree of heterogeneity and variability within and among studies; therefore, the AUROC performances of the ML and severity of illness score models cannot be pooled, and the results cannot be generalized. Every I2 value is >97.7%; most of the 95% CIs of AUROC estimates from various studies did not overlap within the forest plot, suggesting considerable variation among AUROC estimates for model types. The CI for AUROC and the statistical significance of the difference in model performance were inconsistently reported within studies. The high heterogeneity came from the diverse study population and practice location, age of inclusion, primary pathology, medical management leading to the ICU admission, and time prediction window. The heterogenous data management (granularity, frequency of data input, data management, number of predicting variables, prediction timeframe, time series analysis, and training set imbalance) affected model development. It may have resulted in bias, primarily in studies where it has not been addressed (Table 2). Generally, authors reported the ML algorithms with predictive power superior to the clinical scoring system (Table 3); the number of ML models with inferior performance not reported is unknown, which raises the concern of reporting bias. The classification measures of performance were inconsistently reported and required a predefined probability threshold; therefore, models showed different sensitivity and specificity based on the chosen threshold. The variations in the prevalence of the studied outcome secondary to imbalanced data sets make the interpretation of the accuracy difficult. The models’ calibration cannot be interpreted because of limited reporting. The external validation process that is necessary to establish generalization was lacking in 65% (13/20) of studies (Table 2). The limited and variable performance metrics reported precludes a comprehensive model performance comparison among studies. The decision curve analysis and model interpretability (explainability) that are necessary to promote transparency and understanding of the model’s predictive reasoning was addressed in 25% (5/20) of studies. Results of the clinical performance of ML mortality prediction models as alternatives to the severity of illness score are scarce.

    The reviewed studies inconsistently and incompletely captured the descriptive characteristics and other method parameters for ML-based predictive model development. Therefore, we cannot fully assess the superiority or inferiority of ML-based ICU mortality prediction compared with traditional models; however, we recognize the advantage that flexibility in model design offers in the ICU setting.

    Study Limitations

    This review included studies that were retrospective analyses of data sets with known outcome distributions and incorporated the results of interventions. It is unclear which models were developed exclusively for research purposes; hence, they were not validated. We evaluated studies that compared ML-based mortality prediction models with the severity of illness score–based models, although these models relied on different development statistical methods, variable collection times, and outcome measurement methodologies (SOFA).

    The comparison between the artificial intelligence (AI) and severity of illness score models relies only on AUROC values as measures of calibration, discrimination, and classification are not uniformly reported. The random-effects meta-analysis was limited to externally validated models. Owing to the level of heterogeneity, the performance results for most AI and severity of illness score models could not be pooled. The authors recognize that 25% (5/20) of the articles were published between 2004 and 2015 before the TRIPOD (Transparent Reporting of a multivariable prediction model for Individual Prognosis or Diagnosis) recommendations for model development and reporting [18]; thus, they were not aligned with the guidelines.

    The reviewers assessed the models’ ROB and applicability and were aware of the risk of reporting and publication bias favoring the ML models. However, the high heterogeneity among studies prevents an unambiguous interpretation of the funnel plot.

    Conclusions and Recommendations

    The results of our analysis show that the reporting methodology is incomplete, nonadherent to the current recommendations, and consistent with previous observations [16,50]. The lack of consistent reporting of the measures of the reliability calibration (Brier score and calibration curve of reliability deviation), discrimination, and classification of the probabilistic estimates on external data makes the comparative effectiveness of risk prediction models challenging and has been noted by other authors [43].

    Predictive models of mortality can substantially increase patient safety, and by incorporating subtle changes in organ functions that affect outcomes, these models support the early recognition and diagnosis of patients who are deteriorating, thus providing clinicians with additional time to intervene. The heterogeneity of the classification models that was revealed in detail in this review underlines the importance of recognizing the models’ ability for temporal and geographical generalization or proper adaptation to previously unseen data [51]. These concepts apply to both models; similar to the ML models, severity of illness score requires periodical updates and customizations to reflect changes in medical care and regional case pathology over time [6].

    Our findings lead to the following recommendations for model developers:

    1. State whether the developed ML models are intended for clinical practice
    2. If models are intended for clinical applications, provide full transparency of the clinical setting from which the data are acquired and all the model development steps; validate the models externally to ensure generalizability
    3. If intended for clinical practice, report models’ performance metrics, which include measures of discrimination, calibration, and classification, and attach explainer models to facilitate interpretability

    Before using ML and/or severity of illness score models as decision support systems to guide clinical practice, we make the following recommendations for clinicians:

    1. Be cognizant of the similarities or discrepancies between the cohort used for model development and the local practice population, the practice setting, the model’s ability to function prospectively, and the models’ lead times
    2. Acquire knowledge of the model’s performance during testing in the local practice
    3. Ensure that the model is periodically updated to changes in patient characteristics and/or clinical variables and adjusted to new clinical practices and therapeutics
    4. Confirm that the models’ data are monitored and validated and that the model’s performance is periodically updated
    5. When both the severity of illness score and ML models are available, determine one model’s superiority and clinical reliability versus the other through randomized controlled trials
    6. When ML models guide clinical practice, ensure that the model makes the correct recommendation for theright reasonsand consult the explainer model
    7. Identify clinical performance metrics that evaluate the impact of the AI tool on the quality of care, efficiency, productivity, and patient outcomes and account for variability in practice

    AI developers must search for and clinicians must be cognizant of the unintended consequences of AI tools; both must understand human–AI tool interactions. Healthcare organization administrators must be aware of the safety, privacy, causality, and ethical challenges when adopting AI tools and recognize the Food and Drug Administration guiding principles for AI/ML development [52].

    Acknowledgments

    The statistical analysis was funded by the Department of Anesthesiology at Northwestern University.

    The authors would like to thank Dr Shaun Grannis for reviewing the manuscript and funding the publication cost.

    Conflicts of Interest

    None declared.

    ‎
    Multimedia Appendix 1

    R scripts for meta-analysis.

    PDF File (Adobe PDF File), 59 KB
    ‎
    Multimedia Appendix 2

    Forest plots for neural network, decision tree, support vector machine, and Ensemble-based models.

    PDF File (Adobe PDF File), 257 KB
    1. Pirracchio R, Petersen ML, Carone M, Rigon MR, Chevret S, van der Laan MJ. Mortality prediction in intensive care units with the Super ICU Learner Algorithm (SICULA): a population-based study. Lancet Respir Med 2015 Jan;3(1):42-52 [FREE Full text] [CrossRef] [Medline]
    2. Salluh JI, Soares M. ICU severity of illness scores: APACHE, SAPS and MPM. Curr Opin Crit Care 2014 Oct;20(5):557-565. [CrossRef] [Medline]
    3. Keegan MT, Gajic O, Afessa B. Comparison of APACHE III, APACHE IV, SAPS 3, and MPM0III and influence of resuscitation status on model performance. Chest 2012 Oct;142(4):851-858 [FREE Full text] [CrossRef] [Medline]
    4. Breslow MJ, Badawi O. Severity scoring in the critically ill: part 1--interpretation and accuracy of outcome prediction scoring systems. Chest 2012 Jan;141(1):245-252. [CrossRef] [Medline]
    5. Sarkar R, Martin C, Mattie H, Gichoya JW, Stone DJ, Celi LA. Performance of intensive care unit severity scoring systems across different ethnicities. medRxiv (forthcoming) 2021 Jan 20 [FREE Full text] [CrossRef] [Medline]
    6. Strand K, Flaatten H. Severity scoring in the ICU: a review. Acta Anaesthesiol Scand 2008 Apr;52(4):467-478. [CrossRef] [Medline]
    7. Meiring C, Dixit A, Harris S, MacCallum NS, Brealey DA, Watkinson PJ, et al. Optimal intensive care outcome prediction over time using machine learning. PLoS One 2018 Nov 14;13(11):e0206862 [FREE Full text] [CrossRef] [Medline]
    8. Keegan MT, Gajic O, Afessa B. Severity of illness scoring systems in the intensive care unit. Crit Care Med 2011 Jan;39(1):163-169. [CrossRef] [Medline]
    9. Delahanty RJ, Kaufman D, Jones SS. Development and evaluation of an automated machine learning algorithm for in-hospital mortality risk adjustment among critical care patients. Crit Care Med 2018 Jun;46(6):e481-e488. [CrossRef] [Medline]
    10. Balkan B, Essay P, Subbian V. Evaluating ICU clinical severity scoring systems and machine learning applications: APACHE IV/IVa case study. In: Proceedings of the 40th Annual International Conference of the IEEE Engineering in Medicine and Biology Society. 2018 Presented at: EMBC '18; Honolulu, HI, USA; July 18-21, 2018 p. 4073-4076. [CrossRef]
    11. Minne L, Abu-Hanna A, de Jonge E. Evaluation of SOFA-based models for predicting mortality in the ICU: a systematic review. Crit Care 2008;12(6):R161 [FREE Full text] [CrossRef] [Medline]
    12. de Hond AA, Leeuwenberg AM, Hooft L, Kant IM, Nijman SW, van Os HJ, et al. Guidelines and quality criteria for artificial intelligence-based prediction models in healthcare: a scoping review. NPJ Digit Med 2022 Jan 10;5(1):2 [FREE Full text] [CrossRef] [Medline]
    13. Goldstein BA, Navar AM, Carter RE. Moving beyond regression techniques in cardiovascular risk prediction: applying machine learning to address analytic challenges. Eur Heart J 2017 Jun 14;38(23):1805-1814 [FREE Full text] [CrossRef] [Medline]
    14. Breiman L. Statistical modeling: the two cultures (with comments and a rejoinder by the author). Statist Sci 2001 Aug 1;16(3):199-231. [CrossRef]
    15. Hassanipour S, Ghaem H, Arab-Zozani M, Seif M, Fararouei M, Abdzadeh E, et al. Comparison of artificial neural network and logistic regression models for prediction of outcomes in trauma patients: a systematic review and meta-analysis. Injury 2019 Feb;50(2):244-250. [CrossRef] [Medline]
    16. Andaur Navarro CL, Damen JA, Takada T, Nijman SW, Dhiman P, Ma J, et al. Risk of bias in studies on prediction models developed using supervised machine learning techniques: systematic review. BMJ 2021 Oct 20;375:n2281 [FREE Full text] [CrossRef] [Medline]
    17. Bozkurt S, Cahan EM, Seneviratne MG, Sun R, Lossio-Ventura JA, Ioannidis JP, et al. Reporting of demographic data and representativeness in machine learning models using electronic health records. J Am Med Inform Assoc 2020 Dec 09;27(12):1878-1884 [FREE Full text] [CrossRef] [Medline]
    18. Luo W, Phung D, Tran T, Gupta S, Rana S, Karmakar C, et al. Guidelines for developing and reporting machine learning predictive models in biomedical research: a multidisciplinary view. J Med Internet Res 2016 Dec 16;18(12):e323 [FREE Full text] [CrossRef] [Medline]
    19. Stevens LM, Mortazavi BJ, Deo RC, Curtis L, Kao DP. Recommendations for reporting machine learning analyses in clinical research. Circ Cardiovasc Qual Outcomes 2020 Oct;13(10):e006556 [FREE Full text] [CrossRef] [Medline]
    20. Page MJ, McKenzie JE, Bossuyt PM, Boutron I, Hoffmann TC, Mulrow CD, et al. The PRISMA 2020 statement: an updated guideline for reporting systematic reviews. Int J Surg 2021 Apr;88:105906. [CrossRef] [Medline]
    21. Debray TP, Damen JA, Snell KI, Ensor J, Hooft L, Reitsma JB, et al. A guide to systematic review and meta-analysis of prediction model performance. BMJ 2017 Jan 05;356:i6460. [CrossRef] [Medline]
    22. Moons KG, de Groot JA, Bouwmeester W, Vergouwe Y, Mallett S, Altman DG, et al. Critical appraisal and data extraction for systematic reviews of prediction modelling studies: the CHARMS checklist. PLoS Med 2014 Oct;11(10):e1001744 [FREE Full text] [CrossRef] [Medline]
    23. Moons KG, Wolff RF, Riley RD, Whiting PF, Westwood M, Collins GS, et al. PROBAST: a tool to assess risk of bias and applicability of prediction model studies: explanation and elaboration. Ann Intern Med 2019 Jan 01;170(1):W1-33 [FREE Full text] [CrossRef] [Medline]
    24. Nielsen AB, Thorsen-Meyer HC, Belling K, Nielsen AP, Thomas CE, Chmura PJ, et al. Survival prediction in intensive-care units based on aggregation of long-term disease history and acute physiology: a retrospective study of the Danish National Patient Registry and electronic patient records. Lancet Digit Health 2019 Jun;1(2):e78-e89 [FREE Full text] [CrossRef] [Medline]
    25. Nimgaonkar A, Karnad DR, Sudarshan S, Ohno-Machado L, Kohane I. Prediction of mortality in an Indian intensive care unit. Comparison between APACHE II and artificial neural networks. Intensive Care Med 2004 Feb;30(2):248-253. [CrossRef] [Medline]
    26. Xia J, Pan S, Zhu M, Cai G, Yan M, Su Q, et al. A long short-term memory ensemble approach for improving the outcome prediction in intensive care unit. Comput Math Methods Med 2019 Nov 3;2019:8152713 [FREE Full text] [CrossRef] [Medline]
    27. Purushotham S, Meng C, Che Z, Liu Y. Benchmarking deep learning models on large healthcare datasets. J Biomed Inform 2018 Jul;83:112-134 [FREE Full text] [CrossRef] [Medline]
    28. Nanayakkara S, Fogarty S, Tremeer M, Ross K, Richards B, Bergmeir C, et al. Characterising risk of in-hospital mortality following cardiac arrest using machine learning: a retrospective international registry study. PLoS Med 2018 Nov;15(11):e1002709 [FREE Full text] [CrossRef] [Medline]
    29. Meyer A, Zverinski D, Pfahringer B, Kempfert J, Kuehne T, Sündermann SH, et al. Machine learning for real-time prediction of complications in critical care: a retrospective study. Lancet Respir Med 2018 Dec;6(12):905-914. [CrossRef] [Medline]
    30. Lin K, Hu Y, Kong G. Predicting in-hospital mortality of patients with acute kidney injury in the ICU using random forest model. Int J Med Inform 2019 May;125:55-61. [CrossRef] [Medline]
    31. Krishnan GS, Kamath SS. A supervised learning approach for ICU mortality prediction based on unstructured electrocardiogram text reports. In: Proceedings of the 23rd International Conference on Applications of Natural Language to Information Systems. 2018 May 22 Presented at: NLDB '18; June 13-15, 2018; Paris, France p. 126-134. [CrossRef]
    32. Kang MW, Kim J, Kim DK, Oh K, Joo KW, Kim YS, et al. Machine learning algorithm to predict mortality in patients undergoing continuous renal replacement therapy. Crit Care 2020 Feb 06;24(1):42 [FREE Full text] [CrossRef] [Medline]
    33. Johnson AE, Kramer AA, Clifford GD. A new severity of illness scale using a subset of Acute Physiology And Chronic Health Evaluation data elements shows comparable predictive accuracy. Crit Care Med 2013 Jul;41(7):1711-1718. [CrossRef] [Medline]
    34. Holmgren G, Andersson P, Jakobsson A, Frigyesi A. Artificial neural networks improve and simplify intensive care mortality prognostication: a national cohort study of 217,289 first-time intensive care unit admissions. J Intensive Care 2019;7:44 [FREE Full text] [CrossRef] [Medline]
    35. García-Gallo JE, Fonseca-Ruiz NJ, Celi LA, Duitama-Muñoz JF. A machine learning-based model for 1-year mortality prediction in patients admitted to an Intensive Care Unit with a diagnosis of sepsis. Med Intensiva (Engl Ed) 2020 Apr;44(3):160-170. [CrossRef] [Medline]
    36. El-Rashidy N, El-Sappagh S, Abuhmed T, Abdelrazek S, El-Bakry HM. Intensive care unit mortality prediction: an improved patient-specific stacking ensemble model. IEEE Access 2020;8:133541-133564. [CrossRef]
    37. Silva A, Cortez P, Santos MF, Gomes L, Neves J. Mortality assessment in intensive care units via adverse events using artificial neural networks. Artif Intell Med 2006 Mar;36(3):223-234. [CrossRef] [Medline]
    38. Caicedo-Torres W, Gutierrez J. ISeeU: visually interpretable deep learning for mortality prediction inside the ICU. J Biomed Inform 2019 Oct;98:103269 [FREE Full text] [CrossRef] [Medline]
    39. Deshmukh F, Merchant SS. Explainable machine learning model for predicting GI bleed mortality in the intensive care unit. Am J Gastroenterol 2020 Oct;115(10):1657-1668. [CrossRef] [Medline]
    40. Ryan L, Lam C, Mataraso S, Allen A, Green-Saxena A, Pellegrini E, et al. Mortality prediction model for the triage of COVID-19, pneumonia, and mechanically ventilated ICU patients: a retrospective study. Ann Med Surg (Lond) 2020 Nov;59:207-216 [FREE Full text] [CrossRef] [Medline]
    41. Mayaud L, Lai PS, Clifford GD, Tarassenko L, Celi LA, Annane D. Dynamic data during hypotensive episode improves mortality predictions among patients with sepsis and hypotension. Crit Care Med 2013 Apr;41(4):954-962 [FREE Full text] [CrossRef] [Medline]
    42. Wolff RF, Moons KG, Riley RD, Whiting PF, Westwood M, Collins GS, PROBAST Group†. PROBAST: a tool to assess the risk of bias and applicability of prediction model studies. Ann Intern Med 2019 Jan 01;170(1):51-58 [FREE Full text] [CrossRef] [Medline]
    43. Huang C, Li S, Caraballo C, Masoudi FA, Rumsfeld JS, Spertus JA, et al. Performance metrics for the comparative analysis of clinical risk prediction models employing machine learning. Circ Cardiovasc Qual Outcomes 2021 Oct;14(10):e007526. [CrossRef] [Medline]
    44. Christodoulou E, Ma J, Collins GS, Steyerberg EW, Verbakel JY, Van Calster B. A systematic review shows no performance benefit of machine learning over logistic regression for clinical prediction models. J Clin Epidemiol 2019 Jun;110:12-22. [CrossRef] [Medline]
    45. Song X, Liu X, Liu F, Wang C. Comparison of machine learning and logistic regression models in predicting acute kidney injury: a systematic review and meta-analysis. Int J Med Inform 2021 Jul;151:104484 [FREE Full text] [CrossRef] [Medline]
    46. Sufriyana H, Husnayain A, Chen YL, Kuo CY, Singh O, Yeh TY, et al. Comparison of multivariable logistic regression and other machine learning algorithms for prognostic prediction studies in pregnancy care: systematic review and meta-analysis. JMIR Med Inform 2020 Nov 17;8(11):e16503 [FREE Full text] [CrossRef] [Medline]
    47. Louie KS, Seigneurin A, Cathcart P, Sasieni P. Do prostate cancer risk models improve the predictive accuracy of PSA screening? A meta-analysis. Ann Oncol 2015 May;26(5):848-864 [FREE Full text] [CrossRef] [Medline]
    48. Hosmer DW, Lemeshow S. Assessing the fit of the model: 5.2.4. area under the receiver operating characteristic curve. In: Shewhart WA, Wilks SS, editors. Applied Logistic Regression. 2nd Edition. Hoboken, NJ, US: John Wiley & Sons; Sep 13, 2000:143-202.
    49. Higgins JP, Thomas J, Chandler J, Cumpston M, Li T, Page MJ, et al. Cochrane Handbook for Systematic Review of Interventions. London, UK: The Cochrane Collaboration; 2021.
    50. Brnabic A, Hess LM. Systematic literature review of machine learning methods used in the analysis of real-world data for patient-provider decision making. BMC Med Inform Decis Mak 2021 Feb 15;21(1):54 [FREE Full text] [CrossRef] [Medline]
    51. Futoma J, Simons M, Doshi-Velez F, Kamaleswaran R. Generalization in clinical prediction models: the blessing and curse of measurement indicator variables. Crit Care Explor 2021 Jul;3(7):e0453 [FREE Full text] [CrossRef] [Medline]
    52. Good Machine Learning Practice for Medical Device Development: Guiding Principles. U.S. Food & Drug Administration. 2021.   URL: https:/​/www.​fda.gov/​medical-devices/​software-medical-device-samd/​good-machine-learning-practice-medical-device-development-guiding-principles [accessed 2022-05-08]

    ‎
    AI: artificial intelligence
    APACHE: Acute Physiology and Chronic Health Evaluation
    AUROC: area under the receiver operating curve
    CHARMS: Critical Appraisal and Data Extraction for Systematic Reviews of Prediction Modeling Studies
    ICU: intensive care unit
    ML: machine learning
    NN: neural network
    PRISMA: Preferred Reporting Items for Systematic Reviews and Meta-Analyses
    PROBAST: Prediction model Risk of Bias Assessment Tool
    ROB: risk of bias
    SAPS: Simplified Acute Physiology Score
    SOFA: Sequential Organ Failure Assessment
    SVM: support vector machine
    TRIPOD: Transparent Reporting of a multivariable prediction model for Individual Prognosis or Diagnosis

    Edited by C Lovis; submitted 29.11.21; peer-reviewed by P Giabbanelli, R Gore; comments to author 13.02.22; revised version received 24.04.22; accepted 25.04.22; published 31.05.22

    Copyright

    ©Cristina Barboi, Andreas Tzavelis, Lutfiyya NaQiyba Muhammad. Originally published in JMIR Medical Informatics (https://medinform.jmir.org), 31.05.2022.

    This is an open-access article distributed under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work, first published in JMIR Medical Informatics, is properly cited. The complete bibliographic information, a link to the original publication on https://medinform.jmir.org/, as well as this copyright and license information must be included.