Bayesian inference for spatio-temporal stochastic transmission of plant disease in the presence of roguing: A case study to characterise the dispersal of Flavescence dorée.

Estimating the distance at which pathogens disperse from one season to the next is crucial for designing efficient control strategies for invasive plant pathogens and a major milestone in the reduction of pesticide use in agriculture. However, we still lack such estimates for many diseases, especially for insect-vectored pathogens, such as Flavescence dorée (FD). FD is a quarantine disease threatening European vineyards. Its management is based on mandatory insecticide treatments and the removal of infected plants identified during annual surveys. This paper introduces a general statistical framework to model the epidemiological dynamics of FD in a mechanistic manner that can take into account missing hosts in surveyed fields (resulting from infected plant removals). We parameterized the model using Markov chain Monte Carlo (MCMC) and data augmentation from surveillance data gathered in Bordeaux vineyards. The data mainly consist of two snapshot maps of the infectious status of all the plants in three adjacent fields during two consecutive years. We demonstrate that heavy-tailed dispersal kernels best fit the spread of FD and that on average, 50% (resp. 80%) of new infection occurs within 10.5 m (resp. 22.2 m) of the source plant. These values are in agreement with estimates of the flying capacity of Scaphoideus titanus, the leafhopper vector of FD, reported in the literature using mark-capture techniques. Simulations of simple removal scenarios using the fitted model suggest that cryptic infection hampered FD management. Future efforts should explore whether strategies relying on reactive host removal can improve FD management.

A data driven change-point epidemic model for assessing the impact of large gathering and subsequent movement control order on COVID-19 spread in Malaysia.

The second wave of COVID-19 in Malaysia is largely attributed to a four-day mass gathering held in Sri Petaling from February 27, 2020, which contributed to an exponential rise of COVID-19 cases in the country. Starting from March 18, 2020, the Malaysian government introduced four consecutive phases of a Movement Control Order (MCO) to stem the spread of COVID-19. The MCO was implemented through various non-pharmaceutical interventions (NPIs). The reported number of cases reached its peak by the first week of April and then started to reduce, hence proving the effectiveness of the MCO. To gain a quantitative understanding of the effect of MCO on the dynamics of COVID-19, this paper develops a class of mathematical models to capture the disease spread before and after MCO implementation in Malaysia. A heterogeneous variant of the Susceptible-Exposed-Infected-Recovered (SEIR) model is developed with additional compartments for asymptomatic transmission. Further, a change-point is incorporated to model disease dynamics before and after intervention which is inferred based on data. Related statistical analyses for inference are developed in a Bayesian framework and are able to provide quantitative assessments of (1) the impact of the Sri Petaling gathering, and (2) the extent of decreasing transmission during the MCO period. The analysis here also quantitatively demonstrates how quickly transmission rates fall under effective NPI implementation within a short time period. The models and methodology used provided important insights into the nature of local transmissions to decision makers in the Ministry of Health, Malaysia.

Latent likelihood ratio tests for assessing spatial kernels in epidemic models.

One of the most important issues in the critical assessment of spatio-temporal stochastic models for epidemics is the selection of the transmission kernel used to represent the relationship between infectious challenge and spatial separation of infected and susceptible hosts. As the design of control strategies is often based on an assessment of the distance over which transmission can realistically occur and estimation of this distance is very sensitive to the choice of kernel function, it is important that models used to inform control strategies can be scrutinised in the light of observation in order to elicit possible evidence against the selected kernel function. While a range of approaches to model criticism is in existence, the field remains one in which the need for further research is recognised. In this paper, building on earlier contributions by the authors, we introduce a new approach to assessing the validity of spatial kernels-the latent likelihood ratio tests-which use likelihood-based discrepancy variables that can be used to compare the fit of competing models, and compare the capacity of this approach to detect model mis-specification with that of tests based on the use of infection-link residuals. We demonstrate that the new approach can be used to formulate tests with greater power than infection-link residuals to detect kernel mis-specification particularly when the degree of mis-specification is modest. This new tests avoid the use of a fully Bayesian approach which may introduce undesirable complications related to computational complexity and prior sensitivity.

Model diagnostics and refinement for phylodynamic models.

