Megan D. Higgs

Individual heterogeneity in reproductive rates and cost of reproduction in a long‐lived vertebrate

Individual variation in reproductive success is a key feature of evolution, but also has important implications for predicting population responses to variable environments. Although such individual variation in reproductive outcomes has been reported in numerous studies, most analyses to date have not considered whether these realized differences were due to latent individual heterogeneity in reproduction or merely random chance causing different outcomes among like individuals. Furthermore, latent heterogeneity in fitness components might be expressed differently in contrasted environmental conditions, an issue that has only rarely been investigated. Here, we assessed (i) the potential existence of latent individual heterogeneity and (ii) the nature of its expression (fixed vs. variable) in a population of female Weddell seals (Leptonychotes weddellii), using a hierarchical modeling approach on a 30-year mark–recapture data set consisting of 954 individual encounter histories. We found strong support for the existence of latent individual heterogeneity in the population, with “robust” individuals expected to produce twice as many pups as “frail” individuals. Moreover, the expression of individual heterogeneity appeared consistent, with only mild evidence that it might be amplified when environmental conditions are severe. Finally, the explicit modeling of individual heterogeneity allowed us to detect a substantial cost of reproduction that was not evidenced when the heterogeneity was ignored.

A clipped latent variable model for spatially correlated ordered categorical data

We propose a model for a point-referenced spatially correlated ordered categorical response and methodology for inference. Models and methods for spatially correlated continuous response data are widespread, but models for spatially correlated categorical data, and especially ordered multi-category data, are less developed. Bayesian models and methodology have been proposed for the analysis of independent and clustered ordered categorical data, and also for binary and count point-referenced spatial data. We combine and extend these methods to describe a Bayesian model for point-referenced (as opposed to lattice) spatially correlated ordered categorical data. We include simulation results and show that our model offers superior predictive performance as compared to a non-spatial cumulative probit model and a more standard Bayesian generalized linear spatial model. We demonstrate the usefulness of our model in a real-world example to predict ordered categories describing stream health within the state of Maryland.

Use of posterior predictive checks as an inferential tool for investigating individual heterogeneity in animal population vital rates.

The investigation of individual heterogeneity in vital rates has recently received growing attention among population ecologists. Individual heterogeneity in wild animal populations has been accounted for and quantified by including individually varying effects in models for mark–recapture data, but the real need for underlying individual effects to account for observed levels of individual variation has recently been questioned by the work of Tuljapurkar et al. (Ecology Letters, 12, 93, 2009) on dynamic heterogeneity. Model-selection approaches based on information criteria or Bayes factors have been used to address this question. Here, we suggest that, in addition to model-selection, model-checking methods can provide additional important insights to tackle this issue, as they allow one to evaluate a models misfit in terms of ecologically meaningful measures. Specifically, we propose the use of posterior predictive checks to explicitly assess discrepancies between a model and the data, and we explain how to incorporate model checking into the inferential process used to assess the practical implications of ignoring individual heterogeneity. Posterior predictive checking is a straightforward and flexible approach for performing model checks in a Bayesian framework that is based on comparisons of observed data to model-generated replications of the data, where parameter uncertainty is incorporated through use of the posterior distribution. If discrepancy measures are chosen carefully and are relevant to the scientific context, posterior predictive checks can provide important information allowing for more efficient model refinement. We illustrate this approach using analyses of vital rates with long-term mark–recapture data for Weddell seals and emphasize its utility for identifying shortfalls or successes of a model at representing a biological process or pattern of interest. We show how posterior predictive checks can be used to strengthen inferences in ecological studies. We demonstrate the application of this method on analyses dealing with the question of individual reproductive heterogeneity in a population of Antarctic pinnipeds.

A multi-scale assessment of animal aggregation patterns to understand increasing pathogen seroprevalence

Understanding how animal density is related to pathogen transmission is important to develop effective disease control strategies, but requires measuring density at a scale relevant to transmission. However, this is not straightforward or well-studied among large mammals with group sizes that range several orders of magnitude or aggregation patterns that vary across space and time. To address this issue, we examined spatial variation in elk (Cervus canadensis) aggregation patterns and brucellosis across 10 regions in the Greater Yellowstone Area where previous studies suggest the disease may be increasing. We hypothesized that rates of increasing brucellosis would be better related to the frequency of large groups than mean group size or population density, but we examined whether other measures of density would also explain rising seroprevalence. To do this, we measured wintering elk density and group size across multiple spatial and temporal scales from monthly aerial surveys. We used Bayesian hierarchical models and 20 years of serologic data to estimate rates of increase in brucellosis within the 10 regions, and to examine the linear relationships between these estimated rates of increase and multiple measures of aggregation. Brucellosis seroprevalence increased over time in eight regions (one region showed an estimated increase from 0.015 in 1991 to 0.26 in 2011), and these rates of increase were positively related to all measures of aggregation. The relationships were weaker when the analysis was restricted to areas where brucellosis was present for at least two years, potentially because aggregation was related to disease-establishment within a population. Our findings suggest that (1) group size did not explain brucellosis increases any better than population density and (2) some elk populations may have high densities with small groups or lower densities with large groups, but brucellosis is likely to increase in either scenario. In this case, any one control method such as reducing population density or group size may not be sufficient to reduce transmission. This study highlights the importance of examining the density-transmission relationship at multiple scales and across populations before broadly applying disease control strategies.

