Inference in MCMC step selection models

Biometrics. 2020 Jun;76(2):438-447. doi: 10.1111/biom.13170. Epub 2019 Dec 9.

Abstract

Habitat selection models are used in ecology to link the spatial distribution of animals to environmental covariates and identify preferred habitats. The most widely used models of this type, resource selection functions, aim to capture the steady-state distribution of space use of the animal, but they assume independence between the observed locations of an animal. This is unrealistic when location data display temporal autocorrelation. The alternative approach of step selection functions embed habitat selection in a model of animal movement, to account for the autocorrelation. However, inferences from step selection functions depend on the underlying movement model, and they do not readily predict steady-state space use. We suggest an analogy between parameter updates and target distributions in Markov chain Monte Carlo (MCMC) algorithms, and step selection and steady-state distributions in movement ecology, leading to a step selection model with an explicit steady-state distribution. In this framework, we explain how maximum likelihood estimation can be used for simultaneous inference about movement and habitat selection. We describe the local Gibbs sampler, a novel rejection-free MCMC scheme, use it as the basis of a flexible class of animal movement models, and derive its likelihood function for several important special cases. In a simulation study, we verify that maximum likelihood estimation can recover all model parameters. We illustrate the application of the method with data from a zebra.

Keywords: MCMC step selection; Markov chain Monte Carlo; animal movement; local Gibbs sampler; resource selection function; step selection function; utilization distribution.

Publication types

  • Research Support, Non-U.S. Gov't

MeSH terms

  • Algorithms
  • Animal Migration
  • Animals
  • Biometry
  • Computer Simulation
  • Ecology / statistics & numerical data
  • Ecosystem*
  • Equidae
  • Likelihood Functions
  • Markov Chains
  • Models, Biological
  • Monte Carlo Method
  • Population Dynamics / statistics & numerical data