r - 散点图的非参数分位数回归曲线

Question

我用 IV=time 创建了一个散点图(多组 GRP) , DV=concentration .我想添加分位数回归曲线 (0.025,0.05,0.5,0.95,0.975)到我的情节。

顺便说一句，这就是我为创建散点图所做的:

attach(E)  ## E is the name I gave to my data
## Change Group to factor so that may work with levels in the legend
Group<-as.character(Group)
Group<-as.factor(Group)

## Make the colored scatter-plot
mycolors = c('red','orange','green','cornflowerblue')
plot(Time,Concentration,main="Template",xlab="Time",ylab="Concentration",pch=18,col=mycolors[Group])

## This also works identically
## with(E,plot(Time,Concentration,col=mycolors[Group],main="Template",xlab="Time",ylab="Concentration",pch=18))

## Use identify to identify each point by group number (to check)
## identify(Time,Concentration,col=mycolors[Group],labels=Group)
## Press Esc or press Stop to stop identify function

## Create legend
## Use locator(n=1,type="o") to find the point to align top left of legend box
legend('topright',legend=levels(Group),col=mycolors,pch=18,title='Group')

因为我在这里创建的数据是我的大数据的一个小子集，所以它可能看起来可以近似为矩形双曲线。但我还不想在我的自变量和因变量之间调用数学关系。

我想 nlrq来自包裹 quantreg可能是答案，但是当我不知道变量之间的关系时，我不明白如何使用该函数。

我从一篇科学文章中找到了这张图，我想做完全相同的图:


再次感谢您的帮助!

更新

Test.csv
有人指出我的样本数据不可重现。这是我的数据示例。

library(evd)
qcbvnonpar(p=c(0.025,0.05,0.5,0.95,0.975),cbind(TAD,DV),epmar=T,plot=F,add=T)

我也试过qcbvnonpar::evd，但曲线似乎不太平滑。

Answer 1

也许看看 quantreg:::rqss 来平滑样条曲线和分位数回归。
抱歉，示例数据不太好:

set.seed(1234)
period <- 100
x <- 1:100
y <- sin(2*pi*x/period) + runif(length(x),-1,1)


require(quantreg)
mod <- rqss(y ~ qss(x))
mod2 <- rqss(y ~ qss(x), tau=0.75)
mod3 <- rqss(y ~ qss(x), tau=0.25)
plot(x, y)
lines(x[-1], mod$coef[1] + mod$coef[-1], col = 'red')
lines(x[-1], mod2$coef[1] + mod2$coef[-1], col = 'green')
lines(x[-1], mod3$coef[1] + mod3$coef[-1], col = 'green')