t300=c(1,2,2,0,0,2,3,0)
t301=c(1,2,2,0,1,1,2,1)
t302=c(1,2,2,0,2,0,1,2)
t310=c(2,1,1,1,0,2,3,0)
t311=c(2,1,1,1,1,1,2,1)
t312=c(2,1,1,1,2,0,1,2)
t320=c(3,0,0,2,0,2,3,0)
t321=c(3,0,0,2,1,1,2,1)
t322=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(t300 ~ i:k + j, poisson)
fit$dev
fit=glm(t301 ~ i:k + j, poisson)
fit$dev
fit=glm(t302 ~ i:k + j, poisson)
fit$dev
fit=glm(t310 ~ i:k + j, poisson)
fit$dev
fit=glm(t311 ~ i:k + j, poisson)
fit$dev
fit=glm(t312 ~ i:k + j, poisson)
fit$dev
fit=glm(t320 ~ i:k + j, poisson)
fit$dev
fit=glm(t321 ~ i:k + j, poisson)
fit$dev
fit=glm(t322 ~ i:k + j, poisson)
fit$dev