Rによるクロス集計と独立性の検定2
ボーンシュテット&ノーキ「社会統計学」(Bohnstedt & Knoke 1988=1990:88)に紹介されている1984年GSS調査の「信仰する宗教と支持政党」データを用いる(引用元データは百分率クロス集計表なので、サンプルサイズから度数を復元している)。
宗教 | 民主党 | 支持政党なし | 共和党 | 合計 |
---|---|---|---|---|
プロテスタント | 329 | 301 | 284 | 914 |
カトリック | 166 | 142 | 65 | 373 |
ユダヤ教 | 9 | 10 | 7 | 26 |
無宗教 | 29 | 63 | 12 | 104 |
その他 | 10 | 8 | 0 | 18 |
合計 | 543 | 524 | 368 | 1435 |
1984年GSS調査 |
このデータから信仰する宗教と支持政党にどのような傾向を読み取れるか。
- George W. Bohnstedt & David Knoke,1988,Statistcs Fo Social Data Analysys 2nd.edithion(=1990,海野道郎・中村隆訳 社会統計学-社会調査のためのデータ分析入門 ハーベスト社)
スクリプト
convert.crosstable.to.weighted.data関数(socialStatisticsBasic.Rで読み込まれる自作関数)
- convert.crosstable.to.weighted.data(matrix)
-
度数クロス表を重み付きデータに変換する。
- matrix=クロス集計表(matrix型)
convert.weighted.data.to.original関数(socialStatisticsBasic.Rで読み込まれる自作関数)
- convert.weighted.data.to.original(data.frame, freq.name)
-
重み付きデータから元データを復元する。
- data.frame=重み付きデータ(data.frame型)
- freq.name="Freq" 度数変数の項目(変数)名
出力結果
クロス集計表
比率は$row.ratio(行方向の比率)を見る。
χ2検定
- $chisq.test
-
- Peason→「Pearsonのカイ2乗:漸近有意確率(両側)」
- Fisher→「Fhisherの直接法:正確な有意確率(両側)」
- クロス表が2*2ではないのでYates修正は出力されない。
残差分析
符号に着目。プロテスタントは共和支持が相対的に強く、カトリックは民主支持が強い。無宗教は「支持政党なし」。
効果量Cramer's V
- Cramer's V
-
Cramer's Vの目安 効果量V 効果の目安 0.5 大 0.3 中 0.1 小 0 なし
考察例
信仰する宗教によって支持政党に影響するかどうかをクロス集計表により検証を行った。
クロス集計表はχ2検定(Peason)により1%水準で有意であり(χ2(8)=62.80, p<.01, V=0.15)、信仰する宗教によって支持政党が異なる傾向があると言える。
残差分析より、プロテスタントは共和支持が相対的に強く、カトリックは民主支持が強い。また無宗教は支持政党を持たない傾向が強い。