Was sind die roten Linien im Quantil-Regressionsdiagramm (Quantreg-Paket)?

9

Mit plot.rqdem quantregPaket in R können wir die Koeffizientenschätzungsverteilung darstellen und so etwas erhalten: Geben Sie hier die Bildbeschreibung ein

Was sind die gepunkteten roten Linien? Ausgiebiges Googeln hat ergeben, dass die mittlere der Durchschnitt aller 99 Schätzungen ist, aber immer noch nichts über die gepunktete rote Linie weiß.


quelle
6
Ausgiebiges googeln hatte es falsch gemacht: Dies sind die OLS-Schätzungen und sein SE
user603
@ user603 (i) Haben Sie eine Referenz dafür? (ii) Meinst du das CI? Ist das 95% oder 99%?
1
@AndreSilva: Danke, dass du diese Frage in den Vordergrund gestellt hast. Vielen Dank auch für die Bearbeitung der Diagramme in meiner Antwort.
user603

Antworten:

8

Hier ist eine Reproduktion der Engel-Regression aus der Quantreg- Vignette . Um die Darstellung besser lesbar zu machen, konzentriere ich mich auf den zweiten (Steigungs-) Parameter:

library(quantreg)
data(engel)
xx   <- engel$income-mean(engel$income)
fit1 <- summary(rq(engel$foodexp~xx,tau=1:9/10))
plot(fit1, parm=2)

Ich kann jetzt die üblichen OLS-Schätzungen und das (5-95)% -Konfidenzintervall in Grün überzeichnen:

fit2 <- lm(foodexp~xx, data=engel)
abline(h=summary(fit2)$coef[2], 
   col="green", lwd=3, lty=2)                                           #$
abline(h=summary(fit2)$coef[2] + qt(0.95,fit2$df)*summary(fit2)$coef[4], 
   col="green", lwd=3, lty=2)                                           #$
abline(h=summary(fit2)$coef[2] + qt(0.05,fit2$df)*summary(fit2)$coef[4], 
       col="green", lwd=3, lty=2)

und sie stimmen mit den gepunkteten roten Linien auf dem ursprünglichen Plot überein.

Geben Sie hier die Bildbeschreibung ein

user603
quelle