#pdf("kuehl_03_01.pdf") density <- rep(1:5,each=3) density <- 10*density grain <- c(12.2,11.4,12.4, 16.0,15.5,16.5, 18.6,20.2,18.2, 17.6,19.3,17.1, 18.0,16.4,16.6) grain1 <- data.frame(density,grain) attach(grain1) densityf <- factor(density,levels=seq(10,50,10),ordered=TRUE) grain1 tapply(grain,density,mean) tapply(grain,density,mean) grain.mod <- aov(grain~densityf) summary.lm(grain.mod) summary(grain.mod,split=list(densityf=list(linear=1, quadratic=2, cubic=3, quartic=4))) #dev.off()