# NOT RUN {
data(dune, dune.env)
dune.dist <- vegdist(dune)
dune.ano <- with(dune.env, anosim(dune.dist, Management))
df <- fortify(dune.ano)
ggplot(df, aes(x = Class, y = Rank)) +
geom_boxplot(notch = FALSE, varwidth = TRUE)
# }
Run the code above in your browser using DataLab