Podani-family and Baselga-family decompositions of the Jaccard and S<U+00F8>rensen dissimilarity coefficients and their quantitative forms (Ruzicka and percentage difference) into replacement and richness difference components, for species presence-absence or abundance data, as described in Legendre (2014).
beta.div.comp(mat, coef = "J", quant = FALSE, save.abc = FALSE)Community composition data (data.frame or matrix).
Family of coefficients to be computed.
"J" or "Jaccard": Podani family, Jaccard-based indices.
"BS" <U+2013> Baselga family, S<U+00F8>rensen-based indices.
"BJ": Baselga family, Jaccard-based indices.
"N": Podani & Schmera (2011) relativized nestedness index.
The quantitative form of the S<U+00F8>rensen dissimilarity is the percentage difference index. The quantitative form of the Jaccard dissimilarity is the Ruzicka index.
If TRUE, compute the quantitative forms of replacement,
nestedness and D. If FALSE, compute the presence-absence forms of the
coefficients.
If TRUE, save the matrices of parameters a, b and c
used in presence-absence calculations.
A list containing the following results:
repl: Replacement matrix, class = dist.
rich: Richness/abundance difference or nestedness matrix (class
dist). With options "BJ", "BS" and "N", rich contains
nestedness indices. With option "N", the repl[i,j] and rich[i,j] values do
not add up to D[i,j].
D: Dissimilarity matrix (classdist).
part: Beta diversity partitioning vector:
BDtotal (total beta diversity) = sum(D.ij)/(n*(n-1)) (Legendre & De C<U+00E1>ceres 2013). This is equal to sum(d.ij^2)/(n*(n-1)) where d.ij = sqrt(D.ij). The dissimilarities are square-rooted because the Jaccard, S<U+00F8>rensen, Ruzicka and percentage difference indices are not Euclidean.
Repl = Total replacement diversity.
RichDiff|Nes = Total richness difference diversity (or nestedness).
Repl/BDtotal = Total replacement diversity/Total beta diversity.
RichDiff/BDtotal = Total richness difference diversity (or nestedness)/Total beta diversity.
note: Name of the dissimilarity coefficient.
The Jaccard and S<U+00F8>rensen dissimilarity coefficients and their quantitative forms, the Ruzicka and percentage difference indices, all have upper bounds (Dmax) of 1. Hence, when all sites contain a different set of species with no species in common, the maximum value that BDtotal can take is 0.5. See Legendre & De Caceres (2013, p. 958), section Maximum value of BD. This differs form the values produced by function beta.div(): with methods "hellinger", "chord" and "profiles", which have maximum values of sqrt(2), BDtotal has a maximum value of 1 for these dissimilarities.
For species presence-absence data, the dissimilarity coefficients are Jaccard = (b+c)/(a+b+c) and S<U+00F8>rensen = (b+c)/(2*a+b+c) with the usual a,b,c notation. For species abundance data, the dissimilarity coefficients are the Ruzicka index = (B+C)/(A+B+C) and Odum<U+2019>s percentage difference = (B+C)/(2A+B+C) (aka Bray-Curtis in some packages), where
B = sum of abundances at site 1 minus A,
C = sum of abundances at site 2 minus A.
The binary
(quant=FALSE) and quantitative (quant=TRUE) forms of the S and
J indices return the same values when computed for presence-absence data.
Baselga, A. (2010) Partitioning the turnover and nestedness components of beta diversity. Global Ecology and Biogeography, 19, 134<U+2013>143.
Baselga, A. (2012) The relationship between species replacement, dissimilarity derived from nestedness, and nestedness. Global Ecology and Biogeography, 21, 1223<U+2013>1232.
Baselga, A. (2013) Separating the two components of abundance-based dissimilarity: balanced changes in abundance vs. abundance gradients. Methods in Ecology and Evolution, 4, 552<U+2013>557.
Carvalho, J.C., Cardoso, P., Borges, P.A.V., Schmera, D. & Podani, J. (2013) Measuring fractions of beta diversity and their relationships to nestedness: a theoretical and empirical comparison of novel approaches. Oikos, 122, 825<U+2013>834.
Legendre, P. 2014. Interpreting the replacement and richness difference components of beta diversity. Global Ecology and Biogeography, 23, 1324-1334.
Legendre, P. and M. De C<U+00E1>ceres. 2013. Beta diversity as the variance of community data: dissimilarity coefficients and partitioning. Ecology Letters 16: 951-963.
Podani, J., Ricotta, C. & Schmera, D. (2013) A general framework for analyzing beta diversity, nestedness and related community-level phenomena based on abundance data. Ecological Complexity, 15, 52-61.
Podani, J. & Schmera, D. 2011. A new conceptual and methodological framework for exploring and explaining pattern in presence-absence data. Oikos, 120, 1625<U+2013>1638.
# NOT RUN {
if(require(ade4, quietly = TRUE)){
data(doubs)
fish.sp = doubs$fish[-8,] # Fish data; site 8 is removed because no fish were caught
# Compute and partition a matrix of Jaccard indices (presence-absence data)
out1 = beta.div.comp(fish.sp, coef="J", quant=FALSE)
out1$part
# Compute and partition a matrix of percentage difference indices
# (quantitative form of Sorensen index)
out2 = beta.div.comp(fish.sp, coef="S", quant=TRUE)
out2$part
# In paragraph Value, see the description of the 5 elements of vector part.
# Is the fish beta diversity dominated by replacement or richness/abundance difference?
}
# }
Run the code above in your browser using DataLab