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),
sum(residuals(sv7,type="pearson")^2),
deviance(sv7),1-pchisq(deviance(sv7),df.residual(sv7)))
tab6=c(6,df.residual(sv6),
sum(residuals(sv6,type="pearson")^2),
deviance(sv6),1-pchisq(deviance(sv6),df.residual(sv6)))
tab5=c(5,df.residual(sv5),
sum(residuals(sv5,type="pearson")^2),
deviance(sv5),1-pchisq(deviance(sv5),df.residual(sv5)))
tab4=c(4,df.residual(sv4),
sum(residuals(sv4,type="pearson")^2),
deviance(sv4),1-pchisq(deviance(sv4),df.residual(sv4)))
tab1=c(1,df.residual(sv1),
sum(residuals(sv1,type="pearson")^2),
deviance(sv1),1-pchisq(deviance(sv1),df.residual(sv1)))
tab2=c(2,df.residual(sv2),
sum(residuals(sv2,type="pearson")^2),
deviance(sv2),1-pchisq(deviance(sv2),df.residual(sv2)))
tab3=c(3,df.residual(sv3),
sum(residuals(sv3,type="pearson")^2),
deviance(sv3),1-pchisq(deviance(sv3),df.residual(sv3)))
tab0=c(0,df.residual(sv0),
sum(residuals(sv0,type="pearson")^2),
deviance(sv0),1-pchisq(deviance(sv0),df.residual(sv0)))

t(matrix(c(tab7,tab6,tab5,tab4,tab1,tab2,tab3,tab0),5,8))
anova(sv0,sv6)
qchisq(0.95,2)
anova(sv3,sv6)
qchisq(0.95,1)
anova(sv3,sv7)