oocyte <- read.table( "oocyte.dat", header=T ) oocyte$asper <- 100 * asin( sqrt( oocyte$gv / oocyte$toocytes ) ) oocyte$h89 <- ordered( oocyte$h89 ) oocyte$fors <- ordered( oocyte$fors ) oocyte$date <- ordered( oocyte$date ) oocyte.fit <- aov( asper ~ h89 * fors, oocyte ) oocyte.afit <- aov( asper ~ h89 + fors, oocyte )