将不确定性波段添加到散点图中的平滑样条中

Question

我有一个带有光滑样条线的散点图

a<-rep(1:50,len=500)
b<-sample(0:5000,500)
c<-round(seq(0,600,len=500))
data_frame<-as.data.frame(cbind(a,b,c))
names(data_frame)<-c("ID","toxin_level","days_to_event")

plot(data_frame$days_to_event,data_frame$toxin_level, xlim=c(600,0),xlab="days before the event",ylab="Toxin level",type="p")
abline(v=0,col="red")

x <- data_frame$days_to_event
y <- data_frame$toxin_level

fit.sp = smooth.spline(y ~ x, nknots=20)
lines(fit.sp, col="blue")

这是生成的图

​

​

我想知道是否有可能以某种方式将置信区间添加到这条曲线上？

Answer 1

更新：使用scale_x_reverse更精确地匹配您的图形...

library(ggplot2)

ggplot(data_frame, aes(x = days_to_event, y = toxin_level)) + geom_point() +
  geom_vline(xintercept = 0, color = "red") + scale_x_reverse() + 
  xlab("Days before the event") + ylab("Toxin Level") + 
  geom_smooth(method = lm, se = TRUE)

这就给出了：

​

​

或者更符合你的问题：

ggplot(data_frame, aes(x = days_to_event, y = toxin_level)) + geom_point(shape = 1) +
  geom_vline(xintercept = 0, color = "red") + scale_x_reverse() + 
  xlab("Days before the event") + ylab("Toxin Level") + 
  geom_smooth(method = lm, se = TRUE, color = "blue", fill = "lightblue") + 
  theme_bw()

​

​