ggplot errorbar位置多因素问题

Question

我一直在尝试在多级直方图数据上绘制标准误差条，这与以下非常相似：

mth<-rep(c("June","July","August"),length.out=15)
yr<-rep(c("1999","2000","2005","2009","2010"),each=3)
X<-rnorm(15,mean=200,sd=100)
lng<-rep(c(30,31,31),length.out=15)
vrnc<-rnorm(15,mean=740,sd=300)
df<-data.frame(mth,yr,vrnc,lng,X)
dfi<-dim(df)[1]
for(i in 1:dfi){
df$X.se[i]<-sqrt(df$vrnc[i]/df$lng[i])
}

我尝试使用this thread中描述的stat_summary方法放置错误栏。我使用stat_summary的代码如下所示：

ggplot(df,aes(x=yr,y=X,fill=mth))+
    stat_summary(fun.y=mean, geom="bar",position=position_dodge(1)) + 
    stat_summary(fun.ymin=min,fun.ymax=max,geom="errorbar",
        color="grey40",position=position_dodge(1), width=.2) +
    geom_errorbar(data=df,aes(ymin=X-X.se,ymax=X+X.se,position="dodge",width=.2))+
    scale_fill_discrete("mth")

我也尝试过使用此代码：

ggplot(df,aes(x=yr,y=X,fill=mth))+
    geom_bar(stat="identity",position="dodge")+
    geom_errorbar(data=df,aes(ymin=X-X.se,ymax=X+X.se,position="dodge",width=.2))

但结果总是在每年中期出现错误条。

我想知道我做错了什么。每个错误栏应位于数据中引用的月份/年之上。

Answer 1

最近不得不这样做。类似的东西：

ggplot(df,aes(x=yr,y=X,fill=mth))+
  geom_bar(stat="identity", position = position_dodge()) +
  geom_errorbar(aes(ymin=X-X.se, ymax=X+X.se), position = position_dodge(.9), width=.2)