2010/05/27

カイ二乗検定

カイ二乗テストです.

まずデータを R に読み込ませます.ここでは data001 という名前にします.
> data001<-matrix(c(7,5,21,38,15,14,22,63,67,46,6,2,31,28,14,9,17,27,84,48),4,5,byrow=T) でクロス集計のデータ(この場合,5 × 4 のデータです)を R に与えます.byrow=T で横に読み込みます(試しに byrow=F とすればこの意味がわかります).c() のデータを matrix 形式で,data001 の名称で読み込ませるという意味になります. このあと > data001
と入力すれば,data001 の状態がわかります.

> chisq.test(data001)

と入力して,ピアソンのカイ二乗検定の結果をえます.カイ二乗値と自由度,p 値を求めることができます.

残差分析等は出力されませんので注意.これについては群馬大の青木先生のところ(http://aoki2.si.gunma-u.ac.jp/R/my-chisq-test.html)に表記あり.また上記についても同志社大学の金先生の資料(http://mjin.doshisha.ac.jp/R/11.pdf)をネット上で見つけることが可能です.