Столбики погрешностей неправильно отображаются на гистограмме

У меня возникли проблемы с правильным отображением моих баров ошибок. Вот мой фрейм данных:

> d
         Date Plate.1 Plate.2 Plate.3 Plate.4 Plate.5 Plate.6 Plate.7 Plate.8 Sum Average        SE Treatment
1  2014-10-15       2       2       0       5       8      11      11       0  39   4.875 1.6304852         1
2  2014-11-12       5       4      11       6       7       2       9       6  50   6.250 0.9955257         1
3  2014-12-11       1       0       0       0       0       0       0       1   2   0.250 0.1636634         1
4  2015-02-11       0       0       4       0       2       1       0       1   8   1.000 0.5000000         1
5  2015-03-09       0       0       0       0       0       0       0       0   0   0.000 0.0000000         1
6  2014-10-15      30      22      24      19      10      16      19      42 182  22.750 3.4369318         2
7  2014-11-12       8       5       4       9      12      23      14      10  85  10.625 2.1207942         2
8  2014-12-11       0       5      21      19       2       9       4       0  60   7.500 2.9215945         2
9  2015-02-11      16      11       5       8       5      17       0       0  62   7.750 2.3126207         2
10 2015-03-09       0       0       0       1       0       0       0       0   1   0.125 0.1250000         2

Я использую этот код для создания сгруппированного гистограммы:

ggplot(d, aes(x=Date, y=Average, fill = Treatment)) +
  geom_bar(position="dodge", stat="identity") +
  geom_errorbar(aes(ymin=Average-SE, ymax=Average+SE),
                width=.2,                    
                position = position_dodge(.9))

По какой-то причине полосы ошибок отображаются неправильно, и график выглядит так:

введите здесь описание изображения

Надеюсь, решение простое, но оно ускользает от меня. Заранее благодарю за любую помощь!


r bar-chart ggplot2 standard-error
person opalfruits    schedule 02.05.2015    source источник

Ответы (1)


arrow_upward
4
arrow_downward

Обновлено в свете вашего комментария - для этого типа графика вы можете присвоить значение position_dodge переменной и использовать его для всех последующих графиков. Я нашел значение 20 методом проб и ошибок; если другие данные, с которыми вы работаете, сильно различаются в числовом выражении, вам может потребоваться точная настройка.

d <- data.frame(Date=c("2014-10-15", "2014-11-12", "2014-12-11",
                "2015-02-11", "2015-03-09", "2014-10-15",
                "2014-11-12", "2014-12-11", "2015-02-11",
                "2015-03-09"),
                Average=c(4.875, 6.250, 0.25, 1, 0,
                22.75, 10.625, 7.5, 7.75, 0.125),
                Treatment=c(rep(1, 5), rep(2, 5)),
                SE=c(1.6304852, 0.9955257, 0.1636634, 0.5, 0,
                3.4369318, 2.1207942, 2.9215945, 2.3126207, 0.125))
d$Date <- as.Date(d$Date)
d$Treatment <- factor(d$Treatment)
library(ggplot2)
pd1 <- 20
ggplot(d, aes(x=Date, y=Average, fill=Treatment)) + 
    geom_bar(width=pd1,
             position=position_dodge(pd1),
             stat="identity") +
    geom_errorbar(aes(ymin=Average-SE, ymax=Average+SE),
                  width=pd1 / 2,                    
                  position=position_dodge(pd1))

давая:

введите здесь описание изображения

Я думаю, что проблема, с которой вы столкнулись, заключается в том, что использование position="dodge" позволяет ggplot2 вычислить «оптимальное» уклонение. Однако вам нужно явно увидеть/указать значение, чтобы передать его в качестве аргумента geom_errorbar.

Установка width в geom_bar на то же значение, что и position_dodge в обоих geom_bar и geom_errorbar (здесь pd1 <- 20), преодолевает это.

person dardisco    schedule 02.05.2015