cnt=c(716,79,207,25,819,67,186,22) ii=c(1,1,1,1,2,2,2,2)
kk=c(1,2,1,2,1,2,1,2) jj=c(1,1,2,2,1,1,2,2) II=factor(ii)
JJ=factor(jj) KK=factor(kk) sv7 <- glm(cnt ~
II:JJ+II:KK+JJ:KK,family = poisson) sv6 <- glm(cnt ~
II:JJ+II:KK,family = poisson) sv5 <- glm(cnt ~
II:JJ+JJ:KK,family = poisson) sv4 <- glm(cnt ~
II:KK+JJ:KK,family = poisson) sv1 <- glm(cnt ~ II+JJ:KK,family
= poisson) sv2 <- glm(cnt ~ JJ+II:KK,family = poisson) sv3 <-
glm(cnt ~ II:JJ+KK,family = poisson) sv0 <- glm(cnt ~
II+JJ+KK,family = poisson)


tab7=c(7,df.residual(sv7), deviance(sv7),
deviance(sv7)-2*df.residual(sv7),
((deviance(sv0)-deviance(sv7))/(deviance(sv0))),
(1-(deviance(sv7)*df.residual(sv0)/
(deviance(sv0)*df.residual(sv7)))) ) 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(1,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(3,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)))) ) tab0=c(0,df.residual(sv0), deviance(sv0),
deviance(sv0)-2*df.residual(sv0),
((deviance(sv0)-deviance(sv0))/(deviance(sv0))),
(1-(deviance(sv0)*df.residual(sv0)/ (deviance(sv0)*
df.residual(sv0)))) )

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