d <- read.table("Ex12.2.txt", header=T)
str(d)


res1 <- aov(Data~Rep+Variety*Nitrogen,d)
summary(res1)