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

Question

我有需要计算和绘制平均值和 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") 

Answer 1

你能用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)

情节

编辑：

如果有多个变量/行，建议先使用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)