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)
sv <- glm(cnt ~ II:JJ+II:KK,family = poisson)

fitted(sv)
sum(residuals(sv,type="pearson")^2)
deviance(sv)
df.residual(sv)