Nonidentifiability of population size from capture-recapture data with heterogeneous detection probabilities

Heterogeneity in detection probabilities has long been recognized as problematic in mark-recapture studies, and numerous models developed to accommodate its effects. Individual heterogeneity is especially problematic, in that reasonable alternative models may predict essentially identical observations from populations of substantially different sizes. Thus even with very large samples, the analyst will not be able to distinguish among reasonable models of heterogeneity, even though these yield quite distinct inferences about population size. The problem is illustrated with models for closed and open populations.     

Mixture models for estimating the size of a closed population when capture rates vary among individuals

We develop a parameterization of the beta-binomial mixture that provides sensible inferences about the size of a closed population when probabilities of capture or detection vary among individuals. Three classes of mixture models (beta-binomial, logistic-normal, and latent-class) are fitted to recaptures of snowshoe hares for estimating abundance and to counts of bird species for estimating species richness. In both sets of data, rates of detection appear to vary more among individuals (animals or species) than among sampling occasions or locations. The estimates of population size and species richness are sensitive to model-specific assumptions about the latent distribution of individual rates of detection. We demonstrate using simulation experiments that conventional diagnostics for assessing model adequacy, such as deviance, cannot be relied on for selecting classes of mixture models that produce valid inferences about population size. Prior knowledge about sources of individual heterogeneity in detection rates, if available, should be used to help select among classes of mixture models that are to be used for inference.     

Estimation of capture probabilities using generalized estimating equations and mixed effects approaches

Modeling individual heterogeneity in capture probabilities has been one of the most challenging tasks in capture–recapture studies. Heterogeneity in capture probabilities can be modeled as a function of individual covariates, but correlation structure among capture occasions should be taking into account. A proposed generalized estimating equations (GEE) and generalized linear mixed modeling (GLMM) approaches can be used to estimate capture probabilities and population size for capture–recapture closed population models. An example is used for an illustrative application and for comparison with currently used methodology. A simulation study is also conducted to show the performance of the estimation procedures. Our simulation results show that the proposed quasi‐likelihood based on GEE approach provides lower SE than partial likelihood based on either generalized linear models (GLM) or GLMM approaches for estimating population size in a closed capture–recapture experiment. Estimator performance is good if a large proportion of individuals are captured. For cases where only a small proportion of individuals are captured, the estimates become unstable, but the GEE approach outperforms the other methods.     

On the Bayesian estimation of a closed population size in the presence of heterogeneity and model uncertainty

Summary .   We consider the estimation of the size of a closed population, often of interest for wild animal populations, using a capture–recapture study. The estimate of the total population size can be very sensitive to the choice of model used to fit to the data. We consider a Bayesian approach, in which we consider all eight plausible models initially described by Otis et al. (1978, Wildlife Monographs 62, 1–135) within a single framework, including models containing an individual heterogeneity component. We show how we are able to obtain a model-averaged estimate of the total population, incorporating both parameter and model uncertainty. To illustrate the methodology we initially perform a simulation study and analyze two datasets where the population size is known, before considering a real example relating to a population of dolphins off northeast Scotland.     

Upper bound estimators of the population size based on ordinal models for capture‐recapture experiments

Capture‐recapture studies have attracted a lot of attention over the past few decades, especially in applied disciplines where a direct estimate for the size of a population of interest is not available. Epidemiology, ecology, public health, and biodiversity are just a few examples. The estimation of the number of unseen units has been a challenge for theoretical statisticians, and considerable progress has been made in providing lower bound estimators for the population size. In fact, it is well known that consistent estimators for this cannot be provided in the very general case. Considering a case where capture‐recapture studies are summarized by a frequency of frequencies distribution, we derive a simple upper bound of the population size based on the cumulative distribution function. We introduce two estimators of this bound, without any specific parametric assumption on the distribution of the observed frequency counts. The behavior of the proposed estimators is investigated using several benchmark datasets and a large‐scale simulation experiment based on the scheme discussed by Pledger.     

The estimation of population size,migration rates and survival in a stratified population

Summary Estimates of survival, migration rates, and population size are developed for a triple catch marking experiment onn (n>-2) areas with migration among all areas and death in all areas occurring, but no recruitment (birth). This repressents the extension to three sampling times of the method ofChapman andJunge (1956) for estimates in a stratified population. The method is further extented to allow for ‘losses on capture’.     

Estimation of the size of an open population from capture-recapture data using weighted martingale methods

A weighted martingale method, akin to a moving average, is proposed to allow the use of modified closed-population methods in the estimation of the size of a smoothly changing open population when there are frequent capture occasions. We concentrate here on modifications to martingale estimating functions for model Mt, but a wide range of closed-population estimators may be modified in this fashion. The method is motivated by and applied to weekly capture-recapture data from the Mai Po bird sanctuary in Hong Kong. Simulations show that the weighted martingale estimator compared well with the Jolly-Seber estimator when the conditions for the latter to be valid are met, and it performed far better when individuals were allowed to leave and reenter the population. Expressions are derived for the asymptotic bias and variance of the estimator in an appendix.     

On the design of closed recapture experiments

We propose a method to plan the number of occasions of recapture experiments for population size estimation. We do so by fixing the smallest number of capture occasions so that the expected length of the profile confidence interval is less than or equal to a fixed threshold. In some cases, we solve the optimization problem in closed form. For more complex models we use numerical optimization. We detail models assuming homogeneous, time‐varying, subject‐specific capture probabilities, behavioral response to capture, and combining behavioral response with subject‐specific effects. The principle we propose can be extended to plan any other model specification. We formally show the validity of the approach by proving distributional convergence. We illustrate with simulations and challenging examples in epidemiology and ecology. We report that in many cases adding as few as two sampling occasions may substantially reduce the length of confidence intervals.