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!!!.