Recognizing the problems inherent in the assessment of model performance and the trade-off between discrimination and calibration, we also reported Brier scores for the models studied (Tables 2 and 3 of our article1). The Brier score, based on model prediction error, provides an overall estimate of performance. Our data demonstrated a progressive, albeit small, decrease in Brier score (ie, improved model performance) as model complexity increased. The addition of DNR status lowered the Brier scores for both APACHE models and for SAPS 3, suggesting its potential value as a predictor variable. As suggested, we also calculated both the Bayesian information criterion (BIC) and the corrected Akaike information criterion (AIC) for prognostic models with and without inclusion of DNR status.5 For each criterion, and similar to the Brier score analyses, the addition of DNR status was associated with a small improvement in model performance (approximately 3% decrease in AIC and BIC). The small improvements in Brier scores, BIC, and AIC associated with the addition of DNR status were not, however, reflected in statistically significant differences in area under the receiver operating characteristic curve or a consistent directional change in the value of the HLS.