| friend {gnm} | R Documentation |
Cross-classification of the occupation of respondent and that of their closest friend. Data taken from wave 10 (year 2000) of the British Household Panel Survey.
data(friend)
A table of counts, with classifying factors r (respondent's
occupational category; levels 1:31) and c (friend's
occupational category; levels 1:31).
Chan, T.W. and Goldthorpe, J.H. (2004) Is there a status order in contemporary British society: Evidence from the occupational structure of friendship, European Sociological Review, 20, 383–401.
## Not run:
### (Not needed for package testing - same as example(MultHomog))
### Fit an association model with homogeneous row-column effects
rc1 <- gnm(Freq ~ r + c + Diag(r,c) + Nonlin(MultHomog(r, c)),
start = c(coef(lin), init[1:31, 1]), family = poisson,
data = friend)
rc1
### Extend to two-component interaction (using workaround)
original <- as.data.frame(friend)
r2 <- original$r
c2 <- original$c
rc2 <- update(rc1, . ~ . + Nonlin(MultHomog(r2, c2)),
start = c(coef(rc1), init[1:31, 2]))
rc2
## End(Not run)