pdf("rogaine.pdf")

rogaine <- read.fwf("C:\\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

rogaine1.mod <- aov(hair ~ trt*time + Error(subj))

summary(rogaine1.mod)

dev.off()