An empirically adjusted approach to reproductive number estimation for stochastic compartmental models: A case study of two Ebola outbreaks

Bayesian modeling of dynamic behavioral change during an epidemic

For many infectious disease outbreaks, the at-risk population changes their behavior in response to the outbreak severity, causing the transmission dynamics to change in real-time. Behavioral change is often ignored in epidemic modeling efforts, making these models less useful than they could be. We address this by introducing a novel class of data-driven epidemic models which characterize and accurately estimate behavioral change. Our proposed model allows time-varying transmission to be captured by the level of "alarm" in the population, with alarm specified as a function of the past epidemic trajectory. We investigate the estimability of the population alarm across a wide range of scenarios, applying both parametric functions and non-parametric functions using splines and Gaussian processes. The model is set in the data-augmented Bayesian framework to allow estimation on partially observed epidemic data. The benefit and utility of the proposed approach is illustrated through applications to data from real epidemics.

Incorporating infectious duration-dependent transmission into Bayesian epidemic models

Compartmental models are commonly used to describe the spread of infectious diseases by estimating the probabilities of transitions between important disease states. A significant challenge in fitting Bayesian compartmental models lies in the need to estimate the duration of the infectious period, based on limited data providing only symptom onset date or another proxy for the start of infectiousness. Commonly, the exponential distribution is used to describe the infectious duration, an overly simplistic approach, which is not biologically plausible. More flexible distributions can be used, but parameter identifiability and computational cost can worsen for moderately sized or large epidemics. In this article, we present a novel approach, which considers a curve of transmissibility over a fixed infectious duration. The incorporation of infectious duration-dependent (IDD) transmissibility, which decays to zero during the infectious period, is biologically reasonable for many viral infections and fixing the length of the infectious period eases computational complexity in model fitting. Through simulation, we evaluate different functional forms of IDD transmissibility curves and show that the proposed approach offers improved estimation of the time-varying reproductive number. We illustrate the benefit of our approach through a new analysis of the 1995 outbreak of Ebola Virus Disease in the Democratic Republic of the Congo.

Simulating and Forecasting the COVID-19 Spread in a U.S. Metropolitan Region with a Spatial SEIR Model

The global COVID-19 pandemic has taken a heavy toll on health, social, and economic costs since the end of 2019. Predicting the spread of a pandemic is essential to developing effective intervention policies. Since the beginning of this pandemic, many models have been developed to predict its pathways. However, the majority of these models assume homogeneous dynamics over the geographic space, while the pandemic exhibits substantial spatial heterogeneity. In addition, spatial interaction among territorial entities and variations in their magnitude impact the pandemic dynamics. In this study, we used a spatial extension of the SEIR-type epidemiological model to simulate and predict the 4-week number of COVID-19 cases in the Charlotte-Concord-Gastonia Metropolitan Statistical Area (MSA), USA. We incorporated a variety of covariates, including mobility, pharmaceutical, and non-pharmaceutical interventions, demographics, and weather data to improve the model's predictive performance. We predicted the number of COVID-19 cases for up to four weeks in the 10 counties of the studied MSA simultaneously over the time period 29 March 2020 to 13 March 2021, and compared the results with the reported number of cases using the root-mean-squared error (RMSE) metric. Our results highlight the importance of spatial heterogeneity and spatial interactions among locations in COVID-19 pandemic modeling.

Modelling Representative Population Mobility for COVID-19 Spatial Transmission in South Africa

The COVID-19 pandemic starting in the first half of 2020 has changed the lives of everyone across the world. Reduced mobility was essential due to it being the largest impact possible against the spread of the little understood SARS-CoV-2 virus. To understand the spread, a comprehension of human mobility patterns is needed. The use of mobility data in modelling is thus essential to capture the intrinsic spread through the population. It is necessary to determine to what extent mobility data sources convey the same message of mobility within a region. This paper compares different mobility data sources by constructing spatial weight matrices at a variety of spatial resolutions and further compares the results through hierarchical clustering. We consider four methods for constructing spatial weight matrices representing mobility between spatial units, taking into account distance between spatial units as well as spatial covariates. This provides insight for the user into which data provides what type of information and in what situations a particular data source is most useful.

