Performance of the score systems Acute Physiology and Chronic Health Evaluation II and III at an interdisciplinary intensive care unit, after customization
© BioMed Central Ltd on behalf of the copyright holder 2001
Received: 5 June 1998
Accepted: 18 October 2000
Published: 5 January 2001
Mortality predictions calculated using scoring scales are often not accurate in populations other than those in which the scales were developed because of differences in case-mix. The present study investigates the effect of first-level customization, using a logistic regression technique, on discrimination and calibration of the Acute Physiology and Chronic Health Evaluation (APACHE) II and III scales.
Probabilities of hospital death for patients were estimated by applying APACHE II and III and comparing these with observed outcomes. Using the split sample technique, a customized model to predict outcome was developed by logistic regression. The overall goodness-of-fit of the original and the customized models was assessed.
Of 3383 consecutive intensive care unit (ICU) admissions over 3 years, 2795 patients could be analyzed, and were split randomly into development and validation samples. The discriminative powers of APACHE II and III were unchanged by customization (areas under the receiver operating characteristic [ROC] curve 0.82 and 0.85, respectively). Hosmer-Lemeshow goodness-of-fit tests showed good calibration for APACHE II, but insufficient calibration for APACHE III. Customization improved calibration for both models, with a good fit for APACHE III as well. However, fit was different for various subgroups.
The overall goodness-of-fit of APACHE III mortality prediction was improved significantly by customization, but uniformity of fit in different subgroups was not achieved. Therefore, application of the customized model provides no advantage, because differences in case-mix still limit comparisons of quality of care.
KeywordsAcute Physiology and Chronic Health Evaluation customization logistic regression mortality prediction severity of illness
Scoring systems are used in intensive care to control for various case-mix factors in order to compare patient populations. Score-based predictions of mortality in ICU patients may be used for quality assurance and comparison of quality of care [1,2,3]. If a scoring system is intended to be used in a patient population that is different from the original population used in the development of the system (development sample), then it should be validated in this new population [4,5,6,7].
Calibration measures how closely mortality prognosis fits the observed mortality. Poor calibration in a patient sample does not necessarily mean that the quality of care in that particular ICU is better or worse than in the development sample.
Several clinical case-mix as well as nonclinical factors are not accounted for by such scoring systems . The overall fit of a score in a particular patient sample can be improved by customization using logistic regression. This is possible for the whole population in question , but can also be done independently for specific subgroups [10,11]. In the latter case, the customized score can only be used in this subset of patients. If a customized model derived from the whole population is used, then uniformity of fit for the relevant subgroups should still be tested. Knowledge of the influence of subgroups is important, because future changes in case-mix may compromise the improvement achieved by customization.
The aim of the present study was to test the performance of APACHE II and III, after customization of these scales for use in future assessment of quality of care in our unit.
Patients and methods
Over a 3-year period (October 1991-October 1994), 3382 patients were consecutively admitted to the 12-bed interdiciplinary ICU of a 571-bed, university-affiliated community hospital. For the APACHE II analysis, 274 patients who were readmitted to the ICU, 208 patients who were in the ICU for less than 4 h, 16 patients who were admitted for dialysis only, two patients who were younger than 16 years and 87 patients with missing data were excluded. Thus, 2795 patients were included in the analysis. For APACHE III, 79 patients who were admitted to rule out myocardial infarction and 55 cardiosurgical patients were excluded, leaving 2661 for analysis.
Data collection was done according to the criteria and definitions described by the developers of APACHE II and III [12,13]. The data were collected by ward doctors after 4 weeks training in how to use the APACHE system. They had access to a detailed manual, including definitions and procedures. Constant supervision by a documentation assistant included regular comparison of the original with the collected data, and review of completeness. In order to assess reliability of data collection, data from a random sample of 50 patients were recorded by two data collectors independently. Interobserver reliability was analyzed by Kendall's coefficient of concordance and κ statistics. In addition, data collection software, which was provided by APACHE Medical Systems Inc (Washington, DC), automatically checked that the data were plausible. The whole data set was tested using a box-plot technique in order to analyze extreme values seperately. Vital status at hospital discharge was recorded.
The sample was split randomly into a development (n = 1863 for APACHE II and n = 1772 for APACHE III) and a validation sample (n = 932 for APACHE II and n =889 for APACHE III). Development of the original model by logistic regression  led to the following equation:
logit(or) = β0 + Σβi × xi (1)
where β0 is a constant, βi are coefficients, and xi encompasses the various patient factors that are included in the model. The probability of hospital death is calculated as follows:
P(or) = elogit(or)/1 + elogit(or) (2)
In the present study the APACHE II equation was used as indicated by the developers . The APACHE III equation was provided by APACHE Medical Systems Inc, and it has not been published for commercial reasons. In customizing the scales, the original logit was used as the independent variable and hospital death was used as the dependent variable. The new probability of hospital death was calculated as follows:
P(cust) = elogit(cust)/1 + elogit(cust) (3)
and logit(cust) is calculated as follows:
logit(cust) = βc0 + [βc1 × logit(or)] (4)
Customized coefficients used in the calculation of probability of hospital death
Discrimination and calibration were analyzed for the original and customized models. Discriminative power was tested by calculating the areas under the ROC curves , and calibration was calculated using standardized mortality ratio (SMR; observed deaths/expected deaths), with 95% confidence intervals , and using the Hosmer-Lemeshow goodness-of-fit H and C tests . Comparison of the development and validation samples was done by Χ2 test and Mann-Whitney U test, because values were not normally distributed.
Statistical analysis was performed using the SPSS 6.1 software package (SPSS, Chicago, IL, USA). P < 0.05 was considered statistically significant.
Completeness of data was good; excluding just one variable (24-h urine), 94.6% of all necessary data were collected on average for each patient; 24-h urine was available in only 78.1% of patients. Reliability analysis revealed Kendall's coefficients for clinical and laboratory data above 0.9 except for blood gas values (0.878) and 24-h urine (0.870). κ values were low only for diagnosis of renal failure (0.49) and Glasgow Coma Scale score (0.54). Despite that, differences in calculated scores were very low, with Kendall's coefficients above 0.92. Thus, overall reliability of data collection was good.
Demographic and clinical data for APACHE II sample
Age (mean ± SD)
62.5 ± 16.1
61.8 ± 16.3
Direct admission (%)
< 65 years (%)
≥ 65 years (%)
DOS ICU (mean ± SD)
3.95 ± 7.0
3.48 ± 5.4
DOS hospital (mean ± SD)
25.35 ± 22.4
25.26 ± 22.5
Hospital mortality (%)
APACHE II prediction (%)
APACHE III prediction (%)
Discrimination and calibration of APACHE II and III before and after customization
SMR (95% CI)
APACHE II original
APACHE II customized
APACHE III original
APACHE III customized
Discrimination and calibration of APACHE II before and after customization for various subgroups of the validation sample
SMR (95% CI)
< 65 years
≥ 65 years
Discrimination and calibration of APACHE III before and after customization for various subgroups of the validation sample
SMR (95% CI)
< 65 years
≥ 65 years
Customization of APACHE II and III in a large patient population from a single unit led to an improvement in the overall goodness-of-fit of APACHE III, which showed poor calibration in its original version. Despite a similar improvement of fit in several subgroups that were large enough to be tested, good uniformity of fit was not achieved.
These results are comparable with those of a large multi-center study  that analyzed customization of the Mortality Prediction Model. In that study, a second-level customization, in which new coefficients were developed for all single patient factors included in the original model, improved calibration even further. Second-level customization was not attempted in the present patient sample because there were not enough patients for that purpose. Time to collect data in a sufficiently large patient sample in a single unit would probably be so great that real changes in case-mix or ICU treatment might occur during the study, which would confound the results. First-level customization will probably be a more practical method for single units to improve the overall fit of score systems that are to be used for quality assessment.
At present, however, we would not recommend customization routinely. This is because a major problem is still unresolved; although good calibration can be achieved for the whole patient sample, uniformity of fit remains unsatisfactory. This is the case even for APACHE III, which accounts for more case-mix factors, such as diagnostic categories and lead time, than do the other models. Nevertheless, achievement of uniformity is important, because change in case-mix over time will otherwise lead to a loss of accuracy of a customized model. It would be difficult to interpret whether a change in the mortality ratio over time would be due to a change in quality of care or in case-mix.
If a customized model still has a poor fit for a certain subgroup at a specific unit, then customization for this sample can be attempted separately [10,11]. This could be attempted in medical and cardiovascular patients at our unit for APACHE III, because these groups are sufficiently large and because general customization did not lead to a good fit. However, the practicality of such an approach is questionable.
We conclude that the customized models of APACHE II and III did not provide an advantage over the original versions. First-level customization did not solve the problem of differences in case-mix, leading to different score results, which do not necessarily reflect differences in quality of care . Data collection is ongoing to test variation in the original and customized scores over time in a second large sample from our unit.
- Bodmann KF, Ehlers B, Häbel U, Ritschel P, Rühmkorf K: Epidemiologic and prognostic data of 2054 patients in a medical intensive care unit. Dtsch Med Wochenschr. 1997, 122: 919-925.PubMedView ArticleGoogle Scholar
- Knaus WA, Draper EA, Wagner DP, Zimmerman JE: An evaluation of outcome from intensive care in major medical centers. Ann Intern Med. 1986, 104: 410-418.PubMedView ArticleGoogle Scholar
- Shortell SM, Zimmerman JE, Rousseau DM, Gillies RR, Wagner DP, Draper EA, Knaus WA, Duffy J: The performance of intensive care units: does good management make a difference?. Med Care. 1994, 32: 508-525.PubMedView ArticleGoogle Scholar
- Beck DH, Taylor BL, Millar B, Smith GB: Prediction of outcome from intensive care: a prospective cohort study comparing Acute Physiology and Chronic Health Evaluation II and III prognostic systems in a United Kingdom intensive care unit. Crit Care Med. 1997, 25: 9-15. 10.1097/00003246-199701000-00006.PubMedView ArticleGoogle Scholar
- Moreno R, Miranda DR, Fidler V, Van Schilfgaarde R: Evaluation of two outcome prediction models on an independent database. Crit Care Med. 1998, 26: 50-61. 10.1097/00003246-199801000-00016.PubMedView ArticleGoogle Scholar
- Rowan KM, Kerr JH, Major E, McPherson K, Short A, Vessey MP: Intensive Care Society's APACHE II study in Britain and Ireland-II: outcome comparisons of intensive care units after adjustment for case mix by the American APACHE II method. Br Med J. 1993, 307: 977-981.View ArticleGoogle Scholar
- Rowan KM, Kerr JH, Major E, McPherson K, Short A, Vessey MP: Intensive Care Society's APACHE II study in Britain and Ireland-I: variations in case mix of adult admissions to general intensive care units and impact on outcome. Br Med J. 1993, 307: 972-977.View ArticleGoogle Scholar
- Lemeshow S, Le Gall JR: Modeling the severity of illness of ICU patients. A systems update. JAMA. 1994, 272: 1049-1055. 10.1001/jama.272.13.1049.PubMedView ArticleGoogle Scholar
- Moreno R, Apolone G: Impact of different customization strategies in the performance of a general severity score. Crit Care Med. 1997, 25: 2001-2008. 10.1097/00003246-199712000-00017.PubMedView ArticleGoogle Scholar
- Le Gall JR, Lemeshow S, Leleu G, Klar J, Huillard J, Rué M, Teres D, Artigas A: Customized probability models for early severe sepsis in adult intensive care patients. JAMA. 1995, 273: 644-650. 10.1001/jama.273.8.644.PubMedView ArticleGoogle Scholar
- Sarmiento X, Rué M, Guardiola JJ, Toboso JM, Soler M, Artigas A: Assessment of the prognosis of coronary patients. Performance and customization of generic severity indexes. Chest. 1997, 111: 1666-1671.PubMedView ArticleGoogle Scholar
- Knaus WA, Draper EA, Wagner DP, Zimmerman JE: APACHE II: a severity of disease classification system. Crit Care Med. 1985, 13: 818-829.PubMedView ArticleGoogle Scholar
- Knaus WA, Wagner DP, Draper EA, Zimmerman JE, Bergner M, Bastos PG, Sirio CA, Murphy DJ, Lotring T, Damiano A, Harrell FE: The APACHE III prognostic scoring system: risk prediction of hospital mortality for critically ill hospitalized adults. Chest. 1991, 100: 1619-1636.PubMedView ArticleGoogle Scholar
- Hosmer DW, Lemeshow S: Applied Logistic Regression. New York: John Wiley & Sons;. 1989Google Scholar
- Hanley JA, McNeil BJ: The meaning and use of the area under a receiver operating characteristic (ROC) curve. Radiology. 1982, 143: 29-36.PubMedView ArticleGoogle Scholar
- Rapoport J, Teres D, Lemeshow S, Gehlbach S: A method for assessing the clinical performance and cost-effectiveness of intensive care units: a multicenter inception cohort study. Crit Care Med. 1994, 22: 1385-1391.PubMedView ArticleGoogle Scholar
- Lemeshow S, Hosmer DW: A review of goodness of fit statistics for use in the development of logistic regression models. Am J Epidemiol. 1982, 115: 92-106.PubMedGoogle Scholar
- Markgraf R, Deutschinoff G, Pientka L, Scholten T: Comparison of Acute Physiology and Chronic Health Evaluations II and III and Simplified Acute Physiology Score II: a prospective cohort study evaluating these methods to predict outcome in a German interdisciplinary intensive care unit. Crit Care Med. 2000, 28: 26-33. 10.1097/00003246-200001000-00005.PubMedView ArticleGoogle Scholar