运行中值曲线而不是 LOWESS

Question

我目前正在修改我提供的一些 R 代码以满足我的需要。

情况如下:

我们正在绘制约 200 行。然后他们使用 LOWESS 获得最佳拟合曲线。

现在看起来像这样:

lines(lowess(x.lowess, y.lowess), lwd = 3)

其中x.lowess和y.lowess是对应的坐标，各自在一个向量中，如:

> dput(x.lowess)
c(0.268309377138946, 0.511978097193703, 0.785763175906913, 0.974674880219028, ... )
> dput(y.lowess)
c(0.8, 0.5, 0.8, 0.5, ... )

我现在正在寻找运行中值曲线而不是 LOWESS 最佳拟合曲线。

有什么简单的方法/功能可以做到这一点吗？

有关情节的示例，请在 flickr 上查看此内容(抱歉，无法直接上传，我是新来的，不允许:) plot with lowess smoothing curve in red

Answer 1

生成一些示例数据:

set.seed(1001)
x <- runif(1000)
y <- runif(1000)
dat <- data.frame(x,y)

使用 quantreg 包求中位数作为 x 的函数:

library(quantreg)
q1 <- rq(y~x,data=data.frame(x,y))
xvec <- seq(0,1,length=101)
pq <- predict(q1,newdata=data.frame(x=xvec))

绘制基础图形:

plot(x,y,pch=".")
lines(lowess(x,y))
lines(xvec,pq,col=2)

或者使用ggplot2:

library(ggplot2)
theme_set(theme_bw())
qplot(x,y,data=dat,size=I(0.8),alpha=I(0.2))+
    geom_smooth(method="loess")+
    stat_quantile(quantiles=0.5,formula=y~x,colour="red")