分位数回归生成非单调分位数预测,例如Q49> Q50

Question

我希望分位数回归可以预测单调的分位数,即

但是,R中的quantreg包会生成完全没有意义的预测,请参见图:

这有什么理由吗？

以下示例.

library(quantreg)
library(ggplot2)
data(engel)

taus <- seq(0.01,0.99,0.01)
model_qr <- quantreg::rq(foodexp~income,tau=taus,data = engel)

test <- data.frame(income = 200,  foodexp=  300)
result <- data.frame(
            Forecast = as.numeric(predict(model_qr, test)),
            Quantile = taus *100 
            )

ggplot(result, aes(x = Quantile, y = Forecast)) +
    geom_point()