杉野 勇 (SUGINO Isamu)　お茶の水女子大学・人間発達科学専攻・応用社会学コース担当

　本文の標本相関係数の標本抽出分布（一部）を描くスクリプトを紹介しておこう。単回帰分析を深く理解していないと意味が分からないと思うので，解読できなくて良い。シミュレイションなので，実行するたびに少しずつ違う結果が得られる。NやR，n，timesを変えて実行してみると良い。
　timesとnの値を変えて実行したグラフを以下に示す。この場合の標本相関係数の理論的な期待値は約.2986である。

　標本サイズと標本相関係数から，ゼロ仮説「母相関係数=0」のt検定を行う（有意確率を表示する）スクリプトを示す。片側検定と両側検定の両方を表示しているが，両側検定での有意確率は片側検定での有意確率の2倍である。よって，両側検定で有意になる場合，同じ有意水準の片側検定では必ず有意になるが，片側検定で有意になったとしても，同じ有意水準の両側検定で有意になるとは限らない。
nやrに自由な値を設定して試してみよう。rの絶対値はもちろん1を超えてはならない。

　本文の表のように，n，rの幾つかの組合せに対して，95%限界値，有意確率を計算してみよう。本文よりも，限界値に対応する標本相関係数の値など少し追加してある。

　上の結果を解説しよう。n=50，r=.20の場合，検定統計量tは約1.414となり，95%範囲の両端の±2.010635よりずっと中心寄りであり，両端の限界値の外側に位置する5%棄却域には入らない。即ち，「母相関係数はゼロである」とのゼロ仮説を棄却出来ない。両側有意確率は約16.4%である。このサンプルサイズで5%棄却域に入る（＝ゼロ仮説が棄却できる）のは，r=.2787くらいになってからである。

　サンプルサイズと標本相関係数が与えられた時に，信頼区間を計算する方法は，本文脚注にも書いてある通りである。フィッシャーのZ変換を行ったものはそのままZと表記する事が多いが[南風原 2002: 116]，ここでは他の多くの標準正規変量と区別し易いようにあえてLと表記しておく。
　ここで，フィッシャーのZ変換を行ったものが実際に正規分布で近似できるかをシミュレイションで確かめてみよう。上の1で使用した模擬母集団とそこからのランダムサンプリングのスクリプトを利用し，多くの標本相関係数を格納したヴェクトルcorrsまで求められているとし，その続きを示す。まずは，シミュレイションで得られた標本相関係数のヒストグラムを再確認しておこう（シミュレイションなので実行するたびに少しずつ異なる結果になる）。


確かに，かなり正規分布に近付いている様子が分かる。

　本文で「t検定の結果と95%信頼区間，標本相関係数の値は，Rのcor.test(x, y, method=”pearson”)によってすべて出力される。」と書いたので，その出力例を以下に挙げておく。

　本文の，「自由度1,4,6,8,9,12のカイ二乗分布」のグラフを描くスクリプトは以下である。
　次の，棄却域を塗り潰したカイ二乗分布のグラフは下記。自由度と有意水準は変更する事が出来る。
　実際のデータ（演習用データ）から分割表のカイ二乗検定を行うには下記の様にする。欠損値処理が入っているので少し手間がかかっている。ついでにクラメールのVやファイ係数の計算もしておこう。

　分割表の独立性についてのカイ二乗検定において，検定統計量がカイ二乗分布で近似できる条件は，各セルの度数がポワソン分布に従い，そのポワソン分布が正規分布で近似できるということである。本文中の，二つのポワソン分布を重ね描きするスクリプトは下記である。
　少し改変したスクリプトとその結果も例示しておこう。

　社会学では，「カイ二乗検定」と云うと「独立性の検定」と同一視する傾向があるが，これはデータに対して「2変数独立モデル」が適合するかどうかを検定していると考えれば，「適合度検定」の一種であるとも言える。カイ二乗統計量は社会統計学のより発展的な場面で色々と登場してくる事になるが，いずれも「データとモデルのズレ」を表す統計量として統一的に考えるのが良い。ここでは，chisq.test( )関数を使用した最も単純な適合度の検定を紹介しよう。
　少人数の無作為標本に対して，「夫は外で働き，妻は家庭を守るべきである」という考え方について，そう思うかそう思わないか（性別役割分業意識）を，「そう思う／どちらかと言えばそう思う／どちらかと言えばそう思わない／そう思わない」の4件法で尋ねた結果，それぞれ5人，8人，15人，12人であったとする。
　「この結果では偶然ばらついたが，本当はどの回答も等確率で出現する」，つまり同じ比率ずつ分布していると仮定した時，その等比率（もしくは等確率）モデルと実際のデータのズレが偶然の範囲内と言えるかどうかをカイ二乗検定する事が出来るが，これが最も単純な適合度検定である。