Phylodynamic modelling, which studies the joint dynamics of epidemiological and evolutionary processes, has made significant progress in recent years due to increasingly available genomic data and advances in statistical modelling. These advances have greatly improved our understanding of transmission dynamics of many important pathogens. Nevertheless, there remains a lack of effective, targetted diagnostic tools for systematically detecting model mis-specification. Development of such tools is essential for model criticism, refinement, and calibration. The idea of utilising latent residuals for model assessment has already been exploited in general spatio-temporal epidemiological settings. Specifically, by proposing appropriately designed non-centered, re-parameterizations of a given epidemiological process, one can construct latent residuals with known sampling distributions which can be used to quantify evidence of model mis-specification. In this paper, we extend this idea to formulate a novel model-diagnostic framework for phylodynamic models. Using simulated examples, we show that our framework may effectively detect a particular form of mis-specification in a phylodynamic model, particularly in the event of superspreading. We also exemplify our approach by applying the framework to a dataset describing a local foot-and-mouth (FMD) outbreak in the UK, eliciting strong evidence against the assumption of no within-host-diversity in the outbreak. We further demonstrate that our framework can facilitate model calibration in real-life scenarios, by proposing a within-host-diversity model which appears to offer a better fit to data than one that assumes no within-host-diversity of FMD virus.

Evidence-based controls for epidemics using spatio-temporal stochastic models in a Bayesian framework.

The control of highly infectious diseases of agricultural and plantation crops and livestock represents a key challenge in epidemiological and ecological modelling, with implemented control strategies often being controversial. Mathematical models, including the spatio-temporal stochastic models considered here, are playing an increasing role in the design of control as agencies seek to strengthen the evidence on which selected strategies are based. Here, we investigate a general approach to informing the choice of control strategies using spatio-temporal models within the Bayesian framework. We illustrate the approach for the case of strategies based on pre-emptive removal of individual hosts. For an exemplar model, using simulated data and historic data on an epidemic of Asiatic citrus canker in Florida, we assess a range of measures for prioritizing individuals for removal that take account of observations of an emerging epidemic. These measures are based on the potential infection hazard a host poses to susceptible individuals (hazard), the likelihood of infection of a host (risk) and a measure that combines both the hazard and risk (threat). We find that the threat measure typically leads to the most effective control strategies particularly for clustered epidemics when resources are scarce. The extension of the methods to a range of other settings is discussed. A key feature of the approach is the use of functional-model representations of the epidemic model to couple epidemic trajectories under different control strategies. This induces strong positive correlations between the epidemic outcomes under the respective controls, serving to reduce both the variance of the difference in outcomes and, consequently, the need for extensive simulation.

A mechanistic spatio-temporal framework for modelling individual-to-individual transmission-With an application to the 2014-2015 West Africa Ebola outbreak.

In recent years there has been growing availability of individual-level spatio-temporal disease data, particularly due to the use of modern communicating devices with GPS tracking functionality. These detailed data have been proven useful for inferring disease transmission to a more refined level than previously. However, there remains a lack of statistically sound frameworks to model the underlying transmission dynamic in a mechanistic manner. Such a development is particularly crucial for enabling a general epidemic predictive framework at the individual level. In this paper we propose a new statistical framework for mechanistically modelling individual-to-individual disease transmission in a landscape with heterogeneous population density. Our methodology is first tested using simulated datasets, validating our inferential machinery. The methodology is subsequently applied to data that describes a regional Ebola outbreak in Western Africa (2014-2015). Our results show that the methods are able to obtain estimates of key epidemiological parameters that are broadly consistent with the literature, while revealing a significantly shorter distance of transmission. More importantly, in contrast to existing approaches, we are able to perform a more general model prediction that takes into account the susceptible population. Finally, our results show that, given reasonable scenarios, the framework can be an effective surrogate for susceptible-explicit individual models which are often computationally challenging.

A Systematic Bayesian Integration of Epidemiological and Genetic Data.

