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))
summarise(hauser.CR)
hauser.CRdiag <- update(hauser.indep, ~ . + Crossings(Father,Son) + Diag(Father,Son))
summarise(hauser.CRdiag)
Run the code above in your browser using DataLab