#source("c:/Work07/StatDesign07/Programs/R/Examples/Sulphur.R",print.eval=TRUE)# #Does anova for Sulphur data data<-read.table("c:/WorkHome07/StatDesign/DataSets/Sulphur.txt",sep = "",header=T) Nitrogen<-as.character(data[,1]) Sulphur<-as.character(data[,2]) Yield<-data[,3] aovdata <- data.frame(Nitrogen,Sulphur,Yield) #----------This gives the twoway anova table---------------- summary(aov(Yield~Sulphur*Nitrogen,data=aovdata))