Robust qualitative estimation of time-varying contact rates in uncertain epidemics

A data driven change-point epidemic model for assessing the impact of large gathering and subsequent movement control order on COVID-19 spread in Malaysia

The second wave of COVID-19 in Malaysia is largely attributed to a four-day mass gathering held in Sri Petaling from February 27, 2020, which contributed to an exponential rise of COVID-19 cases in the country. Starting from March 18, 2020, the Malaysian government introduced four consecutive phases of a Movement Control Order (MCO) to stem the spread of COVID-19. The MCO was implemented through various non-pharmaceutical interventions (NPIs). The reported number of cases reached its peak by the first week of April and then started to reduce, hence proving the effectiveness of the MCO. To gain a quantitative understanding of the effect of MCO on the dynamics of COVID-19, this paper develops a class of mathematical models to capture the disease spread before and after MCO implementation in Malaysia. A heterogeneous variant of the Susceptible-Exposed-Infected-Recovered (SEIR) model is developed with additional compartments for asymptomatic transmission. Further, a change-point is incorporated to model disease dynamics before and after intervention which is inferred based on data. Related statistical analyses for inference are developed in a Bayesian framework and are able to provide quantitative assessments of (1) the impact of the Sri Petaling gathering, and (2) the extent of decreasing transmission during the MCO period. The analysis here also quantitatively demonstrates how quickly transmission rates fall under effective NPI implementation within a short time period. The models and methodology used provided important insights into the nature of local transmissions to decision makers in the Ministry of Health, Malaysia.

The epidemicity index of recurrent SARS-CoV-2 infections

Several indices can predict the long-term fate of emerging infectious diseases and the effect of their containment measures, including a variety of reproduction numbers (e.g. [Formula: see text]). Other indices evaluate the potential for transient increases of epidemics eventually doomed to disappearance, based on generalized reactivity analysis. They identify conditions for perturbations to a stable disease-free equilibrium ([Formula: see text]) to grow, possibly causing significant damage. Here, we introduce the epidemicity index e0, a threshold-type indicator: if e0 > 0, initial foci may cause infection peaks even if [Formula: see text]. Therefore, effective containment measures should achieve a negative epidemicity index. We use spatially explicit models to rank containment measures for projected evolutions of the ongoing pandemic in Italy. There, we show that, while the effective reproduction number was below one for a sizable timespan, epidemicity remained positive, allowing recurrent infection flare-ups well before the major epidemic rebounding observed in the fall.

A simple criterion to design optimal non-pharmaceutical interventions for mitigating epidemic outbreaks

For mitigating the COVID-19 pandemic, much emphasis is made on implementing non-pharmaceutical interventions to keep the reproduction number below one. However, using that objective ignores that some of these interventions, like bans of public events or lockdowns, must be transitory and as short as possible because of their significant economic and societal costs. Here, we derive a simple and mathematically rigorous criterion for designing optimal transitory non-pharmaceutical interventions for mitigating epidemic outbreaks. We find that reducing the reproduction number below one is sufficient but not necessary. Instead, our criterion prescribes the required reduction in the reproduction number according to the desired maximum of disease prevalence and the maximum decrease of disease transmission that the interventions can achieve. We study the implications of our theoretical results for designing non-pharmaceutical interventions in 16 cities and regions during the COVID-19 pandemic. In particular, we estimate the minimal reduction of each region's contact rate necessary to control the epidemic optimally. Our results contribute to establishing a rigorous methodology to design optimal non-pharmaceutical intervention policies for mitigating epidemic outbreaks.

Accounting for non-stationarity in epidemiology by embedding time-varying parameters in stochastic models

The spread of disease through human populations is complex. The characteristics of disease propagation evolve with time, as a result of a multitude of environmental and anthropic factors, this non-stationarity is a key factor in this huge complexity. In the absence of appropriate external data sources, to correctly describe the disease propagation, we explore a flexible approach, based on stochastic models for the disease dynamics, and on diffusion processes for the parameter dynamics. Using such a diffusion process has the advantage of not requiring a specific mathematical function for the parameter dynamics. Coupled with particle MCMC, this approach allows us to reconstruct the time evolution of some key parameters (average transmission rate for instance). Thus, by capturing the time-varying nature of the different mechanisms involved in disease propagation, the epidemic can be described. Firstly we demonstrate the efficiency of this methodology on a toy model, where the parameters and the observation process are known. Applied then to real datasets, our methodology is able, based solely on simple stochastic models, to reconstruct complex epidemics, such as flu or dengue, over long time periods. Hence we demonstrate that time-varying parameters can improve the accuracy of model performances, and we suggest that our methodology can be used as a first step towards a better understanding of a complex epidemic, in situation where data is limited and/or uncertain.