Microsatellites indicate minimal barriers to mule deer Odocoileus hemionus dispersal across Montana, USA

To better understand the future spread of chronic wasting disease, we conducted a genetic assessment of mule deer Odocoileus hemionus population structure across the state of Montana, USA. Individual based analyses were used to test for population structure in the absence of a priori designations of population membership across the sampling area. Samples from the states of Wyoming, Colorado and Utah were also included in the analysis to provide a geographic context to the levels of population structure observed within Montana. Results showed that mule deer across our entire study region were characterized by weak isolation by distance and a lack of spatial autocorrelation at distances > 10 km. We found evidence for contemporary male bias in dispersal, with female mule deer exhibiting higher mean individual pairwise genetic distance than males. We tested for potential homogenizing effects of past translocations within Montana, but were unable to detect a genetic signature of these events. Our results indicate high levels of connectivity among mule deer populations in Montana and suggest few, if any, detectable barriers to mule deer gene flow or chronic wasting disease transmission.

Inferential consequences of modeling rather than measuring snow accumulation in studies of animal ecology

It is increasingly common for studies of animal ecology to use model-based predictions of environmental variables as explanatory or predictor variables, even though model prediction uncertainty is typically unknown. To demonstrate the potential for misleading inferences when model predictions with error are used in place of direct measurements, we compared snow water equivalent (SWE) and snow depth as predicted by the Snow Data Assimilation System (SNODAS) to field measurements of SWE and snow depth. We examined locations on elk (Cervus canadensis) winter ranges in western Wyoming, because modeled data such as SNODAS output are often used for inferences on elk ecology. Overall, SNODAS predictions tended to overestimate field measurements, prediction uncertainty was high, and the difference between SNODAS predictions and field measurements was greater in snow shadows for both snow variables compared to non-snow shadow areas. We used a simple simulation of snow effects on the probability of an elk being killed by a predator to show that, if SNODAS prediction uncertainty was ignored, we might have mistakenly concluded that SWE was not an important factor in where elk were killed in predatory attacks during the winter. In this simulation, we were interested in the effects of snow at finer scales (< 1 km2) than the resolution of SNODAS. If bias were to decrease when SNODAS predictions are averaged over coarser scales, SNODAS would be applicable to population-level ecology studies. In our study, however, averaging predictions over moderate to broad spatial scales (9-2200 km2) did not reduce the differences between SNODAS predictions and field measurements. This study highlights the need to carefully evaluate two issues when using model output as an explanatory variable in subsequent analysis: (1) the models resolution relative to the scale of the ecological question of interest and (2) the implications of prediction uncertainty on inferences when using model predictions as explanatory or predictor variables.

Discretized and Aggregated: Modeling Dive Depth of Harbor Seals from Ordered Categorical Data with Temporal Autocorrelation

Ordered categorical data are pervasive in environmental and ecological data, and often arise from constraints that require discretizing a continuous variable into ordered categories. A great deal of data have been collected toward the study of marine mammal dive behavior using satellite depth recorders (SDRs), which often discretize a continuous variable such as depth. Additionally, data storage or transmission constraints may also necessitate the aggregation of data over time intervals of a specified length. The categorization and aggregation create a time series of ordered multicategory counts for each animal, which present challenges in terms of statistical modeling and practical interpretation. We describe an intuitive strategy for modeling such aggregated, ordered categorical data allowing for inference regarding the category probabilities and a measure of central tendency on the original scale of the data (e.g., meters), along with incorporation of temporal correlation and overdispersion. The strategy extends covariate-specific cutpoint models for ordinal data. We demonstrate the method in an analysis of SDR dive-depth data collected on harbor seals in Alaska. The primary goal of the analysis is to assess the relationship of covariates, such as time of day, with number of dives and maximum depth of dives. We also predict missing values and introduce novel graphical summaries of the data and results.