r - 分位数回归生成非单调分位数预测，例如 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()