library(statistics,T)
options(echo=T)
aovcat <- aov(Yield ~ Temp + Conc + Cat, catalyst)
aovcat
print.lm(aovcat)
summary(aovcat)
aovall <- aov(Yield ~ Temp * Conc * Cat, catalyst)
summary(aovall)
gunaov <- aov(R ~ M + P/T, gun)
summary(gunaov)
waov1 <- aov( allterms(wafer), wafer)
summary(waov1)
waov2 <- aov( allterms(wafer[,-c(9,10)]), wafer)
summary(waov2)
waov3 <- aov( allterms(wafer,
    response = log(post.dev) ~ .), wafer)
summary(waov3)
wmaov <- aov( allterms(wafer,
  response = cbind(pre.mean, post.mean) ~ .), wafer)
summary(wmaov)

cat("Tire Exp\n")
attach("g.exp")
contrasts(g.exp$Car)
g.aov <- aov(g.response~ Car + Brand )
summary(g.aov)
summary.aov(g.aov, onedf=T)
gm1.exp <- g.exp[-11, ]
detach()
attach("gm1.exp")
gm1.aov <- aov(g.response~ Car * Brand )
summary(gm1.aov)
detach()
cat("BHH 2^4 exp\n")
attach("bhh.exp")
bhh.aov <- aov(conversion ~ catchg*temp*psi*conc)
summary(bhh.aov)
design.apply(bhh.exp, mean, conversion, as.list(1:4))
tapply(conversion, list(catchg, temp, psi, conc), mean)
tapply(conversion, catchg, mean)

detach()

attach("dav.exp")
dav.exp
dav.aov <- aov(daviesresp~react+acidcon+acidamt+reactim+reactem+acidcon:acidamt+acidcon:reactim)
summary(dav.aov)
detach()

cat("glass exp\n")
attach("glass.exp")
glass.aov <- aov(light ~ phosphor*glass)
summary(glass.aov)
detach()


cat("Federer exp\n")
fed.aov <- aov(plants ~ reps + variety * treatment + Error(reps:variety),
	fed.page274)
summary(fed.aov)