adiv (version 1.2)

rare_Rao: Functional Rarefaction for Species Abundance Data

Description

The function Rare_Rao performs distance-based rarefaction curves using species abundance data. It finds the expected functional diversity (if functional distances between species are used) as a function of the sampling effort. Two approaches are available: an analytical solution, a resampling approach.

Usage

rare_Rao(comm, dis, sim = TRUE, resampling = 999, formula = c("QE", "EDI"))

Arguments

comm

a data frame or a matrix with samples as rows, species as columns, and abundance or frequencies as entries.

dis

an object of class dist containing pairwise distances among species. The distance matrix should be squared Euclidean or simply Euclidean (see definition and details in Appendix S1 of Ricotta et al. (2012)). Alternatively, if dis is NULL, the Gini-Simpson index is used.

sim

a logical; if TRUE, the resampling approach is used; if FALSE, the analytical solution is given.

resampling

number of times data are resampled to calculate the mean functional rarefaction curve (used if sim=TRUE)

formula

Either "QE" (default) or "EDI". See details.

Value

If sim = TRUE, the function returns a data frame containing the Expected Rao Quadratic Entropy (column 'ExpRao'), the limits of the 95% Confidence Interval (columns 'LeftIC' and 'RightIC') for each subsample dimension (M) out of the total set of samples (N). If sim = FALSE, the function returns a data frame containing the analytical solution for the Expected Rao Quadratic Entropy (column 'ExpRao') for each subsample dimension (M) out of the total set of samples (N).

Details

If formula = "QE", the definition of the quadratic entropy is: $$QE(\mathbf{p}_i,\mathbf{D})=\sum_{k=1}^S\sum_{l=1}^S p_{k|i}p_{k|j}d_{kl}$$ where \(\mathbf{p}_i=(p_{1|i}, ..., p_{k|i}, ..., p_{S|i})\) is the vector of relative species abundance within sample i; S is the number of species; \(\mathbf{D}=(d_{kl})\) is the matrix of (phylogenetic or functional) dissimilarities among species, and \(d_{kl}\) is the (phylogenetic or functional) dissimilarity between species k and l.

If formula = "EDI", the definition of the quadratic entropy is: $$EDI(\mathbf{p}_i,\mathbf{D})=\sum_{k=1}^S\sum_{l=1}^S p_{k|i}p_{k|j}\frac{d_{kl}^2}{2}$$ EDI stands for the Euclidean Diversity Index of Champely and Chessel (2002) (equation 3 in Pavoine et al. 2004).

References

Ricotta, C., Pavoine, S., Bacaro, G., Acosta, A. (2012) Functional rarefaction for species abundance data. Methods in Ecology and Evolution, 3, 519--525.

Champely, S. and Chessel, D. (2002) Measuring biological diversity using Euclideanmetrics. Environmental and Ecological Statistics, 9, 167--177.

Pavoine, S., Dufour, A.B., Chessel, D. (2004) From dissimilarities among species to dissimilarities among communities: a double principal coordinate analysis. Journal of Theoretical Biology, 228, 523--537.

See Also

See Also QE

Examples

# NOT RUN {
data(aviurba, package="ade4")
distances<-dist.ktab(ktab.list.df(list(aviurba$traits)), type = "N")
# The distances should be squared Euclidean; 
# note that Euclidean distances can be used 
# as they also are squared Euclidean.

abundances<- aviurba$fau
rare_Rao(abundances, distances, sim = TRUE, resampling = 100)
rare_Rao(abundances, distances, sim = FALSE)
# }