Project description:The effective reproduction number, R(t) , plays a key role in the study of infectious diseases, indicating the current average number of new infections caused by an infected individual in an epidemic process. Estimation methods for the time evolution of R(t) , using incidence data, rely on the generation interval distribution, g(τ), which is usually obtained from empirical data or theoretical studies using simple epidemic models. However, for systems that present heterogeneity, either on the host population or in the expression of the disease, there is a lack of data and of a suitable general methodology to obtain g(τ). In this work, we use mathematical models to bridge this gap. We present a general methodology for obtaining explicit expressions of the reproduction numbers and the generation interval distributions, within and between model sub-compartments provided by an arbitrary compartmental model. Additionally, we present the appropriate expressions to evaluate those reproduction numbers using incidence data. To highlight the relevance of such methodology, we apply it to the spread of COVID-19 in municipalities of the state of Rio de Janeiro, Brazil. Using two meta-population models, we estimate the reproduction numbers and the contributions of each municipality in the generation of cases in all others.
Project description:The effective reproduction number, [Formula: see text], is an important epidemiological metric used to assess the state of an epidemic, as well as the effectiveness of public health interventions undertaken in response. When [Formula: see text] is above one, it indicates that new infections are increasing, and thus the epidemic is growing, while an [Formula: see text] is below one indicates that new infections are decreasing, and so the epidemic is under control. There are several established software packages that are readily available to statistically estimate [Formula: see text] using clinical surveillance data. However, there are comparatively few accessible tools for estimating [Formula: see text] from pathogen wastewater concentration, a surveillance data stream that cemented its utility during the COVID-19 pandemic. We present the [Formula: see text] package ern that aims to perform the estimation of the effective reproduction number from real-world wastewater or aggregated clinical surveillance data in a user-friendly way.
Project description:This study aimed to quantify the exposure-lag-response relationship between short-term changes in ambient temperature and absolute humidity and the transmission dynamics of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) in Japan. The prefecture-specific daily time-series of newly confirmed cases, meteorological variables, retail and recreation mobility, and Government Stringency Index were collected for all 47 prefectures of Japan for the study period from 15 February 2020 to 15 October 2022. Generalized conditional Gamma regression models were formulated with distributed lag nonlinear models by adopting the case-time-series design to assess the independent and interactive effects of ambient temperature and absolute humidity on the relative risk (RR) of the time-varying effective reproductive number (Rt). With reference to 17.8 °C, the corresponding cumulative RRs (95% confidence interval) at a mean ambient temperatures of 5.1 °C and 27.9 °C were 1.027 (1.016-1.038) and 0.982 (0.974-0.989), respectively, whereas those at an absolute humidity of 4.2 m/g3 and 20.6 m/g3 were 1.026 (1.017-1.036) and 0.995 (0.985-1.006), respectively, with reference to 10.6 m/g3. Both extremely hot and humid conditions synergistically and slightly reduced the Rt. Our findings provide a better understanding of how meteorological drivers shape the complex heterogeneous dynamics of SARS-CoV-2 in Japan.
Project description:Monitoring the effective reproduction number [Formula: see text] of a rapidly unfolding pandemic in real-time is key to successful mitigation and prevention strategies. However, existing methods based on case numbers, hospital admissions or fatalities suffer from multiple measurement biases and temporal lags due to high test positivity rates or delays in symptom development or administrative reporting. Alternative methods such as web search and social media tracking are less directly indicating epidemic prevalence over time. We instead record age-stratified anonymous contact matrices at a daily resolution using a longitudinal online-offline survey in Hungary during the first two waves of the COVID-19 pandemic. This approach is innovative, cheap, and provides information in near real-time for estimating [Formula: see text] at a daily resolution. Moreover, it allows to complement traditional surveillance systems by signaling periods when official monitoring infrastructures are unreliable due to observational biases.
Project description:Over the last months, cases of SARS-CoV-2 surged repeatedly in many countries but could often be controlled with nonpharmaceutical interventions including social distancing. We analyzed deidentified Global Positioning System (GPS) tracking data from 1.15 to 1.4 million cell phones in Germany per day between March and November 2020 to identify encounters between individuals and statistically evaluate contact behavior. Using graph sampling theory, we estimated the contact index (CX), a metric for number and heterogeneity of contacts. We found that CX, and not the total number of contacts, is an accurate predictor for the effective reproduction number R derived from case numbers. A high correlation between CX and R recorded more than 2 wk later allows assessment of social behavior well before changes in case numbers become detectable. By construction, the CX quantifies the role of superspreading and permits assigning risks to specific contact behavior. We provide a critical CX value beyond which R is expected to rise above 1 and propose to use that value to leverage the social-distancing interventions for the coming months.
Project description:Timely estimation of the transmissibility of a novel pandemic influenza virus was a public health priority in 2009.We extended methods for prospective estimation of the effective reproduction number (Rt) over time in an emerging epidemic to allow for reporting delays and repeated importations. We estimated Rt based on case notifications and hospitalizations associated with laboratory-confirmed pandemic (H1N1) 2009 virus infections in Hong Kong from June through October 2009.Rt declined from around 1.4-1.5 at the start of the local epidemic to around 1.1-1.2 later in the summer, suggesting changes in transmissibility perhaps related to school vacations or seasonality. Estimates of Rt based on hospitalizations of confirmed H1N1 cases closely matched estimates based on case notifications.Real-time monitoring of the effective reproduction number is feasible and can provide useful information to public health authorities for situational awareness and calibration of mitigation strategies.
Project description:Transmission potential of a pathogen, often quantified by the time-varying reproduction number Rt, provides the current pace of infection for a disease and indicates whether an emerging epidemic is under control. In this study, we proposed a novel method, EpiMix, for Rt estimation, wherein we incorporated the impacts of exogenous factors and random effects under a Bayesian regression framework. Using Integrated Nested Laplace Approximation, EpiMix is able to efficiently generate reliable, deterministic Rt estimates. In the simulations and case studies performed, we further demonstrated the method's robustness in low-incidence scenarios, together with other merits, including its flexibility in selecting variables and tolerance of varying reporting rates. All these make EpiMix a potentially useful tool for real-time Rt estimation provided that the serial interval distribution, time series of case counts and external influencing factors are available.
Project description:The time-varying reproduction number (Rt) is an important measure of epidemic transmissibility that directly informs policy decisions and the optimisation of control measures. EpiEstim is a widely used opensource software tool that uses case incidence and the serial interval (SI, time between symptoms in a case and their infector) to estimate Rt in real-time. The incidence and the SI distribution must be provided at the same temporal resolution, which can limit the applicability of EpiEstim and other similar methods, e.g. for contexts where the time window of incidence reporting is longer than the mean SI. In the EpiEstim R package, we implement an expectation-maximisation algorithm to reconstruct daily incidence from temporally aggregated data, from which Rt can then be estimated. We assess the validity of our method using an extensive simulation study and apply it to COVID-19 and influenza data. For all datasets, the influence of intra-weekly variability in reported data was mitigated by using aggregated weekly data. Rt estimated on weekly sliding windows using incidence reconstructed from weekly data was strongly correlated with estimates from the original daily data. The simulation study revealed that Rt was well estimated in all scenarios and regardless of the temporal aggregation of the data. In the presence of weekend effects, Rt estimates from reconstructed data were more successful at recovering the true value of Rt than those obtained from reported daily data. These results show that this novel method allows Rt to be successfully recovered from aggregated data using a simple approach with very few data requirements. Additionally, by removing administrative noise when daily incidence data are reconstructed, the accuracy of Rt estimates can be improved.
Project description:BackgroundDengue, an infectious tropical disease, has recently emerged as one of the most important mosquito-borne viral diseases in the world. We perform a retrospective analysis of the 2011 dengue fever epidemic in Pakistan in order to assess the transmissibility of the disease. We obtain estimates of the basic reproduction number R 0 from epidemic data using different methodologies applied to different epidemic models in order to evaluate the robustness of our estimate.ResultsWe first estimate model parameters by fitting a deterministic ODE vector-host model for the transmission dynamics of single-strain dengue to the epidemic data, using both a basic ordinary least squares (OLS) as well as a generalized least squares (GLS) scheme. Moreover, we perform the same analysis for a direct-transmission ODE model, thereby allowing us to compare our results across different models. In addition, we formulate a direct-transmission stochastic model for the transmission dynamics of dengue and obtain parameter estimates for the stochastic model using Markov chain Monte Carlo (MCMC) methods. In each of the cases we have considered, the estimate for the basic reproduction number R 0 is initially greater than unity leading to an epidemic outbreak. However, control measures implemented several weeks after the initial outbreak successfully reduce R 0 to less than unity, thus resulting in disease elimination. Furthermore, it is observed that there is strong agreement in our estimates for the pre-control value of R 0, both across different methodologies as well across different models. However, there are also significant differences between our estimates for the post-control value of the basic reproduction number across the two different models.ConclusionIn conclusion, we have obtained robust estimates for the value of the basic reproduction number R 0 associated with the 2011 dengue fever epidemic before the implementation of public health control measures. Furthermore, we have shown that there is close agreement between our estimates for the post-control value of R 0 across the different methodologies. Nevertheless, there are also significant differences between the estimates for the post-control value of R 0 across the two different models.
Project description:We compare several popular methods of estimating the basic reproduction number, R0, focusing on the early stages of an epidemic, and assuming weekly reports of new infecteds. We study the situation when data is generated by one of three standard epidemiological compartmental models: SIR, SEIR, and SEAIR; and examine the sensitivity of the estimators to the model structure. As some methods are developed assuming specific epidemiological models, our work adds a study of their performance in both a well-specified (data generating model and method model are the same) and miss-specified (data generating model and method model differ) settings. We also study R0 estimation using Canadian COVID-19 case report data. In this study we focus on examples of influenza and COVID-19, though the general approach is easily extendable to other scenarios. Our simulation study reveals that some estimation methods tend to work better than others, however, no singular best method was clearly detected. In the discussion, we provide recommendations for practitioners based on our results.