построение гистограмм со стандартными ошибками с использованием 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))

Я бы очень хотел преодолеть это!

Спасибо,

Баз


person baz    schedule 16.03.2012    source источник
comment
?barplot2 из пакета gplots -- но и сюжет Google Dynamite   -  person Ben Bolker    schedule 16.03.2012


Ответы (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
  )
)
person Vincent Zoonekynd    schedule 16.03.2012