rm(list = ls())
abt <- read.table(
url("http://stat.unm.edu/~fletcher/LLM/DATA/TAB3-1.DAT"),
#"C:\\E-drive\\Books\\ANREG2\\newdata\\TAB21-4.DAT",
#"C:\\E-drive\\Books\\LOGLIN3\\DATA\\ABORT.DAT",
#"C:\\E-drive\\Books\\LOGLIN3\\DATA\\TAB3-1.DAT",
		  sep="",col.names=c("R","S","A","O","y"))
attach(abt)
abt
#summary(abt)

r=factor(R)
o=factor(O)
s=factor(S)
a=factor(A)

m15 <- glm(y ~ r:s:a+r:s:o+r:o:a+s:o:a ,family = poisson)
m14 <- glm(y ~ r:s:a + r:s:o + r:o:a ,family = poisson)
m13 <- glm(y ~ r:s:a + r:s:o + s:o:a ,family = poisson)
m12 <- glm(y ~ r:s:a + r:o:a + s:o:a ,family = poisson)
m11 <- glm(y ~ r:s:a + r:s:o + o:a   ,family = poisson)
m10 <- glm(y ~ r:s:a + r:o:a + s:o ,family = poisson)
m9 <- glm(y ~ r:s:a + s:o:a + r:o  ,family = poisson)
m8 <- glm(y ~ r:s:a + r:o + s:o + o:a,family = poisson)
m7 <- glm(y ~ r:s:a + r:o + s:o ,family = poisson)
m6 <- glm(y ~ r:s:a + r:o + o:a ,family = poisson)
m5 <- glm(y ~ r:s:a + s:o + o:a ,family = poisson)
m4 <- glm(y ~ r:s:a + r:o   ,family = poisson)
m3 <- glm(y ~ r:s:a + s:o   ,family = poisson)
m2 <- glm(y ~ r:s:a + o:a   ,family = poisson)
m1 <- glm(y ~ r:s:a + o  ,family = poisson)

df=c(m15$df.residual,m14$df.residual,m13$df.residual,
m12$df.residual,m11$df.residual,m10$df.residual,
m9$df.residual,m8$df.residual,
m7$df.residual,m6$df.residual,m5$df.residual,
m4$df.residual,m3$df.residual,m2$df.residual,
m1$df.residual)
G2=c(m15$deviance,m14$deviance,m13$deviance,
m12$deviance,m11$deviance,m10$deviance,m9$deviance,
m8$deviance,m7$deviance,m6$deviance,m5$deviance,
m4$deviance,m3$deviance,m2$deviance,m1$deviance)

A2q=G2-(2*df)
modelm=c(df,G2,A2q)
model=matrix(modelm,15,3,dimnames =
      list(NULL,c("df","G2","A-q")))
model