rm(list = ls())
ct=c(4,42,57,2,7,20,9,4,19,71,7,8,10,31)
s=c(1,1,1,1,1,1,2,2,2,2,2,2,2,2)
a=c(1,1,1,2,2,2,1,1,1,1,2,2,2,2)
h=c(1,3,4,1,3,4,1,2,3,4,1,2,3,4)
S=factor(s)
A=factor(a)
H=factor(h)
m7=glm(ct ~ S:A + S:H + A:H,family=poisson)
summary(m7)
m6=glm(ct ~ S:H + A:H,family=poisson)
summary(m6)
m5=glm(ct ~ S:A + A:H,family=poisson)
summary(m5)
m4=glm(ct ~ S:A + S:H ,family=poisson)
summary(m4)
m3=glm(ct ~ S:A + H,family=poisson)
summary(m3)
m2=glm(ct ~ S:H + A,family=poisson)
summary(m2)
m1=glm(ct ~ S + A:H,family=poisson)
summary(m1)
m0=glm(ct ~ S + A + H,family=poisson)
summary(m0)


tab7=c(7,df.residual(m7),
deviance(m7),1-pchisq(deviance(m7),df.residual(m7)))
tab6=c(6,df.residual(m6),
deviance(m6),1-pchisq(deviance(m6),df.residual(m6)))
tab5=c(5,df.residual(m5),
deviance(m5),1-pchisq(deviance(m5),df.residual(m5)))
tab4=c(4,df.residual(m4),
deviance(m4),1-pchisq(deviance(m4),df.residual(m4)))
tab1=c(1,df.residual(m1),
deviance(m1),1-pchisq(deviance(m1),df.residual(m1)))
tab2=c(2,df.residual(m2),
deviance(m2),1-pchisq(deviance(m2),df.residual(m2)))
tab3=c(3,df.residual(m3),
deviance(m3),1-pchisq(deviance(m3),df.residual(m3)))
tab0=c(0,df.residual(m0),
deviance(m0),1-pchisq(deviance(m0),df.residual(m0)))

t(matrix(c(tab7,tab6,tab5,tab4,tab3,tab2,tab1,tab0),4,8))