cnt=c(350,150,60,112,26,23,19,80)
ii=c(1,1,1,1,2,2,2,2)
jj=c(1,2,1,2,1,2,1,2)
kk=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+JJ:KK,family = poisson)

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