MC(MC)MC: exploring Monte Carlo integration within MCMC for mark–recapture models with individual covariates

Summary Estimating abundance from mark–recapture data is challenging when capture probabilities vary among individuals. Initial solutions to this problem were based on fitting conditional likelihoods and estimating abundance as a derived parameter. More recently, Bayesian methods using full likelihoods have been implemented via reversible jump Markov chain Monte Carlo sampling (RJMCMC) or data augmentation (DA). The latter approach is easily implemented in available software and has been applied to fit models that allow for heterogeneity in both open and closed populations. However, both RJMCMC and DA may be inefficient when modelling large populations. We describe an alternative approach using Monte Carlo (MC) integration to approximate the posterior density within a Markov chain Monte Carlo (MCMC) sampling scheme. We show how this Monte Carlo within MCMC (MCWM) approach may be used to fit a simple, closed population model including a single individual covariate and present results from a simulation study comparing RJMCMC, DA and MCWM. We found that MCWM can provide accurate inference about population size and can be more efficient than both RJMCMC and DA. The efficiency of MCWM can also be improved by using advanced MC methods like antithetic sampling. Finally, we apply MCWM to estimate the abundance of meadow voles (Microtus pennsylvanicus) at the Patuxent Wildlife Research Center in 1982 allowing for capture probabilities to vary as a function body mass.

[1]  Steven C. Amstrup,et al.  Estimation of population size using open capture-recapture models , 2001 .

[2]  R King,et al.  Factors Influencing Soay Sheep Survival: A Bayesian Analysis , 2006, Biometrics.

[3]  Kenneth H. Pollock,et al.  Estimating Transition Probabilities for Stage‐Based Population Projection Matrices Using Capture‐Recapture Data , 1992 .

[4]  J. Andrew Royle,et al.  Parameter-expanded data augmentation for Bayesian analysis of capture–recapture models , 2010, Journal of Ornithology.

[5]  O. Gimenez,et al.  Towards built-in capture–recapture mixed models in program E-SURGE , 2012, Journal of Ornithology.

[6]  K. Pollock.  The use of auxiliary variables in capture-recapture modelling: An overview , 2002 .

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

[8]  J. Andrew Royle,et al.  MODELING ABUNDANCE EFFECTS IN DISTANCE SAMPLING , 2004 .

[9]  J. Alho.  Logistic regression in capture-recapture models. , 1990, Biometrics.

[10]  M. Plummer,et al.  CODA: convergence diagnosis and output analysis for MCMC , 2006 .

[11]  R. Huggins.  On the statistical analysis of capture experiments , 1989 .

[12]  D. Rubin,et al.  Maximum likelihood from incomplete data via the EM - algorithm plus discussions on the paper , 1977 .

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

[14]  C. Crainiceanu,et al.  Semiparametric Regression in Capture–Recapture Modeling , 2006, Biometrics.

[15]  J. Dupuis.  Bayesian estimation of movement and survival probabilities from capture-recapture data , 1995 .

[16]  C. Andrieu,et al.  The pseudo-marginal approach for efficient Monte Carlo computations , 2009, 0903.5480.

[17]  M. Beaumont.  Estimation of population growth or decline in genetically monitored populations. , 2003, Genetics.

[18]  J. Andrew Royle,et al.  Analysis of Multinomial Models With Unknown Index Using Data Augmentation , 2007, Journal of Computational and Graphical Statistics.

[19]  C J Schwarz,et al.  An Extension of the Cormack–Jolly–Seber Model for Continuous Covariates with Application to Microtus pennsylvanicus , 2006, Biometrics.

[20]  William A Link,et al.  A cautionary note on the discrete uniform prior for the binomial N. , 2013, Ecology.

[21]  Richard J. Barker,et al.  Full Open Population Capture–Recapture Models With Individual Covariates , 2010, 1009.3516.

[22]  O Gimenez,et al.  Individual heterogeneity in studies on marked animals using numerical integration: capture-recapture mixed models. , 2010, Ecology.

[23]  Geof H. Givens,et al.  Computational Statistics: Givens/Computational Statistics , 2012 .

[24]  S. Brooks,et al.  On the Bayesian Estimation of a Closed Population Size in the Presence of Heterogeneity and Model Uncertainty , 2007, Biometrics.