Re: [R] Vegan R^2 and tau values for metaMDS

From: Sarah Goslee <>
Date: Mon, 14 Apr 2008 17:00:31 -0400

On Mon, Apr 14, 2008 at 4:47 PM, Gavin Simpson <> wrote:

> Note that the default is to produce a bray-curtis dissimilarity matrix
> from the input species data. As such, I reproduce this dissimilarity
> matrix as arg 1 to cor and then take the Euclidean distances of the
> coordinates on nMDS axes 1:2 (this example is a 2D solution but scales
> to n-dimensions) as the second argument to cor:
> cor(vegdist(dune), dist(sol$points))^2
> In this example, the "R^2" is 0.899222
> If you already have jaccard for your species data, then substitute this
> for 'vegdist(dune)' and sol for whatever your nMDS object is called.

Good point. I forgot to clarify that you should use the same dissimilarity metric for your original data that was used in constructing the ordination. Euclidean distances should be used for the NMDS configuration, though.

> However, this assumes a linear relationship between the original
> distances and the nMDS distances. The relationship need not be linear,
> just monotonic if I recall the details correctly.

If you (the original querent) are interested in the monotonic rather than strictly linear relationship, you could use the "spearman" option to cor().


Sarah Goslee

______________________________________________ mailing list
PLEASE do read the posting guide
and provide commented, minimal, self-contained, reproducible code.
Received on Mon 14 Apr 2008 - 21:08:26 GMT

Archive maintained by Robert King, hosted by the discipline of statistics at the University of Newcastle, Australia.
Archive generated by hypermail 2.2.0, at Tue 15 Apr 2008 - 17:30:28 GMT.

Mailing list information is available at Please read the posting guide before posting to the list.

list of date sections of archive