Bayesian analysis of botanical epidemics using stochastic compartmental models

G. J. Gibson, A. Kleczkowski, C. A. Gilligan

Research output: Contribution to journalArticle

Abstract

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.

Original languageEnglish
Pages (from-to)12120-12124
Number of pages5
JournalProceedings of the National Academy of Sciences
Volume101
Issue number33
DOIs
Publication statusPublished - 17 Aug 2004

Fingerprint

Bayesian analysis
damping-off
Markov chain
biological control
pathogen
fungus
parameter
infection

Keywords

  • Biological control
  • Markov chain methods
  • Plant-pathogen systems

Cite this

@article{aa1ffc5a4e1c4e01aaf40bed386ead8a,
title = "Bayesian analysis of botanical epidemics using stochastic compartmental models",
abstract = "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.",
keywords = "Biological control, Markov chain methods, Plant-pathogen systems",
author = "Gibson, {G. J.} and A. Kleczkowski and Gilligan, {C. A.}",
year = "2004",
month = "8",
day = "17",
doi = "10.1073/pnas.0400829101",
language = "English",
volume = "101",
pages = "12120--12124",
journal = "Proceedings of the National Academy of Sciences",
issn = "0027-8424",
number = "33",

}

Bayesian analysis of botanical epidemics using stochastic compartmental models. / Gibson, G. J.; Kleczkowski, A.; Gilligan, C. A.

In: Proceedings of the National Academy of Sciences, Vol. 101, No. 33, 17.08.2004, p. 12120-12124.

Research output: Contribution to journalArticle

TY - JOUR

T1 - Bayesian analysis of botanical epidemics using stochastic compartmental models

AU - Gibson, G. J.

AU - Kleczkowski, A.

AU - Gilligan, C. A.

PY - 2004/8/17

Y1 - 2004/8/17

N2 - 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.

AB - 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.

KW - Biological control

KW - Markov chain methods

KW - Plant-pathogen systems

UR - http://www.scopus.com/inward/record.url?scp=4344628364&partnerID=8YFLogxK

U2 - 10.1073/pnas.0400829101

DO - 10.1073/pnas.0400829101

M3 - Article

VL - 101

SP - 12120

EP - 12124

JO - Proceedings of the National Academy of Sciences

JF - Proceedings of the National Academy of Sciences

SN - 0027-8424

IS - 33

ER -