カイ二乗と自由度を指定してp値を計算するにはどうすればよいですか?たとえば、カイ二乗= 15、df = 2の正確なp値は何でしょうか?
カイ二乗と自由度を指定してp値を計算するにはどうすればよいですか?たとえば、カイ二乗= 15、df = 2の正確なp値は何でしょうか?
回答:
適用統計では、二乗検定統計は残差の2乗の合計として、または2乗効果の合計から、または対数尤度差から発生します。これらすべてのアプリケーションでの目的は、いくつかのベクトルパラメーターがゼロかどうかをテストすることであり、それが非ゼロであり、カイ二乗統計量は観測された効果の二乗サイズに関連しています。必要なp値はカイ二乗値の右裾確率であり、例ではRで次のようになります。
> pchisq(15, df=2, lower.tail=FALSE)
[1] 0.0005530844
他のdf値または統計値については、明らかにそれらを上記のコードに置き換えるだけです。
Rのすべての累積確率関数は、デフォルトで左テール確率を計算します。ただし、これらにもlower.tail
引数があり、常にこれFALSE
を設定して正しい確率を得ることができます。一部の小学校の教科書で見られるように、を計算するのではなく、これを行うことをお勧めします。
Rは、フォームの濃度又は質量の確率関数のスイートを有しているd*
(例えばdbeta
、dchisq
)、及び形態の分布p*
(例えばpf
、pgamma
)。そこから始めたいと思うかもしれません。
はい、与えられたp値(p)と自由度(df)のカイ2乗値を計算することが可能です。以下はその方法です:
検証のために、最初に、与えられたカイ2乗値= 1.1およびdf = 1のpを計算します。
pchisq(1.1、df = 1、lower.tail = FALSE)#答えはp = 0.2942661です
ここで、pとdfを使用してカイ二乗値を計算することにより、前に戻るために、上記で取得したp = 0.2942661と上記のdf = 1を使用しました。
qchisq(0.2942661、1、lower.tail = FALSE)#最初のソリューションと同様に、答えは1.1です。
したがって、カイ二乗= 15とdf = 2の例を使用すると、ソリューションは次のようになります。
pchisq(15、df = 2、lower.tail = FALSE)#答え:p = 0.0005530844
qchisq(0.0005530844、2、lower.tail = FALSE)#答え:カイ二乗= 15
お役に立てれば!!!
試して、
pchisq(chi,df)
あなたの例では、
pchisq(15,2)
[1] 0.9994469