Many important demographic processes are seasonal, including survival. For many species, mortality risk is significantly higher at certain times of the year than at others, whether because resources are scarce, susceptibility to predators or disease is high, or both. Despite the importance of survival modelling in wildlife sciences, no tools are available to estimate the peak, duration and relative importance of these ‘seasons of mortality’. We present We illustrate the periodic hazard function model and workflow of cyclomort with simulated data. We then estimate mortality seasons of two caribou The
In clades that contain large numbers of undescribed species, DNA barcoding methods can serve as a first pass at identifying species limits. The General Mixed Yule Coalescent (GMYC) approach to species delimitation is accurate under certain conditions that are difficult to verify in the very clades for which the GMYC holds the greatest appeal. To circumvent this challenge, we have developed an R package for assessing how well the statistical model implemented in the GMYC fits empirical data. Our approach uses either a parametric bootstrap or a posterior predictive simulation to evaluate model fit. Computational requirements are modest for our package, and most analyses can be completed within minutes to an hour on a typical laptop, depending on whether a user selects a maximum likelihood or Bayesian framework. Results of simulation testing indicate that our approach is effective for assessing the utility of the GMYC model and suggest that it should be included in the analytical pipeline whenever researchers apply the GMYC to clades with unknown species boundaries.
- NSF-PAR ID:
- 10453760
- Publisher / Repository:
- Wiley-Blackwell
- Date Published:
- Journal Name:
- Methods in Ecology and Evolution
- Volume:
- 12
- Issue:
- 3
- ISSN:
- 2041-210X
- Page Range / eLocation ID:
- p. 487-493
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
For everything there is a season: Analysing periodic mortality patterns with the cyclomort r package
Abstract cyclomort , anr package that estimates the timing, duration and intensity of any number of mortality seasons with reliable confidence intervals. The package includes a model selection approach to determine the number of mortality seasons and to test whether seasons of mortality vary across discrete grouping factors.Rangifer tarandus populations that have strikingly different mortality patterns, including different numbers and timing of mortality peaks, and a marked change in one population over time.cyclomort package was developed to estimate mortality seasons for wildlife, but the package can model any time‐to‐event processes with a periodic component. -
Summary Recent studies have demonstrated that ecological processes that shape community structure and dynamics change along environmental gradients. However, much less is known about how the emergence of the gradients themselves shape the evolution of species that underlie community assembly. In this study, we address how the creation of novel environments leads to community assembly via two nonmutually exclusive processes: immigration and ecological sorting of pre‐adapted clades (ISPC), and recent adaptive diversification (RAD). We study these processes in the context of the elevational gradient created by the uplift of the Central Andes.
We develop a novel approach and method based on the decomposition of species turnover into within‐ and among‐clade components, where clades correspond to lineages that originated before mountain uplift. Effects of ISPC and RAD can be inferred from how components of turnover change with elevation. We test our approach using data from over 500 Andean forest plots.
We found that species turnover between communities at different elevations is dominated by the replacement of clades that originated before the uplift of the Central Andes.
Our results suggest that immigration and sorting of clades pre‐adapted to montane habitats is the primary mechanism shaping tree communities across elevations.
-
spOccupancy: An R package for single‐species, multi‐species, and integrated spatial occupancy models
Abstract Occupancy modelling is a common approach to assess species distribution patterns, while explicitly accounting for false absences in detection–nondetection data. Numerous extensions of the basic single‐species occupancy model exist to model multiple species, spatial autocorrelation and to integrate multiple data types. However, development of specialized and computationally efficient software to incorporate such extensions, especially for large datasets, is scarce or absent.
We introduce the
spOccupancy R package designed to fit single‐species and multi‐species spatially explicit occupancy models. We fit all models within a Bayesian framework using Pólya‐Gamma data augmentation, which results in fast and efficient inference.spOccupancy provides functionality for data integration of multiple single‐species detection–nondetection datasets via a joint likelihood framework. The package leverages Nearest Neighbour Gaussian Processes to account for spatial autocorrelation, which enables spatially explicit occupancy modelling for potentially massive datasets (e.g. 1,000s–100,000s of sites).spOccupancy provides user‐friendly functions for data simulation, model fitting, model validation (by posterior predictive checks), model comparison (using information criteria and k‐fold cross‐validation) and out‐of‐sample prediction. We illustrate the package's functionality via a vignette, simulated data analysis and two bird case studies.The
spOccupancy package provides a user‐friendly platform to fit a variety of single and multi‐species occupancy models, making it straightforward to address detection biases and spatial autocorrelation in species distribution models even for large datasets. -
Abstract Gene flow is increasingly recognized as an important macroevolutionary process. The many mechanisms that contribute to gene flow (e.g. introgression, hybridization, lateral gene transfer) uniquely affect the diversification of dynamics of species, making it important to be able to account for these idiosyncrasies when constructing phylogenetic models. Existing phylogenetic‐network simulators for macroevolution are limited in the ways they model gene flow.
We present
SiPhyNetwork , an R package for simulating phylogenetic networks under a birth–death‐hybridization process.Our package unifies the existing birth–death‐hybridization models while also extending the toolkit for modelling gene flow. This tool can create patterns of reticulation such as hybridization, lateral gene transfer, and introgression.
Specifically, we model different reticulate events by allowing events to either add, remove or keep constant the number of lineages. Additionally, we allow reticulation events to be trait dependent, creating the ability to model the expanse of isolating mechanisms that prevent gene flow. This tool makes it possible for researchers to model many of the complex biological factors associated with gene flow in a phylogenetic context.
-
Abstract Biodiversity studies rely heavily on estimates of species' distributions often obtained through ecological niche modelling. Numerous software packages exist that allow users to model ecological niches using machine learning and statistical methods. However, no existing package with a graphical user interface allows users to perform model calibration and selection based on convex forms such as ellipsoids, which may match fundamental ecological niche shapes better, incorporating tools for exploring, modelling, and evaluating niches and distributions that are intuitive for both novice and proficient users.
Here we describe an
r package, Niche Tool Box (ntbox ), that allows users to conduct all processing steps involved in ecological niche modelling: downloading and curating occurrence data, obtaining and transforming environmental data layers, selecting environmental variables, exploring relationships between geographic and environmental spaces, calibrating and selecting ellipsoid models, evaluating models using binomial and partial ROC tests, assessing extrapolation risk, and performing geographic information system operations via a graphical user interface. A summary of the entire workflow is produced for use as a stand‐alone algorithm or as part of research reports.The method is explained in detail and tested via modelling the threatened feline species
Leopardus wiedii . Georeferenced occurrence data for this species are queried to display both point occurrences and the IUCN extent of occurrence polygon (IUCN, 2007). This information is used to illustrate tools available for accessing, processing and exploring biodiversity data (e.g. number of occurrences and chronology of collecting) and transforming environmental data (e.g. a summary PCA for 19 bioclimatic layers). Visualizations of three‐dimensional ecological niches modelled as minimum volume ellipsoids are developed with ancillary statistics. This niche model is then projected to geographic space, to represent a corresponding potential suitability map.Using
ntbox allows a fast and straightforward means by which to retrieve and manipulate occurrence and environmental data, which can then be implemented in model calibration, projection and evaluation for assessing distributions of species in geographic space and their corresponding environmental combinations.