boral – Bayesian Ordination and Regression Analysis of Multivariate Abundance Data in r

Summary Model-based methods have emerged as a powerful approach for analysing multivariate abundance data in community ecology. Key applications include model-based ordination, modelling the various sources of correlations across species, and making inferences while accounting for these between species correlations. boral (version 0.9.1, licence GPL-2) is an r package available on cran for model-based analysis of multivariate abundance data, with estimation performed using Bayesian Markov chain Monte Carlo methods. A key feature of the boral package is the ability to incorporate latent variables as a parsimonious method of modelling between species correlation. Pure latent variable models offer a model-based approach to unconstrained ordination, for visualizing sites and the indicator species characterizing them on a low-dimensional plot. Correlated response models consist of fitting generalized linear models to each species, while including latent variables to account for residual correlation between species, for example, due to unmeasured covariates.

[1]  David I. Warton,et al.  Finite Mixture of Regression Modeling for High-Dimensional Count and Biomass Data in Ecology , 2013 .

[2]  Martyn Plummer,et al.  JAGS: A program for analysis of Bayesian graphical models using Gibbs sampling , 2003 .

[3]  Steven C. Walker.  Indirect gradient analysis by Markov-chain Monte Carlo , 2015, Plant Ecology.

[4]  Scott D. Foster,et al.  A Poisson–Gamma model for analysis of ecological non-negative continuous data , 2012, Environmental and Ecological Statistics.

[5]  Wim A. Ozinga,et al.  Selecting traits that explain species–environment relationships: a generalized linear mixed model approach , 2013 .

[6]  Aki Vehtari,et al.  Understanding predictive information criteria for Bayesian models , 2013, Statistics and Computing.

[7]  J. Kruskal.  Nonmetric multidimensional scaling: A numerical method , 1964 .

[8]  Sara Taskinen,et al.  Model‐based approaches to unconstrained ordination , 2015 .

[9]  Laura J. Pollock,et al.  Understanding co‐occurrence by modelling species simultaneously with a Joint Species Distribution Model (JSDM) , 2014 .

[10]  M. G. Pittau,et al.  A weakly informative default prior distribution for logistic and other regression models , 2008, 0901.4011.

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

[12]  Nellie Smeenk-Enserink,et al.  Correlations Between Distributions of Hunting Spiders (Lycosidae, Ctenidae) and Environmental Characteristics in a Dune Area , 1974 .

[13]  David I. Warton,et al.  The fourth‐corner solution – using predictive models to understand how species traits interact with the environment , 2014 .

[14]  Francis K. C. Hui,et al.  So Many Variables: Joint Modeling in Community Ecology. , 2015, Trends in ecology & evolution.

[15]  J. Andrew Royle,et al.  Hierarchical Modeling and Inference in Ecology: The Analysis of Data from Populations, Metapopulations and Communities , 2008 .

[16]  Andrew Thomas,et al.  WinBUGS - A Bayesian modelling framework: Concepts, structure, and extensibility , 2000, Stat. Comput..

[17]  David J. Harris.  Generating realistic assemblages with a joint species distribution model , 2015 .

[18]  Yi Wang,et al.  mvabund– an R package for model‐based analysis of multivariate abundance data , 2012 .

[19]  Glenn D. Wylie,et al.  ‘Exciting statistics’: the rapid development and promising future of hierarchical models for population ecology , 2012 .