No Arabic abstract
In this work we demonstrate how to automate parts of the infectious disease-control policy-making process via performing inference in existing epidemiological models. The kind of inference tasks undertaken include computing the posterior distribution over controllable, via direct policy-making choices, simulation model parameters that give rise to acceptable disease progression outcomes. Among other things, we illustrate the use of a probabilistic programming language that automates inference in existing simulators. Neither the full capabilities of this tool for automating inference nor its utility for planning is widely disseminated at the current time. Timely gains in understanding about how such simulation-based models and inference automation tools applied in support of policymaking could lead to less economically damaging policy prescriptions, particularly during the current COVID-19 pandemic.
We consider a continuous-time Markov chain model of SIR disease dynamics with two levels of mixing. For this so-called stochastic households model, we provide two methods for inferring the model parameters---governing within-household transmission, recovery, and between-household transmission---from data of the day upon which each individual became infectious and the household in which each infection occurred, as would be available from first few hundred studies. Each method is a form of Bayesian Markov Chain Monte Carlo that allows us to calculate a joint posterior distribution for all parameters and hence the household reproduction number and the early growth rate of the epidemic. The first method performs exact Bayesian inference using a standard data-augmentation approach; the second performs approximate Bayesian inference based on a likelihood approximation derived from branching processes. These methods are compared for computational efficiency and posteriors from each are compared. The branching process is shown to be an excellent approximation and remains computationally efficient as the amount of data is increased.
The fundamental models of epidemiology describe the progression of an infectious disease through a population using compartmentalized differential equations, but do not incorporate population-level heterogeneity in infection susceptibility. We show that variation strongly influences the rate of infection, while the infection process simultaneously sculpts the susceptibility distribution. These joint dynamics influence the force of infection and are, in turn, influenced by the shape of the initial variability. Intriguingly, we find that certain susceptibility distributions (the exponential and the gamma) are unchanged through the course of the outbreak, and lead naturally to power-law behavior in the force of infection; other distributions often tend towards these eigen-distributions through the process of contagion. The power-law behavior fundamentally alters predictions of the long-term infection rate, and suggests that first-order epidemic models that are parameterized in the exponential-like phase may systematically and significantly over-estimate the final severity of the outbreak.
Knowing COVID-19 epidemiological distributions, such as the time from patient admission to death, is directly relevant to effective primary and secondary care planning, and moreover, the mathematical modelling of the pandemic generally. We determine epidemiological distributions for patients hospitalised with COVID-19 using a large dataset ($N=21{,}000-157{,}000$) from the Brazilian Sistema de Informac{c}~ao de Vigil^ancia Epidemiologica da Gripe database. A joint Bayesian subnational model with partial pooling is used to simultaneously describe the 26 states and one federal district of Brazil, and shows significant variation in the mean of the symptom-onset-to-death time, with ranges between 11.2-17.8 days across the different states, and a mean of 15.2 days for Brazil. We find strong evidence in favour of specific probability density function choices: for example, the gamma distribution gives the best fit for onset-to-death and the generalised log-normal for onset-to-hospital-admission. Our results show that epidemiological distributions have considerable geographical variation, and provide the first estimates of these distributions in a low and middle-income setting. At the subnational level, variation in COVID-19 outcome timings are found to be correlated with poverty, deprivation and segregation levels, and weaker correlation is observed for mean age, wealth and urbanicity.
Several analytical models have been used in this work to describe the evolution of death cases arising from coronavirus (COVID-19). The Death or `D model is a simplified version of the SIR (susceptible-infected-recovered) model, which assumes no recovery over time, and allows for the transmission-dynamics equations to be solved analytically. The D-model can be extended to describe various focuses of infection, which may account for the original pandemic (D1), the lockdown (D2) and other effects (Dn). The evolution of the COVID-19 pandemic in several countries (China, Spain, Italy, France, UK, Iran, USA and Germany) shows a similar behavior in concord with the D-model trend, characterized by a rapid increase of death cases followed by a slow decline, which are affected by the earliness and efficiency of the lockdown effect. These results are in agreement with more accurate calculations using the extended SIR model with a parametrized solution and more sophisticated Monte Carlo grid simulations, which predict similar trends and indicate a common evolution of the pandemic with universal parameters.
We propose a physical theory underlying the temporal evolution of competing virus variants that relies on the existence of (quasi) fixed points capturing the large time scale invariance of the dynamics. To motivate our result we first modify the time-honoured compartmental models of the SIR type to account for the existence of competing variants and then show how their evolution can be naturally re-phrased in terms of flow equations ending at quasi fixed points. As the natural next step we employ (near) scale invariance to organise the time evolution of the competing variants within the effective description of the epidemic Renormalization Group framework. We test the resulting theory against the time evolution of COVID-19 virus variants that validate the theory empirically.