Testäquivalenz von nicht verschachtelten Modellen

12

Angenommen, y ist eine lineare Funktion von und ein Dummy . Meine Hypothese ist, dass selbst wie ein hedonistischer Index eines Vektors anderer Variablen ist, . Ich habe Unterstützung dafür in einer von (dh , , ..., ) auf . Gibt es eine Möglichkeit, die Gleichwertigkeit dieser beiden Modelle zu testen :xddZZ z 1 z 2MANOVAZz1z2znd

Modell 1:y=b0+b1x+b2d+e1

Modell 2:y=g0+ZG+e2

wobei der Spaltenvektor von Parametern ist.G

user3671
quelle

Antworten:

8

Zunächst müssen Sie das Äquivalenzkonzept definieren. Man könnte denken, dass zwei Modelle äquivalent sind, wenn sie fast die gleiche Prognosegenauigkeit liefern (diese wäre für Zeitreihen und Paneldaten relevant), eine andere könnte interessiert sein, wenn die Übereinstimmungen des Modells nahe beieinander liegen . Ersteres ist das Objekt für verschiedene Kreuzvalidierungen (Jack-Knife-Tests in der Regel oder einige Tests außerhalb der Stichprobe, Robs accuracy()tun dies gut), letzteres dient der Minimierung einiger Informationskriterien.

In der Mikroökonometrie ist BIC die Wahl , obwohl Sie auch AIC Betracht ziehen können, wenn Sie mit kleinen Probengrößen arbeiten. Beachten Sie, dass die Auswahl basierend auf der Minimierung des Informationskriteriums auch für verschachtelte Modelle relevant ist.

Eine nette Diskussion in gegeben ist must-have-it - Buchung von Cameron und Trivedi (Kapitel 8.5 bietet eine hervorragende Überprüfung der Methoden), spezifischere theoretische Details sind in Hong und Preston gefunden hier .

Grob gesagt wird empfohlen, aus zwei Modellen das sparsamere zu wählen (mit weniger zu schätzenden Parametern, also mehr Freiheitsgraden). Ein Informationskriterium führt eine spezielle Straffunktion ein, die die Einbeziehung zusätzlicher erklärender Variablen in ein lineares Modell in konzeptioneller Hinsicht einschränkt, ähnlich den durch angepasstes eingeführten Einschränkungen .R2

Möglicherweise sind Sie jedoch nicht nur an der Auswahl des Modells interessiert, das das ausgewählte Informationskriterium minimiert. Das Äquivalenzkonzept impliziert, dass eine Teststatistik formuliert werden sollte. Daher können Sie für Likelihood-Ratio-Tests entweder Cox- oder Voung Tests, Davidson-MacKinnon J- Test. LRJ

Schließlich könnten Sie den Tags zufolge nur an RFunktionen interessiert sein :

library(lmtest)
coxtest(fit1, fit2)
jtest(fit1, fit2)

fit1fit2coxtestLRjtestJ

Dmitrij Celov
quelle
Danke, Dmitrij. Wenn ich das richtig verstehe, sind sowohl coxtest als auch jtest im Wesentlichen modifizierte verschachtelte Tests. Schritt 1: Führen Sie das Modell mit dem kombinierten Pool von Regressoren aus Modell 1 und Modell 2 aus. Schritt 2: Testen Sie Modell1 und Modell2 einzeln als Teilmengen des "Supermodels". Habe ich recht? Gibt es in Bezug auf IC-Messungen auch eine Möglichkeit, die AIC / BIC-Unterschiede zwischen Modell 1 und 2 statistisch zu vergleichen? Hinweis: Ich versuche NICHT, das "beste" Modell auszuwählen, aber Sie haben Recht, ich versuche im Wesentlichen zu testen, ob zwei Modelle die gleiche Passform haben.
user3671
jtestcoxtestLR