library("MASS")
# check variances
tapply(Boston$medv, Boston$rad, var)
# check skewnesss
library("e1071")
tapply(Boston$medv, Boston$rad, skewness)
# Kruskal-Wallis H test :(
kruskal.test(medv~rad, data=Boston)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.