Re: how to calculate "axis variance" in metaMDS, pakage vegan?

Dear all,
What about a r2 statistic in this case? I have seen some people using the r2
of the relationship (a matrix regression) between te original distance and
the final distance obtained with the NMDS.
This r2 is not similar to the stress measure, as we can see above:

library(vegan)
data(varespec)
original.dist<-vegdist(varespec)
r2<-numeric(10)
stress<-numeric(10)
for(i in 1:10){
nmds.resu<-metaMDS(varespec,k=i)
nmds.scores<-scores(nmds.resu)
nmds.dist<-dist(nmds.scores)
r2[i]<-summary(lm(original.dist~nmds.dist))[[8]]
stress[i]<-nmds.resu$stress
}
r2     ## r2 initially increase and as more axis are added it begins to
decrease.
stress
plot(r2,stress)

I have read somewhere something like "dont stress about stress", however I
really dont know what is best, stress or r2. People say that this r2 is
something like "how good you recovered the original distances with the NMDS
ordination".
What about this r2 measure?

victor