Как показать среднее и погрешности на гистограмме

Представим, что на иллюстрации нужно показать как результаты нескольких экспериментов, так и их среднее с погрешностями. Один из способов в R таков:

data=c(1,2,3,4,3,2,1) # вводим данные
barplot(data) # рисуем нашу иллюстрацию
abline(h=mean(data)) # добавляем линию средней
abline(h=mean(data)+sd(data)) # добавляем линии верхнего и нижнего стандартного отклонений
abline(h=mean(data)-sd(data))

Вуаля:

Вместо функций mean() и sd(), естественно, можно вставить либо конкретные числа, либо другие меры средней тенденции и погрешностей. Например, подключив пакет library(epicalc), можно визуализировать среднее с доверительным интервалом:

barplot(data)
abline(h=mean(data))
abline(h=c(ci(data)$lower95ci,ci(data)$upper95ci))