例題:
「4 つの群についてある測定をおこなったところ,表 3 のような結果が得られた。代表値に差があるか有意水準 5% で検定しなさい。また,多重比較を行いなさい。」
群 1 | 13 | 10 | 12 | 19 | ||
---|---|---|---|---|---|---|
群 2 | 21 | 26 | 15 | 14 | 21 | |
群 3 | 27 | 28 | 21 | |||
群 4 | 13 | 16 | 19 | 10 | 12 | 19 |
R による解析:
> dat <- c(13, 10, 12, 19, 21, 26, 15, 14, 21, 27, 28, 21, 13, 16, 19, 10, 12, 19) > group <- rep(1:4, c(4, 5, 3, 6)) 前もって用意されている関数を使う場合 > kruskal.test(dat, group) Kruskal-Wallis rank sum test data: dat and group Kruskal-Wallis chi-squared = 9.7471, df = 3, p-value = 0.02084 新たに定義した関数(多重比較も行う)を使う場合 > kruskal.wallis(dat, group) クラスカル・ウォリス検定(plus 多重比較) data: dat ~ group Kruskal-Wallis chi-squared = 9.7, df = 3, p-value = 0.02 多重比較の結果 chi sq. p-value 1:2 3.2557 0.35385 1:3 7.3061 0.06276 1:4 0.2024 0.97719 2:3 1.3676 0.71315 2:4 2.3082 0.51095 3:4 6.2942 0.09814 または,データをリストで指定する。 > kruskal.wallis(list(c(13, 10, 12, 19), c(21, 26, 15, 14, 21), c(27, 28, 21), c(13, 16, 19, 10, 12, 19))) クラスカル・ウォリス検定(plus 多重比較) data: list(c(13, 10, 12, 19), c(21, 26, 15, 14, 21), c(27, 28, 21), c(13, 16, 19, 10, 12, 19)) Kruskal-Wallis chi-squared = 9.7, df = 3, p-value = 0.02 多重比較の結果 chi sq. p-value 1:2 3.2557 0.35385 1:3 7.3061 0.06276 1:4 0.2024 0.97719 2:3 1.3676 0.71315 2:4 2.3082 0.51095 3:4 6.2942 0.09814