Hi > > Dear All > > I would like to produce interaction boxplots and this seems to work: > > par(mfrow=c(2,2)) > A=sample(rnorm(50,50,10)) > B=sample(rnorm(50,100,10)) > Test=merge(A,B,by=0)#by=0 where 0 is the row.names > TreatA=(gl(2,50,100,labels=c("High","Low"))) > TreatB=rep(gl(2,25,50,labels=c("High","Low")),2) > Newdata=data.frame(TreatA,TreatB,Test) > > bwplot(x~TreatA:TreatB,data=Newdata) > > However, I would prefer the X axis labels to be different, such that there > are two label rows (TreatA and TreatB) something like this: > > TreatA High High Low Low > TreatB High Low High Low > > any guidance on achieving this much appreciated.
Combination of axes=FALSE boxplot(x~TreatA:TreatB,data=Newdata, axes=F) and following axis and mtext commands seems to do what you want axis(2) axis(1, at= 1:4,labels=c("H", "H", "L", "L")) Regards Petr > > > Regards > > Tom > > > [[alternative HTML version deleted]] > > ______________________________________________ > 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. ______________________________________________ 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.