Форумы GIS-Lab.info

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

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) 

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

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)