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

art=glm(ct ~ T:A:S + T:R + A:R,family=poisson)
summary(art)
art$fit


art=glm(ct ~ T:A:S + T:A:R,family=poisson)
summary(art)