クラスカル・ウォリス検定     Last modified: Dec 03, 2001

例題

 「4 つの群についてある測定をおこなったところ,表 3 のような結果が得られた。代表値に差があるか有意水準 5% で検定しなさい。また,多重比較を行いなさい。」

表 3.測定値
群 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


・ 手法の解説ページ
・ 直前のページへ戻る  ・ E-mail to Shigenobu AOKI