Learn R Programming

vegan (version 2.0-10)

designdist: Design your own Dissimilarities

Description

You can define your own dissimilarities using terms for shared and total quantities, number of rows and number of columns. The shared and total quantities can be binary, quadratic or minimum terms. In binary terms, the shared component is number of shared species, and totals are numbers of species on sites. The quadratic terms are cross-products and sums of squares, and minimum terms are sums of parallel minima and row totals.

Usage

designdist(x, method = "(A+B-2*J)/(A+B)",
           terms = c("binary", "quadratic", "minimum"), 
           abcd = FALSE, name)

Arguments

x
Input data.
method
Equation for your dissimilarities. This can use terms J for shared quantity, A and B for totals, N for the number of rows (sites) and P for the number of columns (species). The e
terms
How shared and total components are found. For vectors x and y the "quadratic" terms are J = sum(x*y), A = sum(x^2), B = sum(y^2), and "minimum" terms
abcd
Use 2x2 contingency table notation for binary data: $a$ is the number of shared species, $b$ and $c$ are the numbers of species occurring only one of the sites but not in both, and $d$ is the number of species that occur on neither of the sites.
name
The name you want to use for your index. The default is to combine the method equation and terms argument.

Value

  • designdist returns an object of class dist.

encoding

UTF-8

Details

Most popular dissimilarity measures in ecology can be expressed with the help of terms J, A and B, and some also involve matrix dimensions N and P. Some examples you can define in designdist are: lll{ A+B-2*J "quadratic" squared Euclidean A+B-2*J "minimum" Manhattan (A+B-2*J)/(A+B) "minimum" Bray-Curtis (A+B-2*J)/(A+B) "binary" Sørensen{Sorensen} (A+B-2*J)/(A+B-J) "binary" Jaccard (A+B-2*J)/(A+B-J) "minimum" Ružička{Ruzicka} (A+B-2*J)/(A+B-J) "quadratic" (dis)similarity ratio 1-J/sqrt(A*B) "binary" Ochiai 1-J/sqrt(A*B) "quadratic" cosine complement 1-phyper(J-1, A, P-A, B) "binary" Raup-Crick (but see raupcrick) }

The function designdist can implement most dissimilarity indices in vegdist or elsewhere, and it can also be used to implement many other indices, amongst them, most of those described in Legendre & Legendre (2012). It can also be used to implement all indices of beta diversity described in Koleff et al. (2003), but there also is a specific function betadiver for the purpose.

If you want to implement binary dissimilarities based on the 2x2 contingency table notation, you can set abcd = TRUE. In this notation a = J, b = A-J, c = B-J, d = P-A-B+J. This notation is often used instead of the more more tangible default notation for reasons that are opaque to me.

References

Koleff, P., Gaston, K.J. and Lennon, J.J. (2003) Measuring beta diversity for presence--absence data. J. Animal Ecol. 72, 367--382. Legendre, P. and Legendre, L. (2012) Numerical Ecology. 3rd English ed. Elsevier

See Also

vegdist, betadiver, dist.

Examples

Run this code
## Arrhenius dissimilarity: the value of z in the species-area model
## S = c*A^z when combining two sites of equal areas, where S is the
## number of species, A is the area, and c and z are model parameters.
## The A below is not the area (which cancels out), but number of
## species in one of the sites, as defined in designdist().
data(BCI)
dis <- designdist(BCI, "(log(A+B-J)-log(A+B)+log(2))/log(2)")
## This can be used in clustering or ordination...
ordiplot(cmdscale(dis))
## ... or in analysing beta diversity (without gradients)
summary(dis)

Run the code above in your browser using DataLab