rm(list = ls())
sb <- read.table(
url("http://stat.unm.edu/~fletcher/LLM/DATA/Example3-2-4.dat"),
#"C:\\E-drive\\Books\\ANREG2\\newdata\\EX21-3-2.dat",
#"C:\\E-drive\\Books\\LOGLIN3\\DATA\\Example3-2-4.dat",
		  sep="",col.names=c("y","i","j","k"))

attach(sb)
#summary(sb)
I=factor(i)
J=factor(j)
K=factor(k)
m7 <- glm(y ~ I + J+ K+ J:K,family = poisson)
m7s=summary(m7)
vcov(m7)

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