# tables from Table 16.4
t402=c(2,2,1,0,1,0,2,2)
t401=c(2,2,1,0,0,1,3,1)
t411=c(3,1,0,1,0,1,3,1)
t412=c(3,1,0,1,1,0,2,2)
t502=c(3,2,0,0,0,0,3,2)
# reverse sides of Table 16.4
# G^2 same as first 5
t120=c(1,0,2,2,2,2,1,0)
t110=c(0,1,3,1,2,2,1,0)
t111=c(0,1,3,1,3,1,0,1)
t121=c(1,0,2,2,3,1,0,1)
t020=c(0,0,3,2,3,2,0,0)
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(t402 ~ i:k + j, poisson)
fit$dev
fit=glm(t401 ~ i:k + j, poisson)
fit$dev
fit=glm(t411 ~ i:k + j, poisson)
fit$dev
fit=glm(t412 ~ i:k + j, poisson)
fit$dev
fit=glm(t502 ~ i:k + j, poisson)
fit$dev
fit=glm(t120 ~ i:k + j, poisson)
fit$dev
fit=glm(t110 ~ i:k + j, poisson)
fit$dev
fit=glm(t111 ~ i:k + j, poisson)
fit$dev
fit=glm(t121 ~ i:k + j, poisson)
fit$dev
fit=glm(t020 ~ i:k + j, poisson)
fit$dev