Multistage time-to-event models improve survival inference by partitioning mortality processes of tracked organisms.

Advances in tagging technologies are expanding opportunities to estimate survival of fish and wildlife populations. Yet, capture and handling effects could impact survival outcomes and bias inference about natural mortality processes. We developed a multistage time-to-event model that can partition the survival process into sequential phases that reflect the tagged animal experience, including handling and release mortality, post-release recovery mortality, and subsequently, natural mortality. We demonstrate performance of multistage survival models through simulation testing and through fish and bird telemetry case studies. Models are implemented in a Bayesian framework and can accommodate left, right, and interval censorship events. Our results indicate that accurate survival estimates can be achieved with reasonable sample sizes ( n ≈ 100 + ) and that multimodel inference can inform hypotheses about the configuration and length of survival stages needed to adequately describe mortality processes for tracked specimens. While we focus on survival estimation for tagged fish and wildlife populations, multistage time-to-event models could be used to understand other phenomena of interest such as migration, reproduction, or disease events across a range of taxa including plants and insects.

Practical advice on variable selection and reporting using Akaike information criterion.

The various debates around model selection paradigms are important, but in lieu of a consensus, there is a demonstrable need for a deeper appreciation of existing approaches, at least among the end-users of statistics and model selection tools. In the ecological literature, the Akaike information criterion (AIC) dominates model selection practices, and while it is a relatively straightforward concept, there exists what we perceive to be some common misunderstandings around its application. Two specific questions arise with surprising regularity among colleagues and students when interpreting and reporting AIC model tables. The first is related to the issue of 'pretending' variables, and specifically a muddled understanding of what this means. The second is related to p-values and what constitutes statistical support when using AIC. There exists a wealth of technical literature describing AIC and the relationship between p-values and AIC differences. Here, we complement this technical treatment and use simulation to develop some intuition around these important concepts. In doing so we aim to promote better statistical practices when it comes to using, interpreting and reporting models selected when using AIC.

Including a spatial predictive process in band recovery models improves inference for Lincoln estimates of animal abundance.

Abundance estimation is a critical component of conservation planning, particularly for exploited species where managers set regulations to restrict harvest based on current population size. An increasingly common approach for abundance estimation is through integrated population modeling (IPM), which uses multiple data sources in a joint likelihood to estimate abundance and additional demographic parameters. Lincoln estimators are one commonly used IPM component for harvested species, which combine information on the rate and total number of individuals harvested within an integrated band-recovery framework to estimate abundance at large scales. A major assumption of the Lincoln estimator is that banding and recoveries are representative of the whole population, which may be violated if major sources of spatial heterogeneity in survival or harvest rates are not incorporated into the model. We developed an approach to account for spatial variation in harvest rates using a spatial predictive process, which we incorporated into a Lincoln estimator IPM. We simulated data under different configurations of sample sizes, harvest rates, and sources of spatial heterogeneity in harvest rate to assess potential model bias in parameter estimates. We then applied the model to data collected from a field study of wild turkeys (Meleagris gallapavo) to estimate local and statewide abundance in Maine, USA. We found that the band recovery model that incorporated a spatial predictive process consistently provided estimates of adult and juvenile abundance with low bias across a variety of spatial configurations of harvest rate and sampling intensities. When applied to data collected on wild turkeys, a model that did not incorporate spatial heterogeneity underestimated the harvest rate in some subregions. Consistent with simulation results, this led to overestimation of both local and statewide abundance. Our work demonstrates that a spatial predictive process is a viable mechanism to account for spatial variation in harvest rates and limit bias in abundance estimates. This approach could be extended to large-scale band recovery data sets and has applicability for the estimation of population parameters in other ecological models as well.

Improved inferences about landscape connectivity from spatial capture-recapture by integration of a movement model.

Understanding how broad-scale patterns in animal populations emerge from individual-level processes is an enduring challenge in ecology that requires investigation at multiple scales and perspectives. Complementary to this need for diverse approaches is the recent focus on integrated modeling in statistical ecology. Population-level processes represent the core of spatial capture-recapture (SCR), with many methodological extensions that have been motivated by standing ecological theory and data-integration opportunities. The extent to which these recent advances offer inferential improvements can be limited by the data requirements for quantifying individual-level processes. This is especially true for SCR models that use non-Euclidean distance to relax the restrictive assumption that individual space use is stationary and symmetrical to make inferences about landscape connectivity. To meet the challenges of scale and data quality, we propose integrating an explicit movement model with non-Euclidean SCR for joint estimation of a shared cost parameter between individual and population processes. Here, we define a movement kernel for step selection that uses "ecological distance" instead of Euclidean distance to quantify availability for each movement step in terms of landscape cost. We compare performance of our integrated model to that of existing SCR models using realistic animal movement simulations and data collected on black bears. We demonstrate that an integrated approach offers improvements both in terms of bias and precision in estimating the shared cost parameter over models fit to spatial encounters alone. Simulations suggest these gains were only realized when step lengths were small relative to home range size, and estimates of density were insensitive to whether or not an integrated approach was used. By combining the fine spatiotemporal scale of individual movement processes with the estimation of population density in SCR, integrated approaches such as the one we develop here have the potential to unify the fields of movement, population, and landscape ecology and improve our understanding of landscape connectivity.

Spatial capture-recapture with random thinning for unidentified encounters.

