固定効果ロジスティック回帰のRパッケージ


14

RChamberlainの1980推定器を使用して、個々の固定効果(個別インターセプト)を使用してロジットモデルの係数を推定するためのパッケージを探しています。チェンバレンの固定効果ロジット推定器としてよく知られています。

(少なくとも計量経済学で)バイナリの結果パネルデータを扱う場合、これは古典的な見積もりツールですが、CRANに関連するものは何も見つかりません。

どんな手掛かり?


ここに別の試みがあります:stats.stackexchange.com/questions/10141/…–
アレックス

私は同じ状況に対処していますが、ソリューション/パッケージ/コードを見つけましたか?
マリオGS

回答:


12

条件付きロジスティック回帰(チェンバレンの推定量について話しているときに参照したものと仮定します)はサバイバルパッケージclogit()で利用できます。また、条件付きロジットパラメーターを推定するRコードを含むこのページを見つけました。調査のパッケージには、複雑なサンプリングの場合はGLMと生存モデルのラッパー機能の多くを含んでいるが、私は見ていませんでした。

またlogit.mixedZeligパッケージを確認するか、二項リンクを使用した混合効果モデルのメソッドを提供するlme4パッケージを直接使用してください(lmerまたはを参照glmer)。

あなたは見とりましたRでの計量経済学グラントV.ファーンズワースから、?Rに適用された計量経済学の穏やかな概要を提供しているようです(これは私には馴染みがありません)。


1
実際、「条件付きロジット」は非常に曖昧な用語です。一部のコンテキスト(主にパネルデータを処理する場合)は、チェンバレンの推定量に相当しますが、非常にまれです。ほとんどの場合、結果変数が3つ以上の値を取ることができる断面モデルを指します。すべての提案は、実際にこの最後の可能性を考慮したパッケージを参照しています。混合ロジットと同じ:固定効果ロジットではありません。私はすでにファーンズワースの概要を見てきましたが、この推定量について語るのに十分ではありません。とにかく、答えてくれてありがとう!
-Kamixave

「条件付きロジット」は、3つ以上の結果レベルを持つことを意味しませ。一部の関数はその状況にそれを拡張するかもしれませんが、それはポイントではありません。
アニコ

1
ええ、ただし、条件付きロジットモデルは(私が言ったように)2つ以上の値を取ることができます。これは、Chamberlainがパネルデータ専用に設計されているという事実と同様に、Chamberlainのモデルと簡単に区別できます。したがって、これは関連情報です。通常の条件付きロジットの正確な説明はそうではありません(両方の説明には600文字以上かかります)。
-Kamixave

2

を使用してChamberlainsモデルを実行できますglmer。基本的にはREモデルですが、より多くの変数があります:

glmer(y~X + Z + (1|subject), data, model=binomial("probit"))
  • Xは、固定効果モデルを説明する変数です(単純な場合はZの平均です)。
  • Zはあなたの外生変数です
  • 件名は、不均一性の原因となる変数です

これがお役に立てば幸いです。


2
私は、要求された推定器がそれを可能にしている間、不均一性がXとZに直交するように制限すると思います。
アレックス

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