Modelling Representative Population Mobility for COVID-19 Spatial Transmission in South Africa

The COVID-19 pandemic starting in the first half of 2020 has changed the lives of everyone across the world. Reduced mobility was essential due to it being the largest impact possible against the spread of the little understood SARS-CoV-2 virus. To understand the spread, a comprehension of human mobility patterns is needed. The use of mobility data in modelling is thus essential to capture the intrinsic spread through the population. It is necessary to determine to what extent mobility data sources convey the same message of mobility within a region. This paper compares different mobility data sources by constructing spatial weight matrices at a variety of spatial resolutions and further compares the results through hierarchical clustering. We consider four methods for constructing spatial weight matrices representing mobility between spatial units, taking into account distance between spatial units as well as spatial covariates. This provides insight for the user into which data provides what type of information and in what situations a particular data source is most useful.

Bayesian compartmental models and associated reproductive numbers for an infection with multiple transmission modes

Zoonotic visceral leishmaniasis (ZVL) is a serious neglected tropical disease that is endemic in 98 countries. ZVL is primarily transmitted via a sand fly vector. In the United States, it is enzootic in some canine populations; it is transmitted from infectious mother to pup transplacentally, and vector-borne transmission is absent. This absence affords a unique opportunity to study (1) vertical transmission dynamics in dogs and (2) the importance of vertical transmission in maintaining an infectious reservoir in the presence of a vector. In this paper, we present Bayesian compartmental models and reproductive number formulations to examine (1) and (2), providing a mechanism to plan and evaluate interventions in regions where both transmission modes are present. First, we propose an individual-level susceptible, infectious, removed (SIR) model to study the effect of maternal infection status during pregnancy on pup infection progression. We provide evidence that pups born to diagnostically positive mothers during pregnancy are more likely to become diagnostically positive both earlier in life, and at some point during their lifetime, than those born to diagnostically negative mothers. Second, we propose a population-level SIR model to study the impact of a vertically maintained reservoir on propagating infection in a naive canine population through emergent vector transmission using simulation studies. We also present reproductive numbers to quantify contributions of vertically infected and vector-infected dogs to maintaining infection in the population. We show that a vertically maintained canine reservoir can propagate infection in a theoretical naive population in the presence of a vector.

Transmission center and driving factors of hand, foot, and mouth disease in China: A combined analysis

Hand, foot, and mouth disease (HFMD) has become a major public health issue in China. The disease incidence varies substantially over time and across space. To understand the heterogeneity of HFMD transmission, we compare the spatiotemporal dynamics of HFMD in Qinghai and Shanghai by conducting combined analysis of epidemiological, wavelet time series, and mathematical methods to county-level data from 2009 to 2016. We observe hierarchical epidemic waves in Qinghai, emanating from Huangzhong and in Shanghai from Fengxian. Besides population, we also find that the traveling waves are significantly associated with socio-economic and geographical factors. The population mobility also varies between the two regions: long-distance movement in Qinghai and between-neighbor commuting in Shanghai. Our findings provide important evidence for characterizing the heterogeneity of HFMD transmission and for the design and implementation of interventions, such as deploying optimal vaccine and changing local driving factors in the transmission center, to prevent or limit disease spread in these areas.

Modelling the reproductive power function

This paper discusses methods of estimating the reproductive power and the accompanying survival function of communicable events, e.g. infectious disease transmission. The early stage of an outbreak can be described by the infectiousness of the outbreak process, but in later stages of the outbreak, this is complicated by factors such as changing contact patterns and the impact of control measures. It is important to take these factors into account in order to get a good, if approximate, model for an outbreak process. This paper proposes a non-homogeneous birth process and regression model for the reproductive power function, similar to models in discrete survival analysis. A baseline reproductive power function gives a description of the outbreak when covariates are at their baseline values. As an illustration these methods are applied to an avian influenza (H5N1) outbreak among poultry in Thailand.

Bayesian compartmental model for an infectious disease with dynamic states of infection

