RChomog <- gnm(Freq ~ origin + destination + Diag(origin, destination) +
MultHomog(origin, destination), family = poisson,
data = occupationalStatus, ofInterest = "MultHomog",
constrain = "MultHomog.*1")
coefficients(RChomog)
parameters(RChomog)
Run the code above in your browser using DataLab