rm(list = ls())
abt <- read.table(
url("http://stat.unm.edu/~fletcher/LLM/DATA/TAB3-1.DAT"),
#"C:\\E-drive\\Books\\ANREG2\\newdata\\TAB21-4.DAT",
#"C:\\E-drive\\Books\\LOGLIN3\\DATA\\Tab3-1.dat",
		  sep="",col.names=c("R","S","A","O","y"))
attach(abt)
abt
#summary(abt)

r=factor(R)
o=factor(O)
s=factor(S)
a=factor(A)
m11 <- glm(y ~ r:s:a + r:s:o + o:a   ,family = poisson)


m11s=summary(m11)
m11s
anova(m11)

rpearson=(y-m11$fit)/(m11$fit)^(.5)
rstand=rpearson/(1-hatvalues(m11))^(.5)
infv = c(y,m11$fit,hatvalues(m11),rpearson,rstand,
cooks.distance(m11))
inf=matrix(infv,I(m11s$df[1]+m11s$df[2]),6,dimnames =
list(NULL,c("y","yhat","lev","Pearson","Stand.","C")))
inf