Средний алмазный участок на боксплоте в R

Question

У меня есть поле с точкой, представляющей среднее значение:

Это нормально, но я бы хотел показать доверительный интервал вокруг среднего значения, например, следующим образом:

Как мне наложить эти два графика в R?

Answer 1

Функцию, подобную следующей, легко написать.Введите x и y для среднего значения, верхнего предела, нижнего предела и ширины ромба.Вы также можете передать другие параметры, такие как lty, для типа линии и цвета col для сегментов.

diamondCI <- function(x, y, ul, ll, w = ul-ll, ...){
    hw <- w/2
    segments(x-hw, y, x+hw, y, ...)  # horizontal bar
    segments(x-hw, y, x, ul, ...) # left upper diag
    segments(x, ul, x+hw, y, ...) # right upper diag
    segments(x-hw, y, x, ll, ...) # left lower diag
    segments(x, ll, x+hw, y, ...) # right lwoer diag
    }

Попробуйте это с ...

plot(1,1)
diamondCI(1, 1, 1.2, 0.7, col = 'red', lwd = 3)