R:カイ二乗と自由度を指定してp値を計算する


11

カイ二乗と自由度を指定してp値を計算するにはどうすればよいですか?たとえば、カイ二乗= 15、df = 2の正確なp値は何でしょうか?


p値を決定するには、片側か両側かを指定する必要があります。両側の場合はP(| X |> 15)で、片側の下の尾の場合はP(X <-15)および片側の上の尾の場合はP(X.> 15)です。Xは、2自由度のカイ二乗確率変数です。
マイケルR.チェニック2016

両面です。このテストを実行するためのRのコードは何でしょうか?
マット

Rはよく知りません。あなたが呼ぶことができるカイ二乗関数があると思います。Rのほかに、多くの基本的な統計書で利用可能なカイ2乗分布の表から値を取得できます。テーブルの唯一の問題は、15が表示されない可能性があることです。この場合、intrtpolationによって接近できます。
Michael R.Chernick

自分で訂正させてください。カイ二乗は負になりません。片側検定の場合、P(X> 15)が必要です。私の本では、P(X> 13.815)は0.001です。したがって、p値は<0.001です。
Michael R.Chernick

2
(1)カイ2乗検定が両面になることはまれです。アプリケーションは何ですか。(2)(cc @michael)尾に補間しない。積分を推定します。ところで、2 DFのカイ2乗分布は指数分布であるため、これは初等積分です。上部裾確率のを取得する必要があります。en.wikipedia.org/wiki/Chi-squared_distributionを参照してください。e15/2
whuber

回答:


18

適用統計では、二乗検定統計は残差の2乗の合計として、または2乗効果の合計から、または対数尤度差から発生します。これらすべてのアプリケーションでの目的は、いくつかのベクトルパラメーターがゼロかどうかをテストすることであり、それが非ゼロであり、カイ二乗統計量は観測された効果の二乗サイズに関連しています。必要なp値はカイ二乗値の右裾確率であり、例ではRで次のようになります。

> pchisq(15, df=2, lower.tail=FALSE)
[1] 0.0005530844

他のdf値または統計値については、明らかにそれらを上記のコードに置き換えるだけです。

Rのすべての累積確率関数は、デフォルトで左テール確率を計算します。ただし、これらにもlower.tail引数があり、常にこれFALSEを設定して正しい確率を得ることができます。一部の小学校の教科書で見られるように、を計算するのではなく、これを行うことをお勧めします。1p


1
後退することは可能ですか?P値を使用してカイ2乗値とRコードを使用して自由度を導出しますか?
エリック

2
@Eric Rに次のように入力します。help(qchisq)
Gordon Smyth

3

Rは、フォームの濃度又は質量の確率関数のスイートを有しているd*(例えばdbetadchisq)、及び形態の分布p*(例えばpfpgamma)。そこから始めたいと思うかもしれません。


後退することは可能ですか?P値を使用してカイ2乗値とRコードを使用して自由度を導出しますか?
エリック

0

はい、与えられた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の例を使用すると、ソリューションは次のようになります。

解決策:p値を計算する

pchisq(15、df = 2、lower.tail = FALSE)#答え:p = 0.0005530844

カイ二乗値を取得するには、p = 0.0005530844およびdf = 2を使用します

qchisq(0.0005530844、2、lower.tail = FALSE)#答え:カイ二乗= 15

お役に立てれば!!!


-3

試して、

pchisq(chi,df)

あなたの例では、

pchisq(15,2)

[1] 0.9994469

1
これはp値ではないと思います。おそらく累積確率でしょうか?whuberが言ったように、p値はexp(-15/2)であり、0,999469よりはるかに小さいです。
Michael R. Chernick 2017年

7
これは誤りです。においテストにも合格しません。カイ2乗分布の期待値がそのdfであり、分散が2 * df(wikipedia)であることを考慮してください。したがって、15は平均から約9 SDです。ここで何か問題があることがわかります。本質的に、これを使用して、フィットが「良すぎる」片側検定を実行します(ここを参照)。カイ二乗の一般的な使用法では、、pchisq(15, 2, lower.tail=FALSE)または1-pchisq(15, 2)、が必要です0.0005530844
ガン-モニカの回復

(それはそれはない答えとして削除すべきではないので、これは、でも正しくない場合、答えは、言った私は投票しています。ルックスのOKを。)
GUNG -復活モニカ

@gung +1今日は自分の回答を追加しましたが、最初にコメントを読んでいなかったため、同じ明白な解決策が得られました。
ゴードンスミス2017年

後退することは可能ですか?P値を使用してカイ2乗値とRコードを使用して自由度を導出しますか?
エリック
弊社のサイトを使用することにより、あなたは弊社のクッキーポリシーおよびプライバシーポリシーを読み、理解したものとみなされます。
Licensed under cc by-sa 3.0 with attribution required.