data(Hauser79)
# display table
structable(~Father + Son, data=Hauser79)
hauser.indep <- gnm(Freq ~ Father + Son,
data=Hauser79,
family=poisson)
hauser.CR <- update(hauser.indep,
~ . + Crossings(Father,Son))
LRstats(hauser.CR)
hauser.CRdiag <- update(hauser.indep,
~ . + Crossings(Father,Son) + Diag(Father,Son))
LRstats(hauser.CRdiag)
Run the code above in your browser using DataLab