Genetic sequence data on pathogens have great potential to inform inference of their transmission dynamics ultimately leading to better disease control. Where genetic change and disease transmission occur on comparable timescales additional information can be inferred via the joint analysis of such genetic sequence data and epidemiological observations based on clinical symptoms and diagnostic tests. Although recently introduced approaches represent substantial progress, for computational reasons they approximate genuine joint inference of disease dynamics and genetic change in the pathogen population, capturing partially the joint epidemiological-evolutionary dynamics. Improved methods are needed to fully integrate such genetic data with epidemiological observations, for achieving a more robust inference of the transmission tree and other key epidemiological parameters such as latent periods. Here, building on current literature, a novel Bayesian framework is proposed that infers simultaneously and explicitly the transmission tree and unobserved transmitted pathogen sequences. Our framework facilitates the use of realistic likelihood functions and enables systematic and genuine joint inference of the epidemiological-evolutionary process from partially observed outbreaks. Using simulated data it is shown that this approach is able to infer accurately joint epidemiological-evolutionary dynamics, even when pathogen sequences and epidemiological data are incomplete, and when sequences are available for only a fraction of exposures. These results also characterise and quantify the value of incomplete and partial sequence data, which has important implications for sampling design, and demonstrate the abilities of the introduced method to identify multiple clusters within an outbreak. The framework is used to analyse an outbreak of foot-and-mouth disease in the UK, enhancing current understanding of its transmission dynamics and evolutionary process.

Bayesian inference for an emerging arboreal epidemic in the presence of control.

The spread of Huanglongbing through citrus groves is used as a case study for modeling an emerging epidemic in the presence of a control. Specifically, the spread of the disease is modeled as a susceptible-exposed-infectious-detected-removed epidemic, where the exposure and infectious times are not observed, detection times are censored, removal times are known, and the disease is spreading through a heterogeneous host population with trees of different age and susceptibility. We show that it is possible to characterize the disease transmission process under these conditions. Two innovations in our work are (i) accounting for control measures via time dependence of the infectious process and (ii) including seasonal and host age effects in the model of the latent period. By estimating parameters in different subregions of a large commercially cultivated orchard, we establish a temporal pattern of invasion, host age dependence of the dispersal parameters, and a close to linear relationship between primary and secondary infectious rates. The model can be used to simulate Huanglongbing epidemics to assess economic costs and potential benefits of putative control scenarios.

Bayesian analysis for inference of an emerging epidemic: citrus canker in urban landscapes.

Outbreaks of infectious diseases require a rapid response from policy makers. The choice of an adequate level of response relies upon available knowledge of the spatial and temporal parameters governing pathogen spread, affecting, amongst others, the predicted severity of the epidemic. Yet, when a new pathogen is introduced into an alien environment, such information is often lacking or of no use, and epidemiological parameters must be estimated from the first observations of the epidemic. This poses a challenge to epidemiologists: how quickly can the parameters of an emerging disease be estimated? How soon can the future progress of the epidemic be reliably predicted? We investigate these issues using a unique, spatially and temporally resolved dataset for the invasion of a plant disease, Asiatic citrus canker in urban Miami. We use epidemiological models, Bayesian Markov-chain Monte Carlo, and advanced spatial statistical methods to analyse rates and extent of spread of the disease. A rich and complex epidemic behaviour is revealed. The spatial scale of spread is approximately constant over time and can be estimated rapidly with great precision (although the evidence for long-range transmission is inconclusive). In contrast, the rate of infection is characterised by strong monthly fluctuations that we associate with extreme weather events. Uninformed predictions from the early stages of the epidemic, assuming complete ignorance of the future environmental drivers, fail because of the unpredictable variability of the infection rate. Conversely, predictions improve dramatically if we assume prior knowledge of either the main environmental trend, or the main environmental events. A contrast emerges between the high detail attained by modelling in the spatiotemporal description of the epidemic and the bottleneck imposed on epidemic prediction by the limits of meteorological predictability. We argue that identifying such bottlenecks will be a fundamental step in future modelling of weather-driven epidemics.

New model diagnostics for spatio-temporal systems in epidemiology and ecology.

