「統計学関連なんでもあり」の過去ログ--- 046

No.20936　数量化Ⅰ類のアイテム変数間の相関係数について　　【arai】　2014/03/10(Mon) 18:08

いつも参考にさせて頂いております。

Rを用いて，複数の生活習慣項目（歯磨きあり／なし，など）と抑うつスコアとの関連性を数量化Ⅰ類で検討しようとしております。

多重共線性を排除するために，生活習慣項目どうしの相関係数を見たいのですが，アイテム変数どうしの相関係数は，サンプルスコアをもとに計算されるんですよね。

先生作の関数qt1から，相関係数を算出できますでしょうか？そのプログラムは含まれていないように思うのですが，自分で追加する技能がありません。

qt1への追加プログラムをご教示頂ければ大変ありがたいです。よろしくお願いいたします。

No.20937　Re: 数量化Ⅰ類のアイテム変数間の相関係数について　　【青木繁伸】　2014/03/10(Mon) 18:19

戻り値に追加しました
summary 関数で表示できます（要素名 r）

しかし，多重共線性を相関係数から判定してはいけません。「0.8以上で多重共線性ありとする」などは，迷信です。
数量化I類ではなく，ダミー変数を使用する重回帰分析（カテゴリー変数を factor にして，普通の lm）を行い，その状態で VIF などを見るべきでしょう。

No.20943　Re: 数量化Ⅰ類のアイテム変数間の相関係数について　　【arai】　2014/03/11(Tue) 15:27

青木先生

早速ご教示いただき，どうもありがとうございました。
qt1のsummaryから確認することができました。

また，多重共線性の判断についても，ご助言ありがとうございます。
相関係数は迷信なのですね。勉強になります。

VIF を何とか算出しました。手順としては，make.dummy関数を使って自分でカテゴリー変数をダミー変数に変換してからlmを行い，VIFを算出するということで合っていますでしょうか。lmは，カテゴリー変数をダミー変数に変換しなくても使えますが，それだとfactorのままなのでVIFは算出されないですよね？