分位数回归图(quantreg 包)中的红线是什么?

机器算法验证 r 数据可视化 分位数回归
2022-03-26 18:06:18

plot.rq在 R的包中使用quantreg,我们可以绘制系数估计分布,并得到如下结果: 在此处输入图像描述

什么是红色虚线?广泛的谷歌搜索显示,中间的一个是所有 99 个估计值的平均值,但仍然不知道虚线红线。

1个回答

这是来自 quantreg vignette的 Engel 回归的再现。为了使绘图更具可读性,我关注第二个(斜率)参数:

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

我现在可以用绿色绘制通常的 OLS 估计值及其 (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)

它们与原始图上的红色虚线相匹配。

在此处输入图像描述