Forecasting COVID-19 Counts At A Single Hospital: A Hierarchical Bayesian Approach

We consider the problem of forecasting the daily number of hospitalized COVID-19 patients at a single hospital site, in order to help administrators with logistics and planning. We develop several candidate hierarchical Bayesian models which directly capture the count nature of data via a generalized Poisson likelihood, model time-series dependencies via autoregressive and Gaussian process latent processes, and share statistical strength across related sites. We demonstrate our approach on public datasets for 8 hospitals in Massachusetts, U.S.A. and 10 hospitals in the United Kingdom. Further prospective evaluation compares our approach favorably to baselines currently used by stakeholders at 3 related hospitals to forecast 2-week-ahead demand by rescaling state-level forecasts.

[1]  F. Famoye,et al.  Generalized poisson regression model , 1992 .

[2]  Dima Damen,et al.  Recognizing linked events: Searching the space of feasible explanations , 2009, 2009 IEEE Conference on Computer Vision and Pattern Recognition.

[3]  Generalized Poisson random variate generation , 1997 .

[4]  Neil D. Lawrence,et al.  Variational Gaussian Process Dynamical Systems , 2011, NIPS.

[5]  R. Epstein,et al.  A Predictive Model for Patient Census and Ventilator Requirements at Individual Hospitals During the Coronavirus Disease 2019 (COVID-19) Pandemic: A Preliminary Technical Report , 2020, Cureus.

[6]  Carl E. Rasmussen,et al.  Variational Gaussian Process State-Space Models , 2014, NIPS.

[7]  Andrew Gelman,et al.  Multilevel (Hierarchical) Modeling: What It Can and Cannot Do , 2006, Technometrics.

[8]  J. Wong.  Pandemic Surge Models in the Time of Severe Acute Respiratory Syndrome Coronavirus-2: Wrong or Useful? , 2020, Annals of Internal Medicine.

[9]  Carl E. Rasmussen,et al.  Gaussian processes for machine learning , 2005, Adaptive computation and machine learning.

[10]  Uri Shalit,et al.  Development and validation of a machine learning model predicting illness trajectory and hospital utilization of COVID-19 patients: A nationwide study , 2021, J. Am. Medical Informatics Assoc..

[11]  Andrew Gelman,et al.  The No-U-turn sampler: adaptively setting path lengths in Hamiltonian Monte Carlo , 2011, J. Mach. Learn. Res..

[12]  C. Murray.  Forecasting the impact of the first wave of the COVID-19 pandemic on hospital demand and deaths for the USA and European Economic Area countries , 2020, medRxiv.

[13]  John Salvatier,et al.  Probabilistic programming in Python using PyMC3 , 2016, PeerJ Comput. Sci..

[14]  Nicholas P Jewell,et al.  Predictive Mathematical Models of the COVID-19 Pandemic: Underlying Principles and Value of Projections. , 2020, JAMA.

[15]  A. Gelman.  Multilevel (Hierarchical) Modeling , 2006 .

[16]  Ihme COVID-19 Forecasting Team.  Modeling COVID-19 scenarios for the United States , 2020, Nature medicine.

[17]  Neil D. Lawrence,et al.  Gaussian Process Latent Variable Models for Visualisation of High Dimensional Data , 2003, NIPS.