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

