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

# order -- temp:bact:inoc temp:bact bact:inoc temp:inoc inoc bact temp
mqqnorm( bactroom.bill, group = c("BIT","BT","BI","IT","I","B","T"),
   fulltext = T, xlab = "", ylab = "" )
mtext( "(a) half-normal quantiles", 1, 2 )
mtext( "effects", 2, 2 )

effect.plot( bactroom.bill, xlim = c(.5,7.5), xaxt = "n", xlab = "", ylab = "" )
axis( 1, seq(1,7,by=2), c("T","I","IT","BIT") )
axis( 1, seq(2,6,by=2), c("B","BT","BI") )
mtext( "(b) terms", 1, 2 )
mtext( "MS adjusted effects", 2, 2 )