r - R 中 plot 和 qplot 中的不同 lowess 曲线

Question

我正在比较两个图表，每个图表都叠加了一条非参数低(w)ess曲线。问题是曲线看起来非常不同，尽管它们的参数(例如跨度)是相同的。

y<-rnorm(100)
x<-rgamma(100,2,2)
qplot(x,y)+stat_smooth(span=2/3,se=F)+theme_bw()
plot(x,y)
lines(lowess(y~x))

qplot() 生成的图中似乎有更多的曲率。如您所知，检测曲率在回归分析的诊断中非常重要，我担心如果我要使用 ggplot2，我会得出错误的结论。

你能告诉我如何在 ggplot2 中生成相同的曲线吗？

谢谢

Answer 1

或者，您可以使用 loess(..., degree=1)。这会产生与 lowess(...)

非常相似但不完全相同的结果

set.seed(1)    # for reproducibility
y<-rnorm(100)
x<-rgamma(100,2,2)
plot(x,y)
points(x,loess(y~x,data.frame(x,y),degree=1)$fitted,pch=20,col="red")
lines(lowess(y~x))

使用 ggplot

qplot(x,y)+stat_smooth(se=F,degree=1)+
  theme_bw()+
  geom_point(data=as.data.frame(lowess(y~x)),aes(x,y),col="red")