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)
# what does Crossings do?
cr <- with(Hauser79, Crossings(Father, Son))
head(cr)
# Show the codings for varying Crossings levels
matrix(cr[,1], nrow=5)
matrix(cr[,2], nrow=5)
matrix(cr[,3], nrow=5)
matrix(cr[,4], nrow=5)
Run the code above in your browser using DataLab