A multispecies occupancy model for two or more interacting species

Summary Species occurrence is influenced by environmental conditions and the presence of other species. Current approaches for multispecies occupancy modelling are practically limited to two interacting species and often require the assumption of asymmetric interactions. We propose a multispecies occupancy model that can accommodate two or more interacting species. We generalize the single-species occupancy model to two or more interacting species by assuming the latent occupancy state is a multivariate Bernoulli random variable. We propose modelling the probability of each potential latent occupancy state with both a multinomial logit and a multinomial probit model and present details of a Gibbs sampler for the latter. As an example, we model co-occurrence probabilities of bobcat (Lynx rufus), coyote (Canis latrans), grey fox (Urocyon cinereoargenteus) and red fox (Vulpes vulpes) as a function of human disturbance variables throughout 6 Mid-Atlantic states in the eastern United States. We found evidence for pairwise interactions among most species, and the probability of some pairs of species occupying the same site varied along environmental gradients; for example, occupancy probabilities of coyote and grey fox were independent at sites with little human disturbance, but these two species were more likely to occur together at sites with high human disturbance. Ecological communities are composed of multiple interacting species. Our proposed method improves our ability to draw inference from such communities by permitting modelling of detection/non-detection data from an arbitrary number of species, without assuming asymmetric interactions. Additionally, our proposed method permits modelling the probability two or more species occur together as a function of environmental variables. These advancements represent an important improvement in our ability to draw community-level inference from multiple interacting species that are subject to imperfect detection.

[1]  José J. Lahoz-Monfort,et al.  Ignoring Imperfect Detection in Biological Surveys Is Dangerous: A Response to ‘Fitting and Interpreting Occupancy Models' , 2014, PloS one.

[2]  Sheng Li,et al.  Shared resources between giant panda and sympatric wild and domestic mammals , 2015 .

[3]  D. A. Steen,et al.  Snake co-occurrence patterns are best explained by habitat and hypothesized effects of interspecific interactions. , 2014, The Journal of animal ecology.

[4]  Otso Ovaskainen,et al.  Modeling species co-occurrence by multivariate logistic regression generates new hypotheses on fungal interactions. , 2010, Ecology.

[5]  Mevin B. Hooten,et al.  Spatial occupancy models for large data sets , 2013 .

[6]  K. Veblen.  Savanna glade hotspots: Plant community development and synergy with large herbivores , 2012 .

[7]  J. Andrew Royle,et al.  ESTIMATING SITE OCCUPANCY RATES WHEN DETECTION PROBABILITIES ARE LESS THAN ONE , 2002, Ecology.

[8]  Larissa L. Bailey,et al.  Modeling co-occurrence of northern spotted and barred owls: Accounting for detection probability differences , 2009 .

[9]  J. Nichols,et al.  Investigating species co-occurrence patterns when species are detected imperfectly , 2004 .

[10]  G. Wahba,et al.  Multivariate Bernoulli distribution , 2012, 1206.1874.

[11]  Daniel Fortin,et al.  Considering ecological dynamics in resource selection functions. , 2010, The Journal of animal ecology.

[12]  Richard T. Holmes,et al.  Habitat selection by breeding american redstarts in response to a dominant competitor, the least flycatcher , 1988 .

[13]  Christopher C. Wilmers,et al.  Mesopredator spatial and temporal responses to large predators and human development in the Santa Cruz Mountains of California , 2015 .

[14]  Zhihai He,et al.  Volunteer-run cameras as distributed sensors for macrosystem mammal research , 2015, Landscape Ecology.

[15]  A M Latimer,et al.  Hierarchical models facilitate spatial analysis of large data sets: a case study on invasive plant species in the northeastern United States. , 2009, Ecology letters.

[16]  Wilfried Thuiller,et al.  The influence of interspecific interactions on species range expansion rates , 2014, Ecography.

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

[18]  Tavis D. Forrester,et al.  Cats are Rare Where Coyotes Roam , 2015 .

[19]  W. Godsoe,et al.  How do species interactions affect species distribution models , 2012 .

[20]  Christopher A. Barnes,et al.  Completion of the 2006 National Land Cover Database for the conterminous United States. , 2011 .

[21]  T. Virtanen,et al.  Competitive exclusion within the predator community influences the distribution of a threatened prey species. , 2012, Ecology.

[22]  Robert M Dorazio,et al.  A new parameterization for estimating co-occurrence of interacting species. , 2010, Ecological applications : a publication of the Ecological Society of America.

[23]  Robert M. Dorazio,et al.  A Gibbs sampler for Bayesian analysis of site‐occupancy data , 2012 .

[24]  J. Hines,et al.  Two-species occupancy models: a new parameterization applied to co-occurrence of secretive rails. , 2010, Ecological applications : a publication of the Ecological Society of America.

[25]  D. Lindenmayer,et al.  Fitting and Interpreting Occupancy Models , 2013, PloS one.

[26]  Tavis D. Forrester,et al.  A two-species occupancy model accommodating simultaneous spatial and interspecific dependence. , 2016, Ecology.

[27]  Damon B. Lesmeister,et al.  Spatial and Temporal Structure of a Mesocarnivore Guild in Midwestern North America , 2015 .

[28]  S. Chib,et al.  Bayesian analysis of binary and polychotomous response data , 1993 .

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

[30]  R. A. Hill,et al.  Living in a landscape of fear: the impact of predation, resource availability and habitat structure on primate range use , 2014, Animal Behaviour.

[31]  Erik P Willems,et al.  Predator-specific landscapes of fear and resource distribution: effects on spatial range use. , 2009, Ecology.

[32]  R. O’Hara,et al.  A review of Bayesian variable selection methods: what, how and which , 2009 .