　等確率だとしたら全ての選択肢が10人ずつに選ばれる筈であり，実際の結果はそこからカイ二乗統計量にして5.8ずれているが，これは自由度3(=カテゴリ数4 - 1)のχ²分布に近似的に従うと考えられる。自由度3のχ²分布で5.8以上の範囲の面積（確率）は約12.2%ある。つまり本当は（母集団では）等比率で分布しているのにたまたまこの様な無作為標本が得られることは10回に1回より多く起こると考えられ，特に不自然であるとは言えないので，「本当は（母集団では）等比率である」と云うゼロ仮説は棄却出来ない。
　これを，等比率モデルの適合度は悪くないと考えるのが適合度検定である。「モデルがデータに適合している」と云うゼロ仮説が棄却されなかった事をもって「モデルは正しい」と解釈する事は統計の誤用であるとする立場からするとこの適合度検定の考え方には疑問や批判もある。その立場からすると，「等比率モデルは間違っていると自信をもって言えるだけの証拠が無かっただけで，等比率モデルが正しいと言える訳ではなく，精々判断が保留されるだけである」となる。この辺りは研究者によって考え方が違ってくるかも知れない。
　上の例では，適合度検定と独立性の検定の違いがいまいちはっきりしないかも知れないので，少し違った例を示そう。
　上では等比率仮説は棄却されないが，流石に現在の日本ではこの質問に対して4つの選択肢で同じ比率ずつ意見が分布しているとは考えられない。本当の意見分布を2:3:4:5であると仮定して，データと齟齬を来すかどうかを検討しよう。Rスクリプトと結果は以下の通りである。chisq.test( )関数のオプションp=で，それぞれのカテゴリの生起確率を指定する点だけが異なる。

結果，性別役割分業に否定的な人ほど数が多いと云うこのモデルも有意性検定で棄却されず，データとの適合は悪くないと云う結果になった。むしろズレの指標は等比率モデルよりも遥かに小さく，有意確率もずっと高い。いずれのモデルも（そして他にも多くの可能なモデルが）データに適合する（適合しないとは言い切れない）と云う結果になる。そのいずれかが正しいと証明される訳でも無く，カイ二乗統計量や有意確率の値からいずれかが「より正しい」と言えるかどうかにも議論がある。
　適合度検定と考えられるものには，第6章の基礎1-1や第7章基礎1-5の等分散性の検定（各群の母分散が等しいモデルが適合するかどうか），第6章発展1の正規性の検定（変数の分布が正規分布に適合するかどうか）がある。それ以外にも，第11章発展1-1の構造方程式モデリング(SEM; Structural Equation Modelling)では，設定したモデルの適合度がカイ二乗検定される（それ以外に，GFI[Goodness-of-Fit Index]を始めとして幾つかの適合度指標が算出されるが，GFIらについては統計的検定がされる訳では無い）。第12章基礎1-4では，2項ロジスティック回帰分析に対する適合度検定としてHosmer-Lemeshow検定を紹介した。第12章発展1のログリニア・モデル(log-linear model)も「モデルがデータに適合するか」を重視した分析手法であると言える。
　適合度検定の多くは，統計的検定一般がそうである様に，サンプルサイズが大きいと有意になり易い。大標本では「僅かなズレでも敏感に検出してしまう」のである。適合度検定で統計的に有意になると云う事は，「モデルがデータに適合している」と仮定すると（ゼロ仮説），ズレの指標が偶然とは思えない程に大きくなるので，ゼロ仮説は棄却されると云う事であり，つまり「このモデルはデータには適合しない」と云う結論になるので，モデルが当てはまる事を期待していた分析者にとっては嬉しい結果ではない事も多いだろう。

　上の2-1でデータフレイムdata01，その中の変数data01$sex，data01$q1601b，分割表table1などを既に準備済みであるので，それを利用して調整済み標準化残差（ハーバーマン残差）を表示したりモザイクプロットを表示したりしてみよう。