Spatial capture-recapture (SCR) models have increasingly been used as a basis for combining capture-recapture data types with variable levels of individual identity information to estimate population density and other demographic parameters. Recent examples are the unmarked SCR (or spatial count model), where no individual identities are available and spatial mark-resight (SMR) where individual identities are available for only a marked subset of the population. Currently lacking, though, is a model that allows unidentified samples to be combined with identified samples when there are no separate classes of "marked" and "unmarked" individuals and when the two sample types cannot be considered as arising from two independent observation models. This is a common scenario when using noninvasive sampling methods, for example, when analyzing data on identified and unidentified photographs or scats from the same sites.Here we describe a "random thinning" SCR model that utilizes encounters of both known and unknown identity samples using a natural mechanistic dependence between samples arising from a single observation model. Our model was fitted in a Bayesian framework using NIMBLE.We investigate the improvement in parameter estimates by including the unknown identity samples, which was notable (up to 79% more precise) in low-density populations with a low rate of identified encounters. We then applied the random thinning SCR model to a noninvasive genetic sampling study of brown bear (Ursus arctos) density in Oriental Cantabrian Mountains (North Spain).Our model can improve density estimation for noninvasive sampling studies for low-density populations with low rates of individual identification, by making use of available data that might otherwise be discarded.

Spatial proximity moderates genotype uncertainty in genetic tagging studies.

Accelerating declines of an increasing number of animal populations worldwide necessitate methods to reliably and efficiently estimate demographic parameters such as population density and trajectory. Standard methods for estimating demographic parameters from noninvasive genetic samples are inefficient because lower-quality samples cannot be used, and they assume individuals are identified without error. We introduce the genotype spatial partial identity model (gSPIM), which integrates a genetic classification model with a spatial population model to combine both spatial and genetic information, thus reducing genotype uncertainty and increasing the precision of demographic parameter estimates. We apply this model to data from a study of fishers (Pekania pennanti) in which 37% of hair samples were originally discarded because of uncertainty in individual identity. The gSPIM density estimate using all collected samples was 25% more precise than the original density estimate, and the model identified and corrected three errors in the original individual identity assignments. A simulation study demonstrated that our model increased the accuracy and precision of density estimates 63 and 42%, respectively, using three replicated assignments (e.g., PCRs for microsatellites) per genetic sample. Further, the simulations showed that the gSPIM model parameters are identifiable with only one replicated assignment per sample and that accuracy and precision are relatively insensitive to the number of replicated assignments for high-quality samples. Current genotyping protocols devote the majority of resources to replicating and confirming high-quality samples, but when using the gSPIM, genotyping protocols could be more efficient by devoting more resources to low-quality samples.

Do the antipredator strategies of shared prey mediate intraguild predation and mesopredator suppression?

Understanding the conditions that facilitate top predator effects upon mesopredators and prey is critical for predicting where these effects will be significant. Intraguild predation (IGP) and the ecology of fear are hypotheses used to describe the effects of top predators upon mesopredators and prey species, but make different assumptions about organismal space use. The IGP hypothesis predicts that mesopredator resource acquisition and risk are positively correlated, creating a fitness deficit. But if shared prey also avoid a top predator, then mesopredators may not have to choose between risk and reward. Prey life history may be a critical predictor of how shared prey respond to predation and may mediate mesopredator suppression. We used hierarchical models of species distribution and abundance to test expectations of IGP using two separate triangular relationships between a large carnivore, smaller intraguild carnivore, and shared mammalian prey with different life histories. Following IGP, we expected that a larger carnivore would suppress a smaller carnivore if the shared prey species did not spatially avoid the large carnivore at broad scales. If prey were fearful over broad scales, we expected less evidence of mesopredator suppression. We tested these theoretical hypotheses using remote camera detections across a large spatial extent. Lagomorphs did not appear to avoid coyotes, and fox detection probability was lower as coyote abundance increased. In contrast, white-tailed deer appeared to avoid areas of increased wolf use, and coyote detection probability was not reduced at sites where wolves occurred. These findings suggest that mesopredator suppression by larger carnivores may depend upon the behavior of shared prey, specifically the spatial scale at which they perceive risk. We further discuss how extrinsic environmental factors may contribute to mesopredator suppression.

Improving inferences from short-term ecological studies with Bayesian hierarchical modeling: white-headed woodpeckers in managed forests.

Pilot studies are often used to design short-term research projects and long-term ecological monitoring programs, but data are sometimes discarded when they do not match the eventual survey design. Bayesian hierarchical modeling provides a convenient framework for integrating multiple data sources while explicitly separating sample variation into observation and ecological state processes. Such an approach can better estimate state uncertainty and improve inferences from short-term studies in dynamic systems. We used a dynamic multistate occupancy model to estimate the probabilities of occurrence and nesting for white-headed woodpeckers Picoides albolarvatus in recent harvest units within managed forests of northern California, USA. Our objectives were to examine how occupancy states and state transitions were related to forest management practices, and how the probabilities changed over time. Using Gibbs variable selection, we made inferences using multiple model structures and generated model-averaged estimates. Probabilities of white-headed woodpecker occurrence and nesting were high in 2009 and 2010, and the probability that nesting persisted at a site was positively related to the snag density in harvest units. Prior-year nesting resulted in higher probabilities of subsequent occurrence and nesting. We demonstrate the benefit of forest management practices that increase the density of retained snags in harvest units for providing white-headed woodpecker nesting habitat. While including an additional year of data from our pilot study did not drastically alter management recommendations, it changed the interpretation of the mechanism behind the observed dynamics. Bayesian hierarchical modeling has the potential to maximize the utility of studies based on small sample sizes while fully accounting for measurement error and both estimation and model uncertainty, thereby improving the ability of observational data to inform conservation and management strategies.