Population-level proportions of individuals that fall at different points in the spectrum [of disease severity], from asymptomatic infection to severe disease, are often difficult to observe, but estimating these quantities can provide information about the nature and severity of the disease in a particular population. Logistic and multinomial regression techniques are often applied to infectious disease modeling of large populations and are suited to identifying variables associated with a particular disease or disease state. However, they are less appropriate for estimating infection state prevalence over time because they do not naturally accommodate known disease dynamics like duration of time an individual is infectious, heterogeneity in the risk of acquiring infection, and patterns of seasonality. We propose a Bayesian compartmental model to estimate latent infection state prevalence over time that easily incorporates known disease dynamics. We demonstrate how and why a stochastic compartmental model is a better approach for determining infection state proportions than multinomial regression is by using a novel method for estimating Bayes factors for models with high-dimensional parameter spaces. We provide an example using visceral leishmaniasis in Brazil and present an empirically-adjusted reproductive number for the infection.

Assessing inference of the basic reproduction number in an SIR model incorporating a growth-scaling parameter

The standard mass action, which assumes that infectious disease transmission occurs in well-mixed populations, is popular for formulating compartmental epidemic models. Compartmental epidemic models often follow standard mass action for simplicity and to gain insight into transmission dynamics as it often performs well at reproducing disease dynamics in large populations. In this work, we formulate discrete time stochastic susceptible-infected-removed models with linear (standard) and nonlinear mass action structures to mimic varying mixing levels. Using simulations and real epidemic data, we demonstrate the sensitivity of the basic reproduction number to these mathematical structures of the force of infection. Our results suggest the need to consider nonlinear mass action in order to generate more accurate estimates of the basic reproduction number although its uncertainty increases due to the addition of one growth scaling parameter.

Approximate Bayesian computation for spatial SEIR(S) epidemic models

Approximate Bayesia n Computation (ABC) provides an attractive approach to estimation in complex Bayesian inferential problems for which evaluation of the kernel of the posterior distribution is impossible or computationally expensive. These highly parallelizable techniques have been successfully applied to many fields, particularly in cases where more traditional approaches such as Markov chain Monte Carlo (MCMC) are impractical. In this work, we demonstrate the application of approximate Bayesian inference to spatially heterogeneous Susceptible-Exposed-Infectious-Removed (SEIR) stochastic epidemic models. These models have a tractable posterior distribution, however MCMC techniques nevertheless become computationally infeasible for moderately sized problems. We discuss the practical implementation of these techniques via the open source ABSEIR package for R. The performance of ABC relative to traditional MCMC methods in a small problem is explored under simulation, as well as in the spatially heterogeneous context of the 2014 epidemic of Chikungunya in the Americas.

Integrating independent spatio-temporal replications to assess population trends in disease spread

Glaucoma is the second leading cause of blindness in the USA. A visual field test (perimetry) is used to sample and quantitate visual field function in preselected regions in the eye. These regions can be considered a spatial field with replications across independently measured individuals. At return visits, a new set of visual field measurements is obtained producing a subject specific spatio-temporal dataset. We develop a Bayesian hierarchical modeling framework to analyze these spatio-temporal datasets both for individual level spread and as aggregate population level trends. Our model extends previous research utilizing a dimension reduction matrix and individual specific latent variables. Human characteristics are incorporated into the model to help explain glaucoma progression. One beneficial product of our model is smoothed estimates for individuals. We also specify how progression rates are computed for monitoring purposes so that clinicians can track changes and predict forward in time. Copyright © 2016 John Wiley & Sons, Ltd.

The role of interconnectivity in control of an Ebola epidemic

Several West African countries - Liberia, Sierra Leone and Guinea - experienced significant morbidity and mortality during the largest Ebola epidemic to date, from late 2013 through 2015. The extent of the epidemic was fueled by outbreaks in large urban population centers as well as movement of the pathogen between populations. During the epidemic there was no known vaccine or drug, so effective disease control required coordinated efforts that include both standard medical and community practices such as hospitalization, quarantine and safe burials. Due to the high connectivity of the region, control of the epidemic not only depended on internal strategies but also was impacted by neighboring countries. In this paper, we use a deterministic framework to examine the role of movement between two populations in the overall success of practices designed to minimize the extent of Ebola epidemics. We find that it is possible for even small amounts of intermixing between populations to positively impact the control of an epidemic on a more global scale.