Bayesian analysis of botanical epidemics using stochastic compartmental models.
A stochastic model for an epidemic, incorporating susceptible, latent, and infectious states, is developed. The model represents primary and secondary infection rates and a time-varying host susceptibility with applications to a wide range of epidemiological systems. A Markov chain Monte Carlo algorithm is presented that allows the model to be fitted to experimental observations within a Bayesian framework. The approach allows the uncertainty in unobserved aspects of the process to be represented in the parameter posterior densities. The methods are applied to experimental observations of damping-off of radish (Raphanus sativus) caused by the fungal pathogen Rhizoctonia solani, in the presence and absence of the antagonistic fungus Trichoderma viride, a biological control agent that has previously been shown to affect the rate of primary infection by using a maximum-likelihood estimate for a simpler model with no allowance for a latent period. Using the Bayesian analysis, we are able to estimate the latent period from population data, even when there is uncertainty in discriminating infectious from latently infected individuals in data collection. We also show that the inference that T. viride can control primary, but not secondary, infection is robust to inclusion of the latent period in the model, although the absolute values of the parameters change. Some refinements and potential difficulties with the Bayesian approach in this context, when prior information on parameters is lacking, are discussed along with broader applications of the methods to a wide range of epidemiological systems.
Department of Actuarial Mathematics and Statistics, Heriot Watt University, Riccarton, Edinburgh EH14 4AS, United Kingdom. g.j.gibson@ma.hw.ac.uk
G J
GJ
Gibson
A
A
Kleczkowski
C A
CA
Gilligan
D000465:Algorithms; D001499:Bayes Theorem; D008390:Markov Chains; D008954:Models, Biological; D009010:Monte Carlo Method; D010572:Pest Control, Biological; D010935:Plant Diseases / Q000382:microbiology / Q000706:statistics & numerical data*; D031224:Raphanus / Q000382:microbiology; D012232:Rhizoctonia / Q000472:pathogenicity; D013269:Stochastic Processes; D013997:Time Factors; D014242:Trichoderma / Q000502:physiology