「統計学関連なんでもあり」の過去ログ--- 043

No.12466　RでKruskal-Wallis の検定　　【Kitamura】　2010/04/16(Fri) 12:44

RでKruskal-Wallis の検定をしたいのですが，各群のサンプル数が異なると警告「長さが異なります」が出て，結果を導くことができません。
どうすれば解決できますでしょうか。

お手数をおかけ致しますがよろしくお願い致します。

No.12467　Re: RでKruskal-Wallis の検定　　【青木繁伸】　2010/04/16(Fri) 14:13

「長さが異なります」というのは指定が間違えているからでしょう。正確には，「'x' と 'g' は同じ長さを持たねばなりません」といわれるのでは？
> kruskal.test(x=c(1,2,4,3,   2,1,3,   5,4,4,5,6),  g=c(1,1,1,1,   2,2,2,   3,3,3,3))
 以下にエラー kruskal.test.default(x = c(1, 2, 4, 3, 2, 1, 3, 5, 4, 4, 5, 6),  : 
   'x' と 'g' は同じ長さを持たねばなりません 
この例では，文字通り x, g の要素数（長さ）が違います。

当然，サンプルサイズが違っても，できます。先ほどの例は，実は以下の間違いだったとか？
> kruskal.test(x=c(1,2,4,3,   2,1,3,   5,4,4,5,6),  g=c(1,1,1,1,   2,2,2,   3,3,3,3,3))

	Kruskal-Wallis rank sum test

data:  c(1, 2, 4, 3, 2, 1, 3, 5, 4, 4, 5, 6) and c(1, 1, 1, 1, 2, 2, 2, 3, 3, 3, 3, 3) 
Kruskal-Wallis chi-squared = 7.5585, df = 2, p-value = 0.02284

> kruskal.test(x=c(1,2,4,3,   2,1,3,   5,4,4,5,6),  g=rep(1:3, c(4,3,5)))

	Kruskal-Wallis rank sum test

data:  c(1, 2, 4, 3, 2, 1, 3, 5, 4, 4, 5, 6) and rep(1:3, c(4, 3, 5)) 
Kruskal-Wallis chi-squared = 7.5585, df = 2, p-value = 0.02284
「そんなことじゃない」ということでしたら，あなたがやったことが分かるようにコンソールからコピーして提示してください。