我正在尝试绘制一个带有标准错误的简单条形图,这让我发疯。我确实查找了一些示例并得到了这样的结果:
rt5 <- data.frame(rtgrp=c(37.2,38.0,38.3,38.5,38.9),
mort=c(35,11,16,8,4),
se=c(0.08,0.01,0.005,0.01,0.02))
rt5
xvals=with(rt5,
barplot(mort,names.arg=rtgrp,
xlab="PTEMP_R group mean",ylab="%",ylim=c(0,max(mort+10+se))))
我正在尝试完成脚本的最后一行,但已经使用了很长时间:
with(rt5,
arrows(xvals,mort,xvals,mort+se,length=45,angle=90,code=3))
我真的很想克服这个!
谢谢,
巴兹