Similarity Percentages
Discriminating species between two groups using Bray-Curtis dissimilarities
simper(comm, group, permutations = 0, trace = FALSE,
parallel = getOption("mc.cores"), ...)
## S3 method for class 'simper'
summary(object, ordered = TRUE,
digits = max(3,getOption("digits") - 3), ...)comm |
Community data matrix. |
group |
Factor describing the group structure. Must have at least 2 levels. |
permutations |
a list of control values for the permutations
as returned by the function |
trace |
Trace permutations. |
object |
an object returned by |
ordered |
Logical; Should the species be ordered by their average contribution? |
digits |
Number of digits in output. |
parallel |
Number of parallel processes or a predefined socket
cluster. With |
... |
Parameters passed to other functions. In |
Similarity percentage, simper (Clarke 1993) is based
on the decomposition of Bray-Curtis dissimilarity index (see
vegdist, designdist). The contribution
of individual species i to the overall Bray-Curtis dissimilarity
d[jk] is given by
d[ijk] = abs(x[ij]-x[ik])/sum(x[ij]+x[ik])
where x is the abundance of species i in sampling units j and k. The overall index is the sum of the individual contributions over all S species d[jk] = sum(i=1..S) d[ijk].
The simper functions performs pairwise comparisons of groups
of sampling units and finds the contribution of each species to the
average between-group Bray-Curtis dissimilarity. Although the method
is called simper, it really studied dissimilarities instead of
similarities (Clarke 1993).
The function displays most important species for each pair of
groups. These species contribute at least to 70 % of the
differences between groups. The function returns much more
extensive results (including all species) which can be accessed
directly from the result object (see section Value). Function
summary transforms the result to a list of data frames. With
argument ordered = TRUE the data frames also include the
cumulative contributions and are ordered by species contribution.
The results of simper can be very difficult to interpret and
they are often misunderstood even in publications. The method gives
the contribution of each species to overall dissimilarities, but
these are caused by variation in species abundances, and only partly
by differences among groups. Even if you make groups that are
copies of each other, the method will single out species with high
contribution, but these are not contributions to non-existing
between-group differences but to random noise variation in species
abundances. The most abundant species usually have highest
variances, and they have high contributions even when they do not
differ among groups. Permutation tests study the differences among
groups, and they can be used to find out the species for which the
differences among groups is an important component of their
contribution to dissimilarities.
A list of class "simper" with following items:
species |
The species names. |
average |
Species contribution to average between-group dissimilarity. |
overall |
The average between-group dissimilarity. This is the sum of
the item |
sd |
Standard deviation of contribution. |
ratio |
Average to sd ratio. |
ava, avb |
Average abundances per group. |
ord |
An index vector to order vectors by their contribution or
order |
cusum |
Ordered cumulative contribution. These are based on item
|
p |
Permutation p-value. Probability of getting a larger
or equal average contribution in random permutation of the group
factor. These area only available if |
Eduard Szöcs eduardszoecs@gmail.com
Clarke, K.R. 1993. Non-parametric multivariate analyses of changes in community structure. Australian Journal of Ecology, 18, 117–143.
Function meandist shows the average between-group
dissimilarities (as well as the within-group dissimilarities).
data(dune) data(dune.env) (sim <- with(dune.env, simper(dune, Management))) summary(sim)
Please choose more modern alternatives, such as Google Chrome or Mozilla Firefox.