Markov Chain Monte Carlo Methods for Fitting Spatiotemporal Stochastic Models in Plant Epidemiology

Strategies for controlling plant epidemics are investigated by fitting continuous time spatiotemporal stochastic models to data consisting of maps of disease incidence observed at discrete times. Markov chain Monte Carlo methods are used for fitting two such models to data describing the spread of citrus tristeza virus (CTV) in an orchard. The approach overcomes some of the difficulties encountered when fitting stochastic models to infrequent observations of a continuous process. The results of the analysis cast doubt on the effectiveness of a strategy identified from a previous spatial analysis of the CTV data. Extensions of the approaches to more general models and other problems are also considered.     

Statistical analysis of an endemic disease from a capture–recapture experiment

There are a number of statistical techniques for analysing epidemic outbreaks. However, many diseases are endemic within populations and the analysis of such diseases is complicated by changing population demography. Motivated by the spread of cowpox amongst rodent populations, a combined mathematical model for population and disease dynamics is introduced. A Markov chain Monte Carlo algorithm is then constructed to make statistical inference for the model based on data being obtained from a capture–recapture experiment. The statistical analysis is used to identify the key elements in the spread of the cowpox virus.     

Statement of Withdrawal: Statistical analysis of an endemic disease from a capture–recapture experiment

There are a number of statistical techniques for analysing epidemic outbreaks. However, many diseases are endemic within populations and the analysis of such diseases are complicated by changing population demography. Motivated by the spread of cowpox among rodent populations, a combined mathematical model for population and disease dynamics is introduced. An MCMC algorithm is then constructed to make statistical inference for the model based on data being obtained from a capture–recapture experiment. The statistical analysis is used to identify the key elements in the spread of the cowpox virus.     

Predicting infectious disease outbreak risk via migratory waterfowl vectors

The spread of an emerging infectious disease is a major public health threat. Given the uncertainties associated with vector-borne diseases, in terms of vector dynamics and disease transmission, it is critical to develop statistical models that address how and when such an infectious disease could spread throughout a region such as the USA. This paper considers a spatio-temporal statistical model for how an infectious disease could be carried into the USA by migratory waterfowl vectors during their seasonal migration and, ultimately, the risk of transmission of such a disease to domestic fowl. Modeling spatio-temporal data of this type is inherently difficult given the uncertainty associated with observations, complexity of the dynamics, high dimensionality of the underlying process, and the presence of excessive zeros. In particular, the spatio-temporal dynamics of the waterfowl migration are developed by way of a two-tiered functional temporal and spatial dimension reduction procedure that captures spatial and seasonal trends, as well as regional dynamics. Furthermore, the model relates the migration to a population of poultry farms that are known to be susceptible to such diseases, and is one of the possible avenues toward transmission to domestic poultry and humans. The result is a predictive distribution of those counties containing poultry farms that are at the greatest risk of having the infectious disease infiltrate their flocks assuming that the migratory population was infected. The model naturally fits into the hierarchical Bayesian framework.     

Relative risk analysis of dengue cases using convolution extended into spatio-temporal model

Dengue Hemmorage Fever (DHF) cases have become a serious problem every year in tropical countries such as Indonesia. Understanding the dynamic spread of the disease is essential in order to find an effective strategy in controlling its spread. In this study, a convolution (Poisson-lognormal) model that integrates both uncorrelated and correlated random effects was developed. A spatial–temporal convolution model to accomodate both spatial and temporal variations of the disease spread dynamics was considered. The model was applied to the DHF cases in the city of Kendari, Indonesia. DHF data for 10 districts during the period 2007–2010 were collected from the health services. The data of rainfall and population density were obtained from the local offices in Kendari. The numerical experiments indicated that both the rainfall and the population density played an important role in the increasing DHF cases in the city of Kendari. The result suggested that DHF cases mostly occured in January, the wet session with high rainfall, and in Kadia, the densest district in the city. As people in the city have high mobility while dengue mosquitoes tend to stay localized in their area, the best intervention is in January and in the district of Kadia.     

Spatial measurement error in infectious disease models

Individual-level models (ILMs) for infectious disease can be used to model disease spread between individuals while taking into account important covariates. One important covariate in determining the risk of infection transfer can be spatial location. At the same time, measurement error is a concern in many areas of statistical analysis, and infectious disease modelling is no exception. In this paper, we are concerned with the issue of measurement error in the recorded location of individuals when using a simple spatial ILM to model the spread of disease within a population. An ILM that incorporates spatial location random effects is introduced within a hierarchical Bayesian framework. This model is tested upon both simulated data and data from the UK 2001 foot-and-mouth disease epidemic. The ability of the model to successfully identify both the spatial infection kernel and the basic reproduction number (R ₀) of the disease is tested.     

