裏 RjpWiki

Julia ときどき R, Python によるコンピュータプログラム,コンピュータ・サイエンス,統計学

ダメ出し:なぜ oneway.test を使わないのだろうか?

2012年11月30日 | ブログラミング

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

コメント    この記事についてブログを書く
  • X
  • Facebookでシェアする
  • はてなブックマークに追加する
  • LINEでシェアする
« ダメ出し:par(new=TRUE) は... | トップ | ダメ出し:データに NA が含... »
最新の画像もっと見る

コメントを投稿

ブログラミング」カテゴリの最新記事