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\\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)

sv6 <- glm(y ~ r:s:o + o:a ,family = poisson)
sv5 <- glm(y ~ r:o + s:o + o:a ,family = poisson)
sv4 <- glm(y ~ r + s + o:a ,family = poisson)
sv3 <- glm(y ~ r:a + s + o:a ,family = poisson)
sv2 <- glm(y ~ r + s:o + o:a ,family = poisson)
sv1 <- glm(y ~ r:o + s + o:a ,family = poisson)
sv0 <- glm(y ~ r + s + o + a ,family = poisson)


tab6=c(6,df.residual(sv6), deviance(sv6),
deviance(sv6)-2*df.residual(sv6),
((deviance(sv0)-deviance(sv6))/(deviance(sv0))),
(1-(deviance(sv6)*df.residual(sv0)/
(deviance(sv0)* df.residual(sv6))))  )
tab5=c(5,df.residual(sv5), deviance(sv5),
deviance(sv5)-2*df.residual(sv5),
((deviance(sv0)-deviance(sv5))/(deviance(sv0))),
(1-(deviance(sv5)*df.residual(sv0)/
(deviance(sv0)* df.residual(sv5))))  )
tab4=c(4,df.residual(sv4), deviance(sv4),
deviance(sv4)-2*df.residual(sv4),
((deviance(sv0)-deviance(sv4))/(deviance(sv0))),
(1-(deviance(sv4)*df.residual(sv0)/
(deviance(sv0)* df.residual(sv4))))  )
tab1=c(3,df.residual(sv1), deviance(sv1),
deviance(sv1)-2*df.residual(sv1),
((deviance(sv0)-deviance(sv1))/(deviance(sv0))),
(1-(deviance(sv1)*df.residual(sv0)/
(deviance(sv0)* df.residual(sv1))))  )
tab2=c(2,df.residual(sv2), deviance(sv2),
deviance(sv2)-2*df.residual(sv2),
((deviance(sv0)-deviance(sv2))/(deviance(sv0))),
(1-(deviance(sv2)*df.residual(sv0)/
(deviance(sv0)* df.residual(sv2))))  )
tab3=c(1,df.residual(sv3), deviance(sv3),
deviance(sv3)-2*df.residual(sv3),
((deviance(sv0)-deviance(sv3))/(deviance(sv0))),
(1-(deviance(sv3)*df.residual(sv0)/
(deviance(sv0)* df.residual(sv3))))  )

t(matrix(c(tab6,tab5,tab4,tab3,tab2,tab1),6,6))