tmpar <- par(mar=c(4,4,.5,.5),mfrow=c(1,2))

attach(tukey)
tmp1 <- mandel.plot(humid,y,temp,orient="both",
   xlab=c("(a) humid by temp","(b) temp by humid"),
   ylab=c("","expected yield"))
detach()

par(tmpar)