#source("c:/Work07/StatDesign07/Programs/R/Examples/Diet.R",print.eval=TRUE)# #Does anova for Sulphur data data<-read.table("c:/Work07/StatDesign07/DataSets/Diet.txt",sep = "",header=T) Diet<-as.character(data[,1]) Subject<-as.character(data[,2]) Time<-as.character(data[,3]) BP<-data[,4] aovdata <- data.frame(Diet,Subject,Time,BP) #----------This gives the SplitPlot anova table with no tests---------------- summary(aov(BP~Diet*Subject*Time,data=aovdata)) # #----------This gives the SplitPlot anova table with incorrect tests---------------- summary(aov(BP~Diet*Subject+Time*Diet,data=aovdata)) # #----------This gives the correct SplitPlot anova table ------------------------ summary(aov(BP~Diet+Error(Subject/Diet)+Time*Diet,data=aovdata)) # #--------------Oneway on averages------------------------------------------- Diet<-as.character(c(1,1,1,2,2,2,3,3,3,4,4,4)) BPavg<-2*c(129,128,125.5,128,130.5,132,118.5,120.5,123.5,145,144,146) aovdata <- data.frame(Diet,BPavg) summary(aov(BPavg~Diet,data=aovdata)) # #--------------Below the Line---------------------------------------------- Diet<-as.character(c(1,1,1,2,2,2,3,3,3,4,4,4)) BPdiff<-c(-12,-16,-7,-22,-11,-20,-9,-23,-17,-10,-6,-16) aovdata <- data.frame(Diet,BPdiff) summary(aov(BPdiff~Diet,data=aovdata))