occupationalStatus         package:datasets         R Documentation

_O_c_c_u_p_a_t_i_o_n_a_l _S_t_a_t_u_s _o_f _F_a_t_h_e_r_s _a_n_d _t_h_e_i_r _S_o_n_s

_D_e_s_c_r_i_p_t_i_o_n:

     Cross-classification of a sample of British males according to
     each subject's occupational status and his father's occupational
     status.

_U_s_a_g_e:

     occupationalStatus

_F_o_r_m_a_t:

     A 'table' of counts, with classifying factors 'origin' (father's
     occupational status; levels '1:8') and 'destination' (son's
     occupational status; levels '1:8').

_S_o_u_r_c_e:

     Goodman, L. A. (1979) Simple Models for the Analysis of
     Association in Cross-Classifications having Ordered Categories.
     _J. Am. Stat. Assoc._, *74* (367), 537-552.

     The data set has been in package 'gnm' and been provided by the
     package authors.

_E_x_a_m_p_l_e_s:

     require(stats); require(graphics)

     plot(occupationalStatus)

     ##  Fit a uniform association model separating diagonal effects
     Diag <- as.factor(diag(1:8))
     Rscore <- scale(as.numeric(row(occupationalStatus)), scale = FALSE)
     Cscore <- scale(as.numeric(col(occupationalStatus)), scale = FALSE)
     modUnif <- glm(Freq ~ origin + destination + Diag + Rscore:Cscore,
                    family = poisson, data = occupationalStatus)

     summary(modUnif)
     plot(modUnif) # 4 plots, with warning about  h_ii ~= 1

