用置信区间绘制长期平均值

Question

我要使用30年的数据来绘制每日平均值。

这基本上是我想要的图形，但是我认为我没有以最有效的方式处理数据。我只是learned from a previous question，最终遇到了一些数据格式问题，而2月29日的数据非常困难。

我目前正在处理三列日期（date，ddmm，daymo），它们的格式略有不同。

> str(SFdischg)
'data.frame':   11932 obs. of  5 variables:
 $ date : Factor w/ 11932 levels "1/01/1985","1/01/1986",..: 97 4409 8697 9677 10069 10461 10853 11245 11637 489 ...
 $ ddmm : Factor w/ 366 levels "01-Apr","01-Aug",..: 1 13 25 37 49 61 73 85 97 109 ...
 $ year : int  1984 1984 1984 1984 1984 1984 1984 1984 1984 1984 ...
 $ cfs  : int  1500 1430 1500 1850 1810 1830 1850 1880 1970 1980 ...
 $ daymo: Date, format: "2018-04-01" "2018-04-02" "2018-04-03" ...

您可以看到daymo是同一年。我使用SFdischg$ddmm[is.na(as.Date(SFdischg$ddmm, format="%d-%b"))]创建了它，但是这意味着在2月29日获得前几年的数据时会存在NA。我可以忽略此数据，但是整个步骤似乎有些愚蠢和麻烦。

如何跳过所有这些难题以获得我想要的图，我如何使用date而不是daymo来创建长期均值和“当前”年图？我真的只想使用两列数据：date和cfs。

这是我创建情节的方式：

SF_pre <- subset(SFdischg, year <2016)
SF_year <- subset(SFdischg, year == 2016)

ggplot(SF_pre, aes(x=daymo, y=cfs)) +
stat_summary(fun.data = "mean_cl_boot", geom = "smooth", colour = "blue", se = TRUE) +
  geom_line(data=SF_year, aes(x=daymo, y=cfs), size=1) +
  # legend
  annotate("text", x=as.Date(dmy("12/2/2018")), y=3850, label = "30 year mean") +
  annotate("segment", 
           x=as.Date(dmy("1/1/2018")), xend = as.Date(dmy("15/1/2018")), 
           y=3850, yend = 3850, 
           colour = "blue", size = 1) +
  annotate("text", x=as.Date(dmy("26/1/2018")), y=3650, label = "2016") +
  annotate("segment", 
           x=as.Date(dmy("1/1/2018")), xend = as.Date(dmy("15/1/2018")), 
           y=3650, yend = 3650, 
           colour = "black", size = 1) +
  labs(x="date",y="discharge (cfs)") + 
  scale_x_date(date_breaks = "1 month", 
               date_labels = "%b") +
  theme(plot.title = element_text(color="#666666", face="bold", size=24, hjust=0)) +
  theme(axis.title = element_text(color="#666666", face="bold", size=16)) +
  theme_bw()

已编辑以添加数据

  year  cfs      daymo      date   ddmm
1 1984 1500 2018-04-01 1/04/1984 01-Apr
2 1984 1430 2018-04-02 2/04/1984 02-Apr
3 1984 1500 2018-04-03 3/04/1984 03-Apr
4 1984 1850 2018-04-04 4/04/1984 04-Apr
5 1984 1810 2018-04-05 5/04/1984 05-Apr
6 1984 1830 2018-04-06 6/04/1984 06-Apr