rm(list = ls())
tense <- read.table(
url("http://stat.unm.edu/~fletcher/LLM/DATA/Example3-7-1-logistic.DAT"),
#"C:\\E-drive\\Books\\ANREG2\\newdata\\tab20-10a.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)
ts <- glm(T ~ W + M*D,family = binomial)
tsp=summary(ts)
tsp
anova(ts)

fitted(ts)/(1-fitted(ts))