1

Necesito que los cuatros gráficos que salen con el código que les dejo mas abajo salgan en uno solo ajustando la escala lo mejor posible sin perder información y que cada curva sea de un color poniéndole una etiqueta con el nombre del estimado:

Este es el código original

library (vegan)
data(dune)
data(dune.env)
pool <- with(dune.env, specpool(dune, Management))
pool
op <- par(mfrow=c(1,2))
boxplot(specnumber(dune) ~ Management, data = dune.env,
col = "hotpink", border = "cyan3")
boxplot(specnumber(dune)/specpool2vect(pool) ~ Management,
data = dune.env, col = "hotpink", border = "cyan3")
par(op)
data(BCI)
## Accumulation model
pool <- poolaccum(BCI)
summary(pool, display = "chao")
plot(pool)
## Quantitative model
estimateR(BCI[1:5,])

Esto es un ejemplo reproducible, asi que no deberán tener problema a la hora de correr el código!!!.

4

1 respuesta 1

2

Para ver todas las curvas en un solo gráfico, creo que ya nos empieza a convenir pasar todo a ggplot.

library(tidyverse)
library (vegan)

data(BCI)
pool <- poolaccum(BCI)

df <- summary(pool) %>% 
    map(as.data.frame) %>%   # Cada matriz de la lista a un data.frame
    map(rename, V = 2) %>%   # Renombramos la 2 columna para que el nombre sea consistente
    bind_rows(.id = "index") # juntamos todos los df en uno solo

ggplot(df, aes(x=N, color=index)) +
    xlab("Size") +
    ylab("Richeness") +
    geom_line(aes(y=V), size=1) +
    geom_line(aes(y=`2.5%`), size=.5, linetype=4) +
    geom_line(aes(y=`97.5%`), size=.5, linetype=4) +
    theme_minimal()

El punto clave es transformar el objeto poolacum con un poco de gramática de purrr y dplyr en un único data.frame con todos lo índices, para luego explotarlo mediante ggplot

introducir la descripción de la imagen aquí

19
  • Me parece genial la idea aunque me referia a la grafica que dice riqueza poder tener esas 5 curvas en uno solo con la misma escala es posible esto? Commented el 11 feb. 2019 a las 14:17
  • 1
    Para quitar las grillas, agregale al final + theme(panel.grid.major = element_blank(), panel.grid.minor = element_blank()) . Lo otro, no lo sé, no conozco realmente el paquete vegan. Saludos. Commented el 11 feb. 2019 a las 18:07
  • 1
    Lo más simple es pasarle a ggplot un data.frame filtrado sin el índice que deseas, es decir: ggplot(df[df$index!='jack1',], aes(x=N, color=index)) + .... Commented el 12 feb. 2019 a las 22:46
  • 1
    @AsielCabreraGuerrero, puedes agregar, por ejemplo + scale_x_continuous(limits = c(20, 40)) Commented el 15 feb. 2019 a las 13:59
  • 1
    @AsielCabreraGuerrero, revisa que 97.5% debe estar encerrado por "backticks", fijate el código original de mi respuesta Commented el 15 feb. 2019 a las 20:03

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.