product <- read.table( "proddat.dat", header = T )
product$A <- ordered( product$A, c("L","M","H") )
product$B <- ordered( product$B, c("L","M","H") )
product$C <- ordered( product$C, c("L","H") )
product$D <- ordered( product$D, c("L","H") )
product$E <- factor( product$E )

product.fit1 <- aov( y ~ A + B + D + E, product )
product.fit2 <- aov( y ~ B + A*E, product )
prodres <- product
prodres.fit <- lm( y ~ B + E, product )
prodres$y <- resid( prodres.fit ) + lsmean( prodres.fit, factors=NULL )$pred
product.fit  <- aov( y ~ A + B + D + E + A:E, prodres )