The standard stats::kruskal.test module allows to calculate the kruskal-wallis test on a dataset:
>>> data(diamonds) >>> kruskal.test(price
I would have thought you'd be able to do the following:
data(diamonds, package = "ggplot2") library(coin) library(multcomp) kt <- kruskal_test(price ~ clarity, data = diamonds) glht(kt, mcp(clarity = "Tukey"))
But it seems like multcomp doesn't support coin objects (yet?).
multcomp
coin