Add Kruskal-Wallis test

There is also `tiedrank_adj` in `wilcoxon.jl`, which computes both the ranks and the adjustment for ties (`sum(T_i)` below) in a single pass.