[R] plotCI() with ggplot2

From: Sascha Vieweg <saschaview_at_gmail.com>
Date: Tue, 08 Mar 2011 16:34:15 +0100 (CET)


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_at_gmail.com

______________________________________________
R-help_at_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.
Received on Tue 08 Mar 2011 - 15:36:54 GMT

Archive maintained by Robert King, hosted by the discipline of statistics at the University of Newcastle, Australia.
Archive generated by hypermail 2.2.0, at Tue 08 Mar 2011 - 16:50:19 GMT.

Mailing list information is available at https://stat.ethz.ch/mailman/listinfo/r-help. Please read the posting guide before posting to the list.

list of date sections of archive