ggplot errorbar位置多因素问题

ggplot errorbar position multi-factor problems

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

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))

但结果总是在每年年中出现误差线。

我想知道我做错了什么。每个误差线都应位于数据中引用的 month/year 之上。

最近不得不这样做。类似于:

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)