Performance metrics for models designed to predict individualized treatment effect

Objective Measuring the performance of models designed to predict individualized treatment effect is challenging, because the outcomes of two alternative treatments are inherently unobservable in one patient. The C-for-benefit was proposed to measure discriminative ability. We aimed to propose metrics of calibration and overall performance for models predicting treatment effect. Study Design and Setting Similar to the previously proposed C-for-benefit, we defined the observed treatment effect as the difference between outcomes in pairs of matched patients. Thus, we redefined the E-statistics, the logistic loss and the Brier score into metrics for measuring a model's ability to predict treatment effect. In a simulation study, the metric values of deliberately perturbed models were compared to those of the data generating model. To illustrate the performance metrics, different models predicting treatment effect were applied to the data of the Diabetes Prevention Program. Results As desired, performance metric values of perturbed models were consistently worse than those of the optimal model (Eavg-for-benefit[≥]0.070 versus 0.001, E90-for-benefit[≥]0.115 versus 0.002, log-loss-for-benefit[≥]0.757 versus 0.733, Brier-for-benefit[≥]0.215 versus 0.212). Calibration, discriminative ability, and overall performance of three different models were similar in the case study. Conclusion The proposed metrics are useful to assess the calibration and overall performance of models predicting individualized treatment effect, accessible via (https://github.com/CHMMaas/HTEPredictionMetrics).

[1]  E. Steyerberg,et al.  Concerns with the new SYNTAX score – Authors' reply , 2021, The Lancet.

[2]  John B. Wong,et al.  Predictive approaches to heterogeneous treatment effects: a scoping review , 2020, BMC Medical Research Methodology.

[3]  Maarten van Smeden,et al.  Calibration: the Achilles heel of predictive analytics , 2019, BMC Medicine.

[4]  Ewout W Steyerberg,et al.  Models with interactions overestimated heterogeneity of treatment effects and were prone to treatment mistargeting. , 2019, Journal of clinical epidemiology.

[5]  P. Austin,et al.  The Integrated Calibration Index (ICI) and related metrics for quantifying the calibration of logistic regression models , 2019, Statistics in medicine.

[6]  Ewout W Steyerberg,et al.  The proposed 'concordance-statistic for benefit' provided a useful metric when modeling heterogeneous treatment effects. , 2018, Journal of clinical epidemiology.

[7]  S. Athey,et al.  Generalized random forests , 2016, The Annals of Statistics.

[8]  Yvonne Vergouwe,et al.  Estimates of absolute treatment benefit for individual patients required careful modeling of statistical interactions. , 2015, Journal of clinical epidemiology.

[9]  Susan Athey,et al.  Recursive partitioning for heterogeneous causal effects , 2015, Proceedings of the National Academy of Sciences.

[10]  Jeremy B Sussman,et al.  Improving diabetes prevention with benefit based tailored treatment: risk based reanalysis of Diabetes Prevention Program , 2015, BMJ : British Medical Journal.

[11]  P. Bossuyt,et al.  SYNTAX score II , 2013, The Lancet.

[12]  Stef van Buuren,et al.  MICE: Multivariate Imputation by Chained Equations in R , 2011 .

[13]  Gary King,et al.  MatchIt: Nonparametric Preprocessing for Parametric Causal Inference , 2011 .

[14]  Trevor Hastie,et al.  Regularization Paths for Generalized Linear Models via Coordinate Descent. , 2010, Journal of statistical software.

[15]  C.J.H. Mann Clinical Prediction Models: A Practical Approach to Development, Validation and Updating , 2009 .

[16]  E. Steyerberg Clinical Prediction Models , 2008, Statistics for Biology and Health.

[17]  Gary King,et al.  Matching as Nonparametric Preprocessing for Reducing Model Dependence in Parametric Causal Inference , 2007, Political Analysis.

[18]  Sara T Brookes,et al.  Subgroup analyses in randomized trials: risks of subgroup-specific analyses; power and sample size for the interaction test. , 2004, Journal of clinical epidemiology.

[19]  D. Rubin Estimating causal effects of treatments in randomized and nonrandomized studies. , 1974 .

[20]  Irving John Good,et al.  Some Terminology and Notation in Information Theory , 1956 .

[21]  R Core Team,et al.  R: A language and environment for statistical computing. , 2014 .

[22]  G. Brier VERIFICATION OF FORECASTS EXPRESSED IN TERMS OF PROBABILITY , 1950 .