? UNION-NONUNION WAGE NAMELIST ; X = ONE, ED, URB, MINOR, AGE, EXP, REGS, UNEMPR, OCCW, OCCB $ ? CHOOSING EMPLOYED ONLY REJECT ; EMP # 1 $ ? TWO-STAGE AND MLE CALC ; K = COL(X) $ PROBIT ; LHS = UNION ; RHS = X ; HOLD $ MATRIX ; PB = B $ SELECT ; LHS = LRATE ; RHS = X $ MATRIX ; BETA1 = BSR1 $ ? BSR1 = (b1,s1,rho_{u1}) MATRIX ; B1 = PART(BSR1,1,K) $ SELECT ; LHS = LRATE ; RHS = X ; LIMITS = 1 $ MATRIX ; BETA0 = BSR0 $ ? BASR0 = (b0,s0,rho{u2}) MATRIX ; B0 = PART(BSR0,1,K) $ SELECT ; LHS = LRATE ; RH1 = X ; RH2 = X ; MLE ; ALL ; START = BETA1,BETA0 ; TLF = 0.00001 ; TLB = 0.00001 ; TLG = 0.00001 ; MAXIT = 1000 $ CALC ; ULOGL = LOGL $ ? RESTRICTED MODEL: B_1 = B_2 SELECT ; LHS = LRATE ; RH1 = X ; RH2 = X ; MLE ; ALL ; ; RST = K_BP, K_B, K_B, B3, B4, B5, B6 ; START = BETA1, BETA1 ; TLF = 0.0000001 ; TLB = 0.0000001 ; TLG = 0.0000001 ; MAXIT = 1000 $ CALC ; RLOGL = LOGL $ ? COMPUTING THE LR STATISTICS TITLE ; LR test for b1 = b2 $ CALC ; LIST ; LRT = 2*(ULOGL - RLOGL) ; PVAL = 1 - CHI(LRT,K) $