Гистограмма с нормальной кривой и барами ошибок в R

Question

    #import data
    data = diameters$V1
    error = .005 #mm

    #make histogram
    h <- hist(data, breaks = "FD", density = 10,
              col = "lightblue", xlab = "Diameter", main = "Overall") 

    # Make normal curve
    xfit <- seq(min(data), max(data), length = 40) 
    yfit <- dnorm(xfit, mean = mean(data), sd = sd(data)) 
    yfit <- yfit * diff(h$mids[1:2]) * length(data) 

    #Draw normal curve
    lines(xfit, yfit, col = "black", lwd = 2)

Выход:

Expectation:

Можно ли добавить гистограммы ошибок в гистограмму, используя значение +/- error, без каких-либо внешних библиотек?

Answer 1

Вы должны быть в состоянии нарисовать их с помощью функции arrows():

## Create a histogram from random data
> hist(sample(runif(100)))
> arrows(x0 = 0.15, y0 = 11, x1 = 0.15, y1 = 13, code = 3, length = 0.05, angle = 90)

x0 и x1 задают начальную и конечную координаты x (для прямой вертикальной линии сохраните их одинаковыми)

y0 и y1 задают начальную и конечную координаты y, например, длину линии, которую нужно нарисовать.

code = 3 говорит R нарисовать двухстороннюю «стрелку», angle = 90 делает «стрелку» плоской линией, а length = 0.05 указывает, насколько широкими должны быть полосы ошибок.

Подробнее см. ?arrows.