Parameter inference for an individual based model of chytridiomycosis in frogs.

Individual based models (IBMs) and Agent based models (ABMs) have become widely used tools to understand complex biological systems. However, general methods of parameter inference for IBMs are not available. In this paper we show that it is possible to address this problem with a traditional likelihood-based approach, using an example of an IBM developed to describe the spread of chytridiomycosis in a population of frogs as a case study. We show that if the IBM satisfies certain criteria we can find the likelihood (or posterior) analytically, and use standard computational techniques, such as MCMC, for parameter inference.

[1]  J. Longcore,et al.  BATRACHOCHYTRIUM DENDROBATIDIS GEN. ET SP. NOV., A CHYTRID PATHOGENIC TO AMPHIBIANS , 1999 .

[2]  V. Dukic,et al.  Uncertainty in predictions of disease spread and public health responses to bioterrorism and emerging diseases , 2006, Proceedings of the National Academy of Sciences.

[3]  Volker Grimm,et al.  Using pattern-oriented modeling for revealing hidden information: a key for reconciling ecological theory and application , 2003 .

[4]  M. Saunders,et al.  Plant-Provided Food for Carnivorous Insects: a Protective Mutualism and Its Applications , 2009 .

[5]  W. Kendall,et al.  Seeking a second opinion: uncertainty in disease ecology. , 2010, Ecology letters.

[6]  R. Alford,et al.  Experimental Exposures of Boreal Toads (Bufo boreas) to a Pathogenic Chytrid Fungus (Batrachochytrium dendrobatidis) , 2006, EcoHealth.

[7]  G. Huse.  Individual‐based Modeling and Ecology , 2008 .

[8]  R. Alford,et al.  Life-history trade-offs influence disease in changing climates: strategies of an amphibian pathogen. , 2008, Ecology.

[9]  D. Merl,et al.  A Statistical Framework for the Adaptive Management of Epidemiological Interventions , 2009, PloS one.

[10]  Stephen R Carpenter,et al.  Individual-Based Model for Growth of Young-of-the-Year Walleye: A Piece of the Recruitment Puzzle. , 1991, Ecological applications : a publication of the Ecological Society of America.

[11]  R. Freckleton,et al.  The Ecological Detective: Confronting Models with Data , 1999 .

[12]  Perry de Valpine,et al.  BETTER INFERENCES FROM POPULATION-DYNAMICS EXPERIMENTS USING MONTE CARLO STATE-SPACE LIKELIHOOD METHODS , 2003 .

[13]  D. DeAngelis,et al.  New Computer Models Unify Ecological TheoryComputer simulations show that many ecological patterns can be explained by interactions among individual organisms , 1988 .

[14]  David Welch,et al.  Approximate Bayesian computation scheme for parameter inference and model selection in dynamical systems , 2009, Journal of The Royal Society Interface.

[15]  S. T. Buckland,et al.  Hidden process models for animal population dynamics. , 2006, Ecological applications : a publication of the Ecological Society of America.

[16]  James S. Clark,et al.  Hierarchical Modelling for the Environmental Sciences: Statistical Methods and Applications , 2006 .

[17]  C. Briggs,et al.  INVESTIGATING THE POPULATION‐LEVEL EFFECTS OF CHYTRIDIOMYCOSIS: AN EMERGING INFECTIOUS DISEASE OF AMPHIBIANS , 2005 .

[18]  Cheryl J. Briggs,et al.  Enzootic and epizootic dynamics of the chytrid fungal pathogen of amphibians , 2010, Proceedings of the National Academy of Sciences.

[19]  S. Gueron,et al.  The Dynamics of Herds: From Individuals to Aggregations , 1996 .

[20]  R. Speare,et al.  Impact and Dynamics of Disease in Species Threatened by the Amphibian Chytrid Fungus, Batrachochytrium dendrobatidis , 2009, Conservation biology : the journal of the Society for Conservation Biology.

[21]  R. Alford,et al.  The Novel and Endemic Pathogen Hypotheses: Competing Explanations for the Origin of Emerging Infectious Diseases of Wildlife , 2005 .

