cnt=c(16,7,15,34,5,3,1,1,3,8,1,3)
ii=c(1,1,1,1,1,1,2,2,2,2,2,2)
jj=c(1,2,1,2,1,2,1,2,1,2,1,2)
kk=c(1,1,2,2,3,3,1,1,2,2,3,3)
II=factor(ii)
JJ=factor(jj)
KK=factor(kk)
sv <- glm(cnt ~ II+JJ:KK,family = poisson)

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