A cardinal challenge in epidemiological and ecological modelling is to develop effective and easily deployed tools for model assessment. The availability of such methods would greatly improve understanding, prediction and management of disease and ecosystems. Conventional Bayesian model assessment tools such as Bayes factors and the deviance information criterion (DIC) are natural candidates but suffer from important limitations because of their sensitivity and complexity. Posterior predictive checks, which use summary statistics of the observed process simulated from competing models, can provide a measure of model fit but appropriate statistics can be difficult to identify. Here, we develop a novel approach for diagnosing mis-specifications of a general spatio-temporal transmission model by embedding classical ideas within a Bayesian analysis. Specifically, by proposing suitably designed non-centred parametrization schemes, we construct latent residuals whose sampling properties are known given the model specification and which can be used to measure overall fit and to elicit evidence of the nature of mis-specifications of spatial and temporal processes included in the model. This model assessment approach can readily be implemented as an addendum to standard estimation algorithms for sampling from the posterior distributions, for example Markov chain Monte Carlo. The proposed methodology is first tested using simulated data and subsequently applied to data describing the spread of Heracleum mantegazzianum (giant hogweed) across Great Britain over a 30-year period. The proposed methods are compared with alternative techniques including posterior predictive checking and the DIC. Results show that the proposed diagnostic tools are effective in assessing competing stochastic spatio-temporal transmission models and may offer improvements in power to detect model mis-specifications. Moreover, the latent-residual framework introduced here extends readily to a broad range of ecological and epidemiological models.

Non-exponential tolerance to infection in epidemic systems--modeling, inference, and assessment.

The transmission dynamics of infectious diseases have been traditionally described through a time-inhomogeneous Poisson process, thus assuming exponentially distributed levels of disease tolerance following the Sellke construction. Here we focus on a generalization using Weibull individual tolerance thresholds under the susceptible-exposed-infectious-removed class of models which is widely employed in epidemics. Applications with experimental foot-and-mouth disease and historical smallpox data are discussed, and simulation results are presented. Inference is carried out using Markov chain Monte Carlo methods following a Bayesian approach. Model evaluation is performed, where the adequacy of the models is assessed using methodology based on the properties of Bayesian latent residuals, and comparison between 2 candidate models is also considered using a latent likelihood ratio-type test that avoids problems encountered with relevant methods based on Bayes factors.

Applications of percolation theory to fungal spread with synergy.

There is increasing interest in the use of the percolation paradigm to analyse and predict the progress of disease spreading in spatially structured populations of animals and plants. The wider utility of the approach has been limited, however, by several restrictive assumptions, foremost of which is a strict requirement for simple nearest-neighbour transmission, in which the disease history of an individual is influenced only by that of its neighbours. In a recent paper, the percolation paradigm has been generalized to incorporate synergistic interactions in host infectivity and susceptibility, and the impact of these interactions on the invasive dynamics of an epidemic has been demonstrated. In the current paper, we elicit evidence that such synergistic interactions may underlie transmission dynamics in real-world systems by first formulating a model for the spread of a ubiquitous parasitic and saprotrophic fungus through replicated populations of nutrient sites and subsequently fitting and testing the model using data from experimental microcosms. Using Bayesian computational methods for model fitting, we demonstrate that synergistic interactions are necessary to explain the dynamics observed in the replicate experiments. The broader implications of this work in identifying disease-control strategies that deflect epidemics from invasive to non-invasive regimes are discussed.

Modeling the consistency of hypoglycemic symptoms: high variability in diabetes.

BACKGROUND: The aim of the present study was to examine symptoms of hypoglycemia, to develop a method to quantify individual differences in the consistency of symptom reporting, and to investigate which factors affect these differences. METHODS: Participants recorded their symptoms with every episode of hypoglycemia over a 9-12-month period. A novel logistic-type latent variable model was developed to quantify the consistency of each individual's symptom complex and was used to analyze data from 59 subjects (median age, 57.5 years [range, 22-74 years], 65% male, 77% type 1 diabetes) who had experienced 19 or more hypoglycemic episodes. The association between the calculated consistency parameter and age, sex, type and duration of diabetes, and C-peptide and serum angiotensin converting enzyme concentration was examined using a generalized linear model. Analyses were performed under a Bayesian framework, using Markov chain Monte-Carlo methodology. RESULTS: Individuals exhibited substantial differences in between-episode consistency of their symptom reports, with only a small number of individuals exhibiting high levels of consistency. Men were more consistent than women. No other factors affected consistency in patients with normal hypoglycemia awareness. CONCLUSIONS: By using a novel stochastic model as a quantitative tool to compare the consistency of hypoglycemic symptom reporting, much greater intra-individual variability in symptom reporting was identified than has been recognized previously. This is relevant when instructing patients on identification of hypoglycemic symptoms and in interpreting symptomatic responses during experimentally induced hypoglycemia.

