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+JJ:KK,family = poisson)
anova(sv)