Have a look at the examples of geom_errorbar() at the website: http://had.co.nz/ggplot2/geom_errorbar.html
Best regards, Thierry ---------------------------------------------------------------------------- ir. Thierry Onkelinx Instituut voor natuur- en bosonderzoek team Biometrie & Kwaliteitszorg Gaverstraat 4 9500 Geraardsbergen Belgium Research Institute for Nature and Forest team Biometrics & Quality Assurance Gaverstraat 4 9500 Geraardsbergen Belgium tel. + 32 54/436 185 thierry.onkel...@inbo.be www.inbo.be To call in the statistician after the experiment is done may be no more than asking him to perform a post-mortem examination: he may be able to say what the experiment died of. ~ Sir Ronald Aylmer Fisher The plural of anecdote is not data. ~ Roger Brinner The combination of some data and an aching desire for an answer does not ensure that a reasonable answer can be extracted from a given body of data. ~ John Tukey > -----Oorspronkelijk bericht----- > Van: ggpl...@googlegroups.com > [mailto:ggpl...@googlegroups.com] Namens Sascha Vieweg > Verzonden: dinsdag 8 maart 2011 16:34 > Aan: ggpl...@googlegroups.com; r-help@r-project.org > Onderwerp: plotCI() with ggplot2 > > Hello > > Currently, I plot some coefficients with some intervals using > function "plotCI()" (package "gplots") using the following code: > > (m1 <- matrix(0:5, nrow=2, byrow=T, dimnames=list(c("v1", "v2"), > c("lo", "m", "hi")))) > m2 <- m1 + 1 > library(gplots) > plotCI( > x=1:length(m1[, 1]), > pch="", > xlab="v1/v2", > xlim=c(1-.2, length(m1[, 1])+.2), > ylim=c(-.1, 7.1), yaxp=c(0, 7, 7), > ylab="Coefficients and intervals", > xaxt="n" > ) > plotCI( > x=1:length(m1[, 1])-.1, y=m1[, 2], li=m1[, 1], ui=m1[, 3], > pt.bg="white", cex=1, lty=1, type="o", gap=0, > add=T > ) > plotCI( > x=1:length(m2[, 1])+.1, y=m2[, 2], li=m2[, 1], ui=m2[, 3], > pch=10, pt.bg="white", cex=1, lty=1, type="o", gap=0, > add=T > ) > axis(1, at=1:length(m1[, 1]), labels=NA) text(1:length(m1[, > 1]), -.2, c("v1", "v2"), cex=1, xpd=T, > adj=.5, srt=0) > abline(h=seq(0, 7, 1), col=gray(.7), lty="dotted", lwd=.85) > legend(1.9, 7.1, > legend=c(".7", ".5"), > title="Condition:", > cex=.8 > ) > > I like to produce this plot with package "ggplot2", but I > have just discovered it and don't know where to start. It > would be great if somebody could give me a starting point > with code for that particular task. > > Thanks, *S* > > -- > Sascha Vieweg, saschav...@gmail.com > > -- > You received this message because you are subscribed to the > ggplot2 mailing list. > Please provide a reproducible example: http://gist.github.com/270442 > > To post: email ggpl...@googlegroups.com > To unsubscribe: email ggplot2+unsubscr...@googlegroups.com > More options: http://groups.google.com/group/ggplot2 > ______________________________________________ R-help@r-project.org 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.
