我正试图用标准错误绘制一个简单的条形图,这让我抓狂。我确实查了一些例子,结果是:
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))
我真的很想克服这个!
谢谢,
Baz
length
是箭头的大小(错误条的宽度):45比你的情节要大得多。较小的值应该有效。
with(rt5,
arrows(
xvals,mort,xvals,mort+se,
length=.3, angle=90, code=3,
# Change the colour and line width, to see the error bars
col="navy", lwd=5
)
)