library(quantreg)
v<-rnorm(100,0,0.5)
x<-runif(100,-0.1,5)
y=0
for(i in 1:100){
if(x>=0) y<-1+x+x*v
else y<-1+x+x*(-v)
}
plot(x,y,cex=1,type="p",xlab="x",ylab="QY(t/x)")
abline(lm(y~x),col="red")
taus<-c(0.01,0.3,0.5,0.7,0.99)
for(i in 1:length(taus)){
abline(rq(y~x,tau=taus),col="gray")}
这样的话：画出来的图是直线。