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
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 at 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 at googlegroups.com > [mailto:ggplot2 at googlegroups.com] Namens Sascha Vieweg > Verzonden: dinsdag 8 maart 2011 16:34 > Aan: ggplot2 at googlegroups.com; r-help at 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 at 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 at googlegroups.com > To unsubscribe: email ggplot2+unsubscribe at googlegroups.com > More options: http://groups.google.com/group/ggplot2 >
The code to use would be something like limits <- aes(ymax = y + se, ymin = y - se) # set error bar dimensions ggplot(yourdataframe, aes(x = x, y = y, colour = z)) + geom_point() + geom_errorbar(limits) On Tuesday, March 8, 2011 at 9:34 AM, Sascha Vieweg wrote:> 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@gmail.com > > ______________________________________________ > 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. >[[alternative HTML version deleted]]