Individual heterogeneity in studies on marked animals using numerical integration: capture-recapture mixed models.

In conservation and evolutionary ecology, quantifying and accounting for individual heterogeneity in vital rates of open populations is of particular interest. Individual random effects have been used in capture-recapture models, adopting a Bayesian framework with Markov chain Monte Carlo (MCMC) to carry out estimation and inference. As an alternative, we show how numerical integration via the Gauss-Hermite quadrature (GHQ) can be efficiently used to approximate the capture-recapture model likelihood with individual random effects. We compare the performance of the two approaches (MCMC vs. GHQ) and finite mixture models using two examples, including data on European Dippers and Sociable Weavers. Besides relying on standard statistical tools, GHQ was found to be faster than MCMC simulations. Our approach is implemented in program E-SURGE. Overall, capture recapture mixed models (CR2Ms), implemented either via a GHQ approximation or MCMC simulations, have potential important applications in population biology.

[1]  Emmanuelle Cam,et al.  Dispersal and individual quality in a long lived species , 2004 .

[2]  K. Burnham,et al.  Program MARK: survival estimation from populations of marked animals , 1999 .

[3]  F. Vaida,et al.  Conditional Akaike information for mixed-effects models , 2005 .

[4]  D. Stram,et al.  Variance components testing in the longitudinal mixed effects model. , 1994, Biometrics.

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

[6]  D. Bates,et al.  Mixed-Effects Models in S and S-PLUS , 2001 .

[7]  Olivier Gimenez,et al.  State-space modelling of data on marked individuals , 2007 .

[8]  Jean-Michel Gaillard,et al.  Are extra-pair young better than within-pair young? A comparison of survival and dominance in alpine marmot. , 2007, The Journal of animal ecology.

[9]  A. Shapiro.  Towards a unified theory of inequality constrained testing in multivariate analysis , 1988 .

[10]  Thomas Lenormand,et al.  NONPARAMETRIC ESTIMATION OF NATURAL SELECTION ON A QUANTITATIVE TRAIT USING MARK-RECAPTURE DATA , 2006 .

[11]  Emmanuel Lesaffre,et al.  On the effect of the number of quadrature points in a logistic random effects model: an example , 2001 .

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

[13]  N. Breslow,et al.  Approximate inference in generalized linear mixed models , 1993 .

[14]  A. Agresti,et al.  The Use of Mixed Logit Models to Reflect Heterogeneity in Capture‐Recapture Studies , 1999, Biometrics.

[15]  Xihong Lin.  Variance component testing in generalised linear models with random effects , 1997 .

[16]  K. Liang,et al.  Asymptotic Properties of Maximum Likelihood Estimators and Likelihood Ratio Tests under Nonstandard Conditions , 1987 .

[17]  D. Hedeker,et al.  A random-effects ordinal regression model for multilevel analysis. , 1994, Biometrics.

[18]  Mollie E. Brooks,et al.  Generalized linear mixed models: a practical guide for ecology and evolution. , 2009, Trends in ecology & evolution.

[19]  Rita Covas,et al.  Stabilizing selection on body mass in the sociable weaver Philetairus socius , 2002, Proceedings of the Royal Society of London. Series B: Biological Sciences.

[20]  Roger Pradel,et al.  Program E-Surge: A Software Application for Fitting Multievent Models , 2009 .

[21]  M. Conroy,et al.  Analysis and Management of Animal Populations , 2002 .

[22]  J. Norris,et al.  Capture-Recapture Models with Heterogeneity : I . Cormack-Jolly-Seber Model , 2003 .

[23]  F. Vaida,et al.  MODEL SELECTION FOR PENALIZED SPLINE SMOOTHING USING AKAIKE INFORMATION CRITERIA , 2007 .

[24]  Qing Liu,et al.  A note on Gauss—Hermite quadrature , 1994 .

[25]  Olivier Gimenez,et al.  Discussion: Towards a bayesian analysis template? , 2008 .

[26]  David R. Anderson,et al.  Modeling Survival and Testing Biological Hypotheses Using Marked Animals: A Unified Approach with Case Studies , 1992 .

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

[28]  W. Link,et al.  Individual Covariation in Life‐History Traits: Seeing the Trees Despite the Forest , 2002, The American Naturalist.

[29]  Florian Heiss,et al.  Likelihood approximation by numerical integration on sparse grids , 2008 .

[30]  Gary C. White,et al.  EFFECTS OF INDIVIDUAL HETEROGENEITY IN ESTIMATING THE PERSISTENCE OF SMALL POPULATIONS , 1999 .

[31]  J Andrew Royle,et al.  Web-based Supplementary Materials for “ Modeling Individual Effects in the Cormack-Jolly-Seber Model : A State-space Formulation ” , 2010 .

[32]  B J T Morgan,et al.  The Analysis of Ring‐Recovery Data Using Random Effects , 2003, Biometrics.

[33]  Emmanuel Lesaffre,et al.  Latent‐variable models for longitudinal data with bivariate ordinal outcomes , 2007, Statistics in medicine.

[34]  James S. Clark,et al.  HIERARCHICAL BAYES FOR STRUCTURED, VARIABLE POPULATIONS: FROM RECAPTURE DATA TO LIFE‐HISTORY PREDICTION , 2005 .

[35]  M. Evans,et al.  Methods for Approximating Integrals in Statistics with Special Emphasis on Bayesian Integration Problems , 1995 .

[36]  M. Abramowitz,et al.  Handbook of Mathematical Functions With Formulas, Graphs and Mathematical Tables (National Bureau of Standards Applied Mathematics Series No. 55) , 1965 .

[37]  Robert A. Moffitt,et al.  A COMPUTATIONALLY EFFICIENT QUADRATURE PROCEDURE FOR THE ONE-FACTOR MULTINOMIAL PROBIT MODEL , 1982 .