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

attach( bactroom )

ylab <- c("","bill length (mm)")
xlab <- c("(a)","(b)")
names( ylab ) <- names( xlab ) <- levels( run )

for ( i in levels( run ) )
{
   tmp <- run==i
   mplot( as.numeric( temp[tmp] ), bill[tmp], group=as.character( bact[tmp] ),
      type="b", xaxt="n", xlab="", ylab="", ylim = range( bill ) )
   mtext( as.character( unique( temp ) ), 1, 1, at=as.numeric( unique( temp ) ) )
   mtext( paste( xlab[i], i, "run" ), 1, 2 )
   mtext( ylab[i], 2, 2 )
}
   
detach( )

rm( tmp )

par( tmpar )