Abstract

For the analysis of COVID‐19 pandemic data, we propose Bayesian multinomial and Dirichlet‐multinomial autoregressive models for time‐series of counts of patients in mutually exclusive and exhaustive observational categories, defined according to the severity of the patient status and the required treatment. Categories include hospitalized in regular wards (H) and in intensive care units (ICU), together with deceased (D) and recovered (R). These models explicitly formulate assumptions on the transition probabilities between these categories across time, thanks to a flexible formulation based on parameters that a priori follow normal distributions, possibly truncated to incorporate specific hypotheses having an epidemiological interpretation. The posterior distribution of model parameters and the transition matrices are estimated by a Markov chain Monte Carlo algorithm that also provides predictions and allows us to compute the reproduction number Rt. All estimates and predictions are endowed with an accuracy measure obtained thanks to the Bayesian approach. We present results concerning data collected during the first wave of the pandemic in Italy and Lombardy and study the effect of nonpharmaceutical interventions. Suitable discrepancy measures defined to check and compare models show that the Dirichlet‐multinomial model has an adequate fit and provides good predictive performance in particular for H and ICU patients.

Highlights

  • Proposed Dirichlet-Multinomial distribution and parameterizations Parameter estimation, model checking and comparison Application: Italian data Conclusions

  • We propose Multinomial and Dirichlet-Multinomial statistical autoregressive models to jointly analyse the observed time series of COVID-19 count data

  • It allows us to incorporate expert prior information. It accounts for public health non-pharmaceutical interventions (NPI) enforced to reduce the spread of the epidemic

Read more

Summary

Introduction

Categories include: susceptible not previously ill (S); recovered (R); positive cases in quarantine (Q); hospitalized in regular wards (H) and in intensive care units (ICU), together with deceased (D). The approach is related to the SEIR (Susceptible-Exposed-Infected-Recovered) epidemiological model

It allows us to incorporate expert prior information
Model assumptions
The conditional expected value is
Model checking
ICU D
Conclusions
Main References
Current results
Full Text
Published version (Free)

Talk to us

Join us for a 30 min session where you can share your feedback and ask us any queries you have

Schedule a call