大神们,我刚刚接触R软件,主要是想用R软件进行分位数回归问题的解决,根据网上找到的资料进行摸索。
根据网上的教程,输入命令如下
library(quantreg)
data(engel)
fit1 = rq(foodexp ~ income, tau = 0.5, data = engel)
fit1
summary(fit1)
r1 = resid(fit1)
c1 = coef(fit1)
summary(fit1, se = “nid”)
it1 = summary( rq(foodexp ~ income, tau = 2:98/100) )
fit2 = summary( rq(foodexp ~ income, tau = c(0.05,0.25,0.5,0.75,0.95)) )
windows(5,5)
plot(fit1)
windows(5,5)
plot(fit2)
但是输出的图像没有阴影部分(应该是置信带吧),但是第一次输出图像时是有过一次带置信带的,我真的在网上找不到答案,只能求助各位大神了。。。求大神带我飞