如何在ggplot2中为我的条形图添加SE错误栏?

Sté*_*ane 4 r bar-chart ggplot2 errorbar

我用ggplot2制作了一个简单的条形图,比较了2种昆虫的雄性和雌性的平均寿命(年龄).我的代码看起来像这样,"数据集"就是我的数据集......

    gplot(dataset, aes(Species, Age, fill=Sex))+
stat_summary(fun.y = mean, geom = "bar", position = "dodge")+
scale_fill_manual(values = c("Grey25", "Grey"))+
theme(legend.title = element_blank())+
scale_y_continuous(limits = c(0,15))
Run Code Online (Sandbox Code Playgroud)

我尝试使用以下代码手动输入平均值±SE来设置误差条的限制.为简单起见,我们假设物种1的平均值= 10且SE = 0.5.

geom_errorbar(aes(ymin=9.5, ymax=10.5),width=.2,position=position_dodge(.9))
Run Code Online (Sandbox Code Playgroud)

此代码确实有效,但它为我的绘图中的每个条设置了相同的误差条.

如何为我的绘图中的每个条添加等于相应SE的误差线?

我对ggplot和R一般都很新,所以欢迎任何帮助/建议.

Tho*_*s K 15

你不需要添加stat_summary(geom = "errorbar", fun.data = mean_se, position = "dodge")到你的情节:

library(ggplot2)

ggplot(diamonds, aes(cut, price, fill = color)) +
  stat_summary(geom = "bar", fun.y = mean, position = "dodge") +
  stat_summary(geom = "errorbar", fun.data = mean_se, position = "dodge")
Run Code Online (Sandbox Code Playgroud)

如果您希望事先计算这些值,可以这样做:

library(tidyverse)
pdata <- diamonds %>% 
  group_by(cut, color) %>% 
  summarise(new = list(mean_se(price))) %>% 
  unnest(new)


pdata %>% 
  ggplot(aes(cut, y = y, fill = color)) +
  geom_col(position = "dodge") +
  geom_errorbar(aes(ymin = ymin, ymax = ymax), position = "dodge")
Run Code Online (Sandbox Code Playgroud)