Hola Fernando:

Ya que no nos dices con qué función has realizado tu ordenación, la 
manera más sencilla de visualizar los grupos en el espacio NMDS sería 
realizar un diagrama de dispersión usando los factores como parámetros 
gráficos (color, pch, etc).

Por ejemplo:

plot(NMDS$1, NMDS$2, pch=as.numeric(habitat), col=as.numeric(distalidad))

También podrías calcular un convex hull para cada uno de tus grupos y 
representarlo en forma de polígonos sobre el diagrama de ordenación (ver 
función ordihull() en el paquete vegan. Un mini-tutorial sobre esto lo 
tienes en una de las viñetas de dicho programa (sección 2.2)

https://cran.r-project.org/web/packages/vegan/vignettes/intro-vegan.pdf


Un saludo,

Marcelino


El 24/10/2019 a las 15:22, Fernando Archuby escribió:
> Buenos días. Realizo el NMDS para ordenar muestras compuestas por especies
> biológicas a partir de abundancias, a partir de una matriz de distancias de
> diferencia porcentual (aka Bray Curtis). Hasta ahí, no encuentro problemas.
> Por otro lado, las muestras se clasifican de los siguientes modos:
> - distalidad al centro de la bahía (D,M,P)
> - vivas versus muertas, estas últimas identificadas a partir de las valvas
> (L,D)
> - hábitat contaminado versus hábitat no contaminado (P,U)
> De modo que tengo tres factores y en total 12 grupos.
>
> Quisiera visualizar apropiadamente la distribución de estos grupos en el
> ordenamiento multivariado en 2 dimensiones.
> ¿Alguna idea sobre cómo proceder? ¿Algún ejemplo o tutorial para guiarme?
> Muchas gracias,
> Fernando.
>

-- 
Marcelino de la Cruz Rot
Depto. de Biología y Geología
Física y Química Inorgánica
Universidad Rey Juan Carlos
Móstoles España

_______________________________________________
R-help-es mailing list
R-help-es@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-help-es

Responder a