Geographically dependent individual-level models for infectious diseases transmission

Ensemble learning methods of inference for spatially stratified infectious disease systems

Individual level models are a class of mechanistic models that are widely used to infer infectious disease transmission dynamics. These models incorporate individual level covariate information accounting for population heterogeneity and are generally fitted in a Bayesian Markov chain Monte Carlo (MCMC) framework. However, Bayesian MCMC methods of inference are computationally expensive for large data sets. This issue becomes more severe when applied to infectious disease data collected from spatially heterogeneous populations, as the number of covariates increases. In addition, summary statistics over the global population may not capture the true spatio-temporal dynamics of disease transmission. In this study we propose to use ensemble learning methods to predict epidemic generating models instead of time consuming Bayesian MCMC method. We apply these methods to infer disease transmission dynamics over spatially clustered populations, considering the clusters as natural strata instead of a global population. We compare the performance of two tree-based ensemble learning techniques: random forest and gradient boosting. These methods are applied to the 2001 foot-and-mouth disease epidemic in the U.K. and evaluated using simulated data from a clustered population. It is shown that the spatially clustered data can help to predict epidemic generating models more accurately than the global data.

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.

Analyzing COVID-19 data in the Canadian province of Manitoba: A new approach

The basic homogeneous SEIR (susceptible-exposed-infected-removed) model is a commonly used compartmental model for analysing infectious diseases such as influenza and COVID-19. However, in the homogeneous SEIR model, it is assumed that the population of study is homogeneous and, one cannot incorporate individual-level information (e.g., location of infected people, distance between susceptible and infected individuals, vaccination status) which may be important in predicting new disease cases. Recently, a geographically-dependent individual-level model (GD-ILM) within an SEIR framework was developed for when both regional and individual-level spatial data are available. In this paper, we propose to use an SEIR GD-ILM for each health region of Manitoba (central Canadian province) population to analyse the COVID-19 data. As different health regions of the population under study may act differently, we assume that each health region has its own corresponding parameters determined by a homogeneous SEIR model (such as contact rate, latent period, infectious period). A Monte Carlo Expectation Conditional Maximization (MCECM) algorithm is used for inference. Using estimated parameters we predict the infection rate at each health region of Manitoba over time to identify highly risk local geographical areas. Performance of the proposed approach is also evaluated through simulation studies.

Contextual contact tracing based on stochastic compartment modeling and spatial risk assessment

The current situation of COVID-19 highlights the paramount importance of infectious disease surveillance, which necessitates early monitoring for effective response. Policymakers are interested in data insights identifying high-risk areas as well as individuals to be quarantined, especially as the public gets back to their normal routine. We investigate both requirements by the implementation of disease outbreak modeling and exploring its induced dynamic spatial risk in form of risk assessment, along with its real-time integration back into the disease model. This paper implements a contact tracing-based stochastic compartment model as a baseline, to further modify the existing setup to include the spatial risk. This modification of each individual-level contact's intensity to be dependent on its spatial location has been termed as Contextual Contact Tracing. The results highlight that the inclusion of spatial context tends to send more individuals into quarantine which reduces the overall spread of infection. With a simulated example of an induced spatial high-risk, it is highlighted that the new spatio-SIR model can act as a tool to empower the analyst with a capability to explore disease dynamics from a spatial perspective. We conclude that the proposed spatio-SIR tool can be of great help for policymakers to know the consequences of their decision prior to their implementation.

A Practical Introduction to Mechanistic Modeling of Disease Transmission in Veterinary Science

Computer-based disease spread models are frequently used in veterinary science to simulate disease spread. They are used to predict the impacts of the disease, plan and assess surveillance, or control strategies, and provide insights about disease causation by comparing model outputs with real life data. There are many types of disease spread models, and here we present and describe the implementation of a particular type: individual-based models. Our aim is to provide a practical introduction to building individual-based disease spread models. We also introduce code examples with the goal to make these techniques more accessible to those who are new to the field. We describe the important steps in building such models before, during and after the programming stage, including model verification (to ensure that the model does what was intended), validation (to investigate whether the model results reflect the modeled system), and convergence analysis (to ensure models of endemic diseases are stable before outputs are collected). We also describe how sensitivity analysis can be used to assess the potential impact of uncertainty about model parameters. Finally, we provide an overview of some interesting recent developments in the field of disease spread models.

Spatial modeling of individual-level infectious disease transmission: Tuberculosis data in Manitoba, Canada

Geographically dependent individual level models (GD-ILMs) are a class of statistical models that can be used to study the spread of infectious disease through a population in discrete-time in which covariates can be measured both at individual and area levels. The typical ILMs to illustrate spatial data are based on the distance between susceptible and infectious individuals. A key feature of GD-ILMs is that they take into account the spatial location of the individuals in addition to the distance between susceptible and infectious individuals. As a motivation of this article, we consider tuberculosis (TB) data which is an infectious disease which can be transmitted through individuals. It is also known that certain areas/demographics/communities have higher prevalent of TB (see Section 4 for more details). It is also of interest of policy makers to identify those areas with higher infectivity rate of TB for possible preventions. Therefore, we need to analyze this data properly to address those concerns. In this article, the expectation conditional maximization algorithm is proposed for estimating the parameters of GD-ILMs to be able to predict the areas with the highest average infectivity rates of TB. We also evaluate the performance of our proposed approach through some simulations. Our simulation results indicate that the proposed method provides reliable estimates of parameters which confirms accuracy of the infectivity rates.