pdf("rogaine.pdf") rogaine <- read.fwf("http://www.stat.ufl.edu/~winner/data/rogaine2.dat", width=c(8,8,8,8), col.names=c("trt", "subj", "time", "hair")) rogaine2 <- subset(rogaine,rogaine$time>0) rogaine1 <- data.frame(trt=factor(rogaine2$trt), subj=factor(rogaine2$subj), time=factor(rogaine2$time), hair=rogaine2$hair) attach(rogaine1) rogaine1.mod1 <- aov(hair ~ trt + trt/subj + time + trt:time) summary(rogaine1.mod1) rogaine1.mod2 <- aov(hair ~ trt*time + Error(subj)) summary(rogaine1.mod2) dev.off()