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. 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.