R 分位数回归图

我想知道如何将这些信息引入到绘图中，以显示0.025,0.25,0.50,0.75,0.975处的分位数曲线

library(quantreg)
data(engel)
attach(engel)
qs <- c(0.025,0.25,0.50,0.75,0.975)
cubic.rq <- rq(foodexp ~ poly(income, 3),qs, engel)
f <- coef(cubic.rq)

库（量子力学）
数据（恩格尔）
附件（恩格尔）
qs为了比较，我还包括了线性回归线（用于三阶多项式拟合）和数据点

library(tidyverse)
library(quantreg)
data(engel)

qs <- c(0.025,0.25,0.50,0.75,0.975)

ggplot(engel, aes(income, foodexp)) +
  geom_point(size=1, colour="grey70") +
  geom_quantile(quantiles=qs, formula=y ~ poly(x, 3), colour="red") +
  geom_smooth(method='lm', formula=y ~ poly(x,3), colour="blue", 
              se=FALSE, linetype="11") +
  theme_classic()

库（tidyverse）
图书馆（量子力学）
数据（恩格尔）
qs为了比较，我还包括了线性回归线（用于三阶多项式拟合）和数据点

library(tidyverse)
library(quantreg)
data(engel)

qs <- c(0.025,0.25,0.50,0.75,0.975)

ggplot(engel, aes(income, foodexp)) +
  geom_point(size=1, colour="grey70") +
  geom_quantile(quantiles=qs, formula=y ~ poly(x, 3), colour="red") +
  geom_smooth(method='lm', formula=y ~ poly(x,3), colour="blue", 
              se=FALSE, linetype="11") +
  theme_classic()

库（tidyverse）
图书馆（量子力学）
数据（恩格尔）
这可能就是你想要的：这可能就是你想要的：太好了，非常感谢，只是一个问题。您设置method=“lm”，我是否应该将其更改为method=“rq”？
geom_smooth
语句仅用于绘制蓝色的“标准”线性回归线，与分位数回归无关。这只是为了比较。分位数回归线由
geom_quantile
绘制。太好了，非常感谢，只是一个问题。您设置method=“lm”，我是否应该将其更改为method=“rq”？
geom_smooth
语句仅用于绘制蓝色的“标准”线性回归线，与分位数回归无关。这只是为了比较。分位数回归线由
geom_quantile
绘制。