kosugitti
> result.aov <- aov(height ~ class, data = sample)
> summary(result.aov)
Df Sum Sq Mean Sq F value Pr(>F)
class 2 116 58.21 0.706 0.496
Residuals 97 8000 82.47
何故こんなに回りくどいことをやるのだろうか?
しかも,これは「各群の分散が等しい」ことを仮定している。
t.test を使うなら,oneway.test を使うべきであろう。これは,「各群の分散が等しいことを仮定しない」がデフォルトになっている(t.test と平仄が合っている)。
> oneway.test(height ~ class, data=sample)
One-way analysis of means (not assuming equal variances)
data: height and class
F = 0.6541, num df = 2.000, denom df = 64.509, p-value =
0.5233
> oneway.test(height ~ class, data=sample, var.equal=TRUE)
One-way analysis of means
data: height and class
F = 0.7058, num df = 2, denom df = 97, p-value = 0.4962
※コメント投稿者のブログIDはブログ作成者のみに通知されます