如何将剂量反应曲线拟合到分位数回归模型？

Question

我使用 quantreg 包为 0.10-0.90 分位数拟合了线性结果的分位数回归。现在我想将这些响应曲线拟合到这些模型中。我该怎么做？

我在下面添加了一个我的分位数回归模型的假示例。

###Quantile Regression###

#load library
library(quantreg)

#To perform quantile regression for the 0.10, 0.20, .,0.90 quantiles:
qs <- 1:9/10

###Example at time 0
#PM2.5 model  
PM25_0_modelq = rq(outcome~PM25_0+ covar1 + covar2+ covar3 + covar4 +covar5, data = mydataset, tau = qs) 

summary (PM25_0_modelq )    
qr_no2_0_plot<- plot(summary(PM25_0_modelq ), parm="PM25_0")