tmpar <- par( mar=c(3,2.5,0.5,0), mfrow=c(1,2) )

ci.plot( cloning.fit, lsm=cloning.lsm, rdf=df.resid(cloning.fit), sort=T,
   ylim = c(-150,2000), xlab = "", ylab = "" )
mtext( "(a) with common SD", 1, 2 )
tmpx <- seq( 2, 12, by=2 )
abline( h=0, lty=2 )

ci.plot( cloning.fit, lsm=cloning.means, rdf=cloning.means$N-1, sort=T,
   ylim = c(-150,2000), xlab = "", ylab = "" )
mtext( "(b) with group SDs", 1, 2)
tmpx <- seq( 2, 12, by=2 )
abline( h=0, lty=2 )

par( tmpar )