Hi,

I have a data set with two continous variables that I want to plot MEANS (I 
am not intrerested in median values) on a double-y graph. I also have 2 
factors. I want the factor combinations plotted in different panes.

Dummy dataset:

mydata <- data.frame(factor1 = factor(rep(LETTERS[1:3], each = 40)),
                     factor2 = factor(rep(c(1:4), each = 10)),
                 y1 = rnorm(120, mean = rep(c(0, 3, 5), each = 40),
                           sd = rep(c(1, 2, 3), each = 20)),
                 y2 = rnorm(120, mean = rep(c(6, 7, 8), each = 40),
                            sd = rep(c(1, 2, 3), each = 20)))

I have tried plotrix(), but I everything is overlaid. Also, I am pretty sure 
that the means are wrong as I have assumed that the below calculates the 
mean of each factor level separately and not the mean per level factor 1 AND 
factor 2.

Is there a way of doing this in ggplot2?
I have also tried plotmeans() in the sciplot package, but was unsuccessful.

Sincerely
Anna Zakrisson

library(plotrix)
?brkdn.plot
par(family="serif",font=1)
brkdn.plot("y1", "factor1","factor2", data=mydata,
           mct="mean",md="sd",
           main="",
           cex=0.8,
           stagger=NA,
           xlab="factor1",
           ylab = "y1",
           dispbar=TRUE,
           type="b",pch=c(0),lty=1,col=par("fg"))
par(new=TRUE)
brkdn.plot("y2","factor1", "factor2",data=mydata,
           mct="mean",md="sd",
           main="",
           cex=0.8,
           stagger=NA,
           xlab="",
           ylab = "",
           dispbar=TRUE,
           type="b",pch=1,lty=2,col=par("fg"))
axis(4)
mtext("y2",side=4,line=3)
legend("topleft",col=c("black","black"),bty="n", 
lty=c(1:2),legend=c("y1","y2"))


Anna Zakrisson Braeunlich
PhD Student

Department of Ecology Environment and Plant Sciences
Stockholm University
Svante Arrheniusv. 21A
SE-106 91 Stockholm

Lives in Berlin.
For paper mail:
Katzbachstr. 21
D-10965, Berlin - Kreuzberg
Germany/Deutschland

E-mail: anna.zakrisson@su.se
Tel work: +49-(0)3091541281
Mobile: +49-(0)15777374888
LinkedIn: http://se.linkedin.com/pub/anna-zakrisson-braeunlich/33/5a2/51b

><((((º>`•. . • `•. .• `•. . ><((((º>`•. . • `•. .• 
`•. .><((((º>

	[[alternative HTML version deleted]]

