rm(list = ls())
ct=c(43,16,3,6,11,10,9,18,16)
L=c(1,2,3,1,2,3,1,2,3)
V=c(1,1,1,2,2,2,3,3,3)
l=factor(L)
v=factor(V)
ind=glm(ct ~ v + l,family=poisson)
summary(ind)
t=log(ind$fit)
t2=t*t
m11=glm(ct ~ v + l + v:t,family=poisson)
summary(m11)
m12=glm(ct ~ v + l + l:t,family=poisson)
summary(m12)
m13=glm(ct ~ v + l + t2,family=poisson)
summary(m13)

summary(m11)
t