Skip to main content

Describing organ dysfunction in the intensive care unit: a cohort study of 20,000 patients



Multiple organ dysfunction is a common cause of morbidity and mortality in intensive care units (ICUs). Original development of the Sequential Organ Failure Assessment (SOFA) score was not to predict outcome, but to describe temporal changes in organ dysfunction in critically ill patients. Organ dysfunction scoring may be a reasonable surrogate outcome in clinical trials but further exploration of the impact of case mix on the temporal sequence of organ dysfunction is required. Our aim was to compare temporal changes in SOFA scores between hospital survivors and non-survivors.


We performed a population-based observational retrospective cohort study of critically ill patients admitted from January 1, 2004, to December 31, 2013, to 4 multisystem adult intensive care units (ICUs) in Calgary, Canada. The primary outcome was temporal changes in daily SOFA scores during the first 14 days of ICU admission. SOFA scores were modeled between hospital survivors and non-survivors using generalized estimating equations (GEE) and were also stratified by admission SOFA (≤ 11 versus > 11).


The cohort consisted of 20,007 patients with at least one SOFA score and was mostly male (58.2%) with a median age of 59 (interquartile range [IQR] 44–72). Median ICU length of stay was 3.5 (IQR 1.7–7.5) days. ICU and hospital mortality were 18.5% and 25.5%, respectively. Temporal change in SOFA scores varied by survival and admission SOFA score in a complicated relationship. Area under the receiver operating characteristic (ROC) curve using admission SOFA as a predictor of hospital mortality was 0.77. The hospital mortality rate was 5.6% for patients with an admission SOFA of 0–2 and 94.4% with an admission SOFA of 20–24. There was an approximately linear increase in hospital mortality for SOFA scores of 3–19 (range 8.7–84.7%).


Examining the clinical course of organ dysfunction in a large non-selective cohort of patients provides insight into the utility of SOFA. We have demonstrated that hospital outcome is associated with both admission SOFA and the temporal rate of change in SOFA after admission. It is necessary to further explore the impact of additional clinical factors on the clinical course of SOFA with large datasets.


Multiple organ dysfunction syndrome (MODS)—also variably described as multiple organ failure or multisystem organ failure—is common in intensive care unit (ICU) patients [1]. It is often present at the time of ICU admission and/or present at time of ICU death. MODS has been described as a process of sequential and progressive organ dysfunction and/or failure that precedes death, its presence not infrequently influencing discussions and decisions to limit or withdraw life-sustaining interventions [2,3,4]. The pathophysiology of MODS remains poorly elucidated [1, 5,6,7,8]. In the absence of clearly defined pathophysiology for MODS, clinical diagnostic scoring systems have been developed to assist diagnosis, to describe progression, and to assess prognosis [9].

Though multiple scoring systems have been developed, the most common in practical use today is the Sequential Organ Failure Assessment (SOFA) score [10,11,12]. It was developed by consensus, in contrast to other scores that used statistical techniques to develop models (Logistic Organ Dysfunction System, Multiple Organ Dysfunction Score) [13]. The SOFA score, as with other organ dysfunction scoring systems, does not differentiate between chronic organ dysfunction present due to underlying disease, and the effects of acute organ dysfunction related to critical illness [14].

The original article reporting the development of the SOFA score described two main roles for the score: (1) to describe the natural history of organ dysfunction and (2) to assess the effects of new therapies. SOFA scores have been reported in many clinical trials in critical care [15,16,17]. Many studies have examined the use of SOFA to predict mortality in critically ill patients [17,18,19,20,21]. As such, some investigators have suggested that SOFA scores can be used to inform triage decisions during disaster management [22,23,24,25,26]. Some research has suggested modifications to the SOFA score for specific ICU patient populations such as trauma patients, patients with traumatic brain injury, neurocritical care patients, and many other conditions [27,28,29,30]. Few articles have described the clinical course of organ dysfunction using SOFA scores in populations of ICU patients, one of the original purported uses of such a score [18, 20, 31,32,33,34].

The objective of our study was to describe the clinical course of organ dysfunction among a heterogeneous cohort of critically ill patients admitted to multisystem adult ICUs.


The study was approved by the University of Calgary Conjoint Health Research Ethics Board, and a waiver of consent was granted (REB15-2010). This was a population-based observational retrospective cohort study, consisting of all consecutive critically ill adult patients (≥ 18 years) admitted to one of the four multisystem ICUs in Calgary, Canada, from January 1, 2004, to December 31, 2013. Patients were excluded if they were admitted to a coronary care unit, had undergone cardiac surgery, or survived ICU with a length of stay < 48 h after elective surgery. In the acute care hospitals in Calgary, some have surgical step-down or high observation units that admit elective postoperative patients for ‘routine’ monitoring. Therefore, patients admitted to ICU for postoperative monitoring following uncomplicated elective surgery with an ICU length of stay in the ICU < 48 h were excluded [19].


