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.
set.seed(1) data(friend) ## Fit an association model with homogeneous row-column effects rc1 <- gnm(Freq ~ r + c + Diag(r,c) + Nonlin(MultHomog(r, c)), 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))) rc2