ロジスティック回帰のハットマトリックスからの情報


12

私には明らかであり、複数のサイトで、ハットマトリックスの対角線上の値が線形回帰にどのような情報を与えるかをよく説明しています。

ロジスティック回帰モデルの帽子行列は、私にはあまり明確ではありません。それは、線形回帰を適用して帽子行列から得た情報と同一ですか?これは、CVの別のトピックで見つけた帽子行列の定義です(ソース1):

H=VバツバツVバツ1バツV

Xでは予測変数のベクトル、Vは対角行列です。π1π

言い換えれば、観測値のハットマトリックスの特定の値は、共変量空間内の共変量の位置も示すだけで、その観測値の結果値とは何の関係もないということですか?

これはAgrestiの「Categorical data analysis」という本に書かれています。

観測値のレバレッジが大きいほど、フィットに対する潜在的な影響が大きくなります。通常の回帰と同様に、レバレッジは0から1の間にあり、合計はモデルパラメーターの数になります。通常の回帰とは異なり、帽子の値は近似とモデル行列に依存し、極端な予測値を持つポイントは高いレバレッジを持つ必要はありません。

したがって、この定義から、通常の線形回帰で使用するため、使用できないように見えますか?

ソース1:Rのロジスティック回帰のハットマトリックスの計算方法

回答:


13

表記を少し変更して、H = V 1のようにハットマトリックスを記述します。 Vは、一般的な要素が対角対称行列であり、Vは、J=MjのπXJ[1-πXJ]。表すm個のJ同じ共変量の値を有する個体の群として、X=XJ。あなたが得ることができるj個のTHの対角要素(HJとしてハット行列の) H、J=MjのπXjは

H=V12バツバツVバツ1バツV12
Vvj=mjπバツj[1πバツj]mjバツ=バツjjthhj その後の合計時間jは線形回帰のようなパラメータの数を与えます。あなたの質問に今:
hj=mjπバツj[1πバツj]バツjバツVバツ1バツj
hj

π0.1<π<0.9

ここに画像の説明を入力してください

バツjバツVバツ1バツjhj、したがって、この別個の単調な部分が単独で考慮されることはほとんどありません。

このトピックをさらに深く読みたい場合は、ロジスティックハットマトリックスを導出したPregibon(1981)の論文と、Hosmer and Lemeshow(2000)の本をご覧ください。

弊社のサイトを使用することにより、あなたは弊社のクッキーポリシーおよびプライバシーポリシーを読み、理解したものとみなされます。
Licensed under cc by-sa 3.0 with attribution required.