rm(list = ls()) tense <-
read.table(
url("http://stat.unm.edu/~fletcher/LLM/DATA/Example3-7-1.DAT"),
#"C:\\E-drive\\Books\\ANREG2\\newdata\\tab20-10a.dat",
#"C:\\E-drive\\Books\\LOGLIN3\\DATA\\Example3-7-1.dat", 		
sep="",col.names=c("y","Tn","Wt","Ms","Dr")) attach(tense)
tense #summary(tense)

W=factor(Wt) M=factor(Ms) D=factor(Dr) T=factor(Tn) yy=log(y)

m00 = lm(yy ~ T*W*M*D)
m0a=anova(m00)
SS=m0a[,2]

tab=c(4*sqrt(SS[-16]),4*sqrt(SS[-16])/sqrt(sum(1/y)))

matrix(tab,15,2,dimnames = list(NULL,c("Est", "z")))