[22]  Cheryl J Briggs,et al.  Dynamics of an emerging disease drive large-scale amphibian population extinctions , 2010, Proceedings of the National Academy of Sciences.

[23]  V. Grimm.  Ten years of individual-based modelling in ecology: what have we learned and what could we learn in the future? , 1999 .

[24]  Bruce E. Kendall,et al.  ANALYSIS OF SIZE TRAJECTORY DATA USING AN ENERGETIC-BASED GROWTH MODEL , 2005 .

[25]  R. Speare,et al.  Endemic Infection of the Amphibian Chytrid Fungus in a Frog Community Post-Decline , 2004, PLoS biology.

[26]  J. Boyles,et al.  Could localized warm areas inside cold caves reduce mortality of hibernating bats affected by white-nose syndrome? , 2010 .

[27]  C. Briggs,et al.  Quantifying the disease transmission function: effects of density on Batrachochytrium dendrobatidis transmission in the mountain yellow-legged frog Rana muscosa. , 2007, The Journal of animal ecology.

[28]  G. Booth,et al.  BacSim, a simulator for individual-based modelling of bacterial colony growth. , 1998, Microbiology.

[29]  Stephen T. Buckland,et al.  State-space models for the dynamics of wild animal populations , 2004 .

[30]  E L Ionides,et al.  Inference for nonlinear dynamical systems , 2006, Proceedings of the National Academy of Sciences.

[31]  Uta Berger,et al.  Pattern-Oriented Modeling of Agent-Based Complex Systems: Lessons from Ecology , 2005, Science.

[32]  S. Ron.  Predicting the Distribution of the Amphibian Pathogen Batrachochytrium dendrobatidis in the New World 1 , 2005 .

[33]  O. Ovaskainen,et al.  State-space models of individual animal movement. , 2008, Trends in ecology & evolution.

[34]  J W Wimpenny,et al.  Individual-based modelling of biofilms. , 2001, Microbiology.

[35]  Len Thomas,et al.  Monte carlo inference for state-space models of wild animal populations. , 2009, Biometrics.

[36]  D E Green,et al.  Chytridiomycosis causes amphibian mortality associated with population declines in the rain forests of Australia and Central America. , 1998, Proceedings of the National Academy of Sciences of the United States of America.

[37]  Larry B. Crowder,et al.  Variability in survival of larval fish: disentangling components with a generalized individual-based model , 1996 .

[38]  R. Alford,et al.  Emerging infectious disease and the loss of biodiversity in a Neotropical amphibian community. , 2006, Proceedings of the National Academy of Sciences of the United States of America.

[39]  Alan F. Karr,et al.  State-Space and Unobserved Component Models: Theory and Applications , 2006 .

[40]  D. Kirschner,et al.  A methodology for performing global uncertainty and sensitivity analysis in systems biology. , 2008, Journal of theoretical biology.

[41]  Stephen T. Buckland,et al.  Embedding Population Dynamics Models in Inference , 2007, 0708.3796.

[42]  M. Fisher,et al.  Persistence of the emerging pathogen Batrachochytrium dendrobatidis outside the amphibian host greatly increases the probability of host extinction , 2008, Proceedings of the Royal Society B: Biological Sciences.

[43]  A. Ōkubo.  Dynamical aspects of animal grouping: swarms, schools, flocks, and herds. , 1986, Advances in biophysics.

[44]  Richard Speare,et al.  Spread of Chytridiomycosis Has Caused the Rapid Global Decline and Extinction of Frogs , 2007, EcoHealth.

[45]  S. West,et al.  The relationship between parasitoid size and fitness in the field, a study of Achrysocharoides zwoelferi (Hymenoptera: Eulophidae) , 1996 .

[46]  G. McCracken,et al.  Management of the Panzootic White‐Nose Syndrome through Culling of Bats , 2011, Conservation biology : the journal of the Society for Conservation Biology.

[47]  R. Speare,et al.  Life cycle stages of the amphibian chytrid Batrachochytrium dendrobatidis. , 2005, Diseases of aquatic organisms.