pdf("woolfrict.pdf")

wf1 <- read.fwf("C:\\data\\woolfrict.dat",
width=c(8,8,8,8),
col.names=c("run", "trt", "revs", "shrink"))

ftrt <- factor(Trt, levels=1:4)
levels(ftrt) <- c("Untreated", "AC15sec", "AC4min", "AC15min")

woolfrict <- data.frame(ftrt, run=factor(wf1$run), 
revs=factor(wf1$revs,ordered=TRUE), wf1$shrink)

woolfrict.sp <- aov(shrink ~ ftrt*revs + Error(run))

summary(woolfrict.sp)

summary(wolfrict.sp,split=list(revs=list(linear=1, quadratic=2,
cubic=3, quartic=4, quintic=5, sextic=6)))

dev.off()