row <- rep(1:5,each=5) col <- rep(1:5,times=5) trt <- c(5,1,3,2,4,3,4,2,5,1,2,3,4,1,5, 1,2,5,4,3,4,5,1,3,2) yld <- c(59.45,47.28,54.44,50.14,59.45, 55.16,60.89,56.59,60.17,48.71, 44.41,53.72,55.87,47.99,59.45, 42.26,50.14,55.87,58.74,55.87, 60.89,59.45,49.43,59.45,57.31) wheat <- data.frame(row,col,trt,yld) attach(wheat) row <- factor(row) col <- factor(col) trt <- factor(trt) tapply(yld,row,mean) tapply(yld,col,mean) tapply(yld,trt,mean) wheat.mod1 <- aov(yld ~ trt + row + col) summary(wheat.mod1) TukeyHSD(wheat.mod1,"trt")