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)
m6 <- glm(y ~ T:W + T:M:D + W:M:D,family = poisson)
fitted(m6)
c(fitted(m6)[1]/fitted(m6)[9],fitted(m6)[2]/
fitted(m6)[10],fitted(m6)[3]/fitted(m6)[11],
fitted(m6)[4]/fitted(m6)[12],fitted(m6)[5]/
fitted(m6)[13],fitted(m6)[6]/fitted(m6)[14],
fitted(m6)[7]/fitted(m6)[15],fitted(m6)[8]/
fitted(m6)[16],fitted(m6)[1]/fitted(m6)[9])