96 research outputs found
Testing the consistency of wildlife data types before combining them: the case of camera traps and telemetry.
Wildlife data gathered by different monitoring techniques are often combined to estimate animal density. However, methods to check whether different types of data provide consistent information (i.e., can information from one data type be used to predict responses in the other?) before combining them are lacking. We used generalized linear models and generalized linear mixed-effects models to relate camera trap probabilities for marked animals to independent space use from telemetry relocations using 2 years of data for fishers (Pekania pennanti) as a case study. We evaluated (1) camera trap efficacy by estimating how camera detection probabilities are related to nearby telemetry relocations and (2) whether home range utilization density estimated from telemetry data adequately predicts camera detection probabilities, which would indicate consistency of the two data types. The number of telemetry relocations within 250 and 500 m from camera traps predicted detection probability well. For the same number of relocations, females were more likely to be detected during the first year. During the second year, all fishers were more likely to be detected during the fall/winter season. Models predicting camera detection probability and photo counts solely from telemetry utilization density had the best or nearly best Akaike Information Criterion (AIC), suggesting that telemetry and camera traps provide consistent information on space use. Given the same utilization density, males were more likely to be photo-captured due to larger home ranges and higher movement rates. Although methods that combine data types (spatially explicit capture-recapture) make simple assumptions about home range shapes, it is reasonable to conclude that in our case, camera trap data do reflect space use in a manner consistent with telemetry data. However, differences between the 2 years of data suggest that camera efficacy is not fully consistent across ecological conditions and make the case for integrating other sources of space-use data
Small-scale intraspecific life history variation in herbivorous spider mites (Tetranychus pacificus) is associated with host plant cultivar.
Life history variation is a general feature of arthropod systems, but is rarely included in models of field or laboratory data. Most studies assume that local processes occur identically across individuals, ignoring any genetic or phenotypic variation in life history traits. In this study, we tested whether field populations of Pacific spider mites (Tetranychus pacificus) on grapevines (Vitis vinifera) display significant intraspecific life history variation associated with host plant cultivar. To address this question we collected individuals from sympatric vineyard populations where either Zinfandel or Chardonnay were grown. We then conducted a "common garden experiment" of mites on bean plants (Phaseolus lunatus) in the laboratory. Assay populations were sampled non-destructively with digital photography to quantify development times, survival, and reproductive rates. Two classes of models were fit to the data: standard generalized linear mixed models and a time-to-event model, common in survival analysis, that allowed for interval-censored data and hierarchical random effects. We found a significant effect of cultivar on development time in both GLMM and time-to-event analyses, a slight cultivar effect on juvenile survival, and no effect on reproductive rate. There were shorter development times and a trend towards higher juvenile survival in populations from Zinfandel vineyards compared to those from Chardonnay vineyards. Lines of the same species, originating from field populations on different host plant cultivars, expressed different development times and slightly different survival rates when reared on a common host plant in a common environment
Building integral projection models with non-independent vital rates
Population dynamics are functions of several demographic processes including survival, reproduction, somatic growth, and maturation. The rates or probabilities for these processes can vary by time, by location, and by individual. These processes can co‐vary and interact to varying degrees, e.g., an animal can only reproduce when it is in a particular maturation state. Population dynamics models that treat the processes as independent may yield somewhat biased or imprecise parameter estimates, as well as predictions of population abundances or densities. However, commonly used integral projection models (IPMs) typically assume independence across these demographic processes. We examine several approaches for modelling between process dependence in IPMs and include cases where the processes co‐vary as a function of time (temporal variation), co‐vary within each individual (individual heterogeneity), and combinations of these (temporal variation and individual heterogeneity). We compare our methods to conventional IPMs, which treat vital rates independent, using simulations and a case study of Soay sheep (Ovis aries). In particular, our results indicate that correlation between vital rates can moderately affect variability of some population‐level statistics. Therefore, including such dependent structures is generally advisable when fitting IPMs to ascertain whether or not such between vital rate dependencies exist, which in turn can have subsequent impact on population management or life‐history evolution
Sequential Monte Carlo Methods in the nimble and nimbleSMC R Packages
nimble is an R package for constructing algorithms and conducting inference on hierarchical models. The nimble package provides a unique combination of flexible model specification and the ability to program model-generic algorithms. Specifically, the package allows users to code models in the BUGS language, and it allows users to write algorithms that can be applied to any appropriate model. In this paper, we introduce the nimbleSMC R package. nimbleSMC contains algorithms for state-space model analysis using sequential Monte Carlo (SMC) techniques that are built using nimble. We first provide an overview of state-space models and commonly-used SMC algorithms. We then describe how to build a state-space model in nimble and conduct inference using existing SMC algorithms within nimbleSMC. SMC algorithms within nimbleSMC currently include the bootstrap filter, auxiliary particle filter, ensemble Kalman filter, IF2 method of iterated filtering, and a particle Markov chain Monte Carlo (MCMC) sampler. These algorithms can be run in R or compiled into C++ for more efficient execution. Examples of applying SMC algorithms to linear autoregressive models and a stochastic volatility model are provided. Finally, we give an overview of how model-generic algorithms are coded within nimble by providing code for a simple SMC algorithm. This illustrates how users can easily extend nimble's SMC methods in high-level code
Nested Adaptation of MCMC Algorithms
© 2020 International Society for Bayesian Analysis Markov chain Monte Carlo (MCMC) methods are ubiquitous tools for simulation-based inference in many fields but designing and identifying good MCMC samplers is still an open question. This paper introduces a novel MCMC algorithm, namely, Nested Adaptation MCMC. For sampling variables or blocks of variables, we use two levels of adaptation where the inner adaptation optimizes the MCMC performance within each sampler, while the outer adaptation explores the space of valid kernels to find the optimal samplers. We provide a theoretical foundation for our approach. To show the generality and usefulness of the approach, we describe a framework using only standard MCMC samplers as candidate samplers and some adaptation schemes for both inner and outer iterations. In several benchmark problems, we show that our proposed approach substantially outperforms other approaches, including an automatic blocking algorithm, in terms of MCMC efficiency and computational time
- …