Добавление неопределенности полос на гладкой рейкой в ​​диаграмме рассеяния

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") 

Это Результирующий график

Я wondernig, если это возможно как-то добавить к этой кривой доверительные полосы? Я признателен, я бы хотел, чтобы он был прозрачным синим, но любой цвет, включая серый, в порядке.

Answer 1

Обновлен: с помощью scale_x_reverse, чтобы соответствовать вашему графику более точно ...

Как об этом с помощью ggplot2?

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()