怎么用R做多元分位数回归？

求助怎么用R做多元分位数回归，一元的会做，多元的总是出问题，求大神知道如何做多元的分位数回归，求详细过程，不胜感激，五个论坛币送上！！！！！

大神们都在忙吗？还是我的论坛币太少？论坛币可以再加，求大神指教啊

fit1=summary(rq(y~x1+x2+x3,tau=2:98/100,data))
r1=resid(fit1)# 得到残差序列，并赋值为变量r1
c1 = coef(fit1)# 得到模型的系数，并赋值给变量c1
summary(fit1)# 显示分位数回归的模型和系数
plot(summary(fit1))

7秒记忆 发表于 2017-6-3 19:51
fit1=summary(rq(y~x1+x2+x3,tau=2:98/100,data))
r1=resid(fit1)# 得到残差序列，并赋值为变量r1
c1 = c ...

想请教一下，x1作为自变量，x2、x3...作为控制变量，这样控制变量是能发挥作用的吗