想要绘制分组数据的平均值并添加置信限

我需要计算和绘制平均值和 95% 置信限的数据，但我不知道我的代码哪里出了问题。数据按年份分组，我考虑使用 stat_smooth() 来绘制置信区间区域，但我无法将其付诸实践。如果我尝试，我会收到错误：

错误：geom_path：如果您使用点线或虚线、颜色、 尺寸和线型必须在整条线上保持不变

    library(ggplot2)
    library(dplyr) 
dat <- data.frame(cbind(rep(1:10,10),rnorm(100,600,sd=100),rnorm(100,400,sd=50),rnorm(100,700,sd=100)))
colnames(dat) <- c("year","S4","S5","S6")

    df<-data.frame(dat)
    ggplot(df)+
      geom_line(aes(x=year, y=mean(df$S4),colour=year)) +
      geom_line(aes(x=year, y=quantile(df$S4,0.95),colour=year),linetype="dotted") +
      geom_line(aes(x=year, y=quantile(df$S4,0.05),colour=year),linetype="dotted") 

你能用吗stat_summary:

ggplot(df, aes(x = year, y = S4)) +
  stat_summary(geom = "line", fun = mean, linetype = "solid") +
  stat_summary(geom = "ribbon", fun.data= mean_cl_normal, fun.args = list(conf.int=0.95), alpha=.1)

Plot

Edit:

对于多个变量/行，建议首先将数据采用“长”格式pivot_longer，然后绘制。也许可以对不同的变量使用不同的颜色。

df <- dat %>%
  pivot_longer(cols = c(S4, S5, S6), names_to = "variable", values_to = "value")

ggplot(df, aes(x = year, y = value, color = variable)) +
  stat_summary(geom = "line", fun = mean, linetype = "solid") +
  stat_summary(geom = "ribbon", fun.data= mean_cl_normal, fun.args = list(conf.int=0.95), alpha=.1)