[R] plotCI() with ggplot2

[ONKELINX, Thierry]

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.Onkelinx op 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: ggplot2 op googlegroups.com 
> [mailto:ggplot2 op googlegroups.com] Namens Sascha Vieweg
> Verzonden: dinsdag 8 maart 2011 16:34
> Aan: ggplot2 op googlegroups.com; r-help op 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, saschaview op 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 ggplot2 op googlegroups.com
> To unsubscribe: email ggplot2+unsubscribe op googlegroups.com
> More options: http://groups.google.com/group/ggplot2
>