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

什么是红色虚线?广泛的谷歌搜索显示,中间的一个是所有 99 个估计值的平均值,但仍然不知道虚线红线。
plot.rq在 R的包中使用quantreg,我们可以绘制系数估计分布,并得到如下结果:

什么是红色虚线?广泛的谷歌搜索显示,中间的一个是所有 99 个估计值的平均值,但仍然不知道虚线红线。
这是来自 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)
它们与原始图上的红色虚线相匹配。
