第10回授業(12/4)の目標 カイ2乗検定の実習 WEB を用いたカイ2乗検定と、授業で行った検定結果の正誤の確認方法(宿題)
いつもの架空のデータでの演習 岩原データを用いて、カイ二乗検定の演習を行ってみよう。
WEB 宿題の手順1 千野のホームページの講義テキストの中の「心理統計学」第7章の 7.9 節 「統計ソフト SAS を用いたクロス表のカイ二乗値の計算及び検定の実行手順」、をクリックし、
WEB 宿題の手順1(続き) 7.9.1.1 節 「データ例とその入力方法」、を見ながら、 7.9.1.1 節 「データ例とその入力方法」、を見ながら、 terapad を起動し、各自のデータ20対を 上記WEBページの入力例と同方法に従って入力し、 各自の学籍フォルダの直下の data フォルダに、ファイル名 chisq_ex1.txt なる名前をつけて保存する。
データ入力時の注意 データは、ファイルの左端から2桁でサンプル番号を、 半角ブランクを1つあけて、対の最初の数値(うえの例では、性別と仮定したもの)を入力し、 最後に、やはり半角ブランク1つあけて、対の他方の数値(うえの例では、向性と仮定したもの)を入力すること。
WEB 宿題の手順2 7.9.1.4 節「SAS による具体的手順」をクリックし、そこでの手順に従い、 インターネットエクスプローラ (IE) を開き、千野のホームページの ダウンロードコーナーから、当該プログラ ムを右クリックし、「対象をファイルに保存」 を選択することにより、
WEB 宿題の手順2(続き) Pドライブの上記自らの学籍番号フォルダの下の sasprog フ ォルダに保存する。 保存すべきファイル名は、chisq_ex1.sas であるが、上記の手続きで sasprog フォルダに保存する時、自動的にファイル名がつくので、これを入力する必要はない。
SAS の起動とプログラムの修正 SAS を起動し、うえで保存した SAS プログラム chisq_ex1.sas を、プログラムエディタ上に読み込んで、 filename 文を ecip 用に修正する。 修正は、つぎのようにする:
SAS の起動と プログラムの修正(続き) 修正は、ファイルネーム文、すなわち、 filename data 'p:\data\chisq_ex1.txt'; の引用符の中(赤印の項)を p:\psyphy\chino\psycstat \Report\学籍フォルダ\data \chisq_ex1.txt と書き換え、上書き保存する。
プログラムの実行と出力の保存 修正が終わったら、プログラムを実行する。 プログラムが正常終了し、アウトプットウインドウに出力結果が出たら、それを sasout フォルダの下に chisq_ex1.lst なる名前で保存する。
出力結果の保存時の注意事項 なお、保存ウ インドウが現れたら、ファイルの種類を LIST ファイルと指定すること。この 時、その上のファイル名は chisq_ex1 のみとすること。 これにより、sasout に保存される出力結果には、 chisq_ex1.lst という名前がつく。
SAS による出力結果の見方(1) 出力結果を保存し終わったならば、WEB ページの、7.9.1.2 「出力結果の例」の項を見て、授業中に各自が計算した結果と一致しているかを確認する。 まず、各セルの度数 a, b,c, d については、出力結果の最初の方にある分割表の各セルの中の6つの数値のうちの最上部のものがそれに該当する。
SAS による出力結果の見方(2) 標本のカイ2乗値(χ2)については、分割表の直ぐ下に出力される catx と caty の統計量の項の χ2 値又は、3番目に出力される「連続性補正χ 2 乗値」のいずれかの値を見よ。 カイ2乗統計量の自由度 (d.f.) については、データにかかわらず、2×2分割表では、1。
SAS による出力結果の見方(3) カイ2乗検定結果(棄却、採択)については、SAS の出力結果の中の、分割表の直下の catx と caty の統計量の中の最初に出力される「χ 2 乗値」 又は、3番目に出力される「連続性補正χ 2 乗値」のいずれかにおける p-値の値を見て決める必要があり、イエーツの修正の有無により異なるので注意せよ。
SAS による出力結果の見方(4) もし、イエーツの修正が必要なデータの場合は、 「連続性補正χ 2 乗値」に対応する p-値を見ること。
SAS による出力結果の見方(5) いずれの場合も、 もし、p-値が 0.05 より大きければ、帰無仮説すなわち、「2つの属性間に関連はない」という仮説を採択し、 もし、 p-値が 0.05 未満であれば(5%水準で)帰無仮説を棄却する。
WEB 宿題の期限 今回のWEB の宿題は、 平成21年12月11日までに行うこと。