rm(list = ls())
mel <- read.table(
url("http://stat.unm.edu/~fletcher/LLM/DATA/TAB8-4.DAT"),
#"C:\\E-drive\\Books\\LOGLIN3\\DATA\\TAB8-4.DAT",
		  sep="",col.names=c("hh","ii","jj","kk","ct"))
attach(mel)
mel
#summary(mel)
G=factor(hh)
R=factor(ii)
Im=factor(jj)
S=factor(kk)

out=glm(ct~G:R:Im+G:R:S+G:Im:S+R:Im:S,family=poisson)
summary(out)
out$fit
ct