This is a probability index in which rare species (or descriptors) are overweighted. Accordingly, sharing rare species makes pairs of samples (or objects) more similar than sharing common/frequent species (or descriptors). As Legendre & Legendre (1998) put it when describing this index: "... it is less likely for two sites to both contain the same rare species than a more frequent species. In this sense, agreement for a rare species should be given more importance than for a frequent species, when estimating the similarity between sites.".
dis.goodall(comm, p.simi="steinhaus", approach="proportion")
Dataframe or matrix with samples in rows and species in columns.
The partial similarity index to be used in the Goodall index. Partial match to "steinhaus" (raw abundance data), or "gower" (normalized abundance data).
The two approaches to compute Goodall index. Partial match to "proportion" or "chisquare".
Communities are usually composed of many rare and only a few common/frequent species. Although rare, a species may provide a valuable information in the estimation of resemblance between samples. However, the use of raw abundances makes the contribution of rare species to be negligible to the resulting index value. For instance, dropping a common species from the community data table usually will make much larger differences in the dissimilarity matrix than dropping a rare species.
There are a few indices that give more importance to an individual belonging to a rare than to a common/frequent species (see
dis.nness). Notice, however, that this differential weight of species can also be obtained, for instance, by log-transforming or standardizing data (e.g. dividing by maximum within each species) (Melo in preparation). In this sense, an extreme case in which rare and common species have the same weight is in the use of presence/absence data.
The implementation of Goodall index (Goodall 1966) follows Legendre & Legendre (1998, pp. 269-273). They suggest to use the "steinhaus" index for raw species abundance data and the "gower" partial similarity index (
p.simi) for normalized data. The index can be calculated in two ways or approaches ("proportion" or "chisquare") and values produced by them should be very different, although this is mostly due to scale; they are monotonically correlated. Notice Legendre & Legendre (1998) present a similarity version of this index. The one produced by this function is a dissimilarity, computed simply as 1-similarity. See the examples below for some behaviors of the index.
One important issue made clear by Legendre & Legendre (1998, pp. 270-271) is that the dissimilarity value for a pair of sites depnds on other sites present in the data matrix. If changes are made to some sites, the resulting dissimilarity values for all remaining pairs are modified (see example below).
A "dist" object.
Adriano Sanches Melo
Goodall, D.W. 1966. A new similarity index based on probability. Biometrics 22: 882-907.
Legendre, P & L. Legendre. 1998. Numerical Ecology. 2nd ed. Elsevier.
Melo, A.S. (in preparation) Is it possible to improve recovery of multivariate groups by weighting rare species in similarity indices?
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21
library(vegan) a <- c(1, 1, 0) b <- c(2, 1, 0) c <- c(0, 1, 1) d <- c(0, 1, 2) e <- c(0, 1, 3) dat5 <- rbind(a,b,c,d,e) colnames(dat5) <- c("sp1","sp2","sp3") dat5 # Notice the samples in the pair a-b differ from each other exactly in the same # way as samples in the pair c-d. However, a-b shares a rare species (sp1), # whereas c-d shares a frequent species (sp3, which is also present in e). Thus, # the dissimilarity a-b is the same of c-d using Bray-Curtis, but not using # Goodall index: vegdist(dat5, "bray") dis.goodall(dat5) # As the importance of a species for the Goodall index depends on its overall # frequency in the community data, the deletion of a sample changes results: dis.goodall(dat5[-5,])
Loading required package: vegan Loading required package: permute Loading required package: lattice This is vegan 2.5-7 Loading required package: ape Loading required package: picante Loading required package: nlme Loading required package: adespatial Registered S3 methods overwritten by 'adegraphics': method from biplot.dudi ade4 kplot.foucart ade4 kplot.mcoa ade4 kplot.mfa ade4 kplot.pta ade4 kplot.sepan ade4 kplot.statis ade4 scatter.coa ade4 scatter.dudi ade4 scatter.nipals ade4 scatter.pco ade4 score.acm ade4 score.mix ade4 score.pca ade4 screeplot.dudi ade4 Registered S3 method overwritten by 'spdep': method from plot.mst ape Registered S3 methods overwritten by 'adespatial': method from plot.multispati adegraphics print.multispati ade4 summary.multispati ade4 Loading required package: betapart sp1 sp2 sp3 a 1 1 0 b 2 1 0 c 0 1 1 d 0 1 2 e 0 1 3 a b c d b 0.2000000 c 0.5000000 0.6000000 d 0.6000000 0.6666667 0.2000000 e 0.6666667 0.7142857 0.3333333 0.1428571 a b c d b 0.0 c 0.4 0.4 d 0.4 0.4 0.2 e 0.4 0.4 0.3 0.0 a b c b 0.0000000 c 0.3333333 0.3333333 d 0.3333333 0.3333333 0.0000000
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.