An accessible method for implementing hierarchical models with spatio-temporal abundance data

By: , and 



A common goal in ecology and wildlife management is to determine the causes of variation in population dynamics over long periods of time and across large spatial scales. Many assumptions must nevertheless be overcome to make appropriate inference about spatio-temporal variation in population dynamics, such as autocorrelation among data points, excess zeros, and observation error in count data. To address these issues, many scientists and statisticians have recommended the use of Bayesian hierarchical models. Unfortunately, hierarchical statistical models remain somewhat difficult to use because of the necessary quantitative background needed to implement them, or because of the computational demands of using Markov Chain Monte Carlo algorithms to estimate parameters. Fortunately, new tools have recently been developed that make it more feasible for wildlife biologists to fit sophisticated hierarchical Bayesian models (i.e., Integrated Nested Laplace Approximation, ‘INLA’). We present a case study using two important game species in North America, the lesser and greater scaup, to demonstrate how INLA can be used to estimate the parameters in a hierarchical model that decouples observation error from process variation, and accounts for unknown sources of excess zeros as well as spatial and temporal dependence in the data. Ultimately, our goal was to make unbiased inference about spatial variation in population trends over time.
Publication type Article
Publication Subtype Journal Article
Title An accessible method for implementing hierarchical models with spatio-temporal abundance data
Series title PLoS ONE
DOI 10.1371/journal.pone.0049395
Volume 7
Issue 11
Year Published 2012
Language English
Contributing office(s) Colorado Cooperative Fish and Wildlife Research Unit
Description 8 p.
Larger Work Title PLoS ONE
First page 1
Last page 8
Google Analytic Metrics Metrics page
Additional publication details