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

plot(unfactor(cloning$clone),cloning$titer,xaxt="n",
   xlab="",ylab="")
mtext("(a) clone identifier",1,2)
mtext("titer",2,2)
axis(1,seq(1,13,2),levels(cloning$code)[seq(1,13,2)])
axis(1,seq(2,12,2),levels(cloning$code)[seq(2,12,2)])

tmpar1 <- par(pty="s")
mplot(jitter(predict(cloning.fit)),cloning$titer,
   group=as.character(cloning$code),xlab="",ylab="")
mtext("(b) clone mean",1,2)
mtext("titer",2,2)
axis(1,1000,"1000")
abline(0,1,lty=2)

par(tmpar1)
par(tmpar)