Построение графиков со стандартными ошибками с использованием R
Я пытаюсь построить простой барплот со стандартными ошибками, и это сводит меня с ума. Я посмотрел несколько примеров и дошел до этого:
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
1 ответ
Решение
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
)
)