# NOT RUN {
# Esterase data
data(esterase)
# Fit quantiles 0.25 and 0.75
fit <- rq.counts(Count ~ Esterase, tau = 0.5, data = esterase, M = 50)
cbind(fit$fitted.values, predict(fit, type = "response"))
# }
Run the code above in your browser using DataLab