Dataset Information

Incorporating capture heterogeneity in the estimation of autoregressive coefficients of animal population dynamics using capture-recapture data.

ABSTRACT: Population dynamic models combine density dependence and environmental effects. Ignoring sampling uncertainty might lead to biased estimation of the strength of density dependence. This is typically addressed using state-space model approaches, which integrate sampling error and population process estimates. Such models seldom include an explicit link between the sampling procedures and the true abundance, which is common in capture-recapture settings. However, many of the models proposed to estimate abundance in the presence of capture heterogeneity lead to incomplete likelihood functions and cannot be straightforwardly included in state-space models. We assessed the importance of estimating sampling error explicitly by taking an intermediate approach between ignoring uncertainty in abundance estimates and fully specified state-space models for density-dependence estimation based on autoregressive processes. First, we estimated individual capture probabilities based on a heterogeneity model for a closed population, using a conditional multinomial likelihood, followed by a Horvitz-Thompson estimate for abundance. Second, we estimated coefficients of autoregressive models for the log abundance. Inference was performed using the methodology of integrated nested Laplace approximation (INLA). We performed an extensive simulation study to compare our approach with estimates disregarding capture history information, and using R-package VGAM, for different parameter specifications. The methods were then applied to a real data set of gray-sided voles Myodes rufocanus from Northern Norway. We found that density-dependence estimation was improved when explicitly modeling sampling error in scenarios with low process variances, in which differences in coverage reached up to 8% in estimating the coefficients of the autoregressive processes. In this case, the bias also increased assuming a Poisson distribution in the observational model. For high process variances, the differences between methods were small and it appeared less important to model heterogeneity.

SUBMITTER: Nicolau PG

PROVIDER: S-EPMC7713978 | biostudies-literature | 2020 Dec

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Incorporating capture heterogeneity in the estimation of autoregressive coefficients of animal population dynamics using capture-recapture data.

Nicolau Pedro G PG Sørbye Sigrunn H SH Yoccoz Nigel G NG

Ecology and evolution 20200831 23

Population dynamic models combine density dependence and environmental effects. Ignoring sampling uncertainty might lead to biased estimation of the strength of density dependence. This is typically addressed using state-space model approaches, which integrate sampling error and population process estimates. Such models seldom include an explicit link between the sampling procedures and the true abundance, which is common in capture-recapture settings. However, many of the models proposed to est ...[more]

PMID: 33304489

Similar Datasets

Project description:Background and aimsCapture-recapture (CRC) analysis is recommended for estimating the prevalence of problem drug use or people who inject drugs (PWID). We aim to demonstrate how naive application of CRC can lead to highly misleading results, and to suggest how the problems might be overcome.MethodsWe present a case study of estimating the prevalence of PWID in Bristol, UK, applying CRC to lists in contact with three services. We assess: (i) sensitivity of results to different versions of the dominant (treatment) list: specifically, to inclusion of non-incident cases and of those who were referred directly from one of the other services; (ii) the impact of accounting for a novel covariate, housing instability; and (iii) consistency of CRC estimates with drug-related mortality data. We then incorporate formally the drug-related mortality data and lower bounds for prevalence alongside the CRC into a single coherent model.ResultsFive of 11 models fitted the full data equally well but generated widely varying prevalence estimates, from 2740 [95% confidence interval (CI) = 2670, 2840] to 6890 (95% CI = 3740, 17680). Results were highly sensitive to inclusion of non-incident cases, demonstrating the presence of considerable heterogeneity, and were sensitive to a lesser extent to inclusion of direct referrals. A reduced data set including only incident cases and excluding referrals could be fitted by simpler models, and led to much greater consistency in estimates. Accounting for housing stability improved model fit considerably more than did the standard covariates of age and gender. External data provided validation of results and aided model selection, generating a final estimate of the number of PWID in Bristol in 2011 of 2770 [95% credible interval (Cr-I) = 2570, 3110] or 0.9% (95% Cr-I = 0.9, 1.0%) of the population aged 15-64 years.ConclusionsSteps can be taken to reduce bias in capture-recapture analysis, including: careful consideration of data sources, reduction of lists to less heterogeneous subsamples, use of covariates and formal incorporation of external data.

