An r package for analysing survival using continuous‐time open capture–recapture models

Summary Capture–recapture software packages have proven to be very powerful tools for analysing factors affecting survival in wild populations. However, all such packages are limited to discrete-time protocols. Appropriate survival analysis tools are still lacking for data acquired from continuous-time protocols. We have developed a statistical method and propose an r package for analysing such data based on an extension of classical survival analysis models incorporating an inhomogeneous Poisson process for modelling capture histories. First, data were simulated from a continuous-time protocol. These data were used to (i) compare survival estimation biases of discrete- and continuous-time approaches and (ii) investigate the performance and accuracy of our r package for four types of covariates: factors varying between individuals (like sex), in time (like climatic factors), both in time and between individuals (like physical condition) and age (as a categorical factor). Secondly, the r package has been applied to a real data set for survival analysis of cats in the Kerguelen archipelago (regrouping 682 cats over 20 years) as an illustrative example. Results of the simulated data analysis show that the method performs better than its discrete-time counterpart for analysing data acquired from continuous-time protocols. It provides unbiased parameter estimates for all parameters except those that vary both in time and between individuals – which is not surprising, since in our case, these factors were not updated in continuous time (i.e. only upon capture). When applied to the Kerguelen cat data set, the results suggest that survival is lower in juveniles than in adults and subadults, varies between study sites and increases with physical condition, and this latter effect being more important in females than in males. Sex, season, temporal linear trend in survival and the NDVI vegetation index were also tested but were not found to be significant. However, confidence intervals were too large (due to a low recapture rate) for excluding such effects. Further analyses are still needed for rigorous covariate testing in this context. In conclusion, continuous-time approaches – such as that presented in this paper – should be preferred when data acquired from continuous-time protocols is analysed.

[1]  D. Pontier,et al.  Linking genetic diversity and temporal fluctuations in population abundance of the introduced feral cat (Felis silvestris catus) on the Kerguelen archipelago , 2011, Molecular ecology.

[2]  D. Lin,et al.  Inference for capture-recapture experiments in continuous time with variable capture rates , 1996 .

[3]  D. Pontier,et al.  Assessing the performance of NDVI as a proxy for plant biomass using non-linear models: a case study on the Kerguelen archipelago , 2009, Polar Biology.

[4]  D. Pontier,et al.  High variation in multiple paternity of domestic cats (Felis catus L.) in relation to environmental conditions , 1999, Proceedings of the Royal Society of London. Series B: Biological Sciences.

[5]  A Chao,et al.  The applications of capture‐recapture models to epidemiological data , 2001, Statistics in medicine.

[6]  D. Cox.  Regression Models and Life-Tables , 1972 .

[7]  C. Spearman.  The proof and measurement of association between two things. By C. Spearman, 1904. , 1987, The American journal of psychology.

[8]  D. Pontier,et al.  The mating system of feral cats (Felis catus L.) in a sub-Antarctic environment , 2002, Polar Biology.

[9]  O. Gimenez,et al.  Simple Estimation and Test Procedures in Capture–Mark–Recapture Mixed Models , 2012, Biometrics.

[10]  Anne Chao,et al.  Continuous-time capture-recapture models with covariates , 2002 .

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

[12]  N. Pettorelli,et al.  Using the satellite-derived NDVI to assess ecological responses to environmental change. , 2005, Trends in ecology & evolution.

[13]  J. Thioulouse,et al.  The diet of feral cats (Felis catus L.) at five sites on the Grande Terre, Kerguelen archipelago , 2002, Polar Biology.

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

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

[16]  D. Pontier,et al.  Molecular and behavioural analyses reveal male-biased dispersal between social groups of domestic cats1 , 2004 .

[17]  L. Rivest,et al.  Rcapture: Loglinear Models for Capture-Recapture in R , 2007 .

[18]  Carl James Schwarz,et al.  A Multilevel Model for Continuous Time Population Estimation , 2009, Biometrics.

[19]  Jeffrey L. Laake,et al.  RMark : an R Interface for analysis of capture-recapture data with MARK , 2013 .

[20]  O. Gimenez,et al.  Accounting for Sampling Error When Inferring Population Synchrony from Time-Series Data: A Bayesian State-Space Modelling Approach with Applications , 2014, PloS one.

[21]  P. Yip,et al.  Estimating Population Size for a Continuous Time Frailty Model with Covariates in a Capture–Recapture Study , 2007, Biometrics.

[22]  Kai Lorenzen,et al.  Apparent survival estimation from continuous mark–recapture/resighting data , 2013 .

[23]  David L. Borchers,et al.  Continuous‐time spatially explicit capture–recapture models, with an application to a jaguar camera‐trap survey , 2014 .

[24]  D. Gillespie.  A General Method for Numerically Simulating the Stochastic Time Evolution of Coupled Chemical Reactions , 1976 .