t00=c(1,2,2,0,0,2,3,0)
t01=c(1,2,2,0,1,1,2,1)
t02=c(1,2,2,0,2,0,1,2)
t10=c(2,1,1,1,0,2,3,0)
t11=c(2,1,1,1,1,1,2,1)
t12=c(2,1,1,1,2,0,1,2)
t20=c(3,0,0,2,0,2,3,0)
t21=c(3,0,0,2,1,1,2,1)
t22=c(3,0,0,2,2,0,1,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(t00 ~ i:k + j:k, poisson)
fit$dev
fit=glm(t01 ~ i:k + j:k, poisson)
fit$dev
fit=glm(t02 ~ i:k + j:k, poisson)
fit$dev
fit=glm(t10 ~ i:k + j:k, poisson)
fit$dev
fit=glm(t11 ~ i:k + j:k, poisson)
fit$dev
fit=glm(t12 ~ i:k + j:k, poisson)
fit$dev
fit=glm(t20 ~ i:k + j:k, poisson)
fit$dev
fit=glm(t21 ~ i:k + j:k, poisson)
fit$dev
fit=glm(t22 ~ i:k + j:k, poisson)
fit$dev