+\subsubsection{Anova list}\r
+\r
+<<aovlist>>=\r
+Block <- gl(8, 4)\r
+A <- factor(c(0,1,0,1,0,1,0,1,0,1,0,1,0,1,0,1,0,1,0,1,\r
+ 0,1,0,1,0,1,0,1,0,1,0,1))\r
+B <- factor(c(0,0,1,1,0,0,1,1,0,1,0,1,1,0,1,0,0,0,1,1,\r
+ 0,0,1,1,0,0,1,1,0,0,1,1))\r
+C <- factor(c(0,1,1,0,1,0,0,1,0,0,1,1,0,0,1,1,0,1,0,1,\r
+ 1,0,1,0,0,0,1,1,1,1,0,0))\r
+Yield <- c(101, 373, 398, 291, 312, 106, 265, 450, 106, 306, 324, 449,\r
+ 272, 89, 407, 338, 87, 324, 279, 471, 323, 128, 423, 334,\r
+ 131, 103, 445, 437, 324, 361, 302, 272)\r
+aovdat <- data.frame(Block, A, B, C, Yield)\r
+\r
+old <- getOption("contrasts")\r
+options(contrasts = c("contr.helmert", "contr.poly"))\r
+(fit <- aov(Yield ~ A*B*C + Error(Block), data = aovdat))\r
+class(fit)\r
+summary(fit)\r
+options(contrasts = old)\r
+@\r
+\r
+\p\r
+\r
+<<xtableaovlist, results='asis'>>=\r
+xtable(fit)\r
+@\r
+\r
+\r