pdf("E:\\sta4211\\cat_repmeas.pdf") catzy <- read.table("http://www.stat.ufl.edu/~winner/sta4211/cat_repmeas.dat",header=T) attach(catzy) zylktrt <- factor(zylktrt) timeprd <- factor(timeprd) cat_id <- factor(cat_id) tapply(score,zylktrt,mean) tapply(score,timeprd,mean) tapply(score,list(zylktrt,timeprd),mean) interaction.plot(timeprd,zylktrt,score) cat.mod1 <- aov(score ~ zylktrt + zylktrt/cat_id + timeprd + zylktrt:timeprd) summary(cat.mod1) cat.mod2 <- aov(score ~ zylktrt*timeprd + Error(zylktrt/cat_id)) summary(cat.mod2) dev.off()