support for quantile regression

[yuryzablotski]

Describe the solution you'd like
I would like to be able to report results of quantile regression similarly, like it's currently possible for OLS / GLM. A single quantile would be great, all quantiles at once would be terrific!

How could we do it?
Should be very similar to OLS with the use of "tau" argument, where tau is the quantile.

Code example

library(ISLR)
library(quantreg)

lr <- lm(wage ~ jobclass + age, Wage)
qr10 <- rq(wage ~ jobclass + age, Wage, tau = 0.1)
qr50 <- rq(wage ~ jobclass + age, Wage, tau = 0.5)
qr90 <- rq(wage ~ jobclass + age, Wage, tau = 0.9)

qr_all <- rq(wage ~ jobclass + age, Wage,
tau = seq(.05, .95, by = 0.05))

summary(qr_all) %>% plot()