Ich verwende die Quantilregression, um Prädiktoren für das 90. Perzentil meiner Daten zu finden. Ich mache dies in R mit dem quantreg
Paket. Wie kann ich für die Quantilregression bestimmen, die angibt, wie viel Variabilität durch Prädiktorvariablen erklärt wird?
Was ich wirklich wissen möchte: "Jede Methode, mit der ich herausfinden kann, wie viel Variabilität erklärt wird?" Signifikanzniveaus von P - Werten sind in der Ausgabe des Befehls zur Verfügung: summary(rq(formula,tau,data))
. Wie kann ich fit werden?
Antworten:
Koenker und Machado [ 1 ] beschreiben R 1 , ein lokales Maß für die Anpassungsgüte am jeweiligen ( τ ) -Quantil.[1] R1 τ
Let und für das gesamte Modell der Koeffizient Schätzungen sein, und einem eingeschränkten Modell und lassen und sein die entsprechenden Terme.β^(τ) β~(τ) V^ V~ V
Sie definieren das Anpassungsgütekriterium .R1(τ)=1−V^V~
Koenker gibt Code für hier ,V
Wenn wir also für ein Modell mit einem Intercept-Only ( - oder im folgenden Codeausschnitt) und dann einem uneingeschränkten Modell ( ) berechnen, können wir - zumindest theoretisch - etwas wie das übliche .V V~ V^ R2
V0
R1 <- 1-Vhat/V0
Bearbeiten: In Ihrem Fall ist natürlich das zweite Argument, das
f$tau
in dem Aufruf in der zweiten Codezeile eingefügt wird, der vontau
Ihnen verwendete Wert . Der Wert in der ersten Zeile legt lediglich den Standard fest.Das Erklären der Varianz über den Mittelwert ist nicht das, was Sie mit der quantilen Regression tun. Sie sollten also nicht damit rechnen, ein wirklich gleichwertiges Maß zu haben.
Ich denke nicht, dass sich das Konzept von gut auf die Quantilregression übertragen lässt. Sie können wie hier verschiedene mehr oder weniger analoge Größen definieren, aber unabhängig davon, was Sie auswählen, werden Sie nicht die meisten Eigenschaften haben, die reales in der OLS-Regression hat. Sie müssen sich darüber im Klaren sein, welche Eigenschaften Sie benötigen und welche nicht - in einigen Fällen kann es möglich sein, dass eine Kennzahl das tut, was Sie wollen.R2 R2
-
und verwandte Inferenzprozesse für die quantitative Regression,
Journal of the American Statistical Association, 94 : 448, 1296-1310
quelle
tau
beim Aufruf der Funktion nichts angeben . Ich werde in der Post klären.Das von Koenker und Machado (1999) in JASA vorgeschlagene Pseudo- Maß misst die Anpassungsgüte, indem es die Summe der gewichteten Abweichungen für das interessierende Modell mit der gleichen Summe aus einem Modell vergleicht, in dem nur der Achsenabschnitt auftritt. Es wird berechnet alsR2
Dabei ist das angepasste te Quantil für die Beobachtung , und ist der angepasste Wert aus dem reinen Achsenabschnitt Modell.y^i=ατ+βτx τ i y¯=βτ
Hier ist ein Beispiel in R:
Dies könnte wahrscheinlich eleganter durchgeführt werden.
quelle
R_1(\tau) = 1 -
Das letzte Zeichen ist ein Durcheinander. Könnten Sie das überprüfen? Vielleicht haben Sie ein nicht standardmäßiges Zeichen eingefügt, anstatt Tex zu verwenden.