Bayesian inference on age-specific survival for censored and truncated data.

1. Traditional estimation of age-specific survival and mortality rates in vertebrates is limited to individuals with known age. Although this subject has been studied extensively using effective capture-recapture and capture-recovery models, inference remains challenging because of large numbers of incomplete records (i.e. unknown age of many individuals) and because of the inadequate duration of the studies. 2. Here, we present a hierarchical model for capture-recapture/recovery (CRR) data sets with large proportions of unknown times of birth and death. The model uses a Bayesian framework to draw inference on population-level age-specific demographic rates using parametric survival functions and applies this information to reconstruct times of birth and death for individuals with unknown age. 3. We simulated a set of CRR data sets with varying study span and proportions of individuals with known age, and varying recapture and recovery probabilities. We used these data sets to compare our method to a traditional CRR model, which requires knowledge of individual ages. Subsequently, we applied our method to a subset of a long-term CRR data set on Soay sheep. 4. Our results show that this method performs better than the common CRR model when sample sizes are low. Still, our model is sensitive to the choice of priors with low recapture probability and short studies. In such cases, priors that overestimate survival perform better than those that underestimate it. Also, the model was able to estimate accurately ages at death for Soay sheep, with an average error of 0.94 years and to identify differences in mortality rate between sexes. 5. Although many of the problems in the estimation of age-specific survival can be reduced through more efficient sampling schemes, most ecological data sets are still sparse and with a large proportion of missing records. Thus, improved sampling needs still to be combined with statistical models capable of overcoming the unavoidable limitations of any fieldwork. We show that our approach provides reliable estimates of parameters and unknown times of birth and death even with the most incomplete data sets while being flexible enough to accommodate multiple recapture probabilities and covariates.

[1]  R. Bonduriansky,et al.  Sex effects on life span and senescence in the wild when dates of birth and death are unknown. , 2009, Ecology.

[2]  A. Brix.  Bayesian Data Analysis, 2nd edn , 2005 .

[3]  Hal Caswell,et al.  ESTIMATING POPULATION PROJECTION MATRICES FROM MULTI-STAGE MARK-RECAPTURE DATA , 2002 .

[4]  R. Millar,et al.  Comparison of Hierarchical Bayesian Models for Overdispersed Count Data using DIC and Bayes' Factors , 2009, Biometrics.

[5]  L. Kruuk,et al.  Reproductive Senescence in a Long‐Lived Seabird: Rates of Decline in Late‐Life Performance Are Associated with Varying Costs of Early Reproduction , 2008, The American Naturalist.

[6]  Estimating true age-dependence in survival when only adults can be observed: an example with Black-legged Kittiwakes , 2004 .

[7]  Kenneth H. Pollock,et al.  Capture-Recapture Models Allowing for Age-Dependent Survival and Capture Rates , 1981 .

[8]  D. Cox,et al.  Analysis of Survival Data. , 1985 .

[9]  Byron J. T. Morgan,et al.  Factors influencing Soay sheep survival , 2000 .

[10]  Richard J. Barker,et al.  A unified capture-recapture framework , 2008 .

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

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

[13]  T. Clutton‐Brock,et al.  Soay sheep: dynamics and selection in an island population. , 2003 .

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

[15]  Byron J. T. Morgan,et al.  Integrated recovery/recapture data analysis , 1998 .

[16]  Bradley P. Carlin,et al.  Bayesian measures of model complexity and fit , 2002 .

[17]  R. Ricklefs,et al.  Comparison of aging-related mortality among birds and mammals , 2001, Experimental Gerontology.

[18]  Prior distributions for stratified capture-recapture models , 2002 .

[19]  S. Wanless,et al.  Increased adult mortality and reduced breeding success with age in a population of common guillemot Uria aalge using marked birds of unknown age , 2006 .

[20]  Carl J. Schwarz,et al.  Estimating migration rates using tag-recovery data , 1993 .

