tidy vegan specaccum plot for ggplot2

tidy_specaccum.R

#' Convert vegan specaccum results to tidy format
#'
#' @param x result returned by specaccum
#'
#' @return A data.frame in "long format".
#'
tidy_specaccum <- function(x) {
    data.frame(
      site = x$sites,
      richness = x$richness,
      sd = x$sd)
}

library(vegan)
library(ggplot2)
library(magrittr)

data(BCI)

# Fit rarefaction curves with vegan
sac <- specaccum(BCI)

# Tidying
tidy_sac <- tidy_specaccum(sac)

# Plot with ggplot
tidy_sac %>% 
  ggplot() +
  geom_line(aes(site, richness)) +
  geom_linerange(aes(x = site, ymin = richness - 2*sd, ymax = richness + 2*sd)) +
  ylim(0, NA)