tense <- read.table(
url("http://stat.unm.edu/~fletcher/LLM/DATA/Example3-7-1-logistic.DAT"),
#"C:\\E-drive\\Books\\ANREG2\\newdata\\tab20-10.dat",
#"C:\\E-drive\\Books\\LOGLIN3\\DATA\\Example3-7-1-logistic.dat",
		  sep="",col.names=c("High","Low","Wt","Ms","Dr"))
attach(tense)
tense
#summary(tense)

W=factor(Wt)
M=factor(Ms)
D=factor(Dr)
T=cbind(High,Low)
sv7 <- glm(T ~ W:M+W:D+M:D,family = binomial)
sv6 <- glm(T ~ W:M+W:D,family = binomial)
sv5 <- glm(T ~ W:M+M:D,family = binomial)
sv4 <- glm(T ~ W:D+M:D,family = binomial)
sv1 <- glm(T ~ W+M:D,family = binomial)
sv2 <- glm(T ~ M+W:D,family = binomial)
sv3 <- glm(T ~ W:M+D,family = binomial)
sv0 <- glm(T ~ W+M+D,family = binomial)
svd <- glm(T ~ W+M,family = binomial)
svm <- glm(T ~ W+D,family = binomial)
svw <- glm(T ~ M+D,family = binomial)


tab7=c(7,df.residual(sv7),deviance(sv7),
1-pchisq(deviance(sv7),df.residual(sv7)),
-2*df.residual(sv7)+deviance(sv7))
tab6=c(6,df.residual(sv6),deviance(sv6),
1-pchisq(deviance(sv6),df.residual(sv6)),
-2*df.residual(sv6)+deviance(sv6))
tab5=c(5,df.residual(sv5),deviance(sv5),
1-pchisq(deviance(sv5),df.residual(sv5)),
-2*df.residual(sv5)+deviance(sv5))
tab4=c(4,df.residual(sv4),deviance(sv4),
1-pchisq(deviance(sv4),df.residual(sv4)),
-2*df.residual(sv4)+deviance(sv4))
tab1=c(1,df.residual(sv1),deviance(sv1),
1-pchisq(deviance(sv1),df.residual(sv1)),
-2*df.residual(sv1)+deviance(sv1))
tab2=c(2,df.residual(sv2),deviance(sv2),
1-pchisq(deviance(sv2),df.residual(sv2)),
-2*df.residual(sv2)+deviance(sv2))
tab3=c(3,df.residual(sv3),deviance(sv3),
1-pchisq(deviance(sv3),df.residual(sv3)),
-2*df.residual(sv3)+deviance(sv3))
tab0=c(0,df.residual(sv0),deviance(sv0),
1-pchisq(deviance(sv0),df.residual(sv0)),
-2*df.residual(sv0)+deviance(sv0))
tabd=c(2,df.residual(svd),deviance(svd),
1-pchisq(deviance(svd),df.residual(svd)),
-2*df.residual(svd)+deviance(svd))
tabm=c(3,df.residual(svm),deviance(svm),
1-pchisq(deviance(svm),df.residual(svm)),
-2*df.residual(svm)+deviance(svm))
tabw=c(0,df.residual(svw),deviance(svw),
1-pchisq(deviance(svw),df.residual(svw)),
-2*df.residual(svw)+deviance(svw))

t(matrix(c(tab7,tab6,tab5,tab4,tab3,
tab2,tab1,tab0,tabd,tabm,tabw),5,11))

anova(sv0,sv6)
qchisq(0.95,2)
anova(sv3,sv6)
qchisq(0.95,1)
anova(sv3,sv7)