proddes <- read.table( "proddes.dat", header = T )

proddes.add <- aov( sweet ~ A+B+C+D+E+F, proddes )
proddes.int <- aov( sweet ~ A+B+C+D+E+F+A:B+A:C+A:F+B:F+B:C+F:C+A:E, proddes )