[R] plotCI() with ggplot2

Sascha Vieweg saschaview at gmail.com
Tue Mar 8 16:34:15 CET 2011


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



More information about the R-help mailing list