Я изо всех сил пытаюсь выяснить, какую ошибку я совершаю, пытаясь создать этот барплот со стандартными отклонениями. В результате я получаю строки ошибок, неправильно расположенные (по моему мнению, вертикально) в моих столбцах.
Мой набор данных
'data.Nema': 11 obs. of 6 variables:
$ STATION: Factor w/ 3 levels "FARO","ISLA D",..: 1 1 1 1 2 2 2 2 3 3 ...
$ SEASON : Factor w/ 4 levels "Summer15","Winter15",..: 1 3 2 4 1 3 2 4 1
$ N : int 5 3 5 3 5 3 5 3 5 3 ...
$ mean : num 2285 4419 795 1057 3694 ...
$ sd : num 511 1425 216 597 789 ...
$ se : num 256 1008 108 422 395 ...
мой код:
p4 <- ggplot(data=data.Nema, aes(x=STATION, y=mean, fill="Nematoda")) +
geom_bar(stat = "identity",position=position_dodge()) +
facet_wrap(facets = ~data.Nema$SEASON, ncol =3)+
theme_bw()+
theme(panel.border = element_blank(), panel.grid.major = element_blank(), panel.grid.minor = element_blank(),
panel.background = element_blank(), axis.line = element_line(colour = "black")) +
theme(strip.text.x = element_text(size = 8, colour = "black"))+
theme(axis.text=element_text(size = 8), axis.text.x=element_text(angle = 90, hjust = 1), axis.title.x = element_text(size = 8, vjust = 0.1),
axis.title.y = element_text(size = 9, vjust = 0.3), legend.position="bottom",legend.text=element_text(size=8), legend.title = element_text(size = 8)) + ylab(expression(paste("Individual numbers", " 10",cm^-2 ))) +
xlab("" ) +
geom_errorbar(aes(ymin = data.Nema$mean - sd,
ymax = data.Nema$mean + sd,
width = 0.1), position=position_dodge(0.9)) +
facet_wrap(facets = ~data.Nema$SEASON, ncol =3)
plot(p4)
Спасибо за вашу помощь