Пример из книги Мастицкого С.Э. Визуализация данных с помощью 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)
Не знаю, что произошло, то теперь строится так:
Подскажите, как это исправить?
Спасибо!