Диаграмма bar_plot

[VistaSV30]

Добрый день!
Пример из книги Мастицкого С.Э. Визуализация данных с помощью ggplot2 – М.: ДМК Пресс, 2017.

Код: Выделить всё

v <- read.delim("http://files.figshare.com/1360878/Dreissena.txt")
library(ggplot2); library(doBy)

# Создание векторов с минимальными, максимальными и средними значениями показателя
MIN <- summaryBy(Length ~ Month, FUN = min, 
                 data = subset(v, Lake == "Myastro"))
MAX <- summaryBy(Length ~ Month, FUN = max, 
                 data = subset(v, Lake == "Myastro"))
MEAN <- summaryBy(Length ~ Month, FUN = mean, 
                 data = subset(v, Lake == "Myastro"))
SE <- summaryBy(Length ~ Month, FUN = function(x) sd(x)/sqrt(length(x)), 
                 data = subset(v, Lake == "Myastro"))

# Создание фрэйма с данными для диаграммы
dat <- data.frame(MIN[1], MIN[2], MAX[2], MEAN[2], SE[2])
names(dat) <- c("Month", "Min", "Max", "Mean", "SE")

p <- ggplot(data = dat, aes(x = Month, ymin = Mean - SE, ymax = Mean + SE))
p + geom_bar(aes(y = Mean, fill = Month),
             stat = "identity") + geom_errorbar(width = 0.25) 

Сначала диаграмма строилась как надо:


Не знаю, что произошло, то теперь строится так:


Подскажите, как это исправить?
Спасибо!

[gamm]

забыть про ggplot, и спокойно все строить руками, используя plot

[VistaSV30]

Неужели больше ни кто не сталкивался с такой проблемой?

[gamm]

проблем с ggplot выше крыши, но они не имеют никакого отношения к анализу данных, это мурзилки. А для мурзилок существует фотошоп и корелдрав.

[nplatonov]

Код: Выделить всё

ggplot(data = dat, aes(x = Month))+
   geom_bar(aes(y = Mean, fill = Month),stat = "identity") +
   geom_errorbar(aes(ymin = Mean - SE, ymax = Mean + SE), width = 0.25)

[VistaSV30]

Спасибо еще раз!