# ehd7.cf   multiepisode exponential model with up to four jobs

nvar(

    dfile = rrdat.1,     # data file

    isel = le(c2,4),     # select up to four jobs

    ID     [3.0] = c1,    # identification number
    NOJ    [2.0] = c2,    # spell number
    TStart [3.0] = c3,    # starting time
    TFin   [3.0] = c4,    # ending time
    SEX    [2.0] = c5,    # sex (1 men, 2 women)
    TI     [3.0] = c6,    # interview date
    TB     [3.0] = c7,    # birth date
    TE     [3.0] = c8,    # entry into labor market
    TMAR   [3.0] = c9,    # marriage date (0 if no marriage)
    PRES   [3.0] = c10,   # prestige of current job
    PRESN  [3.0] = c11,   # prestige of next job
    EDU    [2.0] = c12,   # highest educational attainment

    # define additional variables

    ORG [1.0] = 0,                  
    DES [1.0] = lt(TFin,TI),
    TFC [3.0] = if eq(NOJ,1) then TStart else pre(TFC),
    TSP [3.0] = TStart - TFC,
    TFP [3.0] = TFin - TFC + 1,

    COHO2  =   ge(TB,468) & le(TB,504),  # cohort 2
    COHO3  =   ge(TB,588) & le(TB,624),  # cohort 3
    LFX    =   TStart - TE,              # labor force experience
    PNOJ   =   NOJ - 1,                  # previous number of jobs
);

edef(           # define multi-episode data

    id = ID,    # ID
    sn = NOJ,   # spell number
    ts = TSP,   # starting time
    tf = TFP,   # ending time
    org = ORG,  # origin state
    des = DES,  # destination state
);
rate( 
    xa(0,1,1) = EDU,COHO2,COHO3,    PRES,      
    xa(0,1,2) = EDU,COHO2,COHO3,LFX,PRES,      
    xa(0,1,3) = EDU,COHO2,COHO3,LFX,PRES,      
    xa(0,1,4) = EDU,COHO2,COHO3,LFX,PRES,      

)   = 2;