計量2x2のクロス集計表の場合はオッズ比と調整済み標準化残差は似たような統計量になる。どちらが好まれるかは分野によって異なるので、適宜、好まれる方を使用するのがよいと思う。社会学ではそれほど見かけることはないが使用されるはずである。オッズ比はなぜか知らないがあまり好まれない傾向にあるように感じる。ちなみに医学ではオッズ比が好まれ、調整済み標準化残差を見かけたことがない。オッズ比は95%信頼区間を併記するのが普通である。調整済み標準化残差は5%基準、1%基準、0.1%基準をアスタリスク等で表記する。95%信頼区間はどちらかが1を越えていないかを一つ一つ確認しなければならないので、表を見るのがめんどくさい。一方で、調整済み標準化残差のアスタリスク表記は、ぱっと見で分かるので便利だと思う。有名な統計パッケージでは、SPSSでは調整済み標準化残差は計算できるが、オッズ比の計算ができない。stataはマクロを入れると計算できたように記憶している。計算が一番簡単で効率よくできるのはRであるのは間違いない。下準備クロス集計表のマトリックスを作成しておく。nrow=2)データセットからも作成しておく。CPS1985$married)性別と結婚ステータスを選んだ。2x2のクロス表が作りたいのでどちらも二値データの変数である。実際の分析としてはあまり意味がないかもしれない。調整済み標準化残差はRのデフォルト機能で出力できる。-4.166099絶対値4.17である。0.259397こちらは絶対値0.26である。調整済み標準化残差は標準正規分布の区間推定値と同様の値をとる。95%信頼区間は-1.96~1.96であり、残差ののように期待値からは外れていることを示す数値はこの95%信頼区間の外でああればよい。したがって、1.96より大きいと5%水準で有意な差があることになる。ように、1%水準は2.58、0.1%水準は3.29、0.01%水準は3.89以上になる。したがって、例1は4.17であったので0.01%水準で有意だったが、例2は1.96に満たないため、有意ではないということになる。標準正規分布(余談)ちなみに1.96というような数字は標準化しているために使用できる基準である。標準正規分布表からこの数字は読み取れる。正規分布の上側確率http://aoki2.si.gunma-u.ac.jp/CGI-BIN/tgxp.htmlリンク先のでは1.96が例に上がっている。1.96のセルは0.0250であるで、これを2倍(両側にするため)すると0.05=5%水準となる。オッズ比オッズ比を出すには、パッケージのインストールが必要である。いくつか存在しているが、一番いろいろできるのはepitoolsパッケージである。いわゆるオッズ比というとこの方法のことである。2019/10/15)無条件最尤推定いわゆるオッズ比というと無条件最尤推定のことである。oddsratio.wald関数を利用する。oddsratio.wald(d1)出力すると以下。3.098556e-05必要なのは$measureのところなのでそこだけ取り出して出力する。#後ろに$measureをつけるそうするとオッズ比と95%区間のみが表示される。クロス表の関連を総当たりで確かめることをよくするが、そういう時に便利だと思う。)このエントリーは統計学勉強会用の下書きである。勉強したことや気になったことをメモするところ。(7)はてなブログをはじめよう!iDESさんは、はてなブログを使っています。ブログを報告するスターをつけました引用をストックしましたストック一覧を見る閉じる引用するにはまずログインしてくださいログイン閉じる引用をストックできませんでした。再度お試しください閉じる限定公開記事のため引用できません。