kernel.sad | R Documentation |
Fits the SAD to community abundance data based on n-dimensional hypervolumes.
kernel.sad(comm, octaves = TRUE, scale = FALSE, raref = 0, runs = 100)
comm |
A 'Hypervolume' or 'HypervolumeList' object necessarily built using function kernel.build. |
octaves |
a boolean indicating whether octaves should be calculated. |
scale |
scale y-axis to sum 1. |
raref |
An integer specifying the number of individuals for rarefaction (individual based). If raref < 1 no rarefaction is made. If raref = 1 rarefaction is made by the minimum abundance among all sites. If raref > 1 rarefaction is made by the abundance indicated. If not specified, default is 0. |
runs |
Number of resampling runs for rarefaction. If not specified, default is 100. |
The Species Abundance Distribution describes the commonness and rarity in ecological systems. It was recently expanded to accomodate phylegenetic and functional differences between species (Matthews et al., subm.). Classes defined as n = 1, 2-3, 4-7, 8-15, .... Rarefaction allows comparison of sites with different total abundances.
A vector or matrix with the different values per class per community.
Matthews et al. (subm.) Phylogenetic and functional dimensions of the species abundance distribution.
## Not run:
comm = rbind(c(1,3,0,5,3), c(3,2,5,1,0))
colnames(comm) = c("SpA", "SpB", "SpC", "SpD", "SpE")
rownames(comm) = c("Site 1", "Site 2")
trait = data.frame(body = c(1,2,3,4,4), beak = c(1,5,4,1,2))
rownames(trait) = colnames(comm)
hv = kernel.build(comm, trait)
kernel.sad(hv, scale = TRUE)
kernel.sad(hv, octaves = FALSE)
kernel.sad(hv, raref = TRUE)
## End(Not run)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.