ANALYSIS OF DATA FROM A SINGLE EPIDEMIC1

A model with nonrandom latent and infectious periods is suggested for epidemics in a large community. This permits a relatively complete statistical analysis of data from the spread of a single epidemic. An attractive feature of such models is the possibility of exploring how the rate of spread of the disease depends on the number of susceptibles and infectives. An application to smallpox data is included.     

Spatio-temporal modeling of infectious disease dynamics

A stochastic model, which is well suited to capture space–time dependence of an infectious disease, was employed in this study to describe the underlying spatial and temporal pattern of measles in Barisal Division, Bangladesh. The model has two components: an endemic component and an epidemic component; weights are used in the epidemic component for better accounting of the disease spread into different geographical regions. We illustrate our findings using a data set of monthly measles counts in the six districts of Barisal, from January 2000 to August 2009, collected from the Expanded Program on Immunization, Bangladesh. The negative binomial model with both the seasonal and autoregressive components was found to be suitable for capturing space–time dependence of measles in Barisal. Analyses were done using general optimization routines, which provided the maximum likelihood estimates with the corresponding standard errors.     

Perfect simulation for Reed-Frost epidemic models

The Reed-Frost epidemic model is a simple stochastic process with parameter q that describes the spread of an infectious disease among a closed population. Given data on the final outcome of an epidemic, it is possible to perform Bayesian inference for q using a simple Gibbs sampler algorithm. In this paper it is illustrated that by choosing latent variables appropriately, certain monotonicity properties hold which facilitate the use of a perfect simulation algorithm. The methods are applied to real data.     

Spatio-temporal analysis of a plant disease in a non-uniform crop: a Monte Carlo approach

Identification of the type of disease pattern and spread in a field is critical in epidemiological investigations of plant diseases. For example, an aggregation pattern of infected plants suggests that, at the time of observation, the pathogen is spreading from a proximal source. Conversely, a random pattern suggests a lack of spread from a proximal source. Most of the existing methods of spatial pattern analysis work with only one variety of plant at each location and with uniform genetic disease susceptibility across the field. Pecan orchards, used in this study, and other orchard crops are usually composed of different varieties with different levels of susceptibility to disease. A new measure is suggested to characterize the spatio-temporal transmission patterns of disease; a Monte Carlo test procedure is proposed to test whether the transmission of disease is random or aggregated. In addition, we propose a mixed-transmission model, which allows us to quantify the degree of aggregation effect.     

On some effects of variability in the westss epidemic model

The Weiss carrier-borne epidemic model is extended to incorporate random initial susceptible and carrier population sizes and general independent distributions for the infectious periods of the carriers. The qualitative effects of these variability factors on the spread of the disease are examined using different partial orderings for distribution functions.     

Epidemics in heterogeneous communities: estimation of R0 and secure vaccination coverage

A stochastic multitype model for the spread of an infectious disease in a community of heterogeneous individuals is analysed. In particular, estimates of R ₀ (the basic reproduction number) and the critical vaccination coverage are derived, where estimation is based on final size data of an outbreak in the community. It is shown that these key parameters cannot be estimated consistently from data; only upper and lower bounds can be estimated. Confidence regions for the upper bounds are derived, thus giving conservative estimates of R ₀ and the fractions necessary to vaccinate.     

A Test to Detect Within-family Infectivity when the Whole Epidemic Process is Observed

An epidemic model for the spread of an infectious disease in a population of families is considered. The score test of the hypothesis that there is no higher infectivity between family members is constructed under the assumption that the epidemic process is observed continuously up to some time t . The score process is a martingale as a function of t and by letting the number of families tend to infinity, a central limit theorem for the process can be proved. The central limit theorem not only justifies a normal approximation of the test statistic—it also suggests a smaller variance estimator than expected.     

Predicting recessions using trends in the yield spread

The yield spread, measured as the difference between long- and short-term interest rates, is widely regarded as one of the strongest predictors of economic recessions. In this paper, we propose an enhanced recession prediction model that incorporates trends in the value of the yield spread. We expect our model to generate stronger recession signals because a steadily declining value of the yield spread typically indicates growing pessimism associated with the reduced future business activity. We capture trends in the yield spread by considering both the level of the yield spread at a lag of 12 months as well as its value at each of the previous two quarters leading up to the forecast origin, and we evaluate its predictive abilities using both logit and artificial neural network models. Our results indicate that models incorporating information from the time series of the yield spread correctly predict future recession periods much better than models only considering the spread value as of the forecast origin. Furthermore, the results are strongest for our artificial neural network model and logistic regression model that includes interaction terms, which we confirm using both a blocked cross-validation technique as well as an expanding estimation window approach.     