Nanometric depth resolution from multi-focal images in microscopy.

We describe a method for tracking the position of small features in three dimensions from images recorded on a standard microscope with an inexpensive attachment between the microscope and the camera. The depth-measurement accuracy of this method is tested experimentally on a wide-field, inverted microscope and is shown to give approximately 8 nm depth resolution, over a specimen depth of approximately 6 µm, when using a 12-bit charge-coupled device (CCD) camera and very bright but unresolved particles. To assess low-flux limitations a theoretical model is used to derive an analytical expression for the minimum variance bound. The approximations used in the analytical treatment are tested using numerical simulations. It is concluded that approximately 14 nm depth resolution is achievable with flux levels available when tracking fluorescent sources in three dimensions in live-cell biology and that the method is suitable for three-dimensional photo-activated localization microscopy resolution. Sub-nanometre resolution could be achieved with photon-counting techniques at high flux levels.

Multilayered 3D Lidar image construction using spatial models in a Bayesian framework.

Standard 3D imaging systems process only a single return at each pixel from an assumed single opaque surface. However, there are situations when the laser return consists of multiple peaks due to the footprint of the beam impinging on a target with surfaces distributed in depth or with semi-transparent surfaces. If all these returns are processed, a more informative multi-layered 3D image is created. We propose a unified theory of pixel processing for Lidar data using a Bayesian approach that incorporates spatial constraints through a Markov Random Field with a Potts prior model. This allows us to model uncertainty about the underlying spatial process. To palliate some inherent deficiencies of this prior model, we also introduce two proposal distributions, one based on spatial mode jumping, the other on a spatial birth/death process. The different parameters of the several returns are estimated using reversible jump Markov chain Monte Carlo (RJMCMC) techniques in combination with an adaptive strategy of delayed rejection to improve the estimates of the parameters.

Optimal observation times in experimental epidemic processes.

This article describes a method for choosing observation times for stochastic processes to maximise the expected information about their parameters. Two commonly used models for epidemiological processes are considered: a simple death process and a susceptible-infected (SI) epidemic process with dual sources for infection spreading within and from outwith the population. The search for the optimal design uses Bayesian computational methods to explore the joint parameter-data-design space, combined with a method known as moment closure to approximate the likelihood to make the acceptance step efficient. For the processes considered, a small number of optimally chosen observations are shown to yield almost as much information as much more intensively observed schemes that are commonly used in epidemiological experiments. Analysis of the simple death process allows a comparison between the full Bayesian approach and locally optimal designs around a point estimate from the prior based on asymptotic results. The robustness of the approach to misspecified priors is demonstrated for the SI epidemic process, for which the computational intractability of the likelihood precludes locally optimal designs. We show that optimal designs derived by the Bayesian approach are similar for observational studies of a single epidemic and for studies involving replicated epidemics in independent subpopulations. Different optima result, however, when the objective is to maximise the gain in information based on informative and non-informative priors: this has implications when an experiment is designed to convince a naïve or sceptical observer rather than consolidate the belief of an informed observer. Some extensions to the methods, including the selection of information criteria and extension to other epidemic processes with transition probabilities, are briefly addressed.

Estimation of multiple transmission rates for epidemics in heterogeneous populations.

