某所に掲載されていたプログラム
以下のようにするとコンパクト
KrusRep <- function(dat) {
x <- ncol(dat)
mat <- matrix(sapply(dat[, -x],
function(z) kruskal.test(z ~ dat[, x])$statistic), ncol=1)
dimnames(mat) <- list(colnames(dat[, -x]), "chi-squared")
return(mat)
}