Extended Hazard Regression Model for Reliability and Survival Analysis

We propose an extended hazard regression model which allows the spread parameter to be dependent on covariates. This allows a broad class of models which includes the most common hazard models, such as the proportional hazards model, the accelerated failure time model and a proportional hazards/accelerated failure time hybrid model with constant spread parameter. Simulations based on sub-classes of this model suggest that maximum likelihood performs well even when only small or moderate-size data sets are available and the censoring pattern is heavy. The methodology provides a broad framework for analysis of reliability and survival data. Two numerical examples illustrate the results.     

Epidemiologic network inference

Statistics and Computing - In many epidemiologic models, a disease is assumed to spread along a contact network. We aim to infer this network, in addition to the epidemiologic model parameters,...     

Modelling the spread in space and time of an airborne plant disease

Summary.  A spatiotemporal model is developed to analyse epidemics of airborne plant diseases which are spread by spores. The observations consist of measurements of the severity of disease at different times, different locations in the horizontal plane and different heights in the vegetal cover. The model describes the joint distribution of the occurrence and the severity of the disease. The three-dimensional dispersal of spores is modelled by combining a horizontal and a vertical dispersal function. Maximum likelihood combined with a parametric bootstrap is suggested to estimate the model parameters and the uncertainty that is attached to them. The spatiotemporal model is used to analyse a yellow rust epidemic in a wheatfield. In the analysis we pay particular attention to the selection and the estimation of the dispersal functions.     

Bayesian estimation for percolation models of disease spread in plant populations

Statistical methods are formulated for fitting and testing percolation-based, spatio-temporal models that are generally applicable to biological or physical processes that evolve in spatially distributed populations. The approach is developed and illustrated in the context of the spread of Rhizoctonia solani, a fungal pathogen, in radish but is readily generalized to other scenarios. The particular model considered represents processes of primary and secondary infection between nearest-neighbour hosts in a lattice, and time-varying susceptibility of the hosts. Bayesian methods for fitting the model to observations of disease spread through space and time in replicate populations are developed. These use Markov chain Monte Carlo methods to overcome the problems associated with partial observation of the process. We also consider how model testing can be achieved by embedding classical methods within the Bayesian analysis. In particular we show how a residual process, with known sampling distribution, can be defined. Model fit is then examined by generating samples from the posterior distribution of the residual process, to which a classical test for consistency with the known distribution is applied, enabling the posterior distribution of the P-value of the test used to be estimated. For the Rhizoctonia-radish system the methods confirm the findings of earlier non-spatial analyses regarding the dynamics of disease transmission and yield new evidence of environmental heterogeneity in the replicate experiments.     

随机波动HJM框架下信用利差模型及实证研究

将随机波动引入到具有违约风险的HJM模型中,基于无套利条件推导得出随机波动HJM模型框架下信用利差的漂移项限制条件,从而构建了随机波动HJM框架下的信用利差模型,并基于波动率结构的适当设定对模型进行马尔科夫简化,进而利用该模型对中国可违约债券市场信用利差的动态特性进行实证分析。结果发现：短期信用利差的动态特性具有显著的随机波动特征,而随机波动HJM框架下的信用利差模型可以很好地刻画中国可违约债券市场信用利差的动态特性。     

Spatial–temporal variation of casuarina spread in Cauvery delta and north eastern zone of Tamil Nadu,India: a spatial autoregressive model

The low forest cover and productivity are the major obstacles for mitigating the demand supply gap of raw material for forest-based industries, which could be fulfilled from a tree outside forest area. Casuarina is a multi-utile, short rotation tree which adapts to all ecosystems. The casuarina wood is predominantly demanded for fuel, construction and paper industries which is mostly preferred by farmers, traders and industries. This study explores the spatial and temporal variability of casuarina spread in mitigating the gap of demand and supply in Tamil Nadu using a spatial autoregressive model. The spread of casuarina was spatially and temporally significant, which was negatively influenced by the gross area irrigated as main and direct effects and positively in an indirect effect. An assured irrigation forces the farmers to choose traditional agricultural crops for their livelihood in their own district. The increase in the price of casuarina would increase the spread of casuarina in both own district and neighbouring districts. The spread of casuarina would augment the supply of raw material for forest-based industries.