[21]  J. D. Lebreton,et al.  Multistate recapture models: Modelling incomplete individual histories , 2002 .

[22]  A. N. Arnason,et al.  The estimation of population size, migration rates and survival in a stratified population , 1973, Researches on Population Ecology.

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

[24]  Benjamin Gompertz,et al.  On the Nature of the Function Expressive of the Law of Human Mortality , 1815 .

[25]  B T Grenfell,et al.  Age, sex, density, winter weather, and population crashes in Soay sheep. , 2001, Science.

[26]  G. Jolly.  EXPLICIT ESTIMATES FROM CAPTURE-RECAPTURE DATA WITH BOTH DEATH AND IMMIGRATION-STOCHASTIC MODEL. , 1965, Biometrika.

[27]  P Besbeas,et al.  Integrating Mark–Recapture–Recovery and Census Data to Estimate Animal Abundance and Demographic Parameters , 2002, Biometrics.

[28]  William A Link,et al.  Modeling Association among Demographic Parameters in Analysis of Open Population Capture–Recapture Data , 2005, Biometrics.

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

[30]  R. Cormack.  Estimates of survival from the sighting of marked animals , 1964 .

[31]  H. Müller,et al.  Demographic window to aging in the wild: constructing life tables and estimating survival functions from marked individuals of unknown age , 2004, Aging cell.

[32]  Richard A. Parker,et al.  WinBUGS for population ecologists: bayesian modeling using markov chain Monte Carlo methods , 2009 .

[33]  T. Clutton‐Brock,et al.  Predictors of reproductive cost in female Soay sheep , 2005 .

[34]  Benjamin Gompertz,et al.  XXIV. On the nature of the function expressive of the law of human mortality, and on a new mode of determining the value of life contingencies. In a letter to Francis Baily, Esq. F. R. S. &c , 1825, Philosophical Transactions of the Royal Society of London.

[35]  Adrian F. M. Smith,et al.  Sampling-Based Approaches to Calculating Marginal Densities , 1990 .

[36]  Marco Festa-Bianchet,et al.  AGE‐SPECIFIC SURVIVAL IN FIVE POPULATIONS OF UNGULATES: EVIDENCE OF SENESCENCE , 1999 .

[37]  James S. Clark,et al.  Models for Ecological Data: An Introduction , 2007 .

[38]  Alan E. Gelfand,et al.  Model choice: A minimum posterior predictive loss approach , 1998, AISTATS.

[39]  I. Nisbet.  Detecting and measuring senescence in wild birds: experience with long-lived seabirds , 2001, Experimental Gerontology.

[40]  David R. Anderson,et al.  Kullback-Leibler information as a basis for strong inference in ecological studies , 2001 .

[41]  Shripad Tuljapurkar,et al.  The Evolutionary Demography of Ecological Change: Linking Trait Variation and Population Growth , 2007, Science.

[42]  C. Robert,et al.  Deviance information criteria for missing data models , 2006 .

[43]  G. Seber.  A NOTE ON THE MULTIPLE-RECAPTURE CENSUS. , 1965, Biometrika.

[44]  M. Crawley,et al.  Predictors of early survival in Soay sheep: cohort-, maternal- and individual-level variation , 2005, Proceedings of the Royal Society B: Biological Sciences.

[45]  N. Aebischer,et al.  Survival of the kittiwake in relation to sex, year, breeding experience and position in the colony , 1990 .

[46]  C. Jessica E. Metcalf,et al.  Using Bayesian inference to understand the allocation of resources between sexual and asexual reproduction , 2009 .

[47]  L. Kruuk,et al.  Estimating the functional form for the density dependence from life history data. , 2008, Ecology.

[48]  Ruth King,et al.  Bayesian model discrimination for multiple strata capture‐recapture data , 2002 .

[49]  W. Siler.  A Competing‐Risk Model for Animal Mortality , 1979 .