tmpar <- par( mfrow = c(1,2), mar = c(3.1,3.1,0,0) )
on.exit( { par( tmpar ); rm( tmpar ) } )

effect.plot( bacteria.rcbd, xlab="", ylab="", xlim = c(0.5,5.5), xaxt = "n",
   width = .2 )
axis( 1, 1:5, c("run","bact","temp","inter","error") )
axis( 1, c(3,5), c("temp","error") )
mtext( "(a) terms", 1, 2 )
mtext( "MS adjusted effects", 2, 2 )

margin.plot( bacteria.rcbd, factors = c("temp","bact"), xlab="", ylab="" )
mtext( "(b) temperature by bacteria", 1, 2 )
mtext( "bill length (mm)", 2, 2 )