Computing the daily reproduction number of COVID-19 by inverting the renewal equation using a variational technique

The COVID-19 pandemic has undergone frequent and rapid changes in its local and global infection rates, driven by governmental measures or the emergence of new viral variants. The reproduction number Rt indicates the average number of cases generated by an infected person at time t and is a key indicator of the spread of an epidemic. A timely estimation of Rt is a crucial tool to enable governmental organizations to adapt quickly to these changes and assess the consequences of their policies. The EpiEstim method is the most widely accepted method for estimating Rt. But it estimates Rt with a significant temporal delay. Here, we propose a method, EpiInvert, that shows good agreement with EpiEstim, but that provides estimates of Rt several days in advance. We show that Rt can be estimated by inverting the renewal equation linking Rt with the observed incidence curve of new cases, it. Our signal-processing approach to this problem yields both Rt and a restored it corrected for the “weekend effect” by applying a deconvolution and denoising procedure. The implementations of the EpiInvert and EpiEstim methods are fully open source and can be run in real time on every country in the world and every US state.

EpiLPS: a fast and flexible Bayesian tool for near real-time estimation of the time-varying reproduction number

In infectious disease epidemiology, the instantaneous reproduction number R(t) is a timevarying metric defined as the average number of secondary infections generated by individuals who are infectious at time t. It is therefore a crucial epidemiological parameter that assists public health decision makers in the management of an epidemic. We present a new Bayesian tool for robust estimation of the time-varying reproduction number. The proposed methodology smooths the epidemic curve and allows to obtain (approximate) point estimates and credible envelopes of R(t) by employing the renewal equation, using Bayesian P-splines coupled with Laplace approximations of the conditional posterior of the spline vector. Two alternative approaches for inference are presented: (1) an approach based on a maximum a posteriori argument for the model hyperparameters, delivering estimates of R(t) in only a few seconds; and (2) an approach based on a MCMC scheme with underlying Langevin dynamics for efficient sampling of the posterior target distribution. Case counts per unit of time are assumed to follow a Negative Binomial distribution to account for potential excess variability in the data that would not be captured by a classic Poisson model. Furthermore, after smoothing the epidemic curve, a “plug-in” estimate of the reproduction number can be obtained from the renewal equation yielding a closed form expression of R(t) as a function of the spline parameters. The approach is extremely fast and free of arbitrary smoothing assumptions. EpiLPS is applied on data of SARS-CoV-1 in Hong-Kong (2003), influenza A H1N1 (2009) in the USA and current SARS-CoV-2 pandemic (2020-2021) for Belgium, Portugal, Denmark and France.Author summaryThe instantaneous reproduction number R(t) is a key metric that provides important insights into an epidemic outbreak. We present a flexible Bayesian approach called EpiLPS (Epidemiological modeling with Laplacian-P-splines) for smooth estimation of the epidemic curve and R(t). Computational speed and absence of arbitrary assumptions on smoothing makes EpiLPS an interesting tool for near real-time estimation of the reproduction number. An R software package is available (https://github.com/oswaldogressani).

The renewal equation in nonlinear approximation

The family of Markov processes are considered in the article. We study the multidimensional renewal equation in nonlinear approximation. The purpose of the work is to find the limit of renewal function.

One Dimensional Reduction of a Renewal Equation for a Measure-Valued Function of Time Describing Population Dynamics

Despite their relevance in mathematical biology, there are, as yet, few general results about the asymptotic behaviour of measure valued solutions of renewal equations on the basis of assumptions concerning the kernel. We characterise, via their kernels, a class of renewal equations whose measure-valued solution can be expressed in terms of the solution of a scalar renewal equation. The asymptotic behaviour of the solution of the scalar renewal equation, is studied via Feller’s classical renewal theorem and, from it, the large time behaviour of the solution of the original renewal equation is derived.

A renewal equation model to assess roles and limitations of contact tracing for disease outbreak control

We propose a deterministic model capturing essential features of contact tracing as part of public health non-pharmaceutical interventions to mitigate an outbreak of an infectious disease. By incorporating a mechanistic formulation of the processes at the individual level, we obtain an integral equation (delayed in calendar time and advanced in time since infection) for the probability that an infected individual is detected and isolated at any point in time. This is then coupled with a renewal equation for the total incidence to form a closed system describing the transmission dynamics involving contact tracing. We define and calculate basic and effective reproduction numbers in terms of pathogen characteristics and contact tracing implementation constraints. When applied to the case of SARS-CoV-2, our results show that only combinations of diagnosis of symptomatic infections and contact tracing that are almost perfect in terms of speed and coverage can attain control, unless additional measures to reduce overall community transmission are in place. Under constraints on the testing or tracing capacity, a temporary interruption of contact tracing may, depending on the overall growth rate and prevalence of the infection, lead to an irreversible loss of control even when the epidemic was previously contained.

A Convex Optimization Solution for the Effective Reproduction NumberRt

COVID-19 is a global infectious disease that has affected millions of people. With new variants emerging with augmented transmission rates, slowing down of vaccine rollouts, and rising new cases threatening sanitary capabilities to the brink of collapse, there is the need to continue studying more effective forms to track its spread. This paper presents a strategy to compute the effective reproduction numberRt. Our method starts with a form of the renewal equation of the birth process specially suitable to computeRt. After showing that one can express it as a linear system, we proceed to solve it, along with appropriate constraints, using convex optimization. We demonstrate the method’s effectiveness using synthetic and real sequences of infections and comparing it with a leading approach.

A renewal equation model to assess roles and limitations of contact tracing for disease outbreak control

We propose a deterministic model capturing essential features of contact tracing as part of public health non-pharmaceutical interventions to mitigate an outbreak of an infectious disease. By incorporating a mechanistic formulation of the processes at the individual level, we obtain an integral equation (delayed in calendar time and advanced in time since infection) for the probability that an infected individual is detected and isolated at any point in time. This is then coupled with a renewal equation for the total incidence to form a closed system describing the transmission dynamics involving contact tracing. We define and calculate basic and effective reproduction numbers in terms of pathogen characteristics and contact tracing implementation constraints. When applied to the case of SARS-CoV-2, our results show that only combinations of diagnosis of symptomatic infections and contact tracing that are almost perfect in terms of speed and coverage can attain control, unless additional measures to reduce overall community transmission are in place. Under constraints on the testing or tracing capacity, a temporary interruption of contact tracing may, depending on the overall growth rate and prevalence of the infection, lead to an irreversible loss of control even when the epidemic was previously contained.