tmpar <- par( mar=c(4.1,4.1,0,0) )
on.exit( { par( tmpar ); rm( tmpar ) } )

effect.plot( bacteria.bill )
mtext( c("temp.bact","bact.inoc","resid"), 1, 2, at = c(4,6,8) )