Dado o seguinte data frame:
dat<-data.frame(
"sitio" = c("a", "a", "a",'a', "b", "b", "b", 'b', "c", "c", "c",'c'),
"amostra"=c(rep(1,4),rep(2,4),rep(3,4)),
"sp1"=c(1,3,3,2,4,2,1,5,3,6,1,5),
"sp2"=c(1,3,3,2,4,2,1,5,3,6,1,5),
"sp3"=c(1,3,3,2,4,2,1,5,3,6,1,5),
"sp4"=c(2,3,4,1,5,3,1,5,5,8,9,1),
"sp5"=c(3,4,3,1,6,7,5,8,3,1,3,2)
)
Quero ver o agrupamento das espécies de acordo ao sitio de coleta, plotando o resultado do metaMDS
, usando ggplot2.
library(vegan)
library(ggplot2)
dmds<-metaMDS(dat[,3:7], distance = "bray", autotransform = FALSE)
mds1 <- dmds$points[,1]
mds2 <- dmds$points[,2]
plt<-cbind(dat,mds1,mds2)
spp<-dmds$species
ggplot(plt,aes(mds1,mds2, shape=sitio,color=sitio))+
geom_point()+
geom_text(aes(label=spp))
Tentei usar geom_text, etiquetando as espécies do MDS, mas como tem um número diferente de linhas, dá o seguinte erro.
Error: Aesthetics must be either length 1 or the same as the data (12): label
Como poderia juntar o resultado do MDS com o data frame original para ter acesso as espécies e aos sitios, conseguindo visualizar tudo em um gráfico?
[]
Seguindo o exemplo dado pelo @Imonferrari, o gráfico fica assim.