library("vegan")
library("ggplot2")
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