Demographic, clinical, severity of illness, diagnostic, and outcomes data were collected prospectively using an electronic medical record (TRACER). The details of this database have been described elsewhere [35, 36]. In summary, TRACER is a longitudinal database incorporating clinical and outcome data derived from electronic health records for adult ICU patients in the Calgary area. Data is received from various source systems, including the Quantitative Sentinel (GE-Marquette Medical) and eCritical MetaVision (iMDsoft) clinical information systems, Millennium (Cerner Millenium Pathnet) and Sunquest (Sunquest Information) laboratory information systems, and Alberta Health Services corporate databases. This data repository exists to provide data for research, quality and process improvement, as well as to enhance the provision of care to the critically ill.

Only the patient’s first admission to the ICU during the study period and the first 14 days of a patient’s ICU stay was considered. Patients were followed from the day of first ICU admission until discharge from ICU, end of the first 14 days of ICU stay, or death, which ever occurred first.

The calculation of the SOFA score requires the determination every 24 h of the most abnormal physiologic, laboratory, or clinical data within each of the six organ systems (central nervous system (CNS), respiratory, hepatic, renal, coagulation, and cardiovascular). Daily calculation of CNS dysfunction was based on the Glasgow Coma Scale (GCS) which was determined daily by the attending physician as either the actual GCS, in the absence of sedating medications, or the best estimate of GCS without the effects of sedation medications if they were being used. A score from 0 (normal) to 4 (most abnormal) was assigned for each organ system for a minimum SOFA score of 0 to a maximum SOFA score of 24. We followed the exact definitions of the original publication and have previously validated our collection including an audit and manual validation [10, 31]. For each patient, the admission SOFA score was calculated based on the first 24 h of data (therefore if a patient was admitted at 02:00 am, the admission SOFA would be based on the 02:00 am to 01:59 am time period), and thereafter, scores were calculated for each subsequent patient day based on the 07:00 am to 06:59 am time frame. With this method, there may be overlap in data used to calculate the admission SOFA score and the second SOFA score. If a patient was admitted between midnight and 6:59 am, a daily SOFA score in addition to the admission SOFA was not calculated for the short time period between the admission time and 6:59 am. Among all analyzed patient days, 14.5% of patient days were missing a daily SOFA score; however, the majority (81.8%) of those days were the last day of the patient’s stay (with patients most often being in ICU less than half a day on their last day) and would therefore be unlikely to alter results significantly. Missing values were not imputed. For each patient, we also calculated the maximum SOFA score which was the sum of the highest individual organ scores irrespective of day collected and the highest SOFA score which was the highest cumulative 24 h score collected during a patient’s ICU stay.

Statistical analysis

