- Open Access
To develop a regional ICU mortality prediction model during the first 24 h of ICU admission utilizing MODS and NEMS with six other independent variables from the Critical Care Information System (CCIS) Ontario, Canada
Journal of Intensive Care volume 4, Article number: 16 (2016)
The Correction to this article has been published in Journal of Intensive Care 2020 8:9
Intensive care unit (ICU) scoring systems or prediction models evolved to meet the desire of clinical and administrative leaders to assess the quality of care provided by their ICUs. The Critical Care Information System (CCIS) is province-wide data information for all Ontario, Canada level 3 and level 2 ICUs collected for this purpose. With the dataset, we developed a multivariable logistic regression ICU mortality prediction model during the first 24 h of ICU admission utilizing the explanatory variables including the two validated scores, Multiple Organs Dysfunctional Score (MODS) and Nine Equivalents Nursing Manpower Use Score (NEMS) followed by the variables age, sex, readmission to the ICU during the same hospital stay, admission diagnosis, source of admission, and the modified Charlson Co-morbidity Index (CCI) collected through the hospital health records.
This study is a single-center retrospective cohort review of 8822 records from the Critical Care Trauma Centre (CCTC) and Medical-Surgical Intensive Care Unit (MSICU) of London Health Sciences Centre (LHSC), Ontario, Canada between 1 Jan 2009 to 30 Nov 2012. Multivariable logistic regression on training dataset (n = 4321) was used to develop the model and validate by bootstrapping method on the testing dataset (n = 4501). Discrimination, calibration, and overall model performance were also assessed.
The predictors significantly associated with ICU mortality included: age (p < 0.001), source of admission (p < 0.0001), ICU admitting diagnosis (p < 0.0001), MODS (p < 0.0001), and NEMS (p < 0.0001). The variables sex and modified CCI were not significantly associated with ICU mortality. The training dataset for the developed model has good discriminating ability between patients with high risk and those with low risk of mortality (c-statistic 0.787). The Hosmer and Lemeshow goodness-of-fit test has a strong correlation between the observed and expected ICU mortality (χ 2 = 5.48; p > 0.31). The overall optimism of the estimation between the training and testing data set ΔAUC = 0.003, indicating a stable prediction model.
This study demonstrates that CCIS data available after the first 24 h of ICU admission at LHSC can be used to create a robust mortality prediction model with acceptable fit statistic and internal validity for valid benchmarking and monitoring ICU performance.
Patients in the intensive care units (ICUs) have a heterogeneous disease processes and illness severity. Scoring systems developed for ICU patients were introduced 34 years ago with the goal of using physiologic data available at ICU admission to predict individual patient outcomes. Although these predictions have little utility for managing individual patients, they do provide a mechanism for assessing ICU performance by comparing the actual outcome in a given population to the expected outcome determined by the prediction algorithms. The scores that assess the disease severity on admission and are used to predict outcome include the Acute Physiology and Chronic Health Evaluation (APACHE) [1–3], the Simplified Acute Physiological Score (SAPS) , and the Mortality Prediction Model (MPM) . The organ dysfunction scores that assess the presence and severity of organ dysfunction include the Multiple Organ Dysfunction Score (MODS)  and Sequential Organ Failure Assessment (SOFA) . The score that assesses nursing workload is the Therapeutic Intervention Scoring System (TISS) , and the Nine Equivalents of Nursing Manpower Use Score (NEMS)  (Additional file 1: Table S1) assesses ICU resource utilization and efficiency. Many of these measurement systems involve resource-intensive data collection.
In 2007, the Critical Care Services Ontario (CCSO), a division of the Ontario Ministry of Health and Long Term Care, developed the Critical Care Information System (CCIS). The purpose of CCIS is to provide the Ministry, Local Health Integration Networks (LHINs) and hospitals with information on bed availability, critical care utilization, and patient outcomes. The CCIS uses a web-based application to collect real-time information on every patient admitted to a critical care unit in Ontario acute care hospitals. Data captured includes, but it is not limited to the following: demographics, admission and discharge details, MODS on admission, daily NEMS, and patient outcomes such as ICU mortality and other outcomes associated with quality of care.
The MODS is an objective scale that quantifies the severity of multiple organ dysfunction for patients admitted to critical care. The score reflects six major organ systems and the specific physiological data associated with each system . A total of 0–4 points are assigned to each system, where a score of 0 is normal and 4 is the most dysfunctional to give a total maximum score of 24. MODS was not designed to predict mortality, but an increasing MODS does correlate with ICU outcome .
The NEMS was developed from the TISS-28 score and is a less complicated and is more widely used to measure resource utilization in critical care . The score is determined based on the need for any of the nine life support interventions. A weighted point is awarded to each of the nine categories to give a maximum score of 56. NEMS has been validated in large cohorts of ICU patients and is easy to use with minimum inter-observer variability . It has been utilized to classify the different levels of ICUs based on nursing workload efficacy as distinguished from the amount of care being provided.
In an effort to help hospitals analyze and interpret their data, CCSO produces and distributes quarterly reports that include a multitude of utilization and quality indicators of which one is ICU mortality. This data is presented in a manner that promotes benchmarking, but there is currently no means of risk adjustment to ensure that units are comparing themselves to centers with similar case mix and illness severity. Review of these reports shows that there are units with direct correlation of higher mean MODS and ICU mortality, but this is not always the case. There are units with similar MODS but differing mortality rates. The objective of this study is to investigate if existing CCIS data collected by the Critical Care Trauma Center (CCTC) and Medical-Surgical Intensive Care Unit (MSICU) of London Health Sciences Centre (LHSC) can be used to develop and validate an acceptable ICU mortality prediction model that might improve current performance measurement reporting.
Study design and patient population
This is a retrospective study of two adult intensive care units at the LHSC, an academic teaching facility, affiliated with The University of Western Ontario. The CCTC is a 30-bed general medical, surgical, trauma, and oncological unit, and the MSICU is a 25-bed unit that specializes in the care of various patient populations including neurosurgical, cardiovascular surgery, and transplantation patients. In both units, the care is provided by multidisciplinary teams of professional health care providers and is directed by physicians that have specialty training in critical care. This study was approved by the Western Health Research Ethic Board on 13 Nov 2013, IRB 00000940.
Between 1 Jan 2009 to 30 Nov 2012, data was prospectively collected on 4784 admissions to the CCTC and 4297 admissions to the MSICU. The primary endpoint of this study is to develop a mortality prediction model utilizing the available data from CCIS during the first 24-h ICU admission.
Relevant data from both units was exported from CCIS. Comorbidities are not included in CCIS but it is important because it may delay diagnosis, influence treatment decision, are related to complications, may influence chances of survival, and can confound analysis . Based on the APACHE II mortality prediction model, comorbidities can reflect diminished physiological reserve; thus, it is important to incorporate past relevant medical/surgical history into a mortality prediction model. All Canadian hospitals submit information to the Canadian Institute of Health Information (CIHI) which is an independent, not-for-profit corporation that aims to contribute to the improvement of the health of Canadians and the health care system by disseminating quality health information. The CIHI utilizes the Charlson Comorbidity Index (CCI) as their measure of comorbidity  which has been shown to be highly associated with 1-year patient mortality and has been widely used in clinical research . To obtain comorbidities for the prediction model, the ICD-10-CA data for patients admitted to either of the critical care units during the time frame of interest was obtained from the LHSC Health Records Department and only type 1 diagnoses, which specifically refer to pre-admission comorbidity was utilized.
The data from CCIS was exported in three parts. The reference dataset (N = 9081) contains demographic, admit/discharge date and time, admitting diagnosis, and the source of admission of the patients. The second part of the data is the MODS score on the day of admission to the ICU, and the third part of the data is the NEMS score on the day of admission to the ICU. The medical record number (MRN) and ICU admission date were used as the common linking variables to merge the MODS and NEMS data with the reference data into one file. Merging the reference dataset with the MODS dataset resulted in 8953 records, followed by merging with the NEMS dataset that resulted in 8924 records. There were a total of 157 (1.73 % of the original dataset) records missing. To obtain the ICD-10-CA data to calculate the modified CCS, the reference dataset was forwarded to the LHSC Records Department. A total of 8898 records were matched from Health Records, of which 183 records (2 % of the original dataset) from the CCIS dataset were not matched with the hospital records. This is likely due to error in the MRN number and/or failure to capture the ICU admission during discharge coding. Programmed SAS codes were created to extract only the type 1 diagnoses and calculation of the modified CCI score. Then, the dataset containing the modified CCI was combined with the final CCIS dataset (N = 8924) resulted in 8822 records for analysis, from which 2.9 % records were not captured from the original reference dataset, Fig. 1.
Covariates associated with ICU mortality
Based on the literature review, the covariates associated with risk of mortality upon admission to the ICU included age, gender, MODS, NEMS, CCI, source of admission, ICU admission diagnosis, and ICU readmission during the same hospital admission. The continuous, nominal, and ordinal data of the covariates are categorized accordingly (Additional file 2: Table S2). The variable threshold used to divide it into the groups for analysis was done using the mean value for each of the variables for survivors and non-survivors as a reference point. Then, through much iteration with SAS 9.3, the developed groups for each of the variables that give the best discriminatory performances and Hosmer and Lemeshow goodness-of-fit were utilized.
Univariate analysis was conducted for all baseline characteristics, and values were compared for survivors and non-survivors. For continuous variables, data are expressed as mean ± SD and comparisons conducted using the Student’s t test. For categorical variables, data are reported as proportions and comparison made using Pearson’s chi-square test. The prediction model for ICU mortality on admission was constructed by dividing the whole dataset into two random subgroups, “training” and “testing” set. To create the two random sample subgroups (training and testing set), the authors used the SAS 9.3 Ranuni function which generates a random number from a continuous uniform distribution with the interval (0,1) in which we used less than equal to 0.5 for group 1 and greater than 0.5 for group 2. Each 15th record will then be assigned to a random number. After assigning a random number to each record, it is then sorted in ascending or descending order of the random number assigned. A stepwise selection algorithm was also used to select from the eight covariates in the raw logistic regression prediction model. If the covariate significance was less than 0.05, it remained in the model; otherwise, the covariate exited the model. The remaining covariates were then fitted to a raw multivariable logistic regression model.
Often a predictive model’s estimate of these measures from the training set tends to overstate the predictive ability of the chosen model in another dataset. The amount of overestimation is referred to as the “optimism” of the estimate. To further obtain a valid criterion for the performance of the model, the logistic regression model then fit to the bootstrap sample and the corresponding value for the AUC was calculated. The fitted model was then applied to the original dataset, and the value of the AUC was recalculated. The differences in the values for the AUC provide an estimate of the optimism. This process is repeated 500 and 1000 times, and the results are averaged to provide a final bootstrap estimate for the optimism of the AUC.
The ability of the model to estimate mortality and agree with the actual outcome within groups of subjects of similar predicted risks by using the Hosmer and Lemeshow goodness-of-fit statistic obtained by grouping the subjects of the prediction model into k categories of percentiles. A good calibration is considered to be consistent with a small χ 2 value for the Hosmer and Lemeshow test statistic.
All analyses were performed using SAS 9.3 (SAS Institute Inc., Cary, NC, USA). All tests presented are two-sided, and a p value <0.05 is considered significant.
Table 1 compares the baseline characteristics for survivors and non-survivors in a cohort of 8822 subjects. There were 5037 (57.10 %) males with mean age 60.67 ± 17.19 and 3785 (42.90 %) females with mean age 61.16 ± 17.66 with no mortality differences between males and females (23.23 % vs. 22.17 %; p = 0.24) but, the mean age was higher for non-survivors than survivors (66.82 vs. 59.13; p < 0.0001). The admission sources with the highest mortality included the wards, emergency department, and other hospital transfers whereas patients admitted to critical care post-operatively had the lowest mortality (p < 0.0001). There was also a statistically significant association between ICU mortality and ICU admitting diagnosis (p < 0.0001) with cardiovascular/cardiac/vascular diseases having the highest mortality and gastrointestinal disease, the lowest. The MODS and NEMS scores were both significantly associated with ICU mortality (p < 0.0001) whereas the modified CCI did not have a statistically significant association between survivors and non-survivors, p = 0.74. There was also no significant difference in mortality for those re-admitted back to the ICU during the same hospital admission compared to those that were not readmitted (9.19 % vs. 9.01 %; p = 0.81.).
The baseline characteristics partitioned between the groups “training” and ‘testing’ sets were similar (Table 2). The total number of subjects in the training group was 4321 (48.98 %) as compared to 4501 (51.02 %) in the testing group. The combination of the patients from the two ICUs for each of the groups was evenly distributed between the two groups. In the training group, there were 2310 (53.46 %) subjects from CCTC and 2011 (46.54 %) from MSICU, while in the testing group, there were 2324 (51.63 %) from CCTC and 2177 (48.37 %) from MSICU.
In the multivariable logistic regression model, the stepwise selection algorithm eliminated the variables sex, p = 0.20 and readmit, p = 0.16. The c-statistic of the reduced model was smaller to that obtained when all explanatory variables were forced in (c = 0.774). The backward elimination algorithm eliminated readmit, p = 0.16 first, and then sex, p = 0.22, and resulted in a final model very similar to that using stepwise selection (c = 0.774). Because there is already a parsimony of variables in comparison to reported models, the two variables eliminated by both algorithms were forced back into the whole model.
The logistic regression analysis of the training group (Table 3) revealed that the model containing the explanatory variables compared to that with the intercept only significantly impacted the predictive ability of the model with the likelihood ratio, χ 2 = 835.98, p < 0.0001. The overall effect of each of the covariates on mortality revealed that all except gender (χ 2 = 0.59; p = 0.44) and CCI (χ 2 = 4.60; p = 0.10) had a significant independent effect on ICU mortality. Categorically, the odds of mortality are much higher for older patients, specifically ages 40–79 that is 2.23 (95 % confidence interval (CI) 1.64, 3.04), and those of patients ≥80 years old was 5.51 (95 % CI 3.87, 7.84) times higher than those of patients age ≤ 39. For patients admitted to the ICUs from the unit/ward (odds ratio (OR) = 4.93; 95 % CI 3.69, 6.59), other hospitals (OR = 3.054; 95 % CI 2.26, 4.12), emergency departments (OR = 2.71; 95 % CI 2.05, 3.59), and specialty units (OR = 2.66; 95 % CI 1.86, 3.81) had higher odds of ICU mortality as compared to patients admitted from the operating room/post-anesthesia unit. Patients admitted with cardiovascular/cardiac/vascular diagnoses had a higher mortality as compared to other etiologies. Higher scores in MODS and NEMS corresponded to increasing ICU mortality (p ≤ 0.0001). There was a weak association with ICU readmission and mortality (OR = 0.742; 95 % CI 0.56, 0.99; p = 0.04).
The discriminatory performance of the training model revealed the AUC was 0.787. This indicates that the model has good ability to distinguish between patients with a high risk of mortality and those with a low risk of mortality . The comparison of the receiver operating curve (ROC) curves for the training dataset and the testing dataset indicated an area difference of 0.026 (0.787–0.761), which reflects a very narrow gap or the optimism between the two curves, suggesting a small degradation in the model’s performance in prospective testing (Fig. 2). To validate this difference, the bootstrap processes were repeated 500 and 1000 times, and the results were averaged to provide an optimism correction for the AUC of 0.003 (AUC range = 0.758–0.790) which indicates that our model does not overpredict (Additional file 3: Table S3). The AUC comparison between this new model with MODS and NEMS alone in the new model revealed AUC = 0.776 and 0.736, respectively, which are lower than the combined scores AUC = 0.787 (Additional file 4: Table S4). Overall, the combination of the two scores in the model gives better discrimination ability between patients with high and low risks for ICU mortality during the first 24 h of ICU admission.
As a measure of calibration from the model, the Hosmer and Lemeshow goodness-of-fit statistics revealed χ 2 = 5.48 (p > 0.31) indicating strong agreement between observed and expected ICU mortality (Table 4).
The final formula equation for our model is
Log [Mortality (at 24 h ICU admission)] = −5.18 + 0.80[age(40–79)] + 1.71[age(>80)] + 0.60[Sex(male = 0 and female = 1) + 0.98[Other source admission] + 0.00[Operating room/post-anesthesia care] + 1.00[ER admission] + 1.12[Hospital-outside or within LHIN] + 1.60[Ward admission]Cardiovascular/Cardiac/Vascular] + 0.00[−0.81[Other diagnosis] − 0.80[Gastrointestinal] − 0.56[Respiratory] − 0.32[Trauma] + 0.002[Neurological] − 0.30[ICU re-admission] − 0.21[CCI(1–3)] + 0.05[CCI(>3)] + 0.0[NEMS(0–22)] + 0.39[NEMS(23–29)] + 1.02[NEMS(≥300] + 1.18[MODS(1–4)] + 1.91[MODS(5–8)] + 2.90[MODS(9–120] + 3.56[MODS(≥130].
Currently, many health care performance measurement systems are based on administrative databases. These systems are often developed to fulfill the needs of funding agencies and support the quality improvement plans of individual hospitals but rarely provide the necessary level of risk adjustment to provide meaningful comparison, over time or across facilities. It is also well known that prognostic research has received limited attention as compared to etiological, diagnostic, and therapeutic research. The development and application of robust prognostic models are essential for valid benchmarking. Prognostic scoring systems have been developed by the critical care specialty in an effort to quantify the severity of illness of a given patient or groups of patient [16–18]. Adjustment for the severity of illness enables one to monitor the performance of an ICU over time and to allow comparison of ICUs in the same or different hospitals. However, the fact that many prognostic models currently exist suggests that the optimum model has not yet been fully established and any of the developed prognostic models will have a limited effective life span [16, 19], due to changes in clinical practice over time and improved health care that can alter the risk of mortality for a given clinical situation. Thus, prognostic models require periodic updating. Major revisions of prognostic models that were published between 2005 and 2007 include APACHE IV (AUC = 0.88, χ 2 = 16.9, p = 0.08) , SAPS 3 (AUC = 0.848, χ 2 = 14.29, p = 0.16) , and MPM0 III (AUC = 0.823, χ 2 = 11.62, p = 0.31) . A recent review evaluated the latest versions of these models and concluded that although they represent great improvement compared to the previous ones, regular updates and local customizations are required . Also, the huge resources burden needed to collect a significantly large amount of data for the variables to generate these scores is daunting. The present study aimed to use available data already collected by our ICUs and other ICUs for a very limited number of variables for the two scores, MODS and NEMS in the province of Ontario, Canada, as mandated by CCIS as well as comorbidities diagnoses collected by our hospital health records to develop a mortality prognostic model. The AUC for our model was 0.787, which is considered acceptable or very good in differentiating between survivors and non-survivors [24, 25]. This model is well calibrated, showing good agreement between predicted and actual outcomes for all risk strata (Hosmer and Lemeshow χ 2 = 5.4761, p = 0.3146) .
To the best of our knowledge, this study is the first to use validated organ dysfunction score, MODS and severity assessment based on nursing workload, and NEMS coupled with five other variables selected. Overall, only two out the seven independent variables, gender and CCI, were not significant in the prediction of ICU mortality. Two retrospective studies contradicted our findings with respect to gender. One study of 24,778 patients admitted to the ICUs across Ontario, Canada in 2001–2002 revealed that females had a higher ICU mortality than males with an adjusted OR = 1.20 (95 % CI 1.10–1.31, p < 0.001) . Another study of 18,757 patients diagnosed with sepsis in 98 ICUs between 2003 and 2006 reported an adjusted OR = 1.11 (95 % CI 1.04–1.19, p < 0.01) . This increased mortality in female patients was not fully understood but could be explained by differences in the presentation of critical illness, decision-making, or unmeasured confounding factors that may contribute to these findings. The other possibility that our patient cohort was reported at a later time period may result in improved ICU access and earlier care of patients using the Critical Care Outreach Team (CCRT) [29–31].
We know that scoring systems used in the ICUs have been introduced and developed over the last 30 years. These models allow an assessment of the severity of disease and provide an estimate of ICU and hospital mortality. The MODS score independently has been used in many clinical studies and it has an excellent discriminating predictor of mortality in ICU patients [32–37]. The MODS in our study is a very strong predictor of mortality in the first 24 h of ICU admission, and it correlated very well with other scores such as the SOFA score and APACHE II score in terms of mortality prediction . However, the complexity of ICU care goes beyond the severity of illness or organ failure, the level of nursing workload, and NEMS as related to the ICU resource utilization that also correlated well with ICU mortality . Many other factors also have been shown to increase risks of in-hospital mortality after admission to the ICU, including increasing age and severity of acute illness, certain pre-existing medical conditions, source of admission, physiological measurements, and biochemical/hematological indices . By utilizing those other covariates, it may not necessarily improve the discrimination ability of predicted model but rather avoid the pitfall of either underpredicting or overpredicting ICU mortality when using only a single covariate in the prediction model. Although the NEMS in our study is not as a strong mortality predictor in the first 24 h of ICU admission as the MODS, with the combined scores in a prediction model, it provided a much better basis for evaluation of treatment results and documentation of the ICUs’ resource needs . Having knowledge of both severity of organ dysfunction and degree of resource utilization will provide a better basis for assessing whether ICU treatment(s) and/or administrative protocol(s) needs to be modified to improve patient care.
With comorbidities, specifically severe chronic organ system insufficiency or immunocompromised, those patients markedly influence outcomes  and this is supported by other outcome prediction scores. In our study, comorbidities were not predictive because we used the pre-admit comorbidities that existed prior to admission to the ICU as opposed to the conditions that were aggravated or developed subsequently. Another limitation was the actual condition captured in the CCI score that was developed on breast cancer patients and not in ICU patients to predict 1-year patient mortality using comorbidity data obtained from hospital chart review [42, 43]. The CCI is a validated weighted score, the weight for each of the co-morbidities may not fully reflect the severity of the disease and it may or may not include specific or unusual illnesses, and therefore, it is not an assessment of the impact of all illnesses on the overall health of the patient. Furthermore, the ICD-10-CA data is abstracted by medical record clerks and not entered by health care providers and can be subjected to errors based on lack of documentation and misinterpretation. Our model fit could be improved by capturing active chronic health status at ICU admission utilizing the APACHE II chronic health points which reflects diminished physiological reserve and markedly influence outcome [2, 41].
Acute diagnosis was not used in earlier prediction models with the exception of the APACHE II to IV scores. It was not until 1993 that MPM II started to include acute diagnosis to the model and SAPS 3 and MPM III followed suit in 2005 and 2007. However, the predictive accuracy over diagnosis showed that the performance of a prediction model can vary in different diagnostic groups . This is in agreement with previous research [45–47], which suggested that prognostic models can underpredict or overpredict mortality in specific patient sub-groups. The admitting diagnosis in our study is classified using broad, system-based categories and did not include specific diagnostic information to allow for comparisons within the generalized diagnoses, between subgroups or between-study populations. Although these broad diagnostic groups include specific diagnoses that are similar based on the system involved, the exact diagnoses within a group can have completely different treatments and outcomes. For example, the “cardiovascular/cardiac/vascular” diagnosis group used in this study has the highest mortality but includes less severe exact diagnoses with lower mortality rates. Patients with abdominal aorta aneurysm carried a much higher mortality than patients with myocardial infarction and cardiac bypass surgery thus skewing the mortality risk. Another limitation, some of the system-based groups such as metabolic/endocrine, genitourinary, musculoskeletal, skin, oncology, hematology, and “other” was collapsed together due to its small number of patients per group. This represented a significant heterogeneous population within a diagnostic group which would be difficult to interpret the statistical prediction accuracy.
Studies carried out in numerous countries indicated that the source of the patient admission is associated with mortality. Patients transferred from the ward within the same hospital showed a greater ICU mortality when compared with those coming from other sources [48–51]. This is in agreement with the present study where patients admitted from the ward had the highest mortality (OR = 4.93, 95 % CI 3.69–6.59, p < 0.0001), post-surgical patients had the lowest.
Those patient readmitted to the ICU did not have significant mortality differences to those patients not readmitted, which is contrary to published literature . This difference could be due to the implementation of the Critical Care Resuscitation Team (CCRT) that may intervene earlier on the wards of patients’ acute illness [53, 54]. Also, various service teams including CCRT are improving end of life (EOL) discussions with those patients previously admitted to the ICU, thus avoids a readmission [55, 56]. Our institution implemented the CCRT service in 2007.
Scoring systems in critical care have evolved to meet the desire of clinical and administrative leaders to assess the quality of care provided by the ICUs. Mortality is a key ICU quality metric and reflects many aspects of ICU care, including use of best practice, accurate diagnosis, and effective and timely therapies. Our model is locally calibrated to two ICUs in London, Ontario, Canada only, and the results may not be generalizable to other critical care units. But collectively, all ICUs in the province of Ontario, Canada gather the same data information; it is logical that a model be developed to benchmark ICU performance and improve the usability of the current reporting system. This study demonstrates that data from the CCIS can be used to create a mortality prediction model with good calibration and discrimination. Inclusion of data to capture active chronic health status and refinement of the acute diagnosis classification could further improve predictive ability of the developed model.
Acute physiology and chronic health evaluation
Area under the curve
Charlson’s Co-Morbidity Index
Critical care information system
Critical care services Ontario
Critical care trauma centre
Canadian Institute of Health Information
International statistical classification of diseases and related health problems – tenth version – Canada
Intensive care unit
Local health information networks
London health sciences centre
Multiple organs dysfunction score
Mortality prediction model
Medical and surgical intensive care unit
Nine equivalent manpower use score
Receiver operating curve
Simplified acute physiological score
Sequential organ failure assessment
Therapeutic intervention scoring system
Knaus WA, Zimmerman JE, Wagner DP, Draper EA, Lawrence DE. APACHE-acute physiology and chronic health evaluation: a physiologically based classification system. Crit Care Med. 1981;9(8):591–7.
Knaus WA, Draper EA, Wagner DP, Zimmerman JE. APACHE II: a severity of disease classification system. Crit Care Med. 1985;13(10):818–29.
Knaus WA, Wagner DP, Draper EA, et al. The APACHE III prognostic system. Risk prediction of hospital mortality for critically ill hospitalized adults. Chest. 1991;100(6):1619–36.
LeGall JR, Loirat P, Alperovitch A. APACHE II—a severity of disease classification system. Crit Care Med. 1986;14(8):754–5.
Lemeshow S, Teres D, Klar J, Avrunin JS, Gehlbach SH, Rapoport J. Mortality Probability Models (MPM II) based on an international cohort of intensive care unit patients. JAMA. 1993;270(20):2478–86.
Marshall JC, Cook DJ, Christou NV, Bernard GR, Sprung CL, Sibbald WJ. Multiple organ dysfunction score: a reliable descriptor of a complex clinical outcome. Crit Care Med. 1995;23(10):1638–52.
Vincent JL, Moreno R, Takala J, et al. The SOFA (Sepsis-related Organ Failure Assessment) score to describe organ dysfunction/failure. On behalf of the Working Group on Sepsis-Related Problems of the European Society of Intensive Care Medicine. Intensive Care Med. 1996;22(7):707–10.
Cullen DJ, Civetta JM, Briggs BA, Ferrara LC. Therapeutic intervention scoring system: a method for quantitative comparison of patient care. Crit Care Med. 1974;2(2):57–60.
Miranda DR, Nap R, de Rijk A, Schaufeli W, Iapichino G. Nursing activities score. Crit Care Med. 2003;31(2):374–82.
Reis MD, Moreno R, Iapichino G. Nine equivalents of nursing manpower use score (NEMS). Intensive Care Med. 1997;23(7):760–5.
Rothen HU, Kung V, Ryser DH, Zurcher R, Regli B. Validation of “nine equivalents of nursing manpower use score” on an independent data sample. Intensive Care Med. 1999;25(6):606–11.
Hall SF. A user’s guide to selecting a comorbidity index for clinical research. J Clin Epidemiol. 2006;59(8):849–55.
Moreno R, Reis MD. Nursing staff in intensive care in Europe: the mismatch between planning and practice. Chest. 1998;113(3):752–8.
Quan H, Li B, Couris CM, et al. Updating and validating the Charlson comorbidity index and score for risk adjustment in hospital discharge abstracts using data from 6 countries. Am J Epidemiol. 2011;173(6):676–82.
Fawdett T. An introduction to ROC analysis. Pattern Recogn Lett. 2006;27:861–74. Ref Type: Generic.
Keegan MT, Gajic O, Afessa B. Severity of illness scoring systems in the intensive care unit. Crit Care Med. 2011;39(1):163–9.
Rubenfeld GD, Angus DC, Pinsky MR, Curtis JR, Connors Jr AF, Bernard GR. Outcomes research in critical care: results of the American Thoracic Society Critical Care Assembly Workshop on Outcomes Research. The Members of the Outcomes Research Workshop. Am J Respir Crit Care Med. 1999;160(1):358–67.
Vincent JL, Moreno R. Clinical review: scoring systems in the critically ill. Crit Care. 2010;14(2):207.
Beck DH, Smith GB, Pappachan JV. The effects of two methods for customising the original SAPS II model for intensive care patients from South England. Anaesthesia. 2002;57(8):785–93.
Zimmerman JE, Kramer AA, McNair DS, Malila FM. Acute Physiology and Chronic Health Evaluation (APACHE) IV: hospital mortality assessment for today’s critically ill patients. Crit Care Med. 2006;34(5):1297–310.
Moreno RP, Metnitz PG, Almeida E, et al. SAPS 3—from evaluation of the patient to evaluation of the intensive care unit. Part 2: Development of a prognostic model for hospital mortality at ICU admission. Intensive Care Med. 2005;31(10):1345–55.
Higgins TL, Teres D, Copes WS, Nathanson BH, Stark M, Kramer AA. Assessing contemporary intensive care unit outcome: an updated Mortality Probability Admission Model (MPM0-III). Crit Care Med. 2007;35(3):827–35.
Salluh JI, Soares M. ICU severity of illness scores: APACHE, SAPS and MPM. Curr Opin Crit Care. 2014;20(5):557–65.
Hosmer DW, Hosmer T, Le CS, Lemeshow S. A comparison of goodness-of-fit tests for the logistic regression model. Stat Med. 1997;16(9):965–80.
Afessa B, Gajic O, Keegan MT. Severity of illness and organ failure assessment in adult intensive care units. Crit Care Clin. 2007;23(3):639–58.
Lemeshow S, Hosmer Jr DW. A review of goodness of fit statistics for use in the development of logistic regression models. Am J Epidemiol. 1982;115(1):92–106.
Fowler RA, Sabur N, Li P, et al. Sex- and age-based differences in the delivery and outcomes of critical care. CMAJ. 2007;177(12):1513–9.
Pietropaoli AP, Glance LG, Oakes D, Fisher SG. Gender differences in mortality in patients with severe sepsis or septic shock. Gend Med. 2010;7(5):422–37.
Bellomo R, Goldsmith D, Uchino S, et al. A prospective before-and-after trial of a medical emergency team. Med J Aust. 2003;179(6):283–7.
Jones D, Bellomo R, Bates S, et al. Long term effect of a medical emergency team on cardiac arrests in a teaching hospital. Crit Care. 2005;9(6):R808–15.
Sebat F, Musthafa AA, Johnson D, et al. Effect of a rapid response system for patients in shock on time to treatment and mortality during 5 years. Crit Care Med. 2007;35(11):2568–75.
Jacobs S, Zuleika M, Mphansa T. The Multiple Organ Dysfunction Score as a descriptor of patient outcome in septic shock compared with two other scoring systems. Crit Care Med. 1999;27(4):741–4.
Vanderschueren S, De WA, Malbrain M, et al. Thrombocytopenia and prognosis in intensive care. Crit Care Med. 2000;28(6):1871–6.
Bota DP, Melot C, Lopes FF, Nguyen Ba V, Vincent JL. The multiple Organ Dysfunction Score (MODS) versus the Sequential Organ Failure Assessment (SOFA) in outcome prediction. Intensive Care Med. 2002;28(11):1619-24.
Khwannimit B. Serial evaluation of the MODS, SOFA and LOD scores to predict ICU mortality in mixed critically ill patients. J Med Assoc Thai. 2008;91(9):1336–42.
Dover M, Tawfick W, Hynes N. Sultan S. Vascular: Evaluation of illness severity scoring systems and risk prediction in vascular intensive care admissions; 2015.
Xiao K, Guo C, Su L, Yan P, Li X, Xie L. Prognostic value of different scoring models in patients with multiple organ dysfunction syndrome associated with acute COPD exacerbation. J Thorac Dis. 2015;7(3):329–36.
Peres BD, Melot C, Lopes FF, Nguyen BV, Vincent JL. The Multiple Organ Dysfunction Score (MODS) versus the Sequential Organ Failure Assessment (SOFA) score in outcome prediction. Intensive Care Med. 2002;28(11):1619–24.
Haagensen R, Jamtli B, Moen H, Stokland O. Experiences with scoring systems SAPS II and NEMS for registration of activities in an intensive care unit. Tidsskr Nor Laegeforen. 2001;121(6):687–90.
Bouch DP, Thompson JP. Severity scoring systems in the critically ill. Contining Educ Anaesthesia Crit Care Pain. 2008;8(5):181–5. Ref Type: Generic.
Wagner DP, Knaus WA, Draper EA. Statistical validation of a severity of illness measure. Am J Public Health. 1983;73(8):878–84.
Charlson ME, Pompei P, Ales KL, MacKenzie CR. A new method of classifying prognostic comorbidity in longitudinal studies: development and validation. J Chronic Dis. 1987;40(5):373–83.
Needham DM, Scales DC, Laupacis A, Pronovost PJ. A systematic review of the Charlson comorbidity index using Canadian administrative databases: a perspective on risk adjustment in critical care research. J Crit Care. 2005;20(1):12–9.
Paul E, Bailey M, Pilcher D. Risk prediction of hospital mortality for adult patients admitted to Australian and New Zealand intensive care units: development and validation of the Australian and New Zealand Risk of Death model. J Crit Care. 2013;28(6):935–41.
Brinkman S, Bakhshi-Raiez F, Abu-Hanna A, et al. External validation of Acute Physiology and Chronic Health Evaluation IV in Dutch intensive care units and comparison with Acute Physiology and Chronic Health Evaluation II and Simplified Acute Physiology Score II. J Crit Care. 2011;26(1):105–8.
Nierman DM, Schechter CB, Cannon LM, Meier DE. Outcome prediction model for very elderly critically ill patients. Crit Care Med. 2001;29(10):1853–9.
Zhu BP, Lemeshow S, Hosmer DW, Klar J, Avrunin J, Teres D. Factors affecting the performance of the models in the Mortality Probability Model II system and strategies of customization: a simulation study. Crit Care Med. 1996;24(1):57–63.
Odetola FO, Rosenberg AL, Davis MM, Clark SJ, Dechert RE, Shanley TP. Do outcomes vary according to the source of admission to the pediatric intensive care unit? Pediatr Crit Care Med. 2008;9(1):20–5.
Escarce JJ, Kelley MA. Admission source to the medical intensive care unit predicts hospital death independent of APACHE II score. JAMA. 1990;264(18):2389–94.
Knaus WA, Wagner DP, Zimmerman JE, Draper EA. Variations in mortality and length of stay in intensive care units. Ann Intern Med. 1993;118(10):753–61.
Hill AD, Vingilis E, Martin CM, Hartford K, Speechley KN. Interhospital transfer of critically ill patients: demographic and outcomes comparison with nontransferred intensive care unit patients. J Crit Care. 2007;22(4):290–5.
Kramer AA, Higgins TL, Zimmerman JE. The association between ICU readmission rate and patient outcomes. Crit Care Med. 2013;41(1):24–33.
Choi S, Lee J, Shin Y, et al. Effects of a medical emergency team follow-up programme on patients discharged from the medical intensive care unit to the general ward: a single-centre experience. J Eval Clin Pract. 2015;16.
Niven DJ, Bastos JF, Stelfox HT. Critical care transition programs and the risk of readmission or death after discharge from an ICU: a systematic review and meta-analysis. Crit Care Med. 2014;42(1):179–87.
Pattison N, O’Gara G, Wigmore T. Negotiating transitions: involvement of critical care outreach teams in end-of-life decision making. Am J Crit Care. 2015;24(3):232–40.
Jones D, Moran J, Winters B, Welch J. The rapid response system and end-of-life care. Curr Opin Crit Care. 2013;19(6):616–23.
We are grateful to Mrs. Kathryn Barber, Health Information Analyst at the LHSC Health Records department matched the CCIS dataset with the hospital records to obtain the ICD-10-CA for calculation of the Charlson’s Co-morbidity Index.
No financial supports were utilized for this study.
The authors declare that they have no competing interests.
RK conceived and designed the study and performed the analysis and interpretation of the data and drafting and revising of the manuscript. FP participated in the analysis and interpretation of the data and revision of the manuscript critically for important intellectual content. AD carried out revising of the manuscript critically for important intellectual content and has given the final approval of the version to be published. All authors read and approved the final manuscript.
Nine equivalents of nursing manpower use score by Miranda DR, Nap R, de RA, Schaufeli W, Iapichino G. Nursing activities score. Crit Care Med 2003;31:374–382. (DOCX 87 kb)
Selection and categorization of the independent variables for the logistic regression model associated with ICU mortality. (DOCX 14 kb)
Bootstrap sampling (500 and 1000 times) to estimate the optimism in the AUC (c-statistics) based on the predicted outcomes from a logistic regression model. (DOCX 13 kb)
Multivariable logistic regression analysis of the training set (N = 4321) for Critical Care Trauma Center (CCTC) and Medical-Surgical Intensive Care Unit (MSICU) at London Health Sciences Center (LHSC) with MODS or NEMS only. (DOCX 14 kb)
About this article
Cite this article
Kao, R., Priestap, F. & Donner, A. To develop a regional ICU mortality prediction model during the first 24 h of ICU admission utilizing MODS and NEMS with six other independent variables from the Critical Care Information System (CCIS) Ontario, Canada. j intensive care 4, 16 (2016) doi:10.1186/s40560-016-0143-6
- Intensive care unit (ICU)
- Logistic regression
- Hosmer and Lemeshow goodness-of-fit test