LaborParticipation            package:cmm            R Documentation

_W_o_m_e_n'_s _l_a_b_o_r _p_a_r_t_i_c_i_p_a_t_i_o_n: _1_9_6_7-_1_9_7_1

_D_e_s_c_r_i_p_t_i_o_n:

     The labor participation (yes/no) of 1583 women was measured in
     five consecutive year, 1967-1971,  leading to a 2 x 2 x 2 x 2 x 2
     table.

     The data are tabulated in Bergsma, Croon, and Hagenaars (2009, p.
     128).

     Section 4.3 in Bergsma, Croon and Hagenaars, 2009

_U_s_a_g_e:

     data(LaborParticipation)

_F_o_r_m_a_t:

     A data frame with 1583 observations on the following variables.

     '_Y_e_a_r_1_9_6_7' Participated in 1967 (factor): 1 = No; 2 = Yes.

     '_Y_e_a_r_1_9_6_8' Participated in 1968 (factor): 1 = No; 2 = Yes.

     '_Y_e_a_r_1_9_6_9' Participated in 1969 (factor): 1 = No; 2 = Yes.

     '_Y_e_a_r_1_9_7_0' Participated in 1970 (factor): 1 = No; 2 = Yes.

     '_Y_e_a_r_1_9_7_1' Participated in 1971 (factor): 1 = No; 2 = Yes.


_S_o_u_r_c_e:

     Heckman & Willis (1977).

_R_e_f_e_r_e_n_c_e_s:

     Bergsma, W. P., Croon, M. A., & Hagenaars, J. A. P. (2009).
     _Marginal models for dependent, clustered, and longitudinal
     categorical data._ New York: Springer.

     Heckman, J. J. & Willis, R. J. (1977). A beta-logistic model for
     the analysis of sequential labor force participation by married
     women. _Journal of Political Economy, 85_, 27-58.

_E_x_a_m_p_l_e_s:

     data(LaborParticipation)
     n <- c(t(ftable(LaborParticipation)))

     ##############################################################################################################
     #Sample kappa values

     #matrix for obtaining transition matrices for consecutive years
     at = MarginalMatrix(c("67","68","69","70","71"),list(c("67","68"),c("68","69"),c("69","70"),c("70","71")),c(2,2,2,2,2));

     coeff = SpecifyCoefficient("CohenKappa",2,4);
     SampleStatistics(n,list(coeff,at),ShowParameters=FALSE)

     ##############################################################################################################
     #Fitting models for kappa
     #matrix for obtaining transition matrices for consecutive years
     at = MarginalMatrix(c("67","68","69","70","71"),list(c("67","68"),c("68","69"),c("69","70"),c("70","71")),c(2,2,2,2,2));

     coeff = SpecifyCoefficient("CohenKappa",2,4);

     bt1 = ConstraintMatrix(c(1),c(),c(4)); #equal kappas
     bt2 = rbind(c(1,-2,1,0),c(0,1,-2,1));  #linear trend for kappas
     model = list(bt1,coeff,at)

     m=MarginalModelFit(n,model,ShowParameters=FALSE,ShowProgress=10)

