rm(list = ls())
tense <- read.table(
url("http://stat.unm.edu/~fletcher/LLM/DATA/Example3-7-1.DAT"),
#"C:\\E-drive\\Books\\ANREG2\\newdata\\tab20-10a.dat",
#"C:\\E-drive\\Books\\LOGLIN3\\DATA\\Example3-7-1.dat",
		  sep="",col.names=c("y","Tn","Wt","Ms","Dr"))
attach(tense)
tense
#summary(tense)

W=factor(Wt)
M=factor(Ms)
D=factor(Dr)
T=factor(Tn)

m8 <- glm(y ~ (T + W + M + D)^4,family = poisson)
anova(m8)