0

Estoy intentando hacer un gráfico triangular a partir de los valores de beta diversidad, de manera que me quede en uno un gráfico triangular para una familia con puntos de colores para cada columna y en otro lo mismo para la otra familia. Ya sea usando la función ggplot, ggtern u otra. Pudieran ayudarme por favor o darme alguna sugerencia. Gracias de antemano

Adjunto parte de la matriz

Regiones    βjtu    βjne    βjac
I-II        0.85    0.05    0.90
I-III       0.78    0.02    0.80
I-IV        0.44    0.24    0.68
I-V         0.77    0.08    0.85
I-VI        0.85    0.02    0.85
I-VII       0.85    0.09    0.81
II-III      0.88    0.03    0.91
II-IV       0.51    0.32    0.85
II-V        0.88    0.04    0.89
II-VI       0.88    0.02    0.91
II-VII      0.91    0.02    0.91

Y el código que estoy probando

library (vegan)
library (betapart)
library (ggtern) 
## Raw data and plotting data(nn) 
m <- betadiver(nn) 
plot(m,col=2) 

## The indices betadiver(help=TRUE) 
## The basic Whittaker index 
d <- betadiver(sipoo, "w") 
## This should be equal to Sorensen index (binary Bray-Curtis in ## vegan) 
range(d - vegdist(nn, binary=TRUE))

Pero no se como darle la escala de colores a cada columna, ni como ajustar la escala que convencionalmente viene entre o y 100, y mis valores están entre o y 1.

3
  • Bienvenido, te sugiero que hagas el recorrido de bienvenida y veas Cómo preguntar para poder mejorar tu pregunta y obtener buenas respuestas. Por favor, edita tu pregunta para agregar el código de lo que has intentado hasta ahora y eventualmente los problemas con que te has encontrado. Saludos Commented el 9 may. 2018 a las 19:00
  • Soy nuevo en R, basicamente lo que he realizado es esto: library (vegan) library (betapart) library (ggtern) ## Raw data and plotting data(nn) m <- betadiver(nn) plot(m,col=2) ## The indices betadiver(help=TRUE) ## The basic Whittaker index d <- betadiver(sipoo, "w") ## This should be equal to Sorensen index (binary Bray-Curtis in ## vegan) range(d - vegdist(nn, binary=TRUE))
    – Asiel
    Commented el 9 may. 2018 a las 19:13
  • Pero no se como darle la escala de colores a cada columna, ni como ajustar la escala que convencionalmente viene entre o y 100, y mis valores están entre o y 1
    – Asiel
    Commented el 9 may. 2018 a las 19:13

1 respuesta 1

0

Aclaraciones

Tomo el data.frame de la pregunta como los datos finales a ser gratificados. Dos aclaraciones:

  • No sé si sobre esos datos estás aplicando otras transformaciones posteriormente. De todas formas el código debería funcionar siempre que tengas tres columnas con valores que van de 0 a 1.

  • En la pregunta dices que es una matriz. Si es así en sentido de R (is.matrix(x)) deberías pasarlo primero a data.frame, de lo contrario ggplot2 y ggtern no van a aceptar el input.

Código

    library(tidyverse)   #Para tribble(), puede no ser necesario si ya tienes los datos.
    library(ggtern)      #Para `coord_tern()`
    tribble(~Regiones,       ~βjtu,       ~βjne,     ~βjac,
           "I-II"  ,        0.85,        0.05,      0.90,
           "I-III" ,        0.78,        0.02,      0.80,
           "I-IV"  ,        0.44,        0.24,      0.68,
           "I-V"   ,        0.77,        0.08,      0.85,
           "I-VI"  ,        0.85,        0.02,      0.85,
           "I-VII" ,        0.85,        0.09,      0.81,
           "II-III",        0.88,        0.03,      0.91,
           "II-IV" ,        0.51,        0.32,      0.85,
           "II-V"  ,        0.88,        0.04,      0.89,
           "II-VI" ,        0.88,        0.02,      0.91,
           "II-VII",        0.91,        0.02,      0.91) -> datos_ternarios

Con esto dejo preparados los datos.

    ggplot(datos_ternarios, aes(x=βjtu, y=βjne, z=βjac, color= Regiones)) + 
    geom_point() +
    coord_tern() 

Con color = Regiones le indico que haga una escala de colores para cada valor único de esa variable. Automáticamente hace la leyenda. También es automático el escalado a 100, aunque tus datos crudos vayan entre 0 y 1.

Resultado

introducir la descripción de la imagen aquí

Este el gráfico básico con los valores por defecto, puedes ajustar más cosas como las líneas de guía, color del fondo, etc. Te recomiendo visitar http://www.ggtern.com para ver opciones de ajustes. Cuando se usa coord_tern() algunas cosas funcionan un poco diferente que en los gráficos de ggplot con coordenadas cartesianas.

Cualquier otra duda puedes ponerla en comentarios o actualizar tu pregunta.

PD: si fueras a usar el gráfico en una publicación el autor de ggtern solicita que incluyas el paquete en la bibliografía de tu paper. Considerando que el software que ofrece es libre y que eso podría ayudar a su carrera creo que sería de buen gusto hacerlo. No es obligatorio.

11
  • Me prece bien, si los graficos están previamente tratados con la funcion betapart del vegan, me parece bien pero la idea es graficar los valores de cada columna, independientemente de la leyenda de las regiones, o sea algo como esto:
    – Asiel
    Commented el 10 may. 2018 a las 17:02
  • library (vegan) library (betapart) library (ggtern) mi xlsx = data.frame(bjtu = runif(100), bjne = runif(100), bjac = runif(100), beta_diversidad = runif(100)) ggtern(data = mi xlsx, aes(x=bjtu, y=bjne, z=bjac)) + Tlab('bjtu') + Llab('bjne')+ Rlab(expression(paste(bjac))) + geom_point(aes(colour = beta diversidad)) + scale_colour_gradientn(colours = rev(terrain.colors(100)))
    – Asiel
    Commented el 10 may. 2018 a las 17:04
  • la cuestion es que va de 0-100, esa ecala no se puede poner a favor de los datos??, porque de otra forma no se ve casi la disperción de los datos para cada componente.
    – Asiel
    Commented el 10 may. 2018 a las 17:05
  • Para explicarme la idea es hacer un gráfico, con un color para cada región, pero que arriba del punto esté el valor de cada columna para cada punto.
    – Asiel
    Commented el 10 may. 2018 a las 17:26
  • Necesito algo com esto, lo tomé de ggtern.com #Variabels set.seed(1) n = 10; #Number of Data Points nv = 0.05 #Vertical Adjustment #Define Data df = data.frame(x = runif(n), y = runif(n), z = runif(n), label = sprintf("L%s",1:n), color = factor(sample(1:3,n,replace=TRUE)))
    – Asiel
    Commented el 10 may. 2018 a las 17:27

Tu Respuesta

By clicking “Publica tu respuesta”, you agree to our terms of service and acknowledge you have read our privacy policy.

¿No es la respuesta que buscas? Examina otras preguntas con la etiqueta o formula tu propia pregunta.