例題:
「12 匹のラットに 3 種類の餌を与えたときの肝臓の重量は表 1 のようであった。餌の種類により肝臓の重量の平均値に差があるといえいるか,有意水準 5% で検定しなさい。また,多重比較を行いなさい。」
A餌 | 3.42 | 3.84 | 3.96 | 3.76 | |
---|---|---|---|---|---|
B餌 | 3.17 | 3.63 | 3.47 | 3.44 | 3.39 |
C餌 | 3.64 | 3.72 | 3.91 |
R による解析:
> dat <- c(3.42, 3.84, 3.96, 3.76, 3.17, 3.63, 3.47, 3.44, 3.39, 3.64, 3.72, 3.91) > group <- rep(1:3, c(4, 5, 3)) > kruskal.wallis(dat, group) # この関数の定義を見る クラスカル・ウォリス検定(plus 多重比較) data: dat ~ group Kruskal-Wallis chi-squared = 5.5, df = 2, p-value = 0.06 多重比較の結果 chi sq. p-value 1:2 4.104274 0.1285 1:3 0.003663 0.9982 2:3 3.702564 0.1570 データをリストで与えることもできる。 > kruskal.wallis(list(c(3.42, 3.84, 3.96, 3.76), c(3.17, 3.63, 3.47, 3.44, 3.39), c(3.64, 3.72, 3.91))) クラスカル・ウォリス検定(plus 多重比較) data: list(c(3.42, 3.84, 3.96, 3.76), c(3.17, 3.63, 3.47, 3.44, 3.39), c(3.64, 3.72, 3.91)) Kruskal-Wallis chi-squared = 5.5, df = 2, p-value = 0.06 多重比較の結果 chi sq. p-value 1:2 4.104274 0.1285 1:3 0.003663 0.9982 2:3 3.702564 0.1570