Các đường màu đỏ trong biểu đồ hồi quy lượng tử (gói quantreg) là gì?


9

Sử dụng plot.rqtrong quantreggói trong R, chúng ta có thể vẽ phân phối ước tính hệ số và nhận được một cái gì đó như thế này: nhập mô tả hình ảnh ở đây

Các đường màu đỏ chấm là gì? Googling mở rộng đã tiết lộ rằng một trong những trung bình là trung bình của tất cả 99 ước tính, nhưng vẫn không biết về đường màu đỏ chấm.


6
Googling mở rộng đã sai: đây là các ước tính OLS và SE của nó
user603

@ user603 (i) Bạn có tài liệu tham khảo cho việc này không? (ii) Ý bạn là CI? Đây là 95% hay 99%?

1
@AndreSilva: cảm ơn vì đã đưa câu hỏi này lên hàng đầu. Ngoài ra, cảm ơn vì đã chỉnh sửa các ô trong câu trả lời của tôi.
user603

Câu trả lời:


8

Đây là một bản tái tạo hồi quy Engel từ họa tiết quantreg . Để làm cho cốt truyện dễ đọc hơn, tôi tập trung vào tham số thứ hai (độ dốc):

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

Bây giờ tôi có thể overplot, màu xanh lục, ước tính OLS thông thường và khoảng tin cậy (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)

và chúng khớp với các đường màu đỏ chấm trên cốt truyện gốc.

nhập mô tả hình ảnh ở đây

Khi sử dụng trang web của chúng tôi, bạn xác nhận rằng bạn đã đọc và hiểu Chính sách cookieChính sách bảo mật của chúng tôi.
Licensed under cc by-sa 3.0 with attribution required.