Choosing priors in Bayesian ecological models by simulating from the prior predictive distribution

Bayesian data analysis is increasingly used in ecology, but prior specification remains focused on choosing non-informative priors (e.g., flat or vague priors). One barrier to choosing more informative priors is that priors must be specified on model parameters (e.g., intercepts, slopes, sigmas), but prior knowledge often exists on the level of the response variable. This is particularly true for common models in ecology, like generalized linear mixed models, which may have a link function and dozens of parameters, each of which needs a prior distribution. We suggest that this difficulty can be overcome by simulating from the prior predictive distribution and visualizing the results on the scale of the response variable. In doing so, some common choices for non-informative priors on parameters can easily be seen to produce biologically impossible values of response variables. Such implications of prior choices are difficult to foresee without visualization. We demonstrate a workflow for prior selection using simulation and visualization with two ecological examples (predator-prey body sizes and spider responses to food competition). This approach is not new, but its adoption by ecologists will help to better incorporate prior information in ecological models, thereby maximizing one of the benefits of Bayesian data analysis.

[1]  Katharine M. Banner,et al.  The use of Bayesian priors in Ecology: The good, the bad and the not great , 2020, Methods in Ecology and Evolution.

[2]  M. Betancourt,et al.  Toward a principled Bayesian workflow in cognitive science. , 2019, Psychological methods.

[3]  Katharine M. Banner,et al.  Evidence of region‐wide bat population decline from long‐term monitoring and Bayesian occupancy models with empirically informed priors , 2019, Ecology and evolution.

[4]  Nathan P. Lemoine,et al.  Moving beyond noninformative priors: why and how to choose weakly informative priors in Bayesian analyses , 2019, Oikos.

[5]  Andrew Gelman,et al.  The experiment is just as important as the likelihood in understanding the prior : A cautionary note on robust cognitive modelling , 2019 .

[6]  Aki Vehtari,et al.  Visualization in Bayesian workflow , 2017, Journal of the Royal Statistical Society: Series A (Statistics in Society).

[7]  J. Wesner,et al.  Predator foraging strategy mediates the effects of predators on local and emigrating prey , 2018 .

[8]  Paul-Christian Bürkner,et al.  brms: An R Package for Bayesian Multilevel Models Using Stan , 2017 .

[9]  Andrew Gelman,et al.  The Prior Can Often Only Be Understood in the Context of the Likelihood , 2017, Entropy.

[10]  Mark Novak,et al.  Bayesian characterization of uncertainty in species interaction strengths , 2017, Oecologia.

[11]  Christian Robert,et al.  Statistical Rethinking , 2017 .

[12]  Justin C. Touchon,et al.  The mismatch between current statistical practice and doctoral training in ecology , 2016 .

[13]  Willem E. Frankenhuis,et al.  Bayesian Models of , 2016 .

[14]  Jerry W. Warmbold Effects of fish on aquatic and terrestrial ecosystems , 2016 .

[15]  Richard McElreath,et al.  Statistical Rethinking: A Bayesian Course with Examples in R and Stan , 2015 .

[16]  A. Gelman,et al.  Stan , 2015 .

[17]  A. Gelman,et al.  Beyond subjective and objective in statistics , 2015, 1508.05453.

[18]  Mevin B. Hooten,et al.  Bayesian Models: A Statistical Primer for Ecologists , 2015 .

[19]  Mevin B. Hooten,et al.  A guide to Bayesian model selection for ecologists , 2015 .

[20]  Peter A Vesk,et al.  The neglected tool in the Bayesian ecologist's shed: a case study testing informative priors' effect on model accuracy , 2014, Ecology and evolution.

[21]  F. Korner‐Nievergelt,et al.  Normal Linear Models , 2015 .

[22]  Fränzi Korner-Nievergelt,et al.  Bayesian Data Analysis in Ecology Using Linear Models with R, BUGS, and Stan : Including Comparisons to Frequentist Statistics , 2015 .

[23]  R Core Team,et al.  R: A language and environment for statistical computing. , 2014 .

[24]  Nicholas K Dulvy,et al.  Ecosystem ecology: size-based constraints on the pyramids of life. , 2013, Trends in ecology & evolution.

[25]  David B. Dunson,et al.  Bayesian data analysis, third edition , 2013 .

[26]  James D. Stamey,et al.  Hidden Dangers of Specifying Noninformative Priors , 2012 .

[27]  Brian Dennis,et al.  Paths to statistical fluency for ecologists , 2010 .

[28]  Kerrie L. Mengersen,et al.  Elicitator: An expert elicitation tool for regression in ecology , 2010, Environ. Model. Softw..

[29]  Andrew Gelman,et al.  Why We (Usually) Don't Have to Worry About Multiple Comparisons , 2009, 0907.2478.

[30]  Werner Ulrich,et al.  Consumer-resource body-size relationships in natural food webs. , 2006, Ecology.

[31]  M. McCarthy,et al.  Profiting from prior information in Bayesian analyses of ecological data , 2005 .

[32]  Werner Ulrich,et al.  BODY SIZES OF CONSUMERS AND THEIR RESOURCES , 2005 .

[33]  Aaron M. Ellison,et al.  Bayesian inference in ecology , 2004 .

[34]  A. Gelman Prior distributions for variance parameters in hierarchical models (comment on article by Browne and Draper) , 2004 .

[35]  R. Christensen,et al.  A New Perspective on Priors for Generalized Linear Models , 1996 .

[36]  Mervyn Thomas,et al.  A Novel Bayesian Approach to Assessing Impacts of Rain Forest Logging , 1996 .

[37]  Wayne S. Smith,et al.  Interactive Elicitation of Opinion for a Normal Linear Model , 1980 .

[38]  Ward Edwards,et al.  Bayesian statistical inference for psychological research. , 1963 .

[39]  D. Lindley The Use of Prior Probability Distributions in Statistical Inference and Decisions , 1961 .