4

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?

[Seria algo assim, mas com os nomes das espécies]

Seguindo o exemplo dado pelo @Imonferrari, o gráfico fica assim.

Gráfico final seguindo o exemplo de Imonferrari

9
  • Como seria seu gráfico final? Você quer postar sp1/sp2... nas cordenadas de mds1 e mds2 em spp? Commented 8/09/2020 às 20:27
  • Quero que os nomes das espécies apareçam nos mesmos pontos dos 'sitios' 'a', 'b' e 'c'.
    – GEO
    Commented 8/09/2020 às 23:24
  • esse ultimo plot é o resultado que você espera ou o resultado que você está gerando? Commented 9/09/2020 às 0:20
  • É o que espero, é só um exemplo, esse fiz usando geom_text(aes(label="spp")).
    – GEO
    Commented 9/09/2020 às 0:54
  • Certo, mas olha as espécies seriam as sp1, sp2 né? É que posso te ajudar com o ggplot mas não conheço o tipo de dado, aí gostaria de entender mais de como você quer o gráfico final Commented 9/09/2020 às 1:03

1 Resposta 1

1

Você pode pivotar:

library(tidyr)
library(dplyr)

teste <- plt %>% pivot_longer(
  cols = starts_with('sp'),
  names_to = 'species',
  values_to = 'values'
)

E depois plotar:

ggplot(teste, aes(mds1, mds2, shape = sitio, color = sitio))+
  geom_point()+
  geom_text(label = teste$species)
6
  • Isso, pensei em algo assim, colocar todas as spp numa coluna, mas não achei o jeito. Agora, que faço porque no meu data frame todas as espécies começam com letras diferentes e tenho 55.
    – GEO
    Commented 9/09/2020 às 1:44
  • As colunas seguem uma ordem? Caso elas sigam você pode fazer um slicing, onde está colsdentro do pivot_longer você colocaria algo: cols = colnames(plt)[3:7] nesse nosso data.frame de exemplo retornará : '"sp1" "sp2" "sp3" "sp4" "sp5"'. Commented 9/09/2020 às 9:30
  • Esse exemplo responde a sua dúvida? Commented 9/09/2020 às 13:41
  • Responde sim, obrigado. Rodei com seu exemplo no o df original, mas o gráfico não fica bom. Vou coloca-lo para que você o veja.
    – GEO
    Commented 9/09/2020 às 17:06
  • Tem como disponibilizar o dataset que você está usando completo? O que você quer que melhore no gráfico? Commented 9/09/2020 às 17:28

Você deve fazer log-in para responder a esta pergunta.

Esta não é a resposta que você está procurando? Pesquise outras perguntas com a tag .