我正在尝试在ggplot2中使用stat_summary将误差线添加到折线图中，但是当我对图形进行分面时，它将不起作用
我的资料：
date week year location imidacloprid block wickhami virescens sexta
1 15-May 1 2015 kinston tp 1 0 0 0
2 15-May 1 2015 kinston gh 1 0 0 0
3 15-May 1 2015 kinston utc 1 0 0 0
4 15-May 1 2015 kinston gh 2 0 0 0
5 15-May 1 2015 kinston utc 2 0 0 0
6 15-May 1 2015 kinston tp 2 0 0 0
'data.frame': 576 obs. of 9 variables:
$ date : Factor w/ 27 levels "1-Jul","12-Jun",..: 4 4 4 4 4 4 4 4 4 4 ...
$ week : Factor w/ 12 levels "1","2","3","4",..: 1 1 1 1 1 1 1 1 1 1 ...
$ year : Factor w/ 2 levels "2015","2016": 1 1 1 1 1 1 1 1 1 1 ...
$ location : Factor w/ 2 levels "kinston","rocky mount": 1 1 1 1 1 1 1 1 1 1 ...
$ imidacloprid: Factor w/ 3 levels "gh","tp","utc": 2 1 3 1 3 2 3 2 1 2 ...
$ block : Factor w/ 4 levels "1","2","3","4": 1 1 1 2 2 2 3 3 3 4 ...
$ wickhami : num 0 0 0 0 0 0 0 0 0 0 ...
$ virescens : num 0 0 0 0 0 0 0 0 0 0 ...
$ sexta : num 0 0 0 0 0 0 0 0 0 0 ...
汇总用于绘图的数据：
wickhami_sum = summarySE(bug_subset_final,
measurevar="wickhami",
groupvars=c("imidacloprid","week","year"))
imidacloprid week year N wickhami sd se ci
1 gh 1 2015 8 0.0000 0.0000000 0.00000000 0.0000000
2 gh 1 2016 8 0.0000 0.0000000 0.00000000 0.0000000
3 gh 2 2015 8 0.0000 0.0000000 0.00000000 0.0000000
4 gh 2 2016 8 0.0000 0.0000000 0.00000000 0.0000000
5 gh 3 2015 8 0.0000 0.0000000 0.00000000 0.0000000
6 gh 3 2016 8 0.1250 0.2314550 0.08183171 0.1935012
7 gh 4 2015 8 0.0000 0.0000000 0.00000000 0.0000000
8 gh 4 2016 8 0.5000 0.4629100 0.16366342 0.3870025
9 gh 5 2015 8 0.5000 0.3779645 0.13363062 0.3159862
下面的代码对我来说没有问题，并生成了两年数据合并后的折线图，并通过
stat_summary生成了误差线ggplot(wickhami_sum, aes(x=week, y=wickhami,linetype=imidacloprid,group=imidacloprid))+
stat_summary(fun.data=mean_se,geom="errorbar",width=.2,color="black",position=position_dodge(0.2))+
stat_summary(fun.y=mean,geom="line",position=position_dodge(0.2))
但是，当我尝试按年份（如下所示）对数据进行分面时，我无法获得stat_summary来生成错误栏并获得以下错误消息
ggplot(wickhami_sum, aes(x=week, y=wickhami,linetype=imidacloprid,group=imidacloprid))+
stat_summary(fun.y=mean,geom="line",position=position_dodge(0.2))+facet_grid(year~.)+
stat_summary(fun.data=mean_se,geom="errorbar",width=.2,color="black",position=position_dodge(0.2))
Warning message:
Removed 72 rows containing missing values (geom_errorbar).
我尝试扩展y轴的范围/界限以包括误差线，但是我仍然收到相同的警告消息，但没有误差线。
我希望使用stat_summary为多面图生成误差线，而不必再次计算标准误差。感谢您理解为什么多面不允许stat_summary正常运行或我未正确执行操作时的任何帮助。
最佳答案
我想这是正在发生的情况：未刻面绘图中每周有两行数据，但是每面刻面面板中每周只有一行数据，导致标准误差计算返回NA。 stat_summary用于未汇总的数据，并在内部进行数据汇总。将bug_subset_final与stat_summary一起使用，或切换到geom_errorbar继续使用wickhami_sum。详细信息如下。
您已经预先汇总了数据，但是stat_summary用于处理原始数据并在内部计算汇总值。在您传递给ggplot的摘要数据帧wickhami_sum中，每周有两行，2015年每星期一行，2016年每星期一行。每周和每年的所有数据均已折叠至通过汇总操作，每周和每年都有一行。
因此，在未刻面的图中，每周都有两行数据可供stat_summary使用。但是在分面图中，它试图从单个观察值计算标准误差，这可能返回NA，因此没有绘制任何东西。即使在未刻面的图中，也将根据每年的两个平均值来计算误差线，这也不是您想要的。
而是继续使用wickhami_sum，但不要使用stat_summary：
geom_errorbar(aes(ymin = wickhami - se, ymax=wickhami + se))
或者，将原始数据（看起来像是
bug_subset_final）与stat_summary一起使用：ggplot(bug_subset_final, aes(x=week, y=wickhami)) +
stat_summary(fun.data=mean_se, geom="errorbar)`.
关于r - 刻面时使用stat_summary为折线图生成误差线的问题，我们在Stack Overflow上找到一个类似的问题：https://stackoverflow.com/questions/47724507/