t200=c(0,2,3,0,1,2,2,0)
t210=c(1,1,2,1,1,2,2,0)
t220=c(2,0,1,2,1,2,2,0)
t201=c(0,2,3,0,2,1,1,1)
t211=c(1,1,2,1,2,1,1,1)
t221=c(2,0,1,2,2,1,1,1)
t202=c(0,2,3,0,3,0,0,2)
t212=c(1,1,2,1,3,0,0,2)
t222=c(2,0,1,2,3,0,0,2)

i=c(1,2,1,2,1,2,1,2)
j=c(1,1,2,2,1,1,2,2)
k=c(1,1,1,1,2,2,2,2)
i=factor(i)
j=factor(j)
k=factor(k)
fit=glm(t200 ~ i:k + j, poisson)
fit$dev
fit=glm(t210 ~ i:k + j, poisson)
fit$dev
fit=glm(t220 ~ i:k + j, poisson)
fit$dev
fit=glm(t201 ~ i:k + j, poisson)
fit$dev
fit=glm(t211 ~ i:k + j, poisson)
fit$dev
fit=glm(t221 ~ i:k + j, poisson)
fit$dev
fit=glm(t202 ~ i:k + j, poisson)
fit$dev
fit=glm(t212 ~ i:k + j, poisson)
fit$dev
fit=glm(t222 ~ i:k + j, poisson)
fit$dev