Ich habe gerade zwei Millionen Regressionen durchgeführt - Integrierte Wahrscheinlichkeit

9

Ich arbeite derzeit daran, eine Methode zu implementieren, die in einem populären Artikel mit dem Titel "Ich habe gerade zwei Millionen Regressionen durchgeführt" verwendet wird. Die Grundidee dahinter ist, dass es bestimmte Fälle gibt, in denen nicht klar ist, welche Steuerelemente in das Modell aufgenommen werden sollen. In einem solchen Fall können Sie nach dem Zufallsprinzip Steuerelemente zeichnen, Millionen verschiedener Regressionen ausführen und dann sehen, wie Ihre interessierende Variable reagiert hat. Wenn es in allen Spezifikationen im Allgemeinen das gleiche Vorzeichen hat, können wir es als robuster betrachten als eine Variable, deren Vorzeichen sich immer ändert.

Das meiste Papier ist sehr klar. Das Papier gewichtet jedoch alle diese unterschiedlichen Regressionen auf folgende Weise: Die integrierte Wahrscheinlichkeit der gegebenen Spezifikation wird durch die Summe aller integrierten Wahrscheinlichkeiten für alle Spezifikationen geteilt.

Das Problem, das ich habe, ist, dass ich nicht sicher bin, wie die integrierte Wahrscheinlichkeit mit den OLS-Regressionen zusammenhängt, die ich ausführen möchte (in Stata). Das Googeln von Themen wie "stata integrierte Wahrscheinlichkeit" war eine Sackgasse, da ich immer wieder auf Dinge wie die logistische Regression mit gemischten Effekten stoße. Ich gebe zu, dass diese Modelle zu komplex sind, als dass ich sie verstehen könnte.

Meine derzeitige Arbeit besteht darin, dass in der Literatur verschiedene Gewichtungsschemata verwendet werden, die ich (irgendwie) verstehe. Beispielsweise ist es möglich, jede Regression basierend auf dem Likelihood-Ratio-Index zu gewichten. Es gibt sogar ein R-Paket, das den lri als Gewichte verwendet. Natürlich möchte ich auch das Original implementieren.

Irgendein Rat?

Papierlink: http://down.cenet.org.cn/upfile/34/2009112141315178.pdf

NikolaiB
quelle
1
Dieser Thread kann einige Ihrer Bedenken ansprechen
Mike Hunter
1
Ich habe einmal eine Funktion in MATLAB geschrieben, die das Ergebnis von Sala-i-Martin nachbildet (die übrigens bei der Modellauswahl nicht wirklich auf dem neuesten Stand der Technik sind), siehe dropbox.com/s/mqa7qvhn7w5pkag/… . Die integrierte Wahrscheinlichkeit (nicht sicher, worauf Sie sich genau beziehen) ist wahrscheinlich nur die potenzierte Log-Wahrscheinlichkeit.
Christoph Hanck
Vielen Dank! Ich beziehe mich auf Gleichung 4 auf Seite 179. Es heißt "Wo die Gewichte proportional zu den (integrierten) Wahrscheinlichkeiten sind"
NikolaiB

Antworten:

1

Für OLS können Sie weiterhin die Wahrscheinlichkeitsfunktion berechnen (die potenzierte Protokollwahrscheinlichkeit, wie Christoph Hanck im Kommentar erwähnt). Es ist nur das gute alte . Stata speichert dies wie nach dem Ausführen einer Regression mitLi=i(2πσ2).5exp(.5(yixiβ)2)e(ll)regress

Dann konstruieren Sie Gewichte als .wi=LijLj

Schließlich konstruieren Sie gewichtete Mittelwerte Ihrer Regressionskoeffizienten unter Verwendung von als Gewichte.wi

Superpronker
quelle