set.seed(1)
### Fit an association model with homogeneous row-column effects
rc1 <- gnm(Freq ~ r + c + Diag(r,c) + MultHomog(r, c),
family = poisson, data = friend)
rc1
if (FALSE) {
### Extend to two-component interaction
rc2 <- update(rc1, . ~ . + MultHomog(r, c, inst = 2),
etastart = rc1$predictors)
rc2
}
Run the code above in your browser using DataLab