求助：r语言分位数回归engel数据穷人与富人的比较时stepfun函数的作用是什么

 

attach(engel)
engel
z = rq(foodexp ~ income, tau=-1) 
z$sol     # 这里包含了每个分位点下的系数估计结果
x.poor = quantile(income, 0.1) # 10%分位点的收入
x.rich = quantile(income, 0.9) # 90%分位点的收入
ps = z$sol[1,]              # 每个分位点的tau值
ps
qs.poor = c( c(1,x.poor) %*% z$sol[4:5,] ) # 10%分位点的收入的消费估计值
qs.rich = c( c(1,x.rich) %*% z$sol[4:5,] )   # 90%分位点的收入的消费估计值
windows(10,5)
par(mfrow=c(1,2))                 # 把绘图区域划分为一行两列
plot(c(ps,ps),c(qs.poor,qs.rich),type="n",     # type=”n”表示初始化图形区域，但不画图
     xlab=expression(tau), ylab="quantile")
plot(stepfun(ps,c(qs.poor[1],qs.poor)), do.points=F,
     add=T)
plot(stepfun(ps,c(qs.poor[1],qs.rich)), do.points=F,
     add=T, col.hor="gra