Apa garis merah dalam plot regresi kuantil (paket kuantreg)?

9

Menggunakan plot.rqdalam quantregpaket di R, kita dapat memplot distribusi estimasi koefisien, dan mendapatkan sesuatu seperti ini: masukkan deskripsi gambar di sini

Apa garis merah putus-putus? Googling ekstensif telah mengungkapkan bahwa yang di tengah adalah rata-rata dari semua 99 perkiraan, tetapi masih tidak tahu tentang garis merah putus-putus.


sumber
6
Googling yang luas salah: ini adalah perkiraan OLS dan pengguna SE
nya603
@ user603 (i) Apakah Anda punya referensi untuk ini? (ii) Maksud Anda CI? Apakah ini 95% atau 99%?
1
@AndreSilva: terima kasih telah mengajukan pertanyaan ini ke depan. Juga, terima kasih telah mengedit plot dalam jawaban saya.
user603

Jawaban:

8

Berikut adalah reproduksi regresi Engel dari sketsa kuantreg . Untuk membuat plot lebih mudah dibaca, saya fokus pada parameter kedua (kemiringan):

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

Sekarang saya dapat melakukan overplot, berwarna hijau, perkiraan OLS biasa dan interval kepercayaan (5-95)%:

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)

dan mereka cocok dengan garis merah putus-putus pada plot asli.

masukkan deskripsi gambar di sini

pengguna603
sumber