Dataset Information

The Influence of Matrix Size on Statistical Properties of Co-Occurrence and Limiting Similarity Null Models.

ABSTRACT: Null models exploring species co-occurrence and trait-based limiting similarity are increasingly used to explore the influence of competition on community assembly; however, assessments of common models have not thoroughly explored the influence of variation in matrix size on error rates, in spite of the fact that studies have explored community matrices that vary considerably in size. To determine how smaller matrices, which are of greatest concern, perform statistically, we generated biologically realistic presence-absence matrices ranging in size from 3-50 species and sites, as well as associated trait matrices. We examined co-occurrence tests using the C-Score statistic and independent swap algorithm. For trait-based limiting similarity null models, we used the mean nearest neighbour trait distance (NN) and the standard deviation of nearest neighbour distances (SDNN) as test statistics, and considered two common randomization algorithms: abundance independent trait shuffling (AITS), and abundance weighted trait shuffling (AWTS). Matrices as small as three × three resulted in acceptable type I error rates (p < 0.05) for both the co-occurrence and trait-based limiting similarity null models when exclusive p-values were used. The commonly used inclusive p-value (? or ?, as opposed to exclusive p-values; < or >) was associated with increased type I error rates, particularly for matrices with fewer than eight species. Type I error rates increased for limiting similarity tests using the AWTS randomization scheme when community matrices contained more than 35 sites; a similar randomization used in null models of phylogenetic dispersion has previously been viewed as robust. Notwithstanding other potential deficiencies related to the use of small matrices to represent communities, the application of both classes of null model should be restricted to matrices with 10 or more species to avoid the possibility of type II errors. Additionally, researchers should restrict the use of the AWTS randomization to matrices with fewer than 35 sites to avoid type I errors when testing for trait-based limiting similarity. The AITS randomization scheme performed better in terms of type I error rates, and therefore may be more appropriate when considering systems for which traits are not clustered by abundance.

SUBMITTER: Lavender TM

PROVIDER: S-EPMC4778770 | biostudies-literature | 2016

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

The Influence of Matrix Size on Statistical Properties of Co-Occurrence and Limiting Similarity Null Models.

Lavender Thomas Michael TM Schamp Brandon S BS Lamb Eric G EG

PloS one 20160304 3

Null models exploring species co-occurrence and trait-based limiting similarity are increasingly used to explore the influence of competition on community assembly; however, assessments of common models have not thoroughly explored the influence of variation in matrix size on error rates, in spite of the fact that studies have explored community matrices that vary considerably in size. To determine how smaller matrices, which are of greatest concern, perform statistically, we generated biologica ...[more]

PMID: 26942941

Similar Datasets

Project description:Drawing on a long history in macroecology, correlation analysis of microbiome datasets is becoming a common practice for identifying relationships or shared ecological niches among bacterial taxa. However, many of the statistical issues that plague such analyses in macroscale communities remain unresolved for microbial communities. Here, we discuss problems in the analysis of microbial species correlations based on presence-absence data. We focus on presence-absence data because this information is more readily obtainable from sequencing studies, especially for whole-genome sequencing, where abundance estimation is still in its infancy. First, we show how Pearson's correlation coefficient (r) and Jaccard's index (J)-two of the most common metrics for correlation analysis of presence-absence data-can contradict each other when applied to a typical microbiome dataset. In our dataset, for example, 14% of species-pairs predicted to be significantly correlated by r were not predicted to be significantly correlated using J, while 37.4% of species-pairs predicted to be significantly correlated by J were not predicted to be significantly correlated using r. Mismatch was particularly common among species-pairs with at least one rare species (<10% prevalence), explaining why r and J might differ more strongly in microbiome datasets, where there are large numbers of rare taxa. Indeed 74% of all species-pairs in our study had at least one rare species. Next, we show how Pearson's correlation coefficient can result in artificial inflation of positive taxon relationships and how this is a particular problem for microbiome studies. We then illustrate how Jaccard's index of similarity (J) can yield improvements over Pearson's correlation coefficient. However, the standard null model for Jaccard's index is flawed, and thus introduces its own set of spurious conclusions. We thus identify a better null model based on a hypergeometric distribution, which appropriately corrects for species prevalence. This model is available from recent statistics literature, and can be used for evaluating the significance of any value of an empirically observed Jaccard's index. The resulting simple, yet effective method for handling correlation analysis of microbial presence-absence datasets provides a robust means of testing and finding relationships and/or shared environmental responses among microbial taxa.

Project description:Co-occurring species often overlap in their use of resources and can interact in complex ways. However, shifts in environmental conditions or resource availability can lead to changes in patterns of species co-occurrence, which may be exacerbated by global escalation of human disturbances to ecosystems, including conservation-directed interventions. We investigated the relative abundance and co-occurrence of two naturally sympatric mammal species following two forms of environmental disturbance: wildfire and introduced predator control. Using 14 years of abundance data from repeat surveys at long-term monitoring sites in south-eastern Australia, we examined the association between a marsupial, the common brushtail possum Trichosurus vulpecula, and a co-occurring native rodent, the bush rat Rattus fuscipes. We asked: In a fox-controlled environment, are the abundances of common brushtail possums and bush rats affected by environmental disturbance and each other's presence? Using Bayesian regression models, we tested hypotheses that the abundance of each species would vary with changes in environmental and disturbance variables, and that the negative association between bush rats and common brushtail possums was stronger than the association between bush rats and disturbance. Our analyses revealed that bush rat abundance varied greatly in relation to environmental and disturbance variables, whereas common brushtail possums showed relatively limited variation in response to the same variables. There was a negative association between common brushtail possums and bush rats, but this association was weaker than the initial decline and subsequent recovery of bush rats in response to wildfires. Using co-occurrence analysis, we can infer negative relationships in abundance between co-occurring species, but to understand the impacts of such associations, and plan appropriate conservation measures, we require more information on interactions between the species and environmental variables. Co-occurrence can be a powerful and novel method to diagnose threats to communities and understand changes in ecosystem dynamics.

Dataset Information

The Influence of Matrix Size on Statistical Properties of Co-Occurrence and Limiting Similarity Null Models.

Publications

The Influence of Matrix Size on Statistical Properties of Co-Occurrence and Limiting Similarity Null Models.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets