skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.

Explore Research Products in the PAR
It may take a few hours for recently added research products to appear in PAR search results.


  • NSF Public Access
  • Search Results
  • Maximum likelihood outperforms binning methods for detecting differences in abundance size spectra across environmental gradients
Title: Maximum likelihood outperforms binning methods for detecting differences in abundance size spectra across environmental gradients
Abstract Individual body size distributions (ISD) within communities are remarkably consistent across habitats and spatiotemporal scales and can be represented by size spectra, which are described by a power law. The focus of size spectra analysis is to estimate the exponent () of the power law. A common application of size spectra studies is to detect anthropogenic pressures.Many methods have been proposed for estimating most of which involve binning the data, counting the abundance within bins, and then fitting an ordinary least squares regression in log–log space. However, recent work has shown that binning procedures return biased estimates of compared to procedures that directly estimate using maximum likelihood estimation (MLE). While it is clear that MLE produces less biased estimates of site‐specificλ's, it is less clear how this bias affects the ability to test for changes inλacross space and time, a common question in the ecological literature.Here, we used simulation to compare the ability of two normalised binning methods (equal logarithmic and log2bins) and MLE to (1) recapture known values of , and (2) recapture parameters in a linear regression measuring the change in across a hypothetical environmental gradient. We also compared the methods using two previously published body size datasets across a natural temperature gradient and an anthropogenic pollution gradient.Maximum likelihood methods always performed better than common binning methods, which demonstrated consistent bias depending on the simulated values of . This bias carried over to the regressions, which were more accurate when was estimated using MLE compared to the binning procedures. Additionally, the variance in estimates using MLE methods is markedly reduced when compared to binning methods.The error induced by binning methods can be of similar magnitudes as the variation previously published in experimental and observational studies, bringing into question the effect sizes of previously published results. However, while the methods produced different regression slope estimates, they were in qualitative agreement on the sign of those slopes (i.e. all negative or all positive). Our results provide further support for the direct estimation of and its relative variation across environmental gradients using MLE over the more common methods of binning.  more » « less
Award ID(s):
2106067
PAR ID:
10627837
Author(s) / Creator(s):
; ; ;
Publisher / Repository:
British Ecological Society
Date Published:
Journal Name:
Journal of Animal Ecology
Volume:
93
Issue:
3
ISSN:
0021-8790
Page Range / eLocation ID:
267 to 280
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. ; ; ; (, Methods in Ecology and Evolution)
    Abstract A fundamental pattern in ecology is that smaller organisms are more abundant than larger organisms. This pattern is known as the individual size distribution (ISD), which is the frequency distribution of all individual body sizes in an ecosystem.The ISD is described by a power law and a major goal of size spectra analyses is to estimate the exponent of the power law,λ. However, while numerous methods have been developed to do this, they have focused almost exclusively on estimatingλfrom single samples.Here, we develop an extension of the truncated Pareto distribution within the probabilistic modelling language Stan. We use it to estimate multipleλs simultaneously in a hierarchical modelling approach.The most important result is the ability to examine hypotheses related to size spectra, including the assessment of fixed and random effects, within a single Bayesian generalized mixed model. While the example here uses size spectra, the technique can also be generalized to any data that follow a power law distribution. 
    more » « less
  2. ; ; ; (, Methods in Ecology and Evolution)
    Abstract Understanding patterns of diversity is central to ecology and conservation, yet estimates of diversity are often biased by imperfect detection. In recent years, multi‐species occupancy models (MSOM) have been developed as a statistical tool to account for species‐specific heterogeneity in detection while estimating true measures of diversity. Although the power of these models has been tested in various ways, their ability to estimate gamma diversity—or true community size,Nis a largely unrecognized feature that needs rigorous evaluation.We use both simulations and an empirical dataset to evaluate the bias, precision, accuracy and coverage of estimates ofNfrom MSOM compared to the widely applied iChao2 non‐parametric estimator. We simulated 5,600 datasets across seven scenarios of varying average occupancy and detectability covariates, as well as varying numbers of sites, replicates and true community size. Additionally, we use a real dataset of surveys over 9 years (where species accumulation reached an asymptote, indicating trueN), to estimateNfrom each annual survey.Simulations showed that both MSOM and iChao2 estimators are generally accurate (i.e. unbiased and precise) except under unideal scenarios where mean species occupancy is low. In such scenarios, MSOM frequently overestimatedN. Across all scenarios, MSOM estimates were less certain than iChao2, but this led to over‐confident iChao2 estimates that showed poor coverage. Results from the real dataset largely confirmed the simulation findings, with MSOM estimates showing greater accuracy and coverage than iChao2.Community ecologists have a wide choice of analytical methods, and both iChao2 and MSOM estimates ofNare substantially preferable to raw species counts. The simplicity of non‐parametric estimators has obvious advantages, but our results show that in many cases, MSOM may provide superior estimates that also account more accurately for uncertainty. Both methods can show strong bias when average occupancy is very low, and practitioners should show caution when using estimates derived from either method under such conditions. 
    more » « less
  3. ; ; (, Functional Ecology)
    Abstract Eye size varies notably across taxa. Much work suggests that this variation is driven by contrasting ecological selective pressures. However, evaluations of the relationship between ecological factors and shifts in eye size have largely occurred at the macroevolutionary scale. Experimental tests in nature are conspicuously absent.Trinidadian killifish,Rivulus hartii, are found across fish communities that differ in predation intensity. We recently showed that increased predation is associated with the evolution of a smaller eye. Here, we test how divergent predatory regimes alter the trajectory of eye size evolution using comparative mark–recapture experiments in multiple streams.We found that increases in eye size are associated with enhanced survival, irrespective of predation intensity. More importantly, eye size is associated with enhanced growth in communities that lack predators, while this trend is absent when predators are present.Such results argue that increased competition for food in sites that lack predators is the key driver of eye size evolution. Aplain language summaryis available for this article. 
    more » « less
  4. ; ; ; ; (, PLOS ONE)
    Dam, Hans G. (Ed.)
    Recent research has revealed the diversity and biomass of life across ecosystems, but how that biomass is distributed across body sizes of all living things remains unclear. We compile the present-day global body size-biomass spectra for the terrestrial, marine, and subterranean realms. To achieve this compilation, we pair existing and updated biomass estimates with previously uncatalogued body size ranges across all free-living biological groups. These data show that many biological groups share similar ranges of body sizes, and no single group dominates size ranges where cumulative biomass is highest. We then propagate biomass and size uncertainties and provide statistical descriptions of body size-biomass spectra across and within major habitat realms. Power laws show exponentially decreasing abundance (exponent -0.9±0.02 S.D.,R2= 0.97) and nearly equal biomass (exponent 0.09±0.01,R2= 0.56) across log size bins, which resemble previous aquatic size spectra results but with greater organismal inclusivity and global coverage. In contrast, a bimodal Gaussian mixture model describes the biomass pattern better (R2= 0.86) and suggests small (~10−15g) and large (~107g) organisms outweigh other sizes by one order magnitude (15 and 65 Gt versus ~1 Gt per log size). The results suggest that the global body size-biomass relationships is bimodal, but substantial one-to-two orders-of-magnitude uncertainty mean that additional data will be needed to clarify whether global-scale universal constraints or local forces shape these patterns. 
    more » « less
  5. (, Freshwater Biology)
    Abstract A positive relationship between body size and trophic position is often assumed in ecology, but efforts to confirm the generality of this relationship for freshwater vertebrates have produced mixed results. Some authors have tested for among‐species relationships, using species‐level estimates of average size and trophic position. Others have used individual‐level body size and trophic position data to test for within‐species relationships. However, no study has yet estimatedprevalence, defined here as the fraction of total standing stock biomass within a given ecosystem that consists of taxa exhibiting positive size versus trophic position relationships.Individual‐level estimates of body size and relative trophic position (inferred from bulk‐tissue nitrogen stable isotopes) were collected for vertebrates in six temperate streams. In each stream, all locally occurring vertebrate taxa were collected and closed population depletion samples were used to obtain standing stock biomass estimates.Ordinary least squares regression was used with the individual‐level data to test for positive relationships between body size and relative trophic position (r‐STPRs). A separater‐STPR model was tested for each vertebrate species collected at the six study sites. Linear mixed‐effects modelling was then used to test for differingr‐STPRs among species.Prevalence of ther‐STPR was calculated for each site by summing the standing stock biomass of all taxa that exhibited a statistically significant, positive STPR at a given site, then dividing this number by the total standing stock biomass of vertebrates at the site.Significant, positiver‐STPRs were detected in 15 of 22 species × site regression models. Prevalence of species with positiver‐STPRs ranged from 91 to 100% of the total vertebrate biomass in western streams and from 45 to 66% in eastern streams.Results confirm that positiver‐STPRs are characteristic of many of the vertebrate taxa considered in this study. Furthermore, species that exhibit positiver‐STPRs comprise a clear majority of the standing stock biomass in five of six study streams. By using individual‐level data to account for prevalence, a more complete understanding of size‐dependent trophic dynamics should be possible in freshwater ecosystems. 
    more » « less
  • Citation Formats
  • MLA
  • APA
  • Chicago
  • BibTeX
  • Save / Share this Record
  • Send to Email