Project description:Estimating population sizes has long been of interest, from the estimation of the human or ecological population size within regions or countries to the hidden number of civilian casualties in a war. Total enumeration of the population, for example, via a census, is often infeasible or simply impractical. However, a series of partial enumerations or observations of the population is often possible. This has led to the ideas of capture-recapture methods, which have been extensively used within ecology to estimate the size of wildlife populations, with an associated measure of uncertainty, and are most effectively applied when there are multiple capture occasions. Capture-recapture ideology can be more widely applied to multiple data-sources, by the linkage of individuals across the multiple lists. This is often referred to as Multiple Systems Estimation (MSE). The MSE approach has been preferred when estimating "capture-shy" or hard-to-reach populations, including those caught up in the criminal justice system; or homeless; or trafficked; or civilian casualties of war. Motivated by a range of public policy applications of MSE, each briefly introduced, we discuss practical problems with potentially substantial methodological implications. They include: "period" definition; "case" definition; when an observed count is not a true count of the population of interest but an upper bound due to mismatched definitions; exact or probabilistic matching of "cases" across different lists; demographic or other information about the "case" which may influence capture-propensities; required permissions to access extant-lists; list-creation by research-teams or interested parties; referrals (if presence on list A results - almost surely - in presence on list B); different mathematical models leading to widely different estimated population sizes; uncertainty in estimation; computational efficiency; external validation; hypothesis-generation; and additional independent external information. Returning to our motivational applications, we focus on whether the uncertainty which qualified their estimates was sufficiently narrow to orient public policy; and, if not, what options were available and/or taken to reduce the uncertainty or to seek external validation. We also consider whether MSE was hypothesis-generating: in the sense of having spawned new lines of inquiry.

Project description:Standard procedures for capture-mark-recapture modelling (CMR) for the study of animal demography include running goodness-of-fit tests on a general starting model. A frequent reason for poor model fit is heterogeneity in local survival among individuals captured for the first time and those already captured or seen on previous occasions. This deviation is technically termed a transience effect. In specific cases, simple, uni-state CMR modeling showing transients may allow researchers to assess the role of these transients on population dynamics. Transient individuals nearly always have a lower local survival probability, which may appear for a number of reasons. In most cases, transients arise due to permanent dispersal, higher mortality, or a combination of both. In the case of higher mortality, transients may be symptomatic of a cost of first reproduction. A few studies working at large spatial scales actually show that transients more often correspond to survival costs of first reproduction rather than to permanent dispersal, bolstering the interpretation of transience as a measure of costs of reproduction, since initial detections are often associated with first breeding attempts. Regardless of their cause, the loss of transients from a local population should lower population growth rate. We review almost 1000 papers using CMR modeling and find that almost 40% of studies fitting the searching criteria (N = 115) detected transients. Nevertheless, few researchers have considered the ecological or evolutionary meaning of the transient phenomenon. Only three studies from the reviewed papers considered transients to be a cost of first reproduction. We also analyze a long-term individual monitoring dataset (1988-2012) on a long-lived bird to quantify transients, and we use a life table response experiment (LTRE) to measure the consequences of transients at a population level. As expected, population growth rate decreased when the environment became harsher while the proportion of transients increased. LTRE analysis showed that population growth can be substantially affected by changes in traits that are variable under environmental stochasticity and deterministic perturbations, such as recruitment, fecundity of experienced individuals, and transient probabilities. This occurred even though sensitivities and elasticities of these parameters were much lower than those for adult survival. The proportion of transients also increased with the strength of density-dependence. These results have implications for ecological and evolutionary studies and may stimulate other researchers to explore the ecological processes behind the occurrence of transients in capture-recapture studies. In population models, the inclusion of a specific state for transients may help to make more reliable predictions for endangered and harvested species.

