[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