Как показать полосы ошибок на точечных графиках с помощью R

Question

Так что у меня есть данные, похожие на это

Value   Number
3   1.5
6   1.67
9   1.7
12  1.6
15  1.7
18  1.8
21  1.9
24  1.98
27  1.98
30  1.8
33  1.84
36  1.5
39  1.7
42  1.9
45  1.9
48  2.0
51  1.21
54  1.4
57  2.34
60  2.5
63  2.1
66  1.77

Как бы я сделал точечный график со стандартными барами ошибок.... Я посмотрел вверх, и похоже, что-то вроде этого

errbar(df$Value, df$Number, yplus, yminus, cap = 0.015,
       xlab= deparse(substitute(x)),
       ylab= deparse(substitute(y)))

Однако я незнаком с yplus, yminus? и разложить?? Есть ли другие способы сделать это? Я попытался использовать ggplot и загрузил его с помощью install.packages(ggplot2), но R повторял, что не может найти его, и попытался с этим кодом

> ggplot(data=dataset,aes(x=df$Value,y=df$Number,colour=Code,linetype=Group,ymin=Mean-SE,ymax=Mean+SE)) 
Error: could not find function "ggplot"
> + geom_line()
Error: could not find function "geom_line"
> + scale_x_continuous(breaks=c(1,2))
Error: could not find function "scale_x_continuous"
> + scale_linetype_manual(values=c(2,1))
Error: could not find function "scale_linetype_manual"
> + geom_point()
Error: could not find function "geom_point"
> + geom_errorbar(width=.1,position='dodge')
Error: could not find function "geom_errorbar"

Кроме того, как я должен идти о графике на 2 точечных диаграммах на одном графике с одинаковыми значениями для оси х, а также отображать полосы ошибок? Спасибо, ребята... любое предложение будет оценено

Answer 1

Вы видите ошибку, потому что ggplot2 пакет не был загружен. Добавление library(ggplot2) чтобы скрипт исправил это. Это, конечно, при условии ggplot2 пакет установлен. Если нет, используйте install.packages("ggplot2") чтобы исправить это.

Answer 2

В качестве альтернативы, вы можете использовать R gplots пакет, который довольно приятен и прост в обращении (на мой взгляд, по крайней мере).

# Required package
library(gplots)

# Sample data
x <- seq(3, 66, 3)
y <- c(1.5, 1.67, 1.7, 1.6,  1.7,  1.8,  1.9,  1.98,  1.98,  1.8,  1.84,  
       1.5,  1.7,  1.9,  1.9,  2.0,  1.21,  1.4,  2.34,  2.5,  2.1,  1.77)
xy <- data.frame(x, y)

# Plotting
plotCI(x, y, 
       uiw = .5, gap = 0, pch = 22, pt.bg = "green")

Конечно, вы должны заменить аргументы uiw а также liw с вашими ошибочными векторами.

Ура, Флориан