tmpar <- par( mar = c(3.1,2.1,.5,.5), mfrow = c(2,2) )

lsd.plot( oocyte.afit, factors = c("fors","h89"), xpos = 0.95, ypos = 10,
   xlim = c(0,1.15), xlab = "", ylab = "" )
mtext( "(a) additive fit", 1, 2 )
mtext( "fors", 1, 0, at = par( "usr" )[2], adj = 1 )

margin.plot( oocyte.afit, factors = c("fors","h89"), ypos = 10,
   xlab = "", ylab = "" )
mtext( "(b) additive margin", 1, 2 )

lsd.plot( oocyte.fit, factors = c("fors","h89"), xpos = 0.95, ypos = 40,
   xlim = c(0,1.15), xlab = "", ylab = "" )
mtext( "(c) interactive fit", 1, 2 )
mtext( "fors", 1, 0, at = par( "usr" )[2], adj = 1 )

margin.plot( oocyte.fit, factors = c("fors","h89"), ypos = 40,
   xlab = "", ylab = "" )
mtext( "(d) interactive margin", 1, 2 )

par( tmpar )