The effective reproductive number: Modeling and prediction with application to the multi-wave Covid-19 pandemic

Classical compartmental models of infectious disease assume that spread occurs through a homogeneous population. This produces poor fits to real data, because individuals vary in their number of epidemiologically-relevant contacts, and hence in their ability to transmit disease. In particular, network theory suggests that super-spreading events tend to happen more often at the beginning of an epidemic, which is inconsistent with the homogeneity assumption. In this paper we argue that a flexible decay shape for the effective reproductive number (R_t) indexed by the susceptible fraction (S_t) is a theory-informed modeling choice, which better captures the progression of disease incidence over human populations. This, in turn, produces better retrospective fits, as well as more accurate prospective predictions of observed epidemic curves. We extend this framework to fit multi-wave epidemics, and to accommodate public health restrictions on mobility. We demonstrate the performance of this model by doing a prediction study over two years of the SARS-CoV2 pandemic.