tmpar <- par(mar=c(3.6,3.1,0,0),mfrow=c(1,2))

attach(ntukey)
tukey.plot(lat,y,geno,xlim=c(95,105),
   xlab="",ylab="")
mtext("(a) lat marginals by geno",1,2)
mtext("height",2,2)
tmp <- tukey.plot(geno,y,lat,
   xlab="",ylab="")
mtext("(b) geno marginals by lat",1,2)
detach()