Description Usage Arguments Details Value Author(s) Examples
Functions for recursively culling community matrices of rare taxa and depauperate samples.
1 | cullMatrix(CommunityMatrix, Rarity = 2, Richness = 2, Silent = FALSE)
|
CommunityMatrix |
a matrix |
Rarity |
a whole number |
Richness |
a whole number |
Silent |
logical |
Takes a community matrix (see presenceMatrix
or abundanceMatrix
) and removes all samples with fewer than a certain number of taxa and all taxa that occur below a certain threshold of samples. The function operates recursively, and will check to see if removing a rare taxon drops a sampe below the input minimum richness and vice-versa. This means that it is possible to eliminate all taxa and samples if the rarity and richness minimums are too high. If the Silent
argument is set to FALSE
the function will throw an error and print a warning if no taxa or samples are left after culling. If Silent
is set to TRUE
the function will simply return NULL
. The latter case is useful if many matrices are being culled as a part of a loop, and you do not want to break the loop with an error.
These functions originally appeared in the R script appendix of Holland, S.M. and A. Zaffos (2011) "Niche conservatism along an onshore-offshore gradient". Paleobiology 37:270-286.
A community matrix with depauperate samples and rare taxa removed.
Andrew A. Zaffos
1 2 3 4 5 6 7 8 | # Download a test dataset of pleistocene bivalves.
# DataPBDB<-downloadPBDB(Taxa="Bivalvia",StartInterval="Pleistocene",StopInterval="Pleistocene")
# Create a community matrix with tectonic plates as "samples".
# CommunityMatrix<-abundanceMatrix(DataPBDB,"geoplate")
# Remove taxa that occur in less than 5 samples and samples with fewer than 25 taxa.
# cullMatrix(CommunityMatrix,Rarity=5,Richness=25,Silent=FALSE)
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.