例題:
「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 <- matrix(c(
13, 1, 10, 1, 12, 1, 19, 1,
21, 2, 26, 2, 15, 2, 14, 2, 21, 2,
27, 3, 28, 3, 21, 3,
13, 4, 16, 4, 19, 4, 10, 4, 12, 4, 19, 4
), ncol=2, byrow=T)
前もって用意されている関数を使う場合
> kruskal.test(dat[,1], dat[,2])
Kruskal-Wallis rank sum test
data: dat[, 1] and dat[, 2]
Kruskal-Wallis chi-squared = 9.7471, df = 3, p-value = 0.02084
新たに定義した関数(多重比較も行う)を使う場合
> kruskal.wallis(dat, 1, 2)
$Result
Statistics(Chi-sq. value) d.f. P value
9.74705115 3.00000000 0.02084332
$Statistics
Group 1 Group 2 Group 3 Group 4
Group 1 0.0000000 3.255721 7.306100 0.2024443
Group 2 3.2557208 0.000000 1.367571 2.3082068
Group 3 7.3061003 1.367571 0.000000 6.2941603
Group 4 0.2024443 2.308207 6.294160 0.0000000
$P.value
Group 1 Group 2 Group 3 Group 4
Group 1 1.00000000 0.3538532 0.06275556 0.97719374
Group 2 0.35385318 1.0000000 0.71315348 0.51095052
Group 3 0.06275556 0.7131535 1.00000000 0.09814354
Group 4 0.97719374 0.5109505 0.09814354 1.00000000