A vector of point processes for modeling interactions between and within species using capture‐recapture data

Capture‐recapture (CR) data and corresponding models have been used extensively to estimate the size of wildlife populations when detection probability is less than 1. When the locations of traps or cameras used to capture or detect individuals are known, spatially‐explicit CR models are used to infer the spatial pattern of the individual locations and population density. Individual locations, referred to as activity centers (ACs), are defined as the locations around which the individuals move. These ACs are typically assumed to be independent, and their spatial pattern is modeled using homogeneous Poisson processes. However, this assumption is often unrealistic, since individuals can interact with each other, either within a species or between different species. In this article, we consider a vector of point processes from the general class of interaction point processes and develop a model for CR data that can account for interactions, in particular repulsions, between and within multiple species. Interaction point processes present a challenge from an inferential perspective because of the intractability of the normalizing constant of the likelihood function, and hence standard Markov chain Monte Carlo procedures to perform Bayesian inference cannot be applied. Therefore, we adopt an inference procedure based on the Monte Carlo Metropolis Hastings algorithm, which scales well when modeling more than one species. Finally, we adopt an inference method for jointly sampling the latent ACs and the population size based on birth and death processes. This approach also allows us to adaptively tune the proposal distribution of new points, which leads to better mixing especially in the case of non‐uniformly distributed traps. We apply the model to a CR data‐set on leopards and tigers collected at the Corbett Tiger Reserve in India. Our findings suggest that between species repulsion is stronger than within species, while tiger population density is higher than leopard population density at the park.

[1]  H. Bar,et al.  A spatial capture–recapture model with attractions between individuals , 2020, Environmetrics.

[2]  Jim Griffin,et al.  A hierarchical dependent Dirichlet process prior for modelling bird migration patterns in the UK , 2020 .

[3]  U. Kumar,et al.  Do conservation strategies that increase tiger populations have consequences for other wild carnivores like leopards? , 2019, Scientific Reports.

[4]  Yanan Fan,et al.  Handbook of Approximate Bayesian Computation , 2018 .

[5]  Sofia Olhede,et al.  Detecting multivariate interactions in spatial point patterns with Gibbs models and variable selection , 2017, 1705.00689.

[6]  David L. Borchers,et al.  Spatial Capture–Recapture Models , 2016 .

[7]  Jorge Mateu,et al.  Analysis of multispecies point patterns by using multivariate log‐Gaussian Cox processes , 2016 .

[8]  Glenna F. Nightingale,et al.  Pairwise Interaction Point Processes for Modelling Bivariate Spatial Point Patterns in the Presence of Interaction Uncertainty , 2015 .

[9]  Brian J. Reich,et al.  A spatial capture‐recapture model for territorial species , 2014, 1405.1976.

[10]  J. Andrew Royle,et al.  Spatial Capture-Recapture , 2013 .

[11]  Faming Liang,et al.  A Monte Carlo Metropolis-Hastings Algorithm for Sampling from Distributions with Intractable Normalizing Constants , 2013, Neural Computation.

[12]  J. Andrew Royle,et al.  Spatially explicit models for inference about density in unmarked or partially marked populations , 2011, 1112.3250.

[13]  P. Wegge,et al.  Do tigers displace leopards? If so, why? , 2010, Ecological Research.

[14]  Christophe Andrieu,et al.  A tutorial on adaptive MCMC , 2008, Stat. Comput..

[15]  J Andrew Royle,et al.  A hierarchical model for spatial capture-recapture data. , 2008, Ecology.

[16]  D L Borchers,et al.  Spatially Explicit Maximum Likelihood Methods for Capture–Recapture Studies , 2008, Biometrics.

[17]  J Andrew Royle,et al.  Analysis of Capture–Recapture Models with Individual Covariates Using Data Augmentation , 2008, Biometrics.

[18]  P. Giordani,et al.  Adaptive Independent Metropolis–Hastings by Fast Estimation of Mixtures of Normals , 2008, 0801.1864.

[19]  J. Møller,et al.  An efficient Markov chain Monte Carlo method for distributions with intractable normalising constants , 2006 .

[20]  M. Efford.  Density estimation in live‐trapping studies , 2004 .

[21]  J. Andrew Royle.  N‐Mixture Models for Estimating Population Size from Spatially Replicated Counts , 2004, Biometrics.

[22]  J. Møller,et al.  Statistical Inference and Simulation for Spatial Point Processes , 2003 .

[23]  Anders Brix,et al.  Space‐time Multi Type Log Gaussian Cox Processes with a View to Modelling Weeds , 2001 .

[24]  W. Kendall,et al.  Perfect simulation using dominating processes on ordered spaces, with application to locally stable point processes , 2000, Advances in Applied Probability.

[25]  J. L. Jensen,et al.  Pseudolikelihood for Exponential Family Models of Spatial Point Processes , 1991 .

[26]  Y. Ogata,et al.  Likelihood Analysis of Spatial Point Patterns , 1984 .

[27]  David R. Anderson,et al.  Statistical Inference from Band Recovery Data: A Handbook , 1978 .

[28]  J. Seidensticker.  On the Ecological Separation between Tigers and Leopards , 1976 .

[29]  D. J. Strauss.  A model for clustering , 1975 .

[30]  H. Akaike.  A new look at the statistical model identification , 1974 .

[31]  H. Robbins.  A Stochastic Approximation Method , 1951 .

[32]  L. Vogt.  Statistics For Spatial Data , 2016 .

[33]  Ute Dreher,et al.  Statistical Mechanics Rigorous Results , 2016 .

[34]  G. Leckie.  Bulletin of the International Statistical Institute , 2013 .

[35]  J. Andrew Royle,et al.  Dealing with incomplete and variable detectability in multi-year, multi-site monitoring of ecological populations , 2012 .

[36]  Kenneth H. Pollock,et al.  Stopover Duration Analysis with Departure Probability Dependent on Unknown Time Since Arrival , 2009 .

[37]  M. Conroy,et al.  Modeling demographic processes in marked populations , 2009 .

[38]  P. Diggle,et al.  Spatiotemporal prediction for log‐Gaussian Cox processes , 2001 .

[39]  C. Geyer,et al.  Simulation Procedures and Likelihood Inference for Spatial Point Processes , 1994 .

[40]  B. Ripley.  Modelling Spatial Patterns , 1977 .

[41]  J. Neyman,et al.  Statistical Approach to Problems of Cosmology , 1958 .