pdf("kuehl_06_04.pdf")


salinity <- rep(seq(0,12,6),each=6)
days <- rep(seq(14,28,7),each=2,times=3)
h2o <- c(2.2,3.3,5.0,5.7,13.2,12.4,
         3.7,4.5,5.9,7.2,9.4,11.0,
         2.8,3.4,4.5,5.9,7.6,8.3)

barley <- data.frame(salinity,days,h2o)
attach(barley)

salinity <- factor(salinity,levels=seq(0,12,6),ordered=T)
days <- factor(days,levels=seq(14,28,7),ordered=T)

tapply(h2o,list(salinity,days),mean)

barley.mod1 <- aov(h2o ~ salinity + days + salinity:days)
summary(barley.mod1)

summary(barley.mod1,split=list(salinity=list(linear=1,quadratic=2)))
summary(barley.mod1,split=list(days=list(linear=1,quadratic=2)))

dev.off()