One of the principal challenges in epidemiological modeling is to parameterize models with realistic estimates for transmission rates in order to analyze strategies for control and to predict disease outcomes. Using a combination of replicated experiments, Bayesian statistical inference, and stochastic modeling, we introduce and illustrate a strategy to estimate transmission parameters for the spread of infection through a two-phase mosaic, comprising favorable and unfavorable hosts. We focus on epidemics with local dispersal and formulate a spatially explicit, stochastic set of transition probabilities using a percolation paradigm for a susceptible-infected (S-I) epidemiological model. The S-I percolation model is further generalized to allow for multiple sources of infection including external inoculum and host-to-host infection. We fit the model using Bayesian inference and Markov chain Monte Carlo simulation to successive snapshots of damping-off disease spreading through replicated plant populations that differ in relative proportions of favorable and unfavorable hosts and with time-varying rates of transmission. Epidemiologically plausible parametric forms for these transmission rates are compared by using the deviance information criterion. Our results show that there are four transmission rates for a two-phase system, corresponding to each combination of infected donor and susceptible recipient. Knowing the number and magnitudes of the transmission rates allows the dominant pathways for transmission in a heterogeneous population to be identified. Finally, we show how failure to allow for multiple transmission rates can overestimate or underestimate the rate of spread of epidemics in heterogeneous environments, which could lead to marked failure or inefficiency of control strategies.

Bayesian analysis of Lidar signals with multiple returns.

Time-Correlated Single Photon Counting and Burst Illumination Laser data can be used for range profiling and target classification. In general, the problem is to analyse the response from a histogram of either photon counts or integrated intensities to assess the number, positions and amplitudes of the reflected returns from object surfaces. The goal of our work is a complete characterisation of the 3D surfaces viewed by the laser imaging system. The authors present a unified theory of pixel processing that is applicable to both approaches based on a Bayesian framework which allows for careful and thorough treatment of all types of uncertainties associated with the data. We use reversible jump Markov chain Monte Carlo (RJMCMC) techniques to evaluate the posterior distribution of the parameters and to explore spaces with different dimensionality. Further, we use a delayed rejection step to allow the generated Markov chain to mix better through the use of different proposal distributions. The approach is demonstrated on simulated and real data, showing that the return parameters can be estimated to a high degree of accuracy. We also show some practical examples from both near and far range depth imaging.

Bayesian inference for the spatio-temporal invasion of alien species.

In this paper we develop a Bayesian approach to parameter estimation in a stochastic spatio-temporal model of the spread of invasive species across a landscape. To date, statistical techniques, such as logistic and autologistic regression, have outstripped stochastic spatio-temporal models in their ability to handle large numbers of covariates. Here we seek to address this problem by making use of a range of covariates describing the bio-geographical features of the landscape. Relative to regression techniques, stochastic spatio-temporal models are more transparent in their representation of biological processes. They also explicitly model temporal change, and therefore do not require the assumption that the species' distribution (or other spatial pattern) has already reached equilibrium as is often the case with standard statistical approaches. In order to illustrate the use of such techniques we apply them to the analysis of data detailing the spread of an invasive plant, Heracleum mantegazzianum, across Britain in the 20th Century using geo-referenced covariate information describing local temperature, elevation and habitat type. The use of Markov chain Monte Carlo sampling within a Bayesian framework facilitates statistical assessments of differences in the suitability of different habitat classes for H. mantegazzianum, and enables predictions of future spread to account for parametric uncertainty and system variability. Our results show that ignoring such covariate information may lead to biased estimates of key processes and implausible predictions of future distributions.

Novel bivariate moment-closure approximations.

Nonlinear stochastic models are typically intractable to analytic solutions and hence, moment-closure schemes are used to provide approximations to these models. Existing closure approximations are often unable to describe transient aspects caused by extinction behaviour in a stochastic process. Recent work has tackled this problem in the univariate case. In this study, we address this problem by introducing novel bivariate moment-closure methods based on mixture distributions. Novel closure approximations are developed, based on the beta-binomial, zero-modified distributions and the log-Normal, designed to capture the behaviour of the stochastic SIS model with varying population size, around the threshold between persistence and extinction of disease. The idea of conditional dependence between variables of interest underlies these mixture approximations. In the first approximation, we assume that the distribution of infectives (I) conditional on population size (N) is governed by the beta-binomial and for the second form, we assume that I is governed by zero-modified beta-binomial distribution where in either case N follows a log-Normal distribution. We analyse the impact of coupling and inter-dependency between population variables on the behaviour of the approximations developed. Thus, the approximations are applied in two situations in the case of the SIS model where: (1) the death rate is independent of disease status; and (2) the death rate is disease-dependent. Comparison with simulation shows that these mixture approximations are able to predict disease extinction behaviour and describe transient aspects of the process.