Open Access
December 2017 Joint Species Distribution Modeling: Dimension Reduction Using Dirichlet Processes
Daniel Taylor-Rodríguez, Kimberly Kaufeld, Erin M. Schliep, James S. Clark, Alan E. Gelfand
Bayesian Anal. 12(4): 939-967 (December 2017). DOI: 10.1214/16-BA1031


Species distribution models are used to evaluate the variables that affect the distribution and abundance of species and to predict biodiversity. Historically, such models have been fitted to each species independently. While independent models can provide useful information regarding distribution and abundance, they ignore the fact that, after accounting for environmental covariates, residual interspecies dependence persists. With stacking of individual models, misleading behaviors, may arise. In particular, individual models often imply too many species per location.

Recently developed joint species distribution models have application to presence–absence, continuous or discrete abundance, abundance with large numbers of zeros, and discrete, ordinal, and compositional data. Here, we deal with the challenge of joint modeling for a large number of species. To appreciate the challenge in the simplest way, with just presence/absence (binary) response and say, S species, we have an S-way contingency table with 2S cell probabilities. Even if S is as small as 100 this is an enormous table, infeasible to work with without some structure to reduce dimension.

We develop a computationally feasible approach to accommodate a large number of species (say order 103) that allows us to: 1) assess the dependence structure across species; 2) identify clusters of species that have similar dependence patterns; and 3) jointly predict species distributions. To do so, we build hierarchical models capturing dependence between species at the first or “data” stage rather than at a second or “mean” stage. We employ the Dirichlet process for clustering in a novel way to reduce dimension in the joint covariance structure. This last step makes computation tractable.

We use Forest Inventory Analysis (FIA) data in the eastern region of the United States to demonstrate our method. It consists of presence–absence measurements for 112 tree species, observed east of the Mississippi. As a proof of concept for our dimension reduction approach, we also include simulations using continuous and binary data.


Download Citation

Daniel Taylor-Rodríguez. Kimberly Kaufeld. Erin M. Schliep. James S. Clark. Alan E. Gelfand. "Joint Species Distribution Modeling: Dimension Reduction Using Dirichlet Processes." Bayesian Anal. 12 (4) 939 - 967, December 2017.


Published: December 2017
First available in Project Euclid: 2 November 2016

zbMATH: 1384.62317
MathSciNet: MR3724974
Digital Object Identifier: 10.1214/16-BA1031

Keywords: abundance , hierarchical model , latent variables , Markov chain Monte Carlo , presence–absence

Vol.12 • No. 4 • December 2017
Back to Top