The standard stats::kruskal.test module allows to calculate the kruskal-wallis test on a dataset: >>> data(diamonds)