[R-sig-eco] envfit function: is it possible to display correlation coefficient?

Irene Adamo i.adamo90 at gmail.com
Tue May 1 15:50:09 CEST 2018


 Dear all,
I have already looked at similar questions about this topic but I still
have few doubts.
when I run:
data(varespec, varechem)
library(MASS)
ord <- metaMDS(varespec)
(fit <- envfit(ord, varechem, perm = 999))

I get this output:
***VECTORS

            NMDS1    NMDS2     r2 Pr(>r)
N        -0.05738 -0.99835 0.2536  0.041 *
P         0.61980  0.78476 0.1938  0.105
K         0.76650  0.64225 0.1810  0.098 .
Ca        0.68524  0.72832 0.4119  0.001 ***

are the numbers under NMDS1-2  the direction of the arrow or they show
to which axis is the variable correlated? if I look at N can I conclude
that N is negatively correlated to NMDS2, while K is positively correlated
to NMDS1? is it possible to get the correlation coefficient of a variable
with the axes? here we have r2 but I don’t understand if it refers to both
of the axis or just one.

thanks a lot for any help!

	[[alternative HTML version deleted]]



More information about the R-sig-ecology mailing list