data(Boston, package="MASS")
# ANOVA
fit <- aov(medv~rad, data=Boston)
# LSD test
library("agricolae")
LSD.test(fit, "rad", p.adj="bonferroni", console=T)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.