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


res1 <- aov(Data~Treat*Sex,d)
summary(res1)