如何在我的 ggplot2 2 因子条形图上放置标准误差条？

Question

我一直在尝试将当前的标准偏差条更改为条形图中的标准误差条。

这是均值+标准差的图：

EGG <- data.frame(type = c("this", "this", "that", "that"),
                  chemcon = c(10,11,12,13),
                  day = c("Monday", "Tuesday", "Monday", "Tuesday"))

EGG
#>    type chemcon     day
#> 1  this      10  Monday
#> 2  this      11 Tuesday
#> 3  that      12  Monday
#> 4  that      13 Tuesday

require(ggplot2)
aa <- aggregate(chemcon ~ day + type, data=EGG, FUN=mean)
bb <- aggregate(chemcon ~ day + type, data=EGG, FUN=sd)
cc <- merge(aa, bb, by=c("day", "type"))
colnames(cc)[3:4] <- c("mean", "sd")

ggplot(cc, aes(x = type, y = mean, fill = day))+
  geom_bar(stat="identity", position= "dodge") + #nb you can just use 'dodge' in barplots
  scale_fill_brewer(palette="Paired")+
  theme_minimal() +
  labs(x="", y="chemcon") +
  theme(panel.background = element_blank(),
        axis.line = element_line(colour = "black"),
        panel.grid=element_blank()) +
  geom_errorbar(aes(ymin = mean-sd,
                    ymax = mean+sd), 
                position = "dodge")

我尝试用“se”替换“sd”的聚合函数（FUN），没有成功，然后我尝试创建“se”也没有成功：

se = sd(Egg$chemcon) / sqrt(length(Egg$chemcon))

问题是因为我必须保留“聚合”函数，因为它在尝试表示双因素条形图时效果最好，但我还没有看到有人像这样在其他任何地方使用标准错误。有人可以帮助我了解我所缺少的吗？

Answer 1

从 plotrix 包中尝试 std.error()，如下所示

library(plotrix)
library(ggplot2)
EGG <- data.frame(type = rep(c("this", "this", "that", "that"),2),
                  chemcon = c(10,11,12,13,14,15,16,17),
                  day = rep(c("Monday", "Tuesday"),4))

EGG

aa <- aggregate(chemcon ~ day + type, data=EGG, FUN=mean)
bb <- aggregate(chemcon ~ day + type, data=EGG, FUN=sd)
ee <- aggregate(chemcon ~ day + type, data=EGG, FUN=std.error)  ##  from library(plotrix)
#ee <- aggregate(chemcon ~ day + type, data=EGG, FUN=sd(x)/sqrt(sum(!is.na(x))) )
# cc <- aggregate(chemcon ~ day + type, data=EGG, FUN = function(x) c(mean = mean(x), sd = sd(x), se = (sd(x)/sqrt(sum(!is.na(x)))) ))
cc <- merge(aa, ee, by=c("day", "type"))
colnames(cc)[3:4] <- c("mean", "se")


ggplot(cc, aes(x = type, y = mean, fill = day))+
  geom_bar(stat="identity", position= "dodge") + #nb you can just use 'dodge' in barplots
  scale_fill_brewer(palette="Paired")+
  theme_minimal() +
  labs(x="", y="chemcon") +
  theme(panel.background = element_blank(),
        axis.line = element_line(colour = "black"),
        panel.grid=element_blank()) +
  geom_errorbar(aes(ymin = mean-se,
                    ymax = mean+se), 
                position = "dodge")

如何在我的 ggplot2 2 因子条形图上放置标准误差条？

how to put standard error bars on my ggplot2 2-factor bar graphic?

r

ggplot2

errorbar