　模擬的に2行2列の分割表を作成し，カイ二乗検定の関数を適用してみよう。correct=というオプションを省略した場合，TRUEにした場合，FALSEにした場合をよく比較して欲しい。2行2列分割表ではcorrect=を省略すると自動的にcorrect=TRUEとなる。行数か列数が2を超えると連続修正は存在しない。
　同じ分割表に対してフィッシャーの正確検定を行ってみる。fisher.test( )関数の出力するodds ratioは注意が必要なので，通常のオッズ比を計算するコマンドも最後に付け加えてある。

　演習用データの学歴（3行）×性別役割分業意識（4列）の分割表にカイ二乗検定や正確検定を行ってみよう。2行2列ではないのでオッズ比は計算されない。

1)の答え：

　有意確率は18.3%程度であり，有意水準を1%にしようが5%にしようが10%にしようが有意にはならない。つまりゼロ仮説は棄却出来ない。片側検定として考えるなら有意水準を半分にすればよいが，それでも5%では有意にならない。

2)の答え：

　両側5%検定で有意になり，ゼロ仮説は棄却される。1)とは標本相関係数の値は同じであるが，標本サイズが2倍以上になっており，弱い関連に対してより敏感に検出出来る様になっている。

3)の答え：

高度に有意となり，母集団では独立であるとのゼロ仮説は棄却される。
独立性が否定されたので，行比率や列比率を計算して，どの様な偏りがあるのかを確認しておく。


4)の答え：
　カイ二乗検定のみならず正確検定も行っておこう。

いずれにしても有意にはならず，独立モデルは棄却出来ない。3)と比較すると，分割表の行数・列数は等しく，したがって自由度も同じである。各セルの度数を一律に三分の一にしているので，行比率の表もほぼ同じである。
　しかしカイ二乗統計量の値は，30.377と9.5046と3倍余り異なり，その結果検定結果は1%有意と10%でも有意でないと云う様に対照的である。サンプルサイズが小さいと，同じ様な関連が「偶然に過ぎないかも知れない」と見做されてしまうのである。

5)の答え：
　調整済み標準化残差（ハーバーマン残差）が標準正規分布に近似的に従う事を利用する。
調整済み標準化残差は，カイ二乗検定の結果を格納したオブジェクトのstdresで取り出せる。

有意水準を両側5%とすると，絶対値が1.96を超えているセルが，有意に独立モデルから乖離しているセルであり，[3行, 1列]が期待値より有意に，かなり多い。[3, 2], [3, 3], [1, 1], [1, 2]は期待値より有意に少ない。つまり，3行1列のセルに著しく集中している事が分かる。
　実はこの分割表は，演習用データの学歴3区分と従業上の地位4区分のものであり，高学歴者ほど常時雇用者が多いと云う結果を示しているのである。

　第2章と第5章を中心として，本書では分割表を作成する為に table( )関数を使用しているが，他によく使用される関数として，xtabs( )があるのでここで補っておく。
　まずは7ケース分のデータ・フレイムdata01を作成してから，table( )とxtabs( )で三元分割表を作成してみる。その後，ftable( )を利用して三元分割表を見易くする。
　コマンドの文法と出力スタイルの両方を見比べてみよう。



　この限りでは，特に重回帰分析（一般線型モデル）などに慣れている人ほど，xtabs( )の方を好むだろう。出力の表記だけをとってもxtabs( )の方がやや親切である。
　実際の調査データ分析では，欠損値（欠測missing）の処理が必須である。この点も違いがあるので比べておこう。




　NAの処理は，xtabs( )の方が何だか面倒である。
　分割表からデータ・フレイムへの変換（逆変換？）を行う関数，as.data.frame.table( )があるが，これを使用する場合には，table( )の結果よりもxtabs( )の結果の方が便利である。上のdata01を使って実際に示してみるので，どう云う事か違いをよく見比べて読み取って欲しい。


　上のas.data.frame.table( )関数の結果をよく見ると，出力されたデータ・フレイムにはFreqと云う変数が新たに加わっている。この例では，それぞれのパタンに該当するのはせいぜい1ケースしかなかったのでFreqは0か1の値しかないが，全く同じ変数値プロファイルを持つケースが複数あれば，Freqもその2以上のケース数を表す。
　逆に，その様なFreq変数（変数名は別でも良い）を持つデータ・フレイムから，その度数によって重み付けを行った分割表を作成する事も xtabs( )によって可能になる。


　これを利用すると，単に度数で重みを付けるだけでなく，サンプルの歪みの補正にしばしば使われる重み付け（ウェイト付け）を行う事も出来る。原理的には，重回帰分析などのウェイト付けと同じである。ウェイト付けを行わない場合と行った場合では，そもそも分布が違ったものになるので，分析結果も異なり得る。