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

No.11036　自由度とＦ値からのP値の計算　　【k2w】　2009/10/07(Wed) 23:57

色々と調べたのですが，どうしても分からないので質問させてください。

2つの自由度と1つのF値からF検定のp値を計算したいのですが，エクセルでもRでも，そのような関数がみつかりません。自由度が10000近い値ですので，統計表を見ても載っていません。

どなたか，どのようにしたら算出できるかを教えていただけるとありがたいです。よろしくお願いします。

No.11037　Re: 自由度とF値からのP値の計算　　【青木繁伸】　2009/10/08(Thu) 00:10

なぜ正解にたどり着けなかったのかが不思議。。。
FDist {stats}	R Documentation
The F Distribution

Description

Density, distribution function, quantile function and random generation for
the F distribution with df1 and df2 degrees of freedom (and optional
non-centrality parameter ncp).

Usage

df(x, df1, df2, ncp, log = FALSE)
pf(q, df1, df2, ncp, lower.tail = TRUE, log.p = FALSE)
qf(p, df1, df2, ncp, lower.tail = TRUE, log.p = FALSE)
rf(n, df1, df2, ncp)
実行例
> pf(1.5, 30, 10000, lower.tail=FALSE)
[1] 0.03891821
エクセルにだってあります
=FDIST(1.5,30,10000)
0.038918211
これでいいんでしょうか？

No.11038　Re: 自由度とF値からのP値の計算　　【k2w】　2009/10/08(Thu) 09:03

青木先生

まさにそれが探していた関数です。どうもありがとうございました！

Rでhelp.search("F test")とすると，いちばん近いものがvar.testだったり，エクセルでF検定とすると，分析ツールを使ったものしか無く，しかもF値からP値を求めるものではなかったりしました．．．。

探すべきものは，F testでは無くて，F distributionだったのですね。よく考えてみれば，その通りでした。大変参考になりました。ありがとうございました。

No.11039　Re: 自由度とF値からのP値の計算　　【にゃんちゅう】　2009/10/08(Thu) 17:46

>自由度が10000近い値ですので，統計表を見ても載っていません。

統計表を見るときは∞を見ます。1000以上なら∞で問題はないはずです。