我一直在尝试在与以下数据非常相似的多级直方图上绘制标准误差线：
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])
}
我尝试使用此线程中描述的 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))
但结果总是在每年年中出现误差线。
我想知道我做错了什么。每个误差线都应位于数据中引用的月份/年份之上。