tab <- do.call(rbind, list(data1, data2, data3, data4)) etype <- rep(c(sd1, sd2, sd3, sd4), length(data1)) b <- barplot(tab, beside=T) arrows(unlist(b), unlist(tab) - etype, unlist(b), unlist(tab) + etype, code=3) ------------------------------------------------------------------- Jacques VESLOT
CNRS UMR 8090 I.B.L (2ème étage) 1 rue du Professeur Calmette B.P. 245 59019 Lille Cedex Tel : 33 (0)3.20.87.10.44 Fax : 33 (0)3.20.87.10.31 http://www-good.ibl.fr ------------------------------------------------------------------- Antje a écrit : > hello, > > I would like to create the following barplot: > > I have 4 different data sets (same length + stddev for each data point) > > data1 > sd1 > data2 > sd2 > data3 > sd3 > data4 > sd4 > > now, I'd like to plot in the following way: > > data1[1],data2[1],data3[1],data4[1] with it's sd-values side-by-side at > one x-axis label (named "position 1") and each bar in different colors. > > data1[2],data2[2],data3[2],data4[2] at the next x-axis label (named > "position 2") with the same color scheme > > and so on over the whole length. > > I managed to plot one set in the following way: > > par(mai=c(1.5,1,1,0.6)) > plotInfo <- barplot(data1, las=2, ylim = c(0,plotMax+1), ylab = > "Percentage") > arrows(plotInfo,data1,plotInfo, data1 + sd1, length=0.1, angle=90) > arrows(plotInfo,data1,plotInfo, data1 - sd1, length=0.1, angle=90) > > could anybody give me a help on this? > > Antje > > ______________________________________________ > R-help@stat.math.ethz.ch mailing list > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide http://www.R-project.org/posting-guide.html > and provide commented, minimal, self-contained, reproducible code. > ______________________________________________ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.