Dispersion-based weighting - a new function in vegan
When analysing multivariate community data the first step is often to transform the abundance data, in order to downweight contributions of high abundance taxa to dissimilarity. There are a lot of transformations available, for example:
- log(x + 1), but see O’Hara & Kotze (2010)
- ln(2x + 1), which is used often with mesocosm data
- and manymany more
However these are all global transformations (applied to every species in the dataset).
In an email discussion with Bob Clarke he pointed me to one of his papers Clarke et al. 2006:
Clarke, K. R., M. G. Chapman, P. J. Somerfield, and H. R. Needham. 2006. “Dispersion-based Weighting of Species Counts in Assemblage Analyses.” Marine Ecology. Progress Series 320: 11–27. klick
I somehow missed this paper and it looks like this hasn’t been implemented in R yet. So I wrote it up for myself and the function has been merged into vegan this week (see Changelog).
The new function is named
dispweight and takes three arguments:
The community data matrix, a vector describing the group structure and the number of permutations to use for the permutation test.
It returns a list of four elements:
- D: The dispersion index
- p: p-value of the permutation test
- weights: The weights applied to species
- transformed: The transformed abundance matrix
Here it take the dune dataset and
Management as grouping variable.
First we take look at the data via NMDS:
Lets run the NMDS on dispersion-weighted abundances:
In this example there is not a big difference, but as they write in their paper:
[…] we should not expect a dispersion-weighting approach to radically alter ordination and test results. Indeed, the improvements shown in the real examples of this paper are sometimes small and subtle […].
I haven’t tested the function very much, so let me know if you find any bugs! Note that the code isn’t optimised for speed and follows largely the paper.
- KR Clarke, MG Chapman, PJ Somerfield, HR Needham, (2006) Dispersion-Based Weighting of Species Counts in Assemblage Analyses. Marine Ecology Progress Series 320 11-27 10.3354/meps320011
- Robert B. O’Hara, D. Johan Kotze, (2010) do Not Log-Transform Count Data. Methods in Ecology And Evolution 1 118-122 10.1111/j.2041-210X.2010.00021.x