Project description:Information about how animal abundance varies across landscapes is needed to inform management action but is costly and time-consuming to obtain; surveys of a single population distributed over a large area can take years to complete. Surveys employing small, spatially replicated sampling units improve efficiency, but statistical estimators rely on assumptions that constrain survey design or become less reasonable as larger areas are sampled. Efficient methods that avoid assumptions about similarity of detectability or density among replicates are therefore appealing. Using simulations and data from >3500 black bears sampled on 73 independent study areas in Ontario, Canada, we (1) quantified bias induced by unmodeled spatial heterogeneity in detectability and density; (2) evaluated novel, design-based estimators of average density across replicate study areas; and (3) evaluated two estimators of the variance of average density across study areas: an analytic estimator that assumed an underlying homogeneous spatial Poisson point process for the distribution of animals' activity centers, and an empirical estimator of variance across study areas. In simulations where detectability varied in space, assuming spatially constant detectability yielded density estimates that were negatively biased by 20% to 30%; estimating local detectability and density from local data and treating study areas as independent, equal replicates when estimating average density across study areas using the design-based estimator yielded unbiased estimates at local and landscape scales. Similarly, detectability of black bears varied among study areas and estimates of bear density at landscape scales were higher when no information was shared across study areas when estimating detectability. This approach also maximized precision (relative SEs of estimates of average black bear density ranged from 7% to 18%) and computational efficiency. In simulations, the analytic variance estimator was robust to threefold variation in local densities but the empirical estimator performed poorly. Conducting multiple, similar SECR surveys and treating them as independent replicates during analyses allowed us to efficiently estimate density at multiple scales and extents while avoiding biases caused by pooling spatially heterogeneous data. This approach enables researchers to address a wide range of ecological or management-related questions and is applicable with most types of SECR data.

Project description:Over the last decade, spatial capture-recapture (SCR) models have become widespread for estimating demographic parameters in ecological studies. However, the underlying assumptions about animal movement and space use are often not realistic. This is a missed opportunity because interesting ecological questions related to animal space use, habitat selection, and behavior cannot be addressed with most SCR models, despite the fact that the data collected in SCR studies - individual animals observed at specific locations and times - can provide a rich source of information about these processes and how they relate to demographic rates. We developed SCR models that integrated more complex movement processes that are typically inferred from telemetry data, including a simple random walk, correlated random walk (i.e., short-term directional persistence), and habitat-driven Langevin diffusion. We demonstrated how to formulate, simulate from, and fit these models with standard SCR data using data-augmented Bayesian analysis methods. We evaluated their performance through a simulation study, in which we varied the detection, movement, and resource selection parameters. We also examined different numbers of sampling occasions and assessed performance gains when including auxiliary location data collected from telemetered individuals. Across all scenarios, the integrated SCR movement models performed well in terms of abundance, detection, and movement parameter estimation. We found little difference in bias for the simple random walk model when reducing the number of sampling occasions from T = 25 to T = 15. We found some bias in movement parameter estimates under several of the correlated random walk scenarios, but incorporating auxiliary location data improved parameter estimates and significantly improved mixing during model fitting. The Langevin movement model was able to recover resource selection parameters from standard SCR data, which is particularly appealing because it explicitly links the individual-level movement process with habitat selection and population density. We focused on closed population models, but the movement models developed here can be extended to open SCR models. The movement process models could also be easily extended to accommodate additional "building blocks" of random walks, such as central tendency (e.g., territoriality) or multiple movement behavior states, thereby providing a flexible and coherent framework for linking animal movement behavior to population dynamics, density, and distribution.

Dataset Information

Incorporating capture heterogeneity in the estimation of autoregressive coefficients of animal population dynamics using capture-recapture data.

Publications

Incorporating capture heterogeneity in the estimation of autoregressive coefficients of animal population dynamics using capture-recapture data.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets