Comparison of methods for predicting COVID-19-related death in the general population using the OpenSAFELY platform.
OpenSAFELY Collaborative None., Williamson EJ., Tazare J., Bhaskaran K., McDonald HI., Walker AJ., Tomlinson L., Wing K., Bacon S., Bates C., Curtis HJ., Forbes HJ., Minassian C., Morton CE., Nightingale E., Mehrkar A., Evans D., Nicholson BD., Leon DA., Inglesby P., MacKenna B., Davies NG., DeVito NJ., Drysdale H., Cockburn J., Hulme WJ., Morley J., Douglas I., Rentsch CT., Mathur R., Wong A., Schultze A., Croker R., Parry J., Hester F., Harper S., Grieve R., Harrison DA., Steyerberg EW., Eggo RM., Diaz-Ordaz K., Keogh R., Evans SJW., Smeeth L., Goldacre B.
BACKGROUND: Obtaining accurate estimates of the risk of COVID-19-related death in the general population is challenging in the context of changing levels of circulating infection. METHODS: We propose a modelling approach to predict 28-day COVID-19-related death which explicitly accounts for COVID-19 infection prevalence using a series of sub-studies from new landmark times incorporating time-updating proxy measures of COVID-19 infection prevalence. This was compared with an approach ignoring infection prevalence. The target population was adults registered at a general practice in England in March 2020. The outcome was 28-day COVID-19-related death. Predictors included demographic characteristics and comorbidities. Three proxies of local infection prevalence were used: model-based estimates, rate of COVID-19-related attendances in emergency care, and rate of suspected COVID-19 cases in primary care. We used data within the TPP SystmOne electronic health record system linked to Office for National Statistics mortality data, using the OpenSAFELY platform, working on behalf of NHS England. Prediction models were developed in case-cohort samples with a 100-day follow-up. Validation was undertaken in 28-day cohorts from the target population. We considered predictive performance (discrimination and calibration) in geographical and temporal subsets of data not used in developing the risk prediction models. Simple models were contrasted to models including a full range of predictors. RESULTS: Prediction models were developed on 11,972,947 individuals, of whom 7999 experienced COVID-19-related death. All models discriminated well between individuals who did and did not experience the outcome, including simple models adjusting only for basic demographics and number of comorbidities: C-statistics 0.92-0.94. However, absolute risk estimates were substantially miscalibrated when infection prevalence was not explicitly modelled. CONCLUSIONS: Our proposed models allow absolute risk estimation in the context of changing infection prevalence but predictive performance is sensitive to the proxy for infection prevalence. Simple models can provide excellent discrimination and may simplify implementation of risk prediction tools.