ABSTRACT
Objective To evaluate the clinical utility of automatable prediction models for increasing goals-of-care discussions among hospitalized patients at the end-of-life.
Materials and Methods We developed three Random Forest models and updated the Modified Hospital One-year Mortality Risk model: alternative models to predict one-year mortality (proxy for EOL status) using admission-time data. Admissions from July 2011-2016 were used for training and those from July 2017-2018 were used for temporal validation. We simulated alerts for admissions in the validation cohort and modelled alternative scenarios where alerts lead to code status orders (CSOs) in the EHR. We linked actual CSOs and calculated the expected risk difference (eRD), the number needed to benefit (NNB) and the net benefit (NB) of each model for the patient-centered outcome of a CSO among EOL hospitalizations.
Results Models had a C-statistic of 0.79-0.86 among unique patients. A CSO was documented during 2599 of 3773 hospitalizations at the EOL (68.9%). At a threshold that identified 10% of eligible admissions, the eRD ranged from 5.4% to 10.7% (NNB 5.4-10.9 alerts). Under usual care, a CSO had a 34% PPV for EOL status. Using this to inform the relative cost of FPs, only two models improved NB over usual care. A RF model with diagnostic predictors had the highest clinical utility by either measure, including in sensitivity analyses.
Discussion Automatable prediction models with acceptable temporal validity differed meaningfully in their expected ability to improve patient-centered outcomes over usual care.
Conclusion Decision-analysis should precede implementation of automated prediction models for improving palliative and EOL care outcomes.
INTRODUCTION
End-of-life (EOL) conversations and shared decision-making between clinical staff and hospitalized patients can improve the quality of EOL care.[1,2] In the hospital setting, these conversations inform goals-of-care (GOC) documentation, particularly code status orders (CSOs), which encode the essential preferences for life-supporting therapy.[3] Hospitalizations are frequent at the EOL, and where the need to plan for future care is matched by the opportunity.[4] However, hospitalized patients with a poor prognosis do not benefit from EOL conversations or GOC documentation as often as they should.[5] Closing this gap is challenged by workload constraints and difficulty in prognostication.[6]
Clinical decision support systems (CDSS) that integrate automated prediction models may help increase the prevalence of GOC discussions by generating computerized alerts for patients with a high risk of mortality.[2,7] The rationale is that physicians, when explicitly alerted to the poor prognosis of a patient in their care, will initiate a discussion about GOC if one is appropriate and has not already occurred.
In the translational pathway of prediction models, an increasingly recognized step is the assessment of clinical utility,[8] which should occur before a prospective evaluation of clinical impact.[9,10] CDSSs, particularly those with machine learning (ML) models, are potentially costly to implement[11,12] and their impact highly subject to local factors,[13] giving reason to assess clinical utility before investing in application. Decision-analytic methods for this assessment using observational data are accessible[10,14–16] but are rarely unused for prediction models prompting palliative and EOL care (PEOLC) interventions, resulting in poor evidence of value.[17,18] A few decision-analyses in this area of research have assessed system-perspective monetary benefit;[19,20] a decision-analytic evaluation of clinical benefits and harms from the perspective of patient-centered quality improvement has remained elusive.
In this study, we evaluated the clinical utility of locally applicable prediction models using a routinely collected measure of GOC discussions, CSOs in the EHR. Our primary objective was to evaluate the clinical utility of a novel ML model compared to a published model[21] and models requiring less types of predictors. In the process, we demonstrate innovative strategies to increase the applicability of simple decision-analytic techniques for assessing the utility of automatable prediction models before implementation.
MATERIALS AND METHODS
This retrospective study comparing prediction models includes methods for the development, validation, and decision-analytic evaluation of prediction models. We conform to the TRIPOD guidelines[22] for reporting prognostic modelling methods, and to the relevant aspects of the CHEERS guidelines[23] for reporting decision-analytic methods. The study took place at an integrated university hospital network with two sites and about 700 acute care beds in the city of Sherbrooke, Quebec, Canada (details in the supplement). IRB approval was obtained prior to data collection (IRB of the CIUSSS de l’Estrie – CHUS #2018-2478).
Source of data and participants
All predictor data in the study was collected from the institutional data warehouse, which combines EHR and administrative data. All adult hospitalizations admitted to a non-psychiatric service between July 1st, 2011 and June 30th, 2018 were included in the overall cohort, except for admissions to rarely admitting specialities (e.g., genetics) or admissions with a legal context (e.g., court-ordered). Mortality records were sourced from the Quebec vital statistics registry and considered complete until June 30th, 2019 (additional details in supplement).
The overall cohort was split temporally, with a training cohort defined as admissions occurring between July 1st, 2011 and June 30th, 2016 inclusively, and a testing cohort defined as hospital admissions occurring between July 1st, 2017 and June 30th, 2018 inclusively. This split was designed to simulate the prospective evaluation of a given model had it been trained with all available data just before midnight on June 30th, 2017 and then applied prospectively for one year at our institution. Hospitalizations that occurred between July 1st, 2016 and June 30th, 2017, inclusively, were excluded to prevent any unrealistic leakage of outcome information between the training and testing cohort.
For the evaluation of clinical utility, our population of interest was all hospitalizations where there was enough time for a GOC discussion to occur and where it was not inappropriate or unnecessary given the information available to a CDSS at the point-of-care. Our selection criteria were hospitalizations with overnight stay that were not admissions in obstetrics or palliative care. The included sample for the evaluation of clinical utility were all such hospitalizations in the testing cohort.
Prediction models
We developed a ML model using the Random Forest (RF) algorithm that includes administrative, demographic, and diagnostic predictors accessible at the time of hospital admission to predict one-year mortality (RF-AdminDemoDx, 244 predictors). As an alternative strategy, we updated the Modified Hospital One-year Mortality Risk (mHOMR) model[21] for local application (nine predictors). In addition, we specified two simplified versions of the RF-AdminDemoDx model: one where no diagnostic variables were included (RF-AdminDemo, twelve predictors) and one where only four variables – age, sex, admission service and admission type – were included (RF-Minimal). All four prediction models were feasible to operationalize with the existing informatics infrastructure, though had different requirements in terms of data access and implementation (Table 1). Data generation processes were investigated to align retrospectively extracted variables with what would be available within a few minutes of hospital admission. The models were trained with the training cohort, their temporal validity evaluated with the testing cohort and their clinical utility evaluated with the CDSS-eligible cohort. Model development, specification and validation is fully described in the supplement.
Perspective
The evaluation of clinical utility was conducted from the perspective of a clinician-led quality improvement team that aims to implement a CDSS to increase the prevalence of GOC discussions for patients with a poor prognosis: a promising initiative[3] for a well-established problem.[5] A necessary component of GOC discussions for hospitalized patients are discussions about code status, the documentation of which had been standardized as CSOs in the institutional EHR since 2015. The documentation of resuscitation preferences for a hospitalized patient with a poor prognosis is a positively-valued, patient-centered outcome in the context of EOL communication[24] and its absence for the same population considered a potentially harmful medical error.[2,25] The main objective of decision analysis was to identify the prediction model that maximized this quality indicator. The secondary objective was evaluating the net benefit[14,15] (NB) of prediction models.
Alternative strategies
We simulated the operation of a CDSS that uses alternative prediction models for triggering an alert. Conceptually, alerts would suggest discussing GOC, including CPR preferences, if appropriate.[7] The system would generate alerts after midnight for eligible patients admitted the previous day having a predicted risk greater or equal to a certain decisional threshold. Since intervention harm was minimal and time constraints were known to limit GOC discussions,[2,6] we considered the proportion of admissions with an alert, P(Alert), to be the most appropriate criteria for setting model thresholds. We set a P(Alert) of 10% as a point of reference and expected thresholds that identified between 5-20% of CDSS-eligible admissions to be an appropriate range for sensitivity analyses.
The alternative strategies under consideration were the four mortality alert rules that resulted from applying to each prediction model a P(Alert)-specific threshold.
Outcome definitions
Electronic CSOs were linked to hospitalizations in the testing cohort after model development and did not have any role in predictive validation. These orders could convey one of three resuscitation preferences: wants all resuscitation measures (Full code), does not want CPR but wants endotracheal intubation if necessary (DNR/Intubation-OK), and does not want CPR or intubation (DNR/DNI). We considered a CSO to have occurred during a hospitalization if at least one was documented between one week before the admission date and the discharge date, inclusively. The extra week was added to associate a hospitalization with any CSOs documented during observation in the ED prior to hospital admission. The main outcome was a hospitalization with a CSO among those for patients at the EOL, which we defined as death within one year of admission.
Decision trees
We modelled two decision trees where alerts led to the desired action (Figure 1). The first was based on the conventional assumptions of decision analysis for prediction models,[26,27] where alerts lead to action, and the absence of an alert leads to inaction. The second was a scenario-appropriate adaptation that allowed assessing model utility relative to a strategy of usual care, where alerts lead to action, and the absence of an alert leads to usual care: either action or inaction depending on what had factually occurred for the alert-negative case. For both trees, the benefit to patients of discussing and documenting GOC[2] was attributed to true positives (TPs). The cost of this action is spending clinical time,[2] which was attributed to false positives (FPs). Our valuation procedure and assumptions are further explained in the supplement.
Abbreviations: CDSS, Clinical decision support system; CSO, code status order; EOL, end-of-life; TP, true positive; FP, false positive; TN, true negative; FN, false negative.
Two decision trees modelling the potential outcomes of each hospitalization in the CDSS-eligible cohort under alternative strategies. The strategy of an alert rule, R, implies that a CDSS is implemented and uses R to generate alerts. In the strategy of usual care, UC, CSOs occur as they factually did between July 2017 and July 2018 in the two hospitals of Sherbrooke, Quebec, Canada. In both trees, an alert always implies a CSO. The difference between the two trees is how outcomes unfold in the absence of an alert. In decision tree 1, no alert results in no action (no CSO). In decision tree 2, no alert results in the action that occurred retrospectively in usual care. The first tree models the conventional scenario of decision curve analysis where a prediction rule aims to reduce intervention-related harm, while the second models the scenario of a CDSS that aims to increase a routine good practice that is constrained by time. A TP outcome occurs when a CSO is documented during a hospitalization for a patient who died within one year of admission (EOL status). A FP outcome occurs when a CSO is documented during a hospitalization for a patient who survives more than a year (“not EOL” status). A FN outcome occurs when no CSO is documented during an EOL hospitalization. A TN occurs when no CSO is documented for a “not EOL” hospitalization. The formulas to calculate the expected probability of each outcome for a given strategy are provided to the right of each terminal node.
To distinguish the effect of model-based predictions from the effect of simply generating alerts, we included a fifth model of uniformly random numbers between zero and one. We did not expect alerts from such a model to cause physicians to act in the same way as the validated prediction models, but it would serve to explicit a side-effect of the assumption that all alerts would cause the desired action of a CSO.
Statistical analysis
We described cohort characteristics stratified by EOL and CSO status. We assessed model discrimination using the C-statistic and its calibration using a calibration plot.[8] To assess construct validity of predictions, we regressed DNR preference against predicted risk in the CDSS-eligible cohort.
Our primary measure of clinical utility was the expected risk difference (eRD) compared to usual care of the main outcome, calculated for each rule as
This metric is based on the intention-to-treat estimator[10] and answers the hypothetical question: if every alert based on rule R had led to the desired action (a CSO), how many more hospitalizations at the EOL (as a proportion of all hospitalizations at the EOL) would have had a CSO?
For contextualizing the eRD, we calculated the number needed to benefit (NNB):
The NNB is the number of alerts needed to increase benefit by one outcome over usual care. It is the reciprocal of the expected benefit among those with an alert: P(no CSO, EOL | Alert). Conceptually, it incorporates both the number needed to screen, 1/P(EOL | Alert), and the number needed to treat, 1/P(no CSO | EOL), as originally described,[16] but it was calculated without assuming conditional independence: the subset of patients at the end-of-life successfully screened with an alert would not necessarily have the same chance of beneficial treatment (the counterfactual outcome in the event of “no CSO”) as the set of all patients at the end-of life.
Our secondary measure of clinical utility was the NB, calculated for each strategy, S, as
In our scenario, model threshold was based on estimated availability of clinical time, not necessarily patient-provider preference for CSO. This made threshold potentially unsuitable to inform the exchange rate, calculated in conventional decision curve analysis as Threshold/(1-Threshold).[27] The exchange rate represents the theoretical ratio between the harm of inappropriate inaction (FN) and the harm of inappropriate action (FP), which can be obtained by various means;[14] the threshold method is a convenient simplification in the absence of other utility estimates in the validation set.[27] We calculated a model-independent exchange rate using observed actions of clinicians (i.e., using their revealed preferences[14]), who implicitly decide under uncertainty between the harm of inaction and the time-cost of action:
Substituting Equation 4 in Equation 3 results in a NB of zero for the default strategy S = Usual care. We plotted decision curves for both decision trees and for both a threshold-based and observed exchange rate.
We performed a two-way sensitivity analysis[14] between P(Alert) and the exchange rate in subgroups of service type and hospital site. We assessed subgroup heterogeneity using a forest plot of the expected relative risk (the ratio of the terms in the eRD) in relevant inpatient populations.
Consistent with the decision-analytic design, no p-valued significance tests were performed between the alternative strategies.[28] We bootstrapped 95% confidence intervals for all estimates.[8,29] To verify the potential influence of including all hospitalizations in the CDSS-eligible cohort rather than sampling unique patients, we repeated analyses using the first, last, and a random hospitalization per patient. All statistical analyses were performed with R version 3.6.3[30] (relevant extensions and details in the supplement).
RESULTS
Sample and model description
The participant flow diagram is presented in Figure 2. Between July 1st, 2011 and June 30th, 2018, there were 175041 hospitalizations for adults in a non-psychiatric service at our institution (93295 patients). After excluding 76 hospitalizations with rare circumstances, the training cohort included 122860 hospitalizations between July 1st, 2011 and June 30th, 2016 (70788 patients) and the testing cohort included 26291 hospitalizations between July 1st, 2017 and June 30th, 2018 (20012 patients). There were 22034 hospitalizations (16490 patients) in the CDSS-eligible cohort. Patient-hospitalization characteristics are presented for the CDSS-eligible cohort in Table 2 (description of other cohorts in the supplement). Prediction models had acceptable temporal validity (Table 3, Figure S1). When sampling over unique patients, the C-statistic ranged from 0.84 to 0.89 in the testing cohort and lowered to 0.79 to 0.86 in the CDSS-eligible cohort. Figure 3 describes EOL process indicators as a function of model-predicted risk; all models had good construct validity for DNR preferences.
Abbreviations: CDSS, clinical decision support system.
Abbreviations: CSO, code status order; DNR, do-not-resuscitate. Figure caption: Binary variables regressed against predicted risk of a given model then plotted along with 95% CI bands using the loess algorithm. One random hospitalization per patient sampled from the CDSS-eligible cohort before applying regression. Panels A and B include 16490 patients, where 2248 died within one year of hospitalization and 5241 had CSO documentation during that hospitalization. Panel C includes the 5241 patient/hospitalizations with CSO documentation, where 3552 preferred a DNR status in the last CSO documented before discharge. Note that a DNR includes both “DNR/Intubation-OK” and “DNR/DNI”: those without a DNR in Panel C prefer a “Full Code”.
Code status orders at the end of life
There were 7648 hospitalizations associated with a CSO in the CDSS-eligible cohort (34.7%). Among these, 2599 (34.0%) were associated with death within one-year of admission; clinicians were observed to document a CSO during one hospitalization at the EOL□for every ≈1.9 hospitalizations not at the EOL (observed exchange rate of 2599 TP:5049 FP). On average, clinicians acted as though the harm of FNs was 1.9x as harmful as a FP.
There were 3773 hospitalizations at the EOL in the CDSS-eligible cohort (17.1%). Among these, a CSO was not documented in 1174 cases, meaning a minimal GOC discussion had not been documented for 31.1% of applicable hospitalizations with overnight stays at the EOL. Compared to hospitalizations at the EOL that did have a CSO, these cases were more likely to be elective (OR 9.6 [95% CI, 7.3-12.5]), in surgical specialities (OR 5.1 [95% CI, 4.4-6.1]), for younger patients (mean age 68.3 vs 76.9 years, [95% CI −9.5 to −7.6]), and of shorter duration (mean length of stay 6.5 vs 12.9 days, [95% CI −7.1 to −5.7]).
Clinical utility
Simulated at a P(Alert) of 10%, each model would have generated on average six alerts per day over one year (Figure 4). At this same level of resource use, the eRD varied between 5.4% and 10.7%, and the NNB between 5.4 and 10.9 alerts (Table 4). The RF-AdminDemoDx model had the highest clinical utility. This model also maximized NB in the decision curves regardless of the decision tree or exchange rate used (Figure 5). When routine clinical actions were considered, only the RF-AdminDemoDx and RF-AdminDemo models could increase value above usual care in the range of desired alert frequency (Figure 5-D).
Abbreviations: CHUS, Centre Hospitalier Universitaire de Sherbrooke; CSO, code status order; EOL, end-of-life.
Panel A: The proportion of CDSS-eligible cases with a predicted risk above the threshold, P(Alert), is plotted against model threshold. A region of interest is outlined where thresholds satisfy reasonable workload demands (5-20% of CDSS-eligible admission). Panel B: The region of interest in Panel A is magnified. Panel C and D: Time-series between July 1st, 2017 and June 30th, 2018 of daily events stratified by hospital site. The frequency of actual admissions is compared with the frequency of simulated alerts at thresholds satisfying a reference P(Alert) of 10%. “All” refers to all admissions in the testing cohort (excluding pediatric and psychiatric admissions). The loess algorithm was used to smooth day-to-day variations using a span of 0.2 for all curves. CHUS-Fleurimont refers to site A and CHUS-Hôtel-Dieu refers to site B in the main text. Non-CDSS-eligible cases at site A include mostly obstetrical admissions.
Abbreviations: CSO, code status order; EOL, end of life; FP, false positive; TP, true positive.
Decision curves to assess net benefit as a function of either Threshold or P(Alert). Net benefit was standardized by dividing by outcome prevalence (one-year mortality). In the top panels, it is assumed that model threshold is informative of preference between a TP and FP outcome: the exchange rate is calculated as Threshold/(1-Threshold). In the bottom panels, it is instead assumed that observed clinical actions are informative of the exchange rate, which is assigned a constant value: nTP/nFP outcomes in usual care. This value is the odds that a CSO in usual care is a true positive, and as a probability p = nTP/(nTP+nFP), is where the strategy of usual care intersects zero in Panel A. In the left-side panels, the first decision tree is applied to generate TP and FP outcomes, where prediction-based alerts lead to CSO (i.e., “action” or “treatment”) and the absence of an alert leads to no CSO (i.e., “no action” or “no treatment”). In this tree, “Alerts for none” and “Alerts for all” implies “CSO for none” and “CSO for all”, respectively. In the right-side panels, the second decision tree is applied, where alerts lead to CSO, and the absence of an alert leads to usual care: whatever action had factually occurred for a given case. In this tree, “Alerts for all” still implies “CSO for all”, but “Alerts for none” implies “Usual care”. The strategy of “Alerts for all” is a distant outlier in the bottom panels, corresponding to a constant standardized net benefit of around −1.5. The strategy of “Alerts for none” overlaps “Usual care” in panels B to D.
Sensitivity analysis
The clinical utility of the RF-AdminDemoDx model remained the highest among models in the two-way sensitivity analysis (Figure 6). Subgroup analysis indicated heterogeneity that could influence implementation, including a smaller benefit for all models at site B (Figure S2-S6). Estimates of clinical utility using different sampling strategies did not change the direction or interpretation of results (Table S6-S8, Figure S7-S9).
The net benefit (NB) was calculated for each strategy, for each tile and for each plot; the strategy with the highest net benefit is indicated for the corresponding combination of P(Alert) and exchange rate in each hospital site (column) and service type (row). Any ties were resolved by selecting the strategy with the highest benefit, P(TP), or randomly selecting if a tie persisted. A higher exchange rate indicates a higher value to the time-cost of discussing GOC with a patient not at the EOL (worried about FPs), while a lower exchange rate indicates a higher value to the harm of omitting a GOC discussion for a patient at the EOL (worried about FNs). The overall exchange rate (dotted line) was calculated using Equation 4 in the full cohort and the subgroup exchange rate (solid line) corresponds to the result of Equation 4 among a given subgroup. The mHOMR, RF-AdminDemo, and RF-Minimal models are not referenced because they were never a strategy with the highest net benefit.
DISCUSSION
Improving patient identification for routine PEOLC interventions is a priority for health care stakeholders aiming to reconcile the default policies of life-sustaining therapy with the static truth that all life comes to an end. We performed an up-to-date review of model validation studies in this area of research and provide both a narrative and tabular synthesis of related studies in the supplement. In recent years, there has been a shift from manual screening tools[32] towards automated trigger tools,[33] with the latter shifting from query-based algorithms[34] towards increasingly flexible, but infrastructure-dependant, prediction models.[21,35–38] A challenge with such models is that their usual learning objective, minimizing the error of mortality prediction, is only indirectly related to the clinical objective of maximizing benefit for a resource-limited PEOLC intervention. In this setting of mismatched expectations, predictiveness does not mean usefulness, making it essential to assess clinical utility and not just predictive accuracy.[8–10,16]
In our review of the literature, we did not find any retrospective study evaluating the clinical utility – both benefits and harms – of automatable prediction models for prompting PEOLC interventions. In contrast, almost all studies reported the C-statistic for mortality, and these were generally above 0.8. The context insensitivity of the C-statistic makes it practical for research, but uninformative for practice: more value-based metrics are required to guide decision-makers.[14,17] Prediction models for prompting a PEOLC intervention had varying use-cases for decision support, including GOC discussion, palliative care referral, outpatient follow-up for advance care planning, or hospice referral. The benefits, harms, and resources associated with these actions differ between each other and between health systems; one curve does not fit all.
Strengths of our study included ensuring that retrospectively accessed data represented real-time data and the use of temporal rather than random splitting for validation: simulating prospective application at the point of care. When validated in similar cohorts (not necessarily target population), all models in our study had similar C-statistics as published models (i.e., above 0.8). However, the eRD and NNB for a patient-centered outcome ranged almost twofold, and only two models had a higher NB than usual care with our scenario-appropriate decision tree.
Others have validated the predictive performance of a model, then described physician opinion about the appropriateness of high-risk predictions for intervening.[35,38,39] While informative of construct validity, appropriateness does not represent a model’s usefulness over alternatives. If a mortality alert rule resulted in alerts for every hospitalization – and only hospitalizations – with a DNR in the CDSS-eligible cohort, its PPV for 1-year mortality would be 43.5% (2 314/5 325) and all cases would be appropriate for hypothetical CSO documentation; yet this rule is useless for improving this outcome because it tells clinicians what they already act upon. A similar situation could result from using a model that is highly influenced by terms like “palliative” and “DNR”,[40] or a model that uses historical palliative care consults to predict future consults.[36] Even if alerts correctly predict mortality or benefit, those who would benefit from usual care anyways might be disproportionately identified. More concerningly, those who do not usually benefit may be further marginalized.[7]
Prediction models are often evaluated in biased conditions[41] and rarely compared against routine clinical decision-making.[42] Clinical utility metrics – like an intention-to-treat estimator,[10] the NB,[15] or the NNB[16] – allow for patient-centered comparisons of prediction models with more appropriate assumptions. They can also detect unexpected differences in potential impact, like the difference in utility between our two sites, before any health system investment and exposure to patients. We demonstrated three innovative strategies to increase the applicability of decision analysis for assessing the utility of automatable prediction models.
First, we did not rely on a link between model threshold and clinical preference for net benefit analysis.[27,43] Instead, we linked the threshold to the desired alert frequency, representing resource use, and used other procedures to value outcomes. In doing so, we overcome a limitation of threshold-based NB analysis, which has been remarked as inappropriate for prediction model use-cases that require considering resource availability in addition to patient benefits and harms.[16] Note that the intent behind NB analysis – if not most decision-analytic methods[14] – is that it be adapted and extended to specific scenarios,[27,44] the motivating principle being precisely that off-the-shelf metrics are not necessarily appropriate for all scenarios and stakeholders.[17]
Second, we extended the original decision tree used for decision curve analysis to allow simulating model-augmented outcomes (e.g., that no alert can still lead to CSO if clinically appropriate), rather than model-determined outcomes (e.g., that no alert will lead to no CSO). We would not want or expect the latter for our use-case. By design, the adapted decision tree results in a more modest estimation of utility, one that accounts for the expected value of routine care: models can only increase benefit if it is there to be increased after applying usual clinical decision-making.
Third, we used empiric rates of TP and FP actions to inform an observed exchange rate. This enabled decision curve analysis while comparing models at the same level of resource use, which was not necessarily the same threshold across models. Among those with a CSO, most patients preferred a DNR when the predicted risk of mortality was above 10-15%, but such a threshold would result in alerts for half of CDSS-eligible admissions. While likely acceptable for patients, who have little to lose and much to gain from a routine GOC discussion, this low threshold could imply unreasonable workloads for clinicians and cause alert fatigue.[12] The observed exchange rate is a simple measure of the benefit-for-time trade-off that limits a good practice with minimal intervention-related harm. It is readily reproducible if practice patterns change over time and we believe insightful about clinical decision-making, noting that physicians may be influenced by an inflated perception of GOC-related cost.[45] This technique could facilitate the clinical utility assessment of other models for improving good practices in time-constrained environments, where utilities cannot be inferred from the model threshold. We used CSOs because they were the only electronic indicator of GOC documentation at our institution, but the same technique could be applied for other standardized indicators of the EOL communication process, like Physician Orders for Life-Sustaining Treatment.[46]
Our reproduction of the mHOMR model did not discriminate one-year mortality as well as in Ontario (c=0.84 vs 0.89),[21] but it was relatively simple to generalize to our institution. We cannot say the same of our ML model, which relies on admission diagnoses in Quebec-local French and would need another free-text mapping to be transportable beyond provincial borders (we report all variable definitions to enable this). However, while the local instance of our ML model is less geographically transportable than mHOMR, it is convincingly more useful for future application at our institution. This finding adds evidence to the recommendation that the pursuit of model generalizability should not be at the expense of local clinical utility.[13]
Our study has several limitations. First, resuscitation preference documentation is an essential but limited measure of EOL communication.[24] We did not measure the quality of the GOC discussions that preceded a CSO nor the concordance of preferences with care received.[5] However, the role of this study was to inform implementation and not substitute a prospective evaluation of clinical impact, where these higher-value patient outcomes should be assessed before long-term adoption.[8] Second, decision analysis requires simplifying assumptions to be practical, like assuming alerts would deterministically lead to action.[10,27] To increase the transparency of these assumptions, we compared prediction models to a random model. Third, due to the COVID-19 pandemic, a repeat validation is likely warranted before local application because models rely on non-causal associations, such as between admission service and death, that may have unexpectedly shifted after systemic reorganization. Finally, although evaluating clinical utility of a prediction model is recommended and provides more value-based metrics than evaluating just predictive performance,[8–10,14–17] more research is required to investigate how well these metrics predict the actual impact of a model-based CDSS. Future studies can refine on decision analysis based on this retroactive feedback, like including model-independent effects from behavioural economic–inspired co-interventions.[47]
CONCLUSION
An evaluation of clinical utility is recommended after validating a prediction model because metrics of model predictiveness are not informative of value. This is particularly important for mortality prediction models having the use case of automatically prompting a PEOLC intervention, like a GOC discussion. Decision-analytic techniques to assess utility along patient-centered outcomes are feasible for quality improvement teams. They can help discriminate value from hype, calibrate expectations, and provide valuable information before CDSS implementation. As an adjunct to model validation, the routine evaluation of clinical utility could increase the value of automated predictive analytics implemented at the point-of-care.
Data Availability
The hospitalization data underlying this article cannot be shared publicly due to regulations to protect patient privacy that are overseen by the IRB. All prediction model specifications, decision-analytic model specifications and the observed exchange rate estimates with various sampling procedures are included in the article and its online supplement. The source code for specific aspects of the study, like the data analysis or visualization source code in R or the two-stage bootstrap source code in C++, is available upon request to the corresponding author.
STATEMENTS
All authors had full access to all the data in the study and accept responsibility to submit for publication.
Data availability
The hospitalization data underlying this article cannot be shared publicly due to regulations to protect patient privacy that are overseen by the IRB. All prediction model specifications, decision-analytic model specifications and the observed exchange rate estimates with various sampling procedures are included in the article and its online supplement. The source code for specific aspects of the study, like the data analysis or visualization source code in R or the two-stage bootstrap source code in C++, is available upon request to the corresponding author: Ryeyan.Taseen{at}USherbrooke.ca
Authors’ contributions
Conceptualization: RT, JF
Data curation: RT
Formal analysis: RT
Funding acquisition: RT, JF
Methodology: RT, JF
Resources: RT, JF
Software: RT
Validation: RT, JF
Visualization: RT
Writing – original draft: RT
Writing - review and editing: RT, JF
Funding/Support
This study was supported by a clinician-investigator training grant (MR1-291226) funded jointly by the Fonds de Recherche du Quebec – Santé and the Ministère de la Santé et des Services Sociaux (Dr Taseen) and a clinician-investigator grant (CC-253453) from the Fonds de Recherche du Quebec – Santé (Dr Ethier).
Role of funding source
The Fonds de Recherche du Quebec – Santé and the Ministère de la Santé et des Services Sociaux had no role in the design and conduct of the study; collection, management, analysis, and interpretation of the data; preparation, review, or approval of the manuscript; and decision to submit the manuscript for publication.
Competing interests
None.
Acknowledgements
We thank Luc Lavoie, associate professor of computer science at the Université de Sherbrooke, for guidance on data management. We thank Paul Farand, MD, MSc, for having led the quality improvement effort that resulted in the migration of paper to electronic code status orders at the Centre Hospitalier Universitaire de Sherbrooke and initially suggesting the idea to investigate automated mortality alerts to prompt GOC discussions.
Footnotes
Abstract headings formatted for readability.