Baseline characteristics were summarized with medians with interquartile ranges (IQR) for continuous data and frequencies with proportions for categorical data. For the primary objective, marginal models for longitudinal data were developed using generalized estimating equations (GEE) to model the population-averaged change in mean daily SOFA score within the first 14 days of ICU stay between survivors and non-survivors. Subgroup analyses were performed using strata defined by admission SOFA (≤ 11 vs. > 11). We stratified admission SOFA at a value of 11 because prior studies have suggested a high mortality for those with an admission SOFA > 11, including Ferreira et al. who showed a predicted mortality for admission SOFA > 11 of 95% [18] and Anami et al. who observed a mortality of 89% for admission SOFA > 11 [37]. We conducted additional subgroup analyses examining temporal changes in practice over time using GEE models comparing patients admitted in 2004–2008 to patients admitted in 2009–2013 and by ICU LOS (< 3 days versus 3 to < 7 days versus ≥7 days). Due to mean daily SOFA scores potentially following a nonlinear trend, a piecewise linear GEE model with a threshold of 8 was fitted to allow different rates of change in SOFA for days 1–7 and days 8–14. A threshold of 8 was used based on observing similar results from (1) a nonlinear GEE model; (2) a GEE model using splines with knots at days 2, 3, 5, 7, 9, 11, and 13; and (3) a GEE model with day entered a categorical variable. An exchangeable within-group correlation structure and a robust estimator of variance were assumed. Hospital mortality percentages by early change in daily SOFA scores were compared using chi-squared tests. Area under the receiver operating characteristic (ROC) curve and precision-recall (PR) curve as well as the Brier score were used to quantify the predictive ability of SOFA scores for hospital mortality. Area under the ROC curves were compared using the method described in DeLong et al. [38]. Selected sensitivity analyses were repeated for ICU mortality. Analyses were conducted in R (version 3.5.1), a language and environment for statistical computing and graphics (The R Project for Statistical Computing, [39]. We used R packages geepack (function: geeglm) [40], doBy (function: esticon) [41], ROCR (functions: performance and prediction) [42], pROC (functions: roc and roc.test) [43], and PRROC (functions: roc.curve and pr.curve) [44, 45].


There were 20,068 adult patients admitted to ICUs in Calgary, Canada, between January 1, 2004, and December 31, 2013, among which 61 (0.3%) patients did not have at least one SOFA score recorded and were excluded from the analysis cohort. All 61 of these patients had an ICU length of stay < 7 h with 55 (90.2%) dying in hospital, of which 53 (96.4%) died in ICU. Demographics of the remaining cohort are shown in Table 1. Of the 20,007 patients with at least 1 measured SOFA score, 42% were female and median age was 59 (IQR 44–72) years. Median ICU and hospital length of stay were 3.5 (IQR 1.7–7.5) days and 14 (IQR 6–30) days, respectively. ICU mortality was 18.4% and hospital mortality was 25.2%.

Table 1 Patient demographics

The median admission SOFA score among hospital survivors and non-survivors were 6 (IQR 3–8) and 10 (IQR 7–13), respectively (Table 1). The most common admission diagnostic category among survivors was respiratory (28.4%) while it was cardiovascular for non-survivors (31.8%). The greatest difference in single organ dysfunction at admission between survivors and non-survivors was in the CNS and cardiovascular systems, followed by respiratory and renal systems. The median maximum SOFA score and highest SOFA were 7 (IQR 4–10) and 6 (IQR 4–9) for survivors and 13 (IQR 9–16) and 11 (IQR 8–15) for non-survivors. The highest SOFA most often occurred on day 1 or 2 for survivors and days 1–3 for non-survivors.

Figure 1 demonstrates mortality for admission, highest, and maximum SOFA scores. Across all 3 measures of organ failure, an approximately linear increase was observed with the exception of the extremes of low and high scores. The hospital mortality rate was 5.6% for patients with an admission SOFA score of 0–2. Not surprisingly, very high scores (20–24) at admission were associated with very low survival (5.6%), but not certainty of death. There was an approximately linear increase in hospital mortality for SOFA scores of 3–19 (range 8.7% to 84.7%). Patients with significant organ dysfunction with a highest SOFA score of 15–19 had a survival of 22.9%. When the CNS component was excluded, a similar trend was observed suggesting overestimation of organ failure due to GCS as suggested in previous work [46] was unlikely (Fig. 2).

Fig. 1
figure 1

Hospital mortality by a admission, b maximum, and c highest SOFA scores. SOFA, Sequential Organ Failure Assessment

Fig. 2
figure 2

Hospital mortality by a admission, b maximum, and c highest SOFA scores excluding CNS component. CNS, central nervous system; SOFA, Sequential Organ Failure Assessment

Figure 3 illustrates the clinical course of SOFA for survivors and non-survivors overall and stratified by admission SOFA, admission period, and ICU LOS. Table 2 shows that among those with admission SOFA > 11 during days 1–7 of ICU stay, rate of change in SOFA per day was − 1.02 (95% CI − 1.06, − 0.98) for survivors (p < 0.001) and − 0.52 (95% CI − 0.57, − 0.47) for non-survivors (p < 0.001) while for those with admission SOFA ≤ 11, it was − 0.47 (95% CI − 0.48, − 0.46) for survivors (p < 0.001) and − 0.13 (95% CI − 0.16, − 0.10) for non-survivors (p < 0.001). In non-survivors, we were unable to detect a change in SOFA score beyond day 7, whereas survivors had a slower but continued progressive improvement in organ function. The clinical course of SOFA for patients admitted in 2004–2008 and those admitted in 2009–2013 showed similar results (Fig. 3c, Table 2). Among those with an ICU LOS ≥ 7 days, similar trends were shown when compared to the overall cohort (Fig. 3d, Table 2). Nonlinear quadratic, spline, and categorical GEE models showed similar trends to the piecewise linear GEE models (Fig. 4). Figure 5 illustrates the clinical course of SOFA for ICU survivors compared to ICU non-survivors.

Fig. 3
figure 3

ad SOFA scores by hospital mortality overall and by subgroups (admission SOFA, period, and LOS). SOFA, Sequential Organ Failure Assessment; LOS, length of stay

Table 2 Change in daily SOFA score for hospital survivors and non-survivors overall and by subgroups
Fig. 4
figure 4

ac SOFA scores by hospital mortality overall using a nonlinear, b spline, and c categorical GEE models. The points represent results from the piecewise linear GEE model. SOFA, Sequential Organ Failure Assessment

Fig. 5
figure 5

ad SOFA scores by ICU mortality overall and by subgroups (admission SOFA, period, and LOS). SOFA, Sequential Organ Failure Assessment; LOS, length of stay

Among patients with at least 2 SOFA scores, 3502 (22.2%) patients had an increase in their SOFA score from admission of ≥ 2 points occurring at any time during their ICU admission. Hospital mortality of patients who experienced this ≥ 2 point increase was significantly higher when compared to those who did not (40.3% vs 17.5%, Χ2 = 811, p < 0.001) (Table 3, Fig. 6, Additional file 1: Table S1). Similarly, ICU mortality of patients who experienced this increase was significantly higher than those who did not (32.1% vs 9.9%, χ2 = 1061, p < 0.001) (Table 4, Fig. 7, Additional file 1: Table S1). Additionally, patients with an increase from days 1 to 5 had a significantly higher hospital and ICU mortality than those who did not (hospital mortality: 40.5% vs. 20.7%, χ2 = 257, p < 0.001; ICU mortality: 31.6% vs. 11.2%, χ2 = 398, p < 0.001).

Table 3 Hospital mortality by early change in daily SOFA scores
Fig. 6
figure 6

Hospital mortality by early change in daily SOFA scores. SOFA, Sequential Organ Failure Assessment

Table 4 ICU mortality by early change in daily SOFA scores
Fig. 7
figure 7

ICU mortality by early change in daily SOFA scores. SOFA, Sequential Organ Failure Assessment

Among all 20,007 patients, area under the ROC curve (area under the PR curve in brackets) using admission SOFA as a predictor of hospital mortality was 0.77 (0.57) while it was 0.79 (0.58) using maximum SOFA and 0.81 (0.63) using highest SOFA (Table 5). For ICU mortality, the values were 0.80 (0.53) for admission SOFA, 0.82 (0.54) for maximum SOFA and 0.85 (0.60) for highest SOFA. Area under the ROC curve (area under the PR curve in brackets) was 0.81 (0.62) using the admission APACHE II score for hospital mortality and 0.83 (0.58) for ICU mortality. Among patients with SOFA scores collected on both days 1 (admission) and 2 of their ICU stay (n = 15,736), area under the ROC curve was 0.73 using admission SOFA alone, while the difference between days 1 and 2 improved prediction to an area under the ROC curve of 0.77 (Z = 13.7, p < 0.001). The difference between the admission SOFA and maximum SOFA as well as the difference between the admission SOFA and highest SOFA were also improved with the addition of the admission SOFA (Z = 13.6, p < 0.001 and Z = 17.5, p < 0.001, respectively). For ICU mortality, area under the ROC curve was 0.77 using admission SOFA alone while the difference between days 1 and 2 improved the area under the ROC curve to 0.81 (Tables 5 and 6). The ROC and PR curves for hospital and ICU mortality are presented in Fig. 8.

Table 5 Area under the ROC and PR curves for predictors of hospital mortality
Table 6 Area under the receiver operating characteristic curves for predictors of ICU mortality
Fig. 8
figure 8

ROC and PR curves for predictors of hospital and ICU mortality. ROC, receiver operating characteristic; PR, precision-recall


Our data demonstrate that MODS does have a clinical course that can be described on a population basis though with significant individual variation in trajectory. It also demonstrates that this course is different in hospital non-survivors and survivors, and varies in the first 7 days and the next 7 days based on the severity of admission organ dysfunction. Although using a piecewise linear regression model was a simplification to ease in interpretation, less restrictive models showed similar trends. Among survivors, there is less severe organ dysfunction at admission, and a more rapid improvement in SOFA score compared to non-survivors. The rapidity of daily improvement in SOFA scores over the first 7 days is greater in the cohort with an admission SOFA > 11 compared to < 11. However, from days 8–14, there is no statistically significant difference in improvement in daily SOFA scores in non-survivors, compared to survivors. Most patients had their highest SOFA score in the first 2 days of ICU admission. This may suggest that if “multiple hits” are responsible for MODS, hits are not associated with a clinical course of failure-resolution-failure [5].

Our results are similar to those of Moreno et al. [19]. Their study included 1449 patients from 40 ICUs over 4 continents and demonstrated that maximum SOFA had a higher discriminative ability than the admission SOFA or the Delta SOFA (area under the ROC curve (AUC) of 0.847 vs 0.772 and 0.742, respectively). The relationship of maximum SOFA to hospital mortality in our study was essentially similar to the relationship observed of maximum SOFA to ICU mortality from the Moreno study (AUC 0.82 vs 0.847). One of the recommendations by Moreno et al. was that SOFA scores be compared not only against ICU outcome but with longer outcomes such as 30-day or hospital mortality. Our results complement the earlier findings of Moreno et al. and respond to what they consider to be limitations in their own observations. There was a less of a relationship observed between delta SOFA and mortality in our study and that of Moreno (AUC 0.58 vs 0.742). These differences might reflect differences in the study cohorts, or potential overlap in the calculation of our admission and day 1 SOFA, or that the relationship of delta SOFA to hospital outcome is not as relevant as the relationship with maximum SOFA. It may also suggest that overall trajectory within ICU is a better descriptor of the relationship with hospital mortality. Bingold et al.’s [20] examination of SOFA scores were limited to scores on admission and at 3 days after admission. While Badawi et al. [34] was a much larger study, our patient cohort differed in illness severity with a mean SOFA score of 7.0 in our cohort and a mean SOFA score of 4.5 in their cohort. Additionally, ICU mortality (18.5% vs 4.8%) and hospital mortality (25.5% vs 8.1%) were significantly higher in our patient cohort, and our data was provided beyond 7 days in ICU. These results are also similar to data from Badreldin and colleagues who suggested that a “daily-mean-SOFA” (an average daily SOFA from day 1 to up to day 6—essentially a slope) was more accurate in describing outcome than admission or Max SOFA scores [32]. Holder et al. [33] showed that early serial organ failure scores up to 5 days improve prediction of ICU mortality; however, ICU and hospital mortality rates were very low (6.4% and 11.2%, respectively) in comparison to our study (18.5% and 25.5%, respectively), and coronary care units were included. In spite of including a different study population, the main results from our current study are remarkably similar to those from an earlier study by our group [31]. The current study included 20,007 patients admitted to ICUs in Calgary, Canada, over 10 years (2004–2013), whereas our earlier study included only 1436 patients admitted to Calgary ICUs over 1 year from May 1, 2000, to April 30, 2001. The main difference between the current study and our previous study is the examination in the current study of how SOFA scores change early during ICU admission and the association of these changes with survival. Taken together with our previous study, we have demonstrated that multiple organ dysfunction does not follow a course of progressive and sequential failure. Both studies have also shown that the clinical course of organ dysfunction differs between hospital survivors and non-survivors, but we were additionally able to show this difference varies by admission SOFA and the temporal change in the first week of the patient’s ICU stay. This information is important as it has potential implications for evaluation of therapeutic interventions in the critically ill as has been suggested in recent critical care literature [15, 47,48,49,50]. Temporal change in organ dysfunction could be considered as a measure for guiding therapy. For example, in a clinical trial, the absence of improvement in organ dysfunction over a number of days might suggest the absence of benefit, whereas the rate of organ improvement as measured by SOFA might determine the duration of therapy. Decisions to terminate or duration of therapy might be particularly relevant for therapies with high acquisition costs, or for therapies with potentially high risk from adverse side effect profiles.

In addition to demonstrating its potential utility as a tool to guide therapy, our study suggests that a certain “fatalism” about the severity of organ dysfunction may not be justified [18, 19, 51, 52]. Patients in our study with SOFA admission, SOFA highest, or SOFA maximum scores in the 15–19 range and even 20 or above did not have a certain outcome of death (e.g., survival 22.9% for admission SOFA of 15–19 and 5.6% for admission SOFA of 20–24). Our current results expand on our previous observations and conclusions. Careful consideration of other factors, such as admitting diagnosis, frailty, and underlying burden of illness, though not part of our study, could reasonably be evaluated in future work and should be considered in the meantime prior to applying SOFA scores as means of “triaging” outcome.

Our study is not without limitations. First is a potential limitation and criticism of the SOFA score. The original score and cutoff points within each organ system were developed by consensus rather than by mathematical modelling where the interval between ordinal scores purports a proportional increase in the risk of death. Our results do not support such criticism. Whether examining admission, highest, or maximum SOFA, there is an approximately linear increase in the risk of death as SOFA rises. The association of organ dysfunction with outcomes is complex. Although we have reported predicted probabilities of mortality associated with different measures of SOFA (i.e., maximum SOFA), we do not consider that simple mathematical models can fully explain the association, nor are we confident that organ failure scores will provide robust estimates for prognosis. Rather, we believe that this information is descriptive and must be placed in the context of a patient’s individual course. Alternate methods of examining the association have also not provided results that can be used with confidence for predicting outcome [53,54,55,56,57,58]. A major limitation of this study is the presence of informative censoring of patients who die and are discharged that is dependent on organ function. This was shown in the analysis separating patients into three groups based on ICU length of stay where the profiles differed between the three groups of patients. Further work needs to be done to more closely examine this strong feature of our data. Another limitation arises from our subgroup analyses. Grouping SOFA scores into two categories (≤ 11 versus > 11) is a form of extreme rounding. Although this was done based on previous studies showing increased mortality with SOFA scores > 11 and to show that patients admitted to ICU with a greater degree of organ dysfunction behave differently than those who have low organ dysfunction, caution is advised as simple dichotomization can result in spurious inferences in interpretation.


Multiple organ dysfunction syndrome is a common problem experienced by many patients admitted to ICUs. The severity of organ dysfunction is usually worst at admission. Death is more common in those with more severe admission SOFA scores. There are clear differences in the rates of improvement in SOFA between survivors and non-survivors with the most rapid improvement in survivors. A lack of improvement in SOFA scores beyond 7 days is the temporal course experienced by most non-survivors of critical illness.



Confidence interval


Central nervous system


Glasgow Coma Scale


Generalized estimating equations


Intensive care unit


Interquartile range


Multiple organ dysfunction syndrome




Receiver operating characteristic


Sequential Organ Failure Assessment


  1. Brealey D, Singer M. Multi-organ dysfunction in the critically ill: epidemiology, pathophysiology and management. J R Coll Physicians Lond. 2000;34(5):424–7.

    CAS  PubMed  Google Scholar 

  2. Baue AE. Multiple, progressive, or sequential systems failure: a syndrome of the 1970s. Arch Surg. 1975;110:779–81.

    Article  CAS  Google Scholar 

  3. Baue AE. Multiple organ failure, multiple organ dysfunction syndrome, and systemic inflammatory response syndrome: why no magic bullets? Arch Surg. 1997;132:703–7.

    Article  CAS  Google Scholar 

  4. Heard S, Fink M. Multiple organ failure syndrome. Part 1: epidemiology, prognosis, and pathophysiology. J Int Care Med. 1991;6:279–94.

    Article  Google Scholar 

  5. Deitch EA. Multiple organ failure: pathophysiology and potential future therapy. Ann Surg. 1992;216:117–34.

    Article  CAS  Google Scholar 

  6. Seely AJE, Christou NV. Multiple organ dysfunction syndrome: exploring the paradigm of complex nonlinear systems. Crit Care Med. 2000;28:2193–000.

    Article  CAS  Google Scholar 

  7. Marshal JC. Critical illness is an iatrogenic disorder. Crit Care Med. 2010;38 (suppl:S582–9.

    Article  Google Scholar 

  8. Matsuda N, Hattori Y. Systemic inflammatory response syndrome (SIRS): molecular pathophysiology and gene therapy. J Pharmacol Sci. 2006;101:189–98.

    Article  CAS  Google Scholar 

  9. Baue AE. MOF, MODS, and SIRS: what is in a name or an acronym? Shock. 2006;26:438–49.

    Article  CAS  Google Scholar 

  10. Vincent J-L, Moreno R, Takala J, Willatts S, De Mendonça A, Bruining H, et al. The SOFA (sepsis-related organ failure assessment) score to describe organ dysfunction/failure. Intensive Care Med. 1996;22:707–10.

    Article  CAS  Google Scholar 

  11. 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:1638–52.

    Article  CAS  Google Scholar 

  12. Le Gall JR, Lkar J, Lemeshow S, Saulnier F, Alberti C, Artigas A, et al. The logistic organ dysfunction system: a new way to assess organ dysfunction in the intensive care unit. ICU scoring group. JAMA. 1996;276:802–10.

    Article  Google Scholar 

  13. Zygun DA, Doig CJ. Measuring organ dysfunction. In: Vincent J-L, editor. Yearbook of intensive care and emergency medicine. Berlin: Springer-Verlag; 2002.

    Google Scholar 

  14. Lee KS, Sheen SS, Jung YJ, Park RW, Lee YJ, Chung WY, et al. Consideration of additional factors in sequential organ failure assessment score. J Crit Care. 2014;29(1):185.e9–185.e12.

    Article  Google Scholar 

  15. de Grooth HJ, Geenen IL, Girbes AR, Vincent JL, Parienti JJ, Oudemans-van Straaten HM. SOFA and mortality endpoints in randomized controlled trials: a systematic review and meta-regression analysis. Crit Care. 2017;21:38.

    Article  Google Scholar 

  16. Marshall J. Organ dysfunction as an outcome measure in clinical trials. Eur J Surg Suppl. 1999;584:62–7.

    Article  Google Scholar 

  17. 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.

    Article  Google Scholar 

  18. Ferreira FL, Bota DP, Bross A, Mélot C, Vincent JL. Serial evaluation of the SOFA score to predict outcome in critically ill patients. JAMA. 2001;286(14):1754–8.

    Article  CAS  Google Scholar 

  19. Moreno R, Vincent JL, Matos R, Mendonça A, Cantraine F, Thijs L, et al. The use of maximum SOFA score to quantify organ dysfunction/failure in intensive care. Results of a prospective, multicenter study. Int Care Med. 1999;25:686–96.

    Article  CAS  Google Scholar 

  20. Bingold TM, Lefering R, Zacharowski K, Meybohm P, Waydhas C, Rosenberger P, et al. Individual organ failure and concomitant risk of mortality differs according to the type of ICU admission to ICU, a retrospective study of SOFA score of 23,795 patients. PLoS One. 2015;10(8):e0134329.

    Article  Google Scholar 

  21. Zygun DA, Laupland KB, Fick GH, Sandham JD, Doig CJ. Limited ability of SOFA and MOD scores to discriminate outcome: a prospective evaluation in 1,436 patients. Can J Anaesth. 2005;52(3):302–8.

    Article  Google Scholar 

  22. Christian MD, Hawryluck L, Wax RS, Cook T, Lazar NM, Herridge MS, et al. Development of a triage protocol for critical care during an influenza pandemic. CMAJ. 2006;175:1377–81.

    Article  Google Scholar 

  23. Khan Z, Hulme J, Sherwood N. An assessment of the validity of SOFA score based triage in H1N1 critically ill patients during an influenza pandemic. Anaesthesia. 2009;64:1283–8.

    Article  CAS  Google Scholar 

  24. Sprung CL, Zimmerman JL, Christian MD, Joynt GM, Hick JL, Taylor B, et al. Recommendations for intensive care unit and hospital preparations for an influenza epidemic or mass disaster. Summary report of the European Society of Intensive Care Medicine’s Task Force for intensive care unit triage during an influenza epidemic or mass disaster. Int Care Med. 2010;36:428–43.

    Article  Google Scholar 

  25. Shahpori R, Stelfox HT, Doig CJ, Boiteau PJ, Zygun DA. Sequential organ failure assessment in H1N1 pandemic planning. Crit Care Med. 2011;39(4):827–32.

    Article  Google Scholar 

  26. Kanter RK. Would triage predictors perform better than first-come, first-served in pandemic ventilator allocation? Chest. 2015;147(1):102–8.

    Article  Google Scholar 

  27. Hutchings L, Watkinson P, Young JD, Willett K. Defining multiple organ failure after major trauma: a comparison of the Denver, Sequential Organ Failure Assessment, and Marshal scoring systems. J Trauma and Acute Care Surg. 2017;82(3):534–41.

    Article  Google Scholar 

  28. Zygun D, Berthiaume L, Laupland K, Kortbeek J, Doig C. SOFA is superior to MOD score for the determination of non-neurologic organ dysfunction in patients with severe traumatic brain injury: a cohort study. Crit Care. 2006;10(4):R115.

    Article  Google Scholar 

  29. Zygun DA, Kortbeek JB, Fick GH, Laupland KB, Doig CJ. Non-neurologic organ dysfunction in severe traumatic brain injury. Crit Care Med. 2005;33(3):654–60.

  30. Zygun DA, Doig CJ, Gupta AK, Whiting G, Nicholas C, Shepherd E, et al. Non-neurological organ dysfunction in neurocritical care. J Crit Care. 2003;18(4):238–44.

    Article  Google Scholar 

  31. Doig CJ, Zygun DA, Fick GH, Laupland KB, Boiteau PJ, Shahpori R, et al. Study of clinical course of organ dysfunction in intensive care. Crit Care Med. 2004;32(2):384–90.

    Article  Google Scholar 

  32. Badreldin A, Elsobky S, Lehmann T, Brehm BB, Doenst T, Hekmat K. Daily-mean-SOFA, a new derivative to increase accuracy of mortality prediction in cardiac surgical intensive care units. Thor Cardiovasc Surg. 2012;60(1):43–50.

    Google Scholar 

  33. Holder AL, Overton E, Lyu P, Kempker JA, Nemati S, Razmi F, et al. Serial daily organ failure assessment beyond ICU day 5 does not independently add precision to ICU risk-of-death prediction. Crit Care Med. 2017;45(12):2014–22.

    Article  Google Scholar 

  34. Badawi O, Liu X, Hassan E, Amelung PJ, Swami S. Evaluation of ICU risk models adapted for use as continuous markers of severity of illness throughout the ICU stay. Crit Care Med. 2018;46(3):361–7.

    Article  Google Scholar 

  35. Kramer AH, Zygun DA, Doig CJ, Zuege DJ. Incidence of neurologic death among patients with brain injury: a cohort study in a Canadian health region. CMAJ. 2013;185(18):E838–45.

    Article  Google Scholar 

  36. Brundin-Mather R, Soo A, Zuege D, Niven DJ, Fiest K, Doig CJ, et al. Secondary EMR data for quality improvement and research: a comparison of manual and electronic data collection from an integrated critical care electronic medical record system. J Crit Care. 2018;47:295–301.

    Article  Google Scholar 

  37. Anami EH, Grion CM, Cardoso LT, Kauss IA, Thomazini MC, Zampa HB, Bonametti AM, Matsuo T. Serial evaluation of SOFA score in a Brazilian teaching hospital. Intensive Crit Care Nurs. 2010;26(20):75–82.

    Article  Google Scholar 

  38. DeLong ER, DeLong DM, Clarke-Pearson DL. Comparing the areas under two or more correlated receiver operating characteristic curves: a nonparametric approach. Biometrics. 1998;44(3):837–45.

    Article  Google Scholar 

  39. R Core Team. R: a language and environment for statistical computing. Vienna: R Foundation for Statistical Computing; 2018.

    Google Scholar 

  40. Højsgaard S, Halekoh U, Yan J. The R package geepack for generalized estimating equations. J Stat Softw. 2005;15(2):1–11.

    Google Scholar 

  41. Højsgaard S, Halekoh U (2018). doBy: Groupwise statistics, LSmeans, linear contrasts, utilities. R package version 4.6-2.

  42. Sing T, Sander O, Beerenwinkel N, Lengauer T. ROCR: visualizing classifier performance in R. Bioinformatics. 2005;21(20):3940–1.

    Article  CAS  Google Scholar 

  43. Robin X, Turck N, Hainard A, Tiberti N, Lisacek F, Sanchez JC, Müller M. pROC: an open-source package for R and S+ to analyze and compare ROC curves. BMC Bioinformatics. 2011;12:77.

    Article  Google Scholar 

  44. Knox DB, Lanspa MJ, Pratt CM, Kuttler KG, Jones JP, Brown SM. Glasgow coma scale dominates the association between admission sequential organ failure assessment score and 30-day mortality in a mixed intensive care unit population. J Crit Care. 2014;29(5):780–5.

    Article  Google Scholar 

  45. Keilwagen J, Grosse I, Grau J. Area under precision-recall curves for weighted and unweighted data. PLoS One. 2014;9(3):e92209.

    Article  Google Scholar 

  46. Grau J, Grosse I, Keilwagen J. PRROC: computing and visualizing precision-recall and receiver operating characteristic curves in R. Bioinformatics. 2015;31(15):2595–7.

    Article  CAS  Google Scholar 

  47. Russell J, Singer J, Bernard G, Wheeler A, Fulkerson W, Hudson L, et al. Changing patterns of organ dysfunction in early human sepsis is related to morbidity. Crit Care Med. 2000;28:3405–11.

    Article  CAS  Google Scholar 

  48. Prentice R. Surrogate endpoints in clinical trials: definitions and operational criteria. Stat Med. 1989;8:831–40.

    Article  Google Scholar 

  49. Rubenfeld GD. Surrogate measures of patient-centered outcomes in critical care. In: Angus DC, Carlet J, editors. Surviving intensive care. Berlin: Springer; 2003. p. 169–80.

    Chapter  Google Scholar 

  50. Vincent JL. Endpoints in sepsis trials: more than just 28-day mortality? Crit Care Med. 2004;32:S209–13.

    Article  Google Scholar 

  51. Kleiman DA, Barie PS. Survival in fully manifest multiple organ dysfunction syndrome. Surg Infect. 2014;15(4):445–9.

    Article  Google Scholar 

  52. Cabre L, Mancebo J, Solsona JF, Saura P, Gich I, Blanch L, et al. Multicenter study of the multiple organ dysfunction syndrome in intensive care units: the usefulness of sequential organ failure assessment scores in decision making. Int Care Med. 2005;31:927–33.

    Article  CAS  Google Scholar 

  53. Houthooft R, Ruyssinck J, van der Herten J, Stijven S, Couckuyt I, Gadeyne B, et al. Predictive modelling of survival and length of stay in critically ill patients using sequential organ failure scores. Artif Intell Med. 2015;63(3):191–207.

    Article  Google Scholar 

  54. Rehman MF, Siddiqui MS. Predicting death and disability, is it really possible? A medical ICU prognostication model study. Crit Care Med. 2014;42:2449–50.

    Article  Google Scholar 

  55. Meadow W, Pohlman A, Reynolds D, Rand L, Correia C, Christoph E, et al. Power and limitations of daily prognostication of death in the medical ICU for outcomes in the following 6 months. Crit Care Med. 2014;42:2387–92.

    Article  Google Scholar 

  56. Hsieh YZ, Su MC, Wang CH, Wang PC. Prediction of survival of ICU patients using computational intelligence. Comp Biol Med. 2014;47:13–9.

    Article  Google Scholar 

  57. Sandri M, Berchialla P, Baldi I, Gregori D, De Blasi RA. Dynamic Bayesian networks to predict sequences of organ failures in patients admitted to ICU. J Biomed Informatics. 2014;48:106–13.

    Article  Google Scholar 

  58. Sekulic AD, Trpkovic SV, Pavlovic AP, Marinkovic OM, Ilic AN. Scoring systems in assessing survival of critically ill patients. Med Sci Monit. 2015;21:2621–9.

    Article  Google Scholar 

Download references


The authors wish to acknowledge Malik Agyemang and Cathy Curr (Alberta Health Services in Alberta, Canada) for their help in data acquisition and extraction.


No funding was received by any authors for this study.

Availability of data and materials

The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request.

Author information

Authors and Affiliations



AS was responsible for the study design, data acquisition/analysis/interpretation, manuscript writing/revisions, and drafting of figures and tables. DJZ was responsible for the data interpretation/analysis and manuscript writing/revisions. GHF was responsible for the data interpretation/analysis and manuscript writing/revisions. DJN was responsible for the data interpretation/analysis and manuscript writing/revisions. LRB was responsible for the data interpretation/analysis and manuscript writing/revisions. HTS was responsible for the data interpretation/analysis and manuscript writing/revisions. CJD was responsible for the study conception and design, data acquisition/analysis/interpretation, and manuscript writing/revisions. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Andrea Soo.

Ethics declarations

Ethics approval and consent to participate

This study was approved by the Conjoint Health Research Ethics Board in Calgary, Alberta (REB approval 15-2010).

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Additional information

This work was performed at the University of Calgary, Calgary, Alberta, Canada

Additional file

Additional file 1:

Table S1. Exact p values for hospital and ICU mortality by early change in daily SOFA scores. (DOCX 19 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Soo, A., Zuege, D.J., Fick, G.H. et al. Describing organ dysfunction in the intensive care unit: a cohort study of 20,000 patients. Crit Care 23, 186 (2019).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: