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)