rm(list = ls())
High=c(245,330,388,100,77,51,28,89,102,67,87,62,
   125,234,233,109,197,90)
Low=c(115,152,153,40,37,19,11,37,35,18,12,13,68,
   91,173,47,82,32)
R=c(1,1,1,1,1,1,2,2,2,2,2,2,3,3,3,3,3,3)
E=c(1,2,3,4,5,6,1,2,3,4,5,6,1,2,3,4,5,6)
r=factor(R)
e=factor(E)
T=cbind(High,Low)
m4=glm(T ~ r + e,family=binomial)
summary(m4)
t=log(m4$fit/(1-m4$fit))
t2=t*t
m5=glm(T ~ r + e + r:t,family=binomial)
summary(m5)
m6=glm(T ~ r + e + e:t,family=binomial)
summary(m6)
m7=glm(T ~ r + e + t2,family=binomial)
summary(m7)