Supplementary results for paper about expected wage differentials
Germany
Men
Reduced form model
              ÉÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍ»
              º FIML Estimates of Bivariate Probit Model    º
              º Maximum Likelihood Estimates                º
              º Dependent variable                OCCFT     º
              º Number of observations             1823     º
              º Iterations completed                 29     º
              º Log likelihood function       -1364.958     º
              º Selection model based on FT                 º
              º Means for vars. 1- 12 are after selection.  º
              ÈÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍͼ

  Variable  Coefficient   Standard Error  z=b/s.e. P[³Z³òz]   Mean of X
  ÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄ
          Index    equation for OCCUPN
  Constant   -2.8503         0.53910       -5.287   0.00000
  SCH1       0.18366         0.17836E-01   10.297   0.00000   10.67
  EX         0.29389E-01     0.18556E-01    1.584   0.11324   20.83
  EXSQ      -0.71615E-03     0.37074E-03   -1.932   0.05340   546.6
  PT        -0.96167E-02     0.84095E-01   -0.114   0.90896  0.3172
  Y89        0.26198         0.12941        2.024   0.04293  0.1463
  Y90        0.16169         0.10181        1.588   0.11226  0.3364
  Y91        0.22937         0.11498        1.995   0.04606  0.2094
  REG1      -0.13670         0.16794       -0.814   0.41564  0.1476
  REG2      -0.30213         0.26068       -1.159   0.24645  0.3524E-01
  REG3      -0.17936         0.14996       -1.196   0.23166  0.4428
  REG4      -0.80981E-01     0.15545       -0.521   0.60239  0.3251
          Index    equation for FT
  Constant   0.22860         0.17372        1.316   0.18820
  EX         0.35474E-01     0.27210E-02   13.037   0.00000   19.33
  HHSIZE     0.15509         0.25271E-01    6.137   0.00000   2.592
  LOCUNEMP  -0.28875E-01     0.15356E-01   -1.880   0.06006   7.731
  Y89       -0.21615         0.11020       -1.961   0.04983  0.1476
  Y90       -0.48807E-01     0.96403E-01   -0.506   0.61266  0.3390
  Y91       -0.15502E-01     0.11761       -0.132   0.89514  0.2074
          Disturbance correlation and std. dev.
  RHO(1,2)  -0.57391         0.29336       -1.956   0.05042

ÉÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍ»
Knowledge workers
ÉÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍ»
º Bivariate Probit Sample Selection Model                               º
º Two stage   least squares regression    Weighting variable = ONE      º
º Dependent variable is SALARY    Mean =    2.43494, S.D. =      0.4610 º
º Model size: Observations =     206, Parameters =  10, Deg.Fr. =   196 º
º Residuals:  Sum of squares=    31.9196     Std.Dev. =         0.40355 º
º Fit:        R-squared = 0.22997, Adjusted R-squared =         0.19462 º
º             Note:  Not using OLS.  R-squared is not bounded in [0,1]  º
º Model test: F[  9,    196] =    6.50,    Prob value =         0.00000 º
º Diagnostic: Log-L =   -100.2417, Restricted(á=0) Log-L =    -132.2841 º
º             Amemiya Pr. Crt.=    0.171, Akaike Info. Crt.=      1.070 º
º Full sample contains        1921 observations.                        º
º Selected sample contains     206 observations.                        º
º Estimated disturbance standard deviation =        0.446780            º
º Estimated correlation with selection equation A =-0.343215            º
º Estimated correlation with selection equation B =-0.434682            º
ÈÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍͼ

  Variable  Coefficient   Standard Error  z=b/s.e. P[³Z³òz]   Mean of X
  ÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄ
  Constant    1.5027          1.2538        1.199   0.23071
  SCH1       0.21008E-01     0.49560E-01    0.424   0.67165
  EX         0.59993E-01     0.18048E-01    3.324   0.00089
  EXSQ      -0.14343E-02     0.47412E-03   -3.025   0.00249
  PT         0.10301         0.69271E-01    1.487   0.13699
  Y89        0.23447         0.11366        2.063   0.03913
  Y90        0.30942         0.90726E-01    3.411   0.00065
  Y91        0.23344         0.10046        2.324   0.02015
  Lambda-A  -0.15334         0.47123       -0.325   0.74487
  Lambda-B  -0.19421         0.21407       -0.907   0.36429

ÉÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍ»
Other workers
ÉÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍ»
º Bivariate Probit Sample Selection Model                               º
º Two stage   least squares regression    Weighting variable = ONE      º
º Dependent variable is SALARY    Mean =    2.05674, S.D. =      0.3731 º
º Model size: Observations =    1376, Parameters =  10, Deg.Fr. =  1366 º
º Residuals:  Sum of squares=    142.067     Std.Dev. =         0.32249 º
º Fit:        R-squared = 0.25216, Adjusted R-squared =         0.24724 º
º             Note:  Not using OLS.  R-squared is not bounded in [0,1]  º
º Model test: F[  9,   1366] =   51.18,    Prob value =         0.00000 º
º Diagnostic: Log-L =   -390.2619, Restricted(á=0) Log-L =    -595.1921 º
º             Amemiya Pr. Crt.=    0.105, Akaike Info. Crt.=      0.582 º
º Full sample contains        1921 observations.                        º
º Selected sample contains    1376 observations.                        º
º Estimated disturbance standard deviation =        0.320797            º
º Estimated correlation with selection equation A = -1.72379            º
º Estimated correlation with selection equation B = 0.109971            º
ÈÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍͼ

  Variable  Coefficient   Standard Error  z=b/s.e. P[³Z³òz]   Mean of X
  ÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄ
  Constant    1.1851         0.20347        5.825   0.00000
  SCH1       0.48162E-02     0.23949E-01    0.201   0.84062
  EX         0.48888E-01     0.61202E-02    7.988   0.00000
  EXSQ      -0.81034E-03     0.12568E-03   -6.447   0.00000
  PT         0.71229E-01     0.28531E-01    2.497   0.01254
  Y89       -0.85005E-02     0.49321E-01   -0.172   0.86316
  Y90        0.67259E-01     0.36986E-01    1.819   0.06899
  Y91        0.13399         0.43505E-01    3.080   0.00207
  Lambda-A  -0.55299         0.25334       -2.183   0.02905
  Lambda-B   0.35278E-01     0.13595        0.259   0.79525

ÉÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍ»
Structural model
              ÉÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍ»
              º FIML Estimates of Bivariate Probit Model    º
              º Maximum Likelihood Estimates                º
              º Dependent variable                OCCFT     º
              º Number of observations             1823     º
              º Iterations completed                 20     º
              º Log likelihood function       -1129.088     º
              º Selection model based on FT                 º
              º Means for vars. 1-  7 are after selection.  º
              ÈÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍͼ

  Variable  Coefficient   Standard Error  z=b/s.e. P[³Z³òz]   Mean of X
  ÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄ
          Index    equation for OCCUPN
  Constant   -3.7378         0.42208       -8.856   0.00000
  EXRELW      4.8879         0.44869       10.894   0.00000  0.5286
  EX         0.92383E-02     0.52877E-02    1.747   0.08062   20.83
  REG1      -0.18572         0.17901       -1.038   0.29949  0.1476
  REG2      -0.81727E-01     0.21190       -0.386   0.69973  0.3524E-01
  REG3      -0.22797         0.16826       -1.355   0.17547  0.4428
  REG4      -0.11417         0.17366       -0.657   0.51091  0.3251
          Index    equation for FT
  Constant   0.69507         0.16451        4.225   0.00002
  EX         0.37972E-01     0.25898E-02   14.662   0.00000   19.33
  HHSIZE     0.86886E-01     0.24222E-01    3.587   0.00033   2.592
  LOCUNEMP  -0.28446E-01     0.14100E-01   -2.017   0.04364   7.731
  Y89       -0.79761         0.10270       -7.767   0.00000  0.1476
  Y90       -0.69781         0.88956E-01   -7.844   0.00000  0.3390
  Y91       -0.16943         0.88689E-01   -1.910   0.05608  0.2074
          Disturbance correlation and std. dev.
  RHO(1,2)  -0.99519         0.78394E-02 -126.947   0.00000


    MGLFX=R5*R6*BBB1(2)
    MGLFX   =   1.12819

Single women
Reduced form model      
	        ÉÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍ»
              º FIML Estimates of Bivariate Probit Model    º
              º Maximum Likelihood Estimates                º
              º Dependent variable                OCCFT     º
              º Number of observations              684     º
              º Iterations completed                 22     º
              º Log likelihood function       -648.0525     º
              º Selection model based on FT                 º
              º Means for vars. 1- 12 are after selection.  º
              ÈÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍͼ

  Variable  Coefficient   Standard Error  z=b/s.e. P[³Z³òz]   Mean of X
  ÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄ
          Index    equation for OCCUPN
  Constant   -2.6427          1.1557       -2.287   0.02221
  SCH1       0.18082         0.42941E-01    4.211   0.00003   10.82
  EX         0.68173E-01     0.35837E-01    1.902   0.05713   15.77
  EXSQ      -0.15764E-02     0.75929E-03   -2.076   0.03789   357.6
  PT        -0.32836E-01     0.13651       -0.241   0.80991  0.5077
  Y89       -0.13163         0.21997       -0.598   0.54957  0.1794
  Y90       -0.14674         0.17974       -0.816   0.41426  0.2976
  Y91        0.32145E-01     0.20613        0.156   0.87607  0.1904
  REG1      -0.23533         0.35369       -0.665   0.50582  0.1225
  REG2      -0.93952         0.45956       -2.044   0.04091  0.5689E-01
  REG3      -0.28601         0.30090       -0.951   0.34185  0.4486
  REG4       0.18733E-01     0.29885        0.063   0.95002  0.3195
          Index    equation for FT
  Constant   0.79664         0.24370        3.269   0.00108
  EX         0.11099E-01     0.41266E-02    2.690   0.00715   14.92
  HHSIZE    -0.11529         0.44360E-01   -2.599   0.00935   1.927
  LOCUNEMP  -0.35469E-01     0.20834E-01   -1.702   0.08867   7.857
  Y89       -0.10376         0.14432       -0.719   0.47217  0.1886
  Y90        0.55148E-02     0.12902        0.043   0.96590  0.2909
  Y91        0.41409E-03     0.15492        0.003   0.99787  0.1842
          Disturbance correlation and std. dev.
  RHO(1,2)  -0.43824         0.65997       -0.664   0.50667

ÉÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍ»
Knowledge workers
ÉÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍ»
º Bivariate Probit Sample Selection Model                               º
º Two stage   least squares regression    Weighting variable = ONE      º
º Dependent variable is SALARY    Mean =    2.13707, S.D. =      0.5337 º
º Model size: Observations =      57, Parameters =  10, Deg.Fr. =    47 º
º Residuals:  Sum of squares=    9.77855     Std.Dev. =         0.45613 º
º Fit:        R-squared = 0.25644, Adjusted R-squared =         0.11405 º
º             Note:  Not using OLS.  R-squared is not bounded in [0,1]  º
º Model test: F[  9,     47] =    1.80,    Prob value =         0.09319 º
º Diagnostic: Log-L =    -30.6380, Restricted(á=0) Log-L =     -44.5804 º
º             Amemiya Pr. Crt.=    0.245, Akaike Info. Crt.=      1.426 º
º Full sample contains         716 observations.                        º
º Selected sample contains      57 observations.                        º
º Estimated disturbance standard deviation =        0.557619            º
º Estimated correlation with selection equation A = 0.782872            º
º Estimated correlation with selection equation B =-0.228495            º
ÈÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍͼ

  Variable  Coefficient   Standard Error  z=b/s.e. P[³Z³òz]   Mean of X
  ÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄ
  Constant  -0.35716          1.3314       -0.268   0.78850
  SCH1       0.10355         0.52672E-01    1.966   0.04931
  EX         0.63281E-01     0.31276E-01    2.023   0.04304
  EXSQ      -0.13179E-02     0.96394E-03   -1.367   0.17156
  PT         0.18906         0.16243        1.164   0.24443
  Y89       -0.32898E-01     0.25306       -0.130   0.89657
  Y90       -0.19421         0.20328       -0.955   0.33937
  Y91        0.14719         0.19002        0.775   0.43857
  Lambda-A   0.43654         0.54500        0.801   0.42313
  Lambda-B  -0.12741         0.41589       -0.306   0.75933

ÉÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍ»
Other workers
ÉÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍ»
º Bivariate Probit Sample Selection Model                               º
º Two stage   least squares regression    Weighting variable = ONE      º
º Dependent variable is SALARY    Mean =    1.85290, S.D. =      0.3507 º
º Model size: Observations =     421, Parameters =  10, Deg.Fr. =   411 º
º Residuals:  Sum of squares=    37.1823     Std.Dev. =         0.30078 º
º Fit:        R-squared = 0.26269, Adjusted R-squared =         0.24655 º
º             Note:  Not using OLS.  R-squared is not bounded in [0,1]  º
º Model test: F[  9,    411] =   16.27,    Prob value =         0.00000 º
º Diagnostic: Log-L =    -86.5315, Restricted(á=0) Log-L =    -155.7420 º
º             Amemiya Pr. Crt.=    0.093, Akaike Info. Crt.=      0.459 º
º Full sample contains         716 observations.                        º
º Selected sample contains     421 observations.                        º
º Estimated disturbance standard deviation =        0.291106            º
º Estimated correlation with selection equation A =-0.534527            º
º Estimated correlation with selection equation B =-0.501873            º
ÈÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍͼ

  Variable  Coefficient   Standard Error  z=b/s.e. P[³Z³òz]   Mean of X
  ÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄ
  Constant   0.96478         0.18018        5.354   0.00000
  SCH1       0.37117E-01     0.14604E-01    2.542   0.01103
  EX         0.45137E-01     0.83659E-02    5.395   0.00000
  EXSQ      -0.77202E-03     0.18695E-03   -4.129   0.00004
  PT         0.57838E-01     0.32956E-01    1.755   0.07926
  Y89        0.93486E-01     0.50393E-01    1.855   0.06358
  Y90        0.18459E-01     0.45495E-01    0.406   0.68493
  Y91        0.12911         0.54243E-01    2.380   0.01730
  Lambda-A  -0.15560         0.15393       -1.011   0.31208
  Lambda-B  -0.14610         0.13335       -1.096   0.27327

ÉÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍ»
Structural model
              ÉÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍ»
              º FIML Estimates of Bivariate Probit Model    º
              º Maximum Likelihood Estimates                º
              º Dependent variable                OCCFT     º
              º Number of observations              684     º
              º Iterations completed                  4     º
              º Log likelihood function       -564.5332     º
              º Selection model based on FT                 º
              º Means for vars. 1-  7 are after selection.  º
              ÈÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍͼ

  Variable  Coefficient   Standard Error  z=b/s.e. P[³Z³òz]   Mean of X
  ÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄ
          Index    equation for OCCUPN
  Constant   0.56505          1.6308        0.346   0.72897
  EXRELW      1.0978         0.42872        2.561   0.01045 -0.5279
  EX        -0.87683E-01     0.11910E-01   -7.362   0.00000   15.77
  REG1       0.15020E-01     0.41554        0.036   0.97117  0.1225
  REG2      -0.14530         0.52476       -0.277   0.78186  0.5689E-01
  REG3       0.10239E-02     0.35834        0.003   0.99772  0.4486
  REG4       0.60949E-01     0.41028        0.149   0.88191  0.3195
          Index    equation for FT
  Constant   0.79152         0.24676        3.208   0.00134
  EX         0.13021E-02     0.43461E-02    0.300   0.76448   14.92
  HHSIZE    -0.51668E-01     0.45168E-01   -1.144   0.25267   1.927
  LOCUNEMP  -0.35190E-01     0.21044E-01   -1.672   0.09449   7.857
  Y89       -0.38997E-01     0.14381       -0.271   0.78626  0.1886
  Y90        0.81608E-03     0.12915        0.006   0.99496  0.2909
  Y91       -0.20900E-03     0.15465       -0.001   0.99892  0.1842
          Disturbance correlation and std. dev.
  RHO(1,2)  -0.47185E-02      2.6442       -0.002   0.99858

    MGLFX=R5*R6*BBB1(2)
    MGLFX   =  0.890220E-01

Married women
Reduced form model
              ÉÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍ»
              º FIML Estimates of Bivariate Probit Model    º
              º Maximum Likelihood Estimates                º
              º Dependent variable                OCCFT     º
              º Number of observations             1471     º
              º Iterations completed                 27     º
              º Log likelihood function       -1095.944     º
              º Selection model based on FT                 º
              º Means for vars. 1- 12 are after selection.  º
              ÈÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍͼ

  Variable  Coefficient   Standard Error  z=b/s.e. P[³Z³òz]   Mean of X
  ÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄ
          Index    equation for OCCUPN
  Constant   -3.9681         0.73940       -5.367   0.00000
  SCH1       0.21704         0.40077E-01    5.416   0.00000   10.41
  EX         0.70256E-01     0.46745E-01    1.503   0.13285   21.79
  EXSQ      -0.18642E-02     0.10416E-02   -1.790   0.07348   565.6
  PT         0.10200         0.15384        0.663   0.50730  0.7264
  Y89        0.42076         0.21854        1.925   0.05419  0.1830
  Y90        0.22411         0.19679        1.139   0.25479  0.3179
  Y91        0.42880         0.20227        2.120   0.03401  0.2062
  REG1      -0.71102         0.41839       -1.699   0.08924  0.1291
  REG2       0.26089         0.47634        0.548   0.58390  0.3083E-01
  REG3      -0.32118         0.36576       -0.878   0.37988  0.4547
  REG4      -0.64602         0.38621       -1.673   0.09438  0.3410
          Index    equation for FT
  Constant   0.80740         0.20386        3.961   0.00007
  EX        -0.12400E-01     0.36318E-02   -3.414   0.00064   22.64
  HHSIZE    -0.30241         0.31790E-01   -9.513   0.00000   3.382
  LOCUNEMP  -0.56948E-02     0.14467E-01   -0.394   0.69384   7.629
  Y89        0.91133E-01     0.98914E-01    0.921   0.35687  0.1924
  Y90        0.23961         0.88940E-01    2.694   0.00706  0.2862
  Y91        0.30904         0.10669        2.897   0.00377  0.1740
          Disturbance correlation and std. dev.
  RHO(1,2)   0.52811         0.31605        1.671   0.09473

ÉÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍ»
Knowledge workers
ÉÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍ»
º Bivariate Probit Sample Selection Model                               º
º Two stage   least squares regression    Weighting variable = ONE      º
º Dependent variable is SALARY    Mean =    2.30160, S.D. =      0.4571 º
º Model size: Observations =      53, Parameters =  10, Deg.Fr. =    43 º
º Residuals:  Sum of squares=    5.48868     Std.Dev. =         0.35727 º
º Fit:        R-squared = 0.37748, Adjusted R-squared =         0.24718 º
º             Note:  Not using OLS.  R-squared is not bounded in [0,1]  º
º Model test: F[  9,     43] =    2.90,    Prob value =         0.00907 º
º Diagnostic: Log-L =    -15.1122, Restricted(á=0) Log-L =     -33.2135 º
º             Amemiya Pr. Crt.=    0.152, Akaike Info. Crt.=      0.948 º
º Full sample contains        1528 observations.                        º
º Selected sample contains      53 observations.                        º
º Estimated disturbance standard deviation =        0.388673            º
º Estimated correlation with selection equation A =-0.630701            º
º Estimated correlation with selection equation B =-0.169391            º
ÈÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍͼ

  Variable  Coefficient   Standard Error  z=b/s.e. P[³Z³òz]   Mean of X
  ÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄ
  Constant    1.8821          1.4573        1.291   0.19653
  SCH1      -0.43552E-02     0.59294E-01   -0.073   0.94145
  EX         0.60867E-01     0.30466E-01    1.998   0.04573
  EXSQ      -0.11795E-02     0.90164E-03   -1.308   0.19083
  PT         0.11283         0.13851        0.815   0.41532
  Y89        0.98189E-01     0.16471        0.596   0.55109
  Y90       -0.19682         0.20704       -0.951   0.34179
  Y91        0.21918         0.22202        0.987   0.32354
  Lambda-A  -0.24514         0.50223       -0.488   0.62548
  Lambda-B  -0.65838E-01     0.17623       -0.374   0.70871

ÉÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍ»
Other workers
ÉÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍ»
º Bivariate Probit Sample Selection Model                               º
º Two stage   least squares regression    Weighting variable = ONE      º
º Dependent variable is SALARY    Mean =    1.84405, S.D. =      0.4060 º
º Model size: Observations =     491, Parameters =  10, Deg.Fr. =   481 º
º Residuals:  Sum of squares=    65.6431     Std.Dev. =         0.36942 º
º Fit:        R-squared = 0.17048, Adjusted R-squared =         0.15496 º
º             Note:  Not using OLS.  R-squared is not bounded in [0,1]  º
º Model test: F[  9,    481] =   10.98,    Prob value =         0.00000 º
º Diagnostic: Log-L =   -202.7008, Restricted(á=0) Log-L =    -253.6384 º
º             Amemiya Pr. Crt.=    0.139, Akaike Info. Crt.=      0.866 º
º Full sample contains        1528 observations.                        º
º Selected sample contains     491 observations.                        º
º Estimated disturbance standard deviation =        0.305924            º
º Estimated correlation with selection equation A =-0.806042            º
º Estimated correlation with selection equation B =-0.389199            º
ÈÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍͼ

  Variable  Coefficient   Standard Error  z=b/s.e. P[³Z³òz]   Mean of X
  ÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄ
  Constant   0.96536         0.25076        3.850   0.00012
  SCH1       0.42614E-01     0.22800E-01    1.869   0.06162
  EX         0.38526E-01     0.10722E-01    3.593   0.00033
  EXSQ      -0.77022E-03     0.22947E-03   -3.357   0.00079
  PT         0.14534         0.38861E-01    3.740   0.00018
  Y89        0.61881E-01     0.59280E-01    1.044   0.29654
  Y90       -0.44435E-01     0.48485E-01   -0.916   0.35942
  Y91        0.13292E-01     0.57545E-01    0.231   0.81733
  Lambda-A  -0.24659         0.22275       -1.107   0.26828
  Lambda-B  -0.11907         0.81446E-01   -1.462   0.14377

ÉÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍ»
Structural model
              ÉÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍ»
              º FIML Estimates of Bivariate Probit Model    º
              º Maximum Likelihood Estimates                º
              º Dependent variable                OCCFT     º
              º Number of observations             1471     º
              º Iterations completed                 11     º
              º Log likelihood function       -964.3127     º
              º Selection model based on FT                 º
              º Means for vars. 1-  7 are after selection.  º
              ÈÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍͼ

  Variable  Coefficient   Standard Error  z=b/s.e. P[³Z³òz]   Mean of X
  ÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄ
          Index    equation for OCCUPN
  Constant   0.99012         0.71458        1.386   0.16587
  EXRELW     -4.6512          1.2707       -3.660   0.00025  0.6600
  EX        -0.20379E-01     0.16007E-01   -1.273   0.20297   21.79
  REG1      -0.46303E-01     0.50907       -0.091   0.92753  0.1291
  REG2       0.78107         0.55932        1.396   0.16258  0.3083E-01
  REG3      -0.74368E-01     0.45386       -0.164   0.86984  0.4547
  REG4      -0.43050         0.52022       -0.828   0.40793  0.3410
          Index    equation for FT
  Constant   0.76140         0.20168        3.775   0.00016
  EX        -0.10251E-01     0.36383E-02   -2.818   0.00484   22.64
  HHSIZE    -0.28353         0.31060E-01   -9.129   0.00000   3.382
  LOCUNEMP  -0.31852E-02     0.14635E-01   -0.218   0.82771   7.629
  Y89       -0.24229E-01     0.99321E-01   -0.244   0.80727  0.1924
  Y90        0.20115         0.88524E-01    2.272   0.02307  0.2862
  Y91        0.25343         0.10814        2.343   0.01911  0.1740
          Disturbance correlation and std. dev.
  RHO(1,2)   0.81966         0.29543        2.774   0.00553

    MGLFX=R5*R6*BBB1(2)
    MGLFX   = -0.250852E-08





United Kingdom
Men

Reduced form model
              ÉÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍ»
              º FIML Estimates of Bivariate Probit Model    º
              º Maximum Likelihood Estimates                º
              º Dependent variable                OCCFT     º
              º Number of observations             1257     º
              º Iterations completed                 17     º
              º Log likelihood function       -986.6273     º
              º Selection model based on FT                 º
              º Means for vars. 1- 17 are after selection.  º
              ÈÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍͼ

  Variable  Coefficient   Standard Error  z=b/s.e. P[³Z³òz]   Mean of X
  ÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄ
          Index    equation for OCCUPN
  Constant   -7.5994          1.6130       -4.711   0.00000
  SCH1       0.53307         0.10318        5.166   0.00000   11.55
  EX         0.75911E-01     0.22796E-01    3.330   0.00087   20.88
  EXSQ      -0.12149E-02     0.51695E-03   -2.350   0.01877   537.0
  PT         0.57466E-01     0.97411E-01    0.590   0.55524  0.2938
  Y89        0.23107         0.12227        1.890   0.05878  0.2812
  Y90        0.89538E-01     0.13218        0.677   0.49814  0.2192
  Y91        0.20601         0.17507        1.177   0.23930  0.2363
  SCOT       0.10776         0.23861        0.452   0.65156  0.7727E-01
  WALES     -0.28969         0.30650       -0.945   0.34458  0.3684E-01
  NTH        0.33693         0.26611        1.266   0.20546  0.6289E-01
  YH         0.10046         0.20847        0.482   0.62988  0.1042
  SW         0.15712         0.14359        1.094   0.27386  0.1015
  EA        -0.45025         0.26253       -1.715   0.08633  0.4043E-01
  WM         0.22251         0.17089        1.302   0.19290  0.9434E-01
  EM         0.61120E-01     0.18090        0.338   0.73546  0.7727E-01
  NW         0.23420         0.20576        1.138   0.25504  0.9704E-01
          Index    equation for FT
  Constant    2.0369         0.27362        7.444   0.00000
  EX        -0.63759E-02     0.44301E-02   -1.439   0.15009   21.02
  HHSIZE     0.12968E-01     0.34037E-01    0.381   0.70322   3.197
  LOCUNEMP  -0.76366E-01     0.22634E-01   -3.374   0.00074   7.014
  Y89       -0.20344         0.14204       -1.432   0.15206  0.2768
  Y90       -0.30873         0.15265       -2.022   0.04313  0.2180
  Y91       -0.33060         0.13159       -2.512   0.01199  0.2482
          Disturbance correlation and std. dev.
  RHO(1,2)  -0.38766          1.5424       -0.251   0.80156

ÉÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍ»

Knowledge workers
ÉÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍ»
º Bivariate Probit Sample Selection Model                               º
º Two stage   least squares regression    Weighting variable = ONE      º
º Dependent variable is SALARY    Mean =    2.61724, S.D. =      0.4565 º
º Model size: Observations =     336, Parameters =  10, Deg.Fr. =   326 º
º Residuals:  Sum of squares=    54.2810     Std.Dev. =         0.40805 º
º Fit:        R-squared = 0.19855, Adjusted R-squared =         0.17643 º
º             Note:  Not using OLS.  R-squared is not bounded in [0,1]  º
º Model test: F[  9,    326] =    8.97,    Prob value =         0.00000 º
º Diagnostic: Log-L =   -170.5098, Restricted(á=0) Log-L =    -212.7706 º
º             Amemiya Pr. Crt.=    0.171, Akaike Info. Crt.=      1.074 º
º Full sample contains        1257 observations.                        º
º Selected sample contains     336 observations.                        º
º Estimated disturbance standard deviation =        0.670151            º
º Estimated correlation with selection equation A = 0.467440            º
º Estimated correlation with selection equation B = -1.28562            º
ÈÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍͼ

  Variable  Coefficient   Standard Error  z=b/s.e. P[³Z³òz]   Mean of X
  ÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄ
  Constant  -0.21018          1.9118       -0.110   0.91246
  SCH1       0.14307         0.11060        1.294   0.19580
  EX         0.73506E-01     0.19339E-01    3.801   0.00014
  EXSQ      -0.14067E-02     0.39696E-03   -3.544   0.00039
  PT         0.22042         0.60104E-01    3.667   0.00025
  Y89        0.73781E-01     0.89378E-01    0.825   0.40909
  Y90        0.20817         0.88808E-01    2.344   0.01907
  Y91        0.25191         0.94776E-01    2.658   0.00786
  Lambda-A   0.31326         0.41178        0.761   0.44682
  Lambda-B  -0.86156         0.54694       -1.575   0.11520

ÉÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍ»
Other workers
ÉÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍ»
º Bivariate Probit Sample Selection Model                               º
º Two stage   least squares regression    Weighting variable = ONE      º
º Dependent variable is SALARY    Mean =    2.20996, S.D. =      0.4332 º
º Model size: Observations =     777, Parameters =  10, Deg.Fr. =   767 º
º Residuals:  Sum of squares=    122.941     Std.Dev. =         0.40036 º
º Fit:        R-squared = 0.14491, Adjusted R-squared =         0.13487 º
º             Note:  Not using OLS.  R-squared is not bounded in [0,1]  º
º Model test: F[  9,    767] =   14.44,    Prob value =         0.00000 º
º Diagnostic: Log-L =   -386.2226, Restricted(á=0) Log-L =    -452.0732 º
º             Amemiya Pr. Crt.=    0.162, Akaike Info. Crt.=      1.020 º
º Full sample contains        1257 observations.                        º
º Selected sample contains     777 observations.                        º
º Estimated disturbance standard deviation =        0.560742            º
º Estimated correlation with selection equation A =-0.358451            º
º Estimated correlation with selection equation B = -1.59623            º
ÈÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍͼ

  Variable  Coefficient   Standard Error  z=b/s.e. P[³Z³òz]   Mean of X
  ÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄ
  Constant   0.54139         0.74310        0.729   0.46627
  SCH1       0.11183         0.63379E-01    1.764   0.07765
  EX         0.37549E-01     0.10819E-01    3.471   0.00052
  EXSQ      -0.56828E-03     0.20612E-03   -2.757   0.00583
  PT         0.10770         0.34738E-01    3.100   0.00193
  Y89       -0.68130E-01     0.56372E-01   -1.209   0.22683
  Y90        0.81228E-01     0.56990E-01    1.425   0.15407
  Y91        0.16367         0.69887E-01    2.342   0.01918
  Lambda-A  -0.20100         0.23808       -0.844   0.39854
  Lambda-B  -0.89507         0.50393       -1.776   0.07570

ÉÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍ»
             
Structural model
		 ÉÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍ»
              º FIML Estimates of Bivariate Probit Model    º
              º Maximum Likelihood Estimates                º
              º Dependent variable                OCCFT     º
              º Number of observations             1257     º
              º Iterations completed                 18     º
              º Log likelihood function       -981.9258     º
              º Selection model based on FT                 º
              º Means for vars. 1- 12 are after selection.  º
              ÈÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍͼ

  Variable  Coefficient   Standard Error  z=b/s.e. P[³Z³òz]   Mean of X
  ÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄ
          Index    equation for OCCUPN
  Constant  -0.37695         0.11099       -3.396   0.00068
  EXRELW      5.7157         0.42736       13.374   0.00000  0.3178E-01
  EX        -0.23126E-01     0.41970E-02   -5.510   0.00000   20.88
  SCOT      -0.10999E-01     0.15177       -0.072   0.94223  0.7727E-01
  WALES     -0.39653         0.26704       -1.485   0.13758  0.3684E-01
  NTH       -0.63872         0.19105       -3.343   0.00083  0.6289E-01
  YH        -0.51433         0.16952       -3.034   0.00241  0.1042
  SW         0.14446         0.14074        1.026   0.30472  0.1015
  EA         0.14326E-01     0.19509        0.073   0.94146  0.4043E-01
  WM        -0.12390         0.15616       -0.793   0.42755  0.9434E-01
  EM        -0.32516         0.15944       -2.039   0.04141  0.7727E-01
  NW        -0.33877E-02     0.13525       -0.025   0.98002  0.9704E-01
          Index    equation for FT
  Constant    1.7586         0.27556        6.382   0.00000
  EX        -0.65621E-02     0.43790E-02   -1.499   0.13399   21.02
  HHSIZE     0.55341E-01     0.33217E-01    1.666   0.09570   3.197
  LOCUNEMP  -0.69327E-01     0.22524E-01   -3.078   0.00208   7.014
  Y89       -0.21539E-01     0.14150       -0.152   0.87901  0.2768
  Y90       -0.73813E-01     0.15048       -0.491   0.62377  0.2180
  Y91       -0.28313         0.13252       -2.136   0.03264  0.2482
          Disturbance correlation and std. dev.
  RHO(1,2)   0.95584         0.87006        1.099   0.27195

ÉÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍ»
    MGLFX=R5*R6*BBB1(2)
    MGLFX   =  0.262677E-05

Single women

Reduced form model
              ÉÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍ»
              º FIML Estimates of Bivariate Probit Model    º
              º Maximum Likelihood Estimates                º
              º Dependent variable                OCCFT     º
              º Number of observations              372     º
              º Iterations completed                  2     º
              º Log likelihood function       -398.0928     º
              º Selection model based on FT                 º
              º Means for vars. 1- 17 are after selection.  º
              ÈÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍͼ

  Variable  Coefficient   Standard Error  z=b/s.e. P[³Z³òz]   Mean of X
  ÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄ
          Index    equation for OCCUPN
  Constant   -1.7096          1.5584       -1.097   0.27263
  SCH1       0.15309         0.69406E-01    2.206   0.02740   11.73
  EX         0.24958E-01     0.35161E-01    0.710   0.47781   17.57
  EXSQ      -0.12469E-02     0.86530E-03   -1.441   0.14958   426.7
  PT        -0.83525E-01     0.18701       -0.447   0.65514  0.6913
  Y89       -0.50018E-01     0.26973       -0.185   0.85289  0.2087
  Y90        0.74639E-01     0.24272        0.308   0.75845  0.2261
  Y91       -0.64723E-01     0.27724       -0.233   0.81541  0.2826
  SCOT      -0.66991E-02     0.41509       -0.016   0.98712  0.8261E-01
  WALES      0.15324         0.49520        0.309   0.75698  0.3478E-01
  NTH        0.73609E-01     0.51087        0.144   0.88543  0.6522E-01
  YH        -0.23444         0.43813       -0.535   0.59259  0.6957E-01
  SW        -0.84153E-01     0.36875       -0.228   0.81948  0.6957E-01
  EA        -0.76888E-01     0.47891       -0.161   0.87245  0.3478E-01
  WM        -0.39860E-02     0.33954       -0.012   0.99063  0.1043
  EM        -0.11324         0.36137       -0.313   0.75401  0.7391E-01
  NW         0.70142E-02     0.33635        0.021   0.98336  0.1261
          Index    equation for FT
  Constant    1.1534         0.36785        3.136   0.00172
  EX        -0.61825E-02     0.62406E-02   -0.991   0.32183   18.15
  HHSIZE    -0.43134E-01     0.53117E-01   -0.812   0.41676   2.702
  LOCUNEMP  -0.36453E-01     0.29870E-01   -1.220   0.22232   7.412
  Y89       -0.11599         0.19891       -0.583   0.55982  0.2097
  Y90       -0.12289         0.19617       -0.626   0.53102  0.2285
  Y91       -0.11019         0.17661       -0.624   0.53266  0.3011
          Disturbance correlation and std. dev.
  RHO(1,2)  -0.69981E-02      3.1577       -0.002   0.99823

ÉÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍ»
Knowledge workers
ÉÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍ»
º Bivariate Probit Sample Selection Model                               º
º Two stage   least squares regression    Weighting variable = ONE      º
º Dependent variable is SALARY    Mean =    2.09876, S.D. =      0.4618 º
º Model size: Observations =     117, Parameters =  10, Deg.Fr. =   107 º
º Residuals:  Sum of squares=    16.0524     Std.Dev. =         0.38733 º
º Fit:        R-squared = 0.29049, Adjusted R-squared =         0.23081 º
º             Note:  Not using OLS.  R-squared is not bounded in [0,1]  º
º Model test: F[  9,    107] =    4.87,    Prob value =         0.00002 º
º Diagnostic: Log-L =    -49.8164, Restricted(á=0) Log-L =     -75.1188 º
º             Amemiya Pr. Crt.=    0.163, Akaike Info. Crt.=      1.023 º
º Full sample contains         372 observations.                        º
º Selected sample contains     117 observations.                        º
º Estimated disturbance standard deviation =         2.09865            º
º Estimated correlation with selection equation A =-0.564983E-01        º
º Estimated correlation with selection equation B = -1.42030            º
ÈÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍͼ

  Variable  Coefficient   Standard Error  z=b/s.e. P[³Z³òz]   Mean of X
  ÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄ
  Constant    2.6776          1.9323        1.386   0.16584
  SCH1       0.21849E-01     0.70637E-01    0.309   0.75708
  EX         0.25511E-01     0.27646E-01    0.923   0.35613
  EXSQ      -0.17178E-03     0.71600E-03   -0.240   0.81039
  PT         0.16446         0.85226E-01    1.930   0.05365
  Y89        0.13820         0.46575        0.297   0.76667
  Y90        0.31113         0.52358        0.594   0.55235
  Y91        0.45710         0.45802        0.998   0.31828
  Lambda-A  -0.11857         0.76224       -0.156   0.87638
  Lambda-B   -2.9807          3.5050       -0.850   0.39509

ÉÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍ»
Other workers
ÉÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍ»
º Bivariate Probit Sample Selection Model                               º
º Two stage   least squares regression    Weighting variable = ONE      º
º Dependent variable is SALARY    Mean =    1.93828, S.D. =      0.5192 º
º Model size: Observations =     113, Parameters =  10, Deg.Fr. =   103 º
º Residuals:  Sum of squares=    20.5825     Std.Dev. =         0.44702 º
º Fit:        R-squared = 0.25212, Adjusted R-squared =         0.18677 º
º             Note:  Not using OLS.  R-squared is not bounded in [0,1]  º
º Model test: F[  9,    103] =    3.86,    Prob value =         0.00031 º
º Diagnostic: Log-L =    -64.1235, Restricted(á=0) Log-L =     -85.7728 º
º             Amemiya Pr. Crt.=    0.218, Akaike Info. Crt.=      1.312 º
º Full sample contains         372 observations.                        º
º Selected sample contains     113 observations.                        º
º Selectivity corrected std. error for the á on EXSQ     is negative.   º
º The set of OLS standard errors is given.                              º
º Estimated disturbance standard deviation =        0.344742            º
º Estimated correlation with selection equation A = -4.26603            º
º Estimated correlation with selection equation B = -8.10781            º
ÈÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍͼ

  Variable  Coefficient   Standard Error  z=b/s.e. P[³Z³òz]   Mean of X
  ÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄ
  Constant   0.59836         0.80873        0.740   0.45937
  SCH1       0.86829E-01     0.95617E-01    0.908   0.36383
  EX         0.10130E-01     0.23651E-01    0.428   0.66841
  EXSQ       0.60949E-03     0.75919E-03    0.803   0.42208
  PT         0.16633         0.13083        1.271   0.20361
  Y89        0.17897E-01     0.13047        0.137   0.89089
  Y90        0.11685         0.16021        0.729   0.46579
  Y91        0.37603         0.12521        3.003   0.00267
  Lambda-A   -1.4707          1.1104       -1.324   0.18535
  Lambda-B   -2.7951          1.0529       -2.655   0.00794

ÉÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍ»
Structural model
              ÉÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍ»
              º FIML Estimates of Bivariate Probit Model    º
              º Maximum Likelihood Estimates                º
              º Dependent variable                OCCFT     º
              º Number of observations              372     º
              º Iterations completed                  3     º
              º Log likelihood function       -343.4179     º
              º Selection model based on FT                 º
              º Means for vars. 1- 12 are after selection.  º
              ÈÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍͼ

  Variable  Coefficient   Standard Error  z=b/s.e. P[³Z³òz]   Mean of X
  ÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄ
          Index    equation for OCCUPN
  Constant    2.5791          1.9643        1.313   0.18919
  EXRELW    -0.97663          1.0343       -0.944   0.34505   1.345
  EX        -0.72669E-01     0.29442E-01   -2.468   0.01358   17.57
  SCOT      -0.21328         0.49068       -0.435   0.66380  0.8261E-01
  WALES      0.26796         0.68975        0.388   0.69766  0.3478E-01
  NTH       -0.19410E-01     0.66355       -0.029   0.97666  0.6522E-01
  YH        -0.15462         0.48719       -0.317   0.75096  0.6957E-01
  SW        -0.15636         0.46313       -0.338   0.73564  0.6957E-01
  EA        -0.44566E-01     0.48000       -0.093   0.92603  0.3478E-01
  WM        -0.81104E-01     0.43460       -0.187   0.85196  0.1043
  EM        -0.15049         0.41781       -0.360   0.71871  0.7391E-01
  NW         0.99790E-01     0.45992        0.217   0.82823  0.1261
          Index    equation for FT
  Constant    1.1528         0.38337        3.007   0.00264
  EX        -0.13928E-01     0.67420E-02   -2.066   0.03884   18.15
  HHSIZE    -0.45316E-01     0.56096E-01   -0.808   0.41919   2.702
  LOCUNEMP  -0.41893E-01     0.30983E-01   -1.352   0.17633   7.412
  Y89       -0.11611         0.20321       -0.571   0.56774  0.2097
  Y90       -0.12304         0.20502       -0.600   0.54843  0.2285
  Y91       -0.11045         0.18207       -0.607   0.54409  0.3011
          Disturbance correlation and std. dev.
  RHO(1,2)  -0.78759E-02      2.9816       -0.003   0.99789

ÉÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍ»
    MGLFX=R5*R6*BBB1(2)
    MGLFX   = -0.187983

Married women

Reduced form model
              ÉÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍ»
              º FIML Estimates of Bivariate Probit Model    º
              º Maximum Likelihood Estimates                º
              º Dependent variable                OCCFT     º
              º Number of observations             1279     º
              º Iterations completed                 30     º
              º Log likelihood function       -1186.984     º
              º Selection model based on FT                 º
              º Means for vars. 1- 17 are after selection.  º
              ÈÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍͼ

  Variable  Coefficient   Standard Error  z=b/s.e. P[³Z³òz]   Mean of X
  ÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄ
          Index    equation for OCCUPN
  Constant   -4.8420         0.60910       -7.949   0.00000
  SCH1       0.35905         0.47346E-01    7.583   0.00000   11.43
  EX         0.17972E-01     0.27862E-01    0.645   0.51890   22.75
  EXSQ      -0.16707E-03     0.62795E-03   -0.266   0.79020   604.4
  PT        -0.49546         0.13022       -3.805   0.00014  0.8262
  Y89       -0.11399E-01     0.14005       -0.081   0.93513  0.2660
  Y90        0.87976E-01     0.13981        0.629   0.52919  0.2647
  Y91        0.41624E-01     0.15533        0.268   0.78872  0.2072
  SCOT      -0.60247E-01     0.20279       -0.297   0.76640  0.9225E-01
  WALES     -0.26626         0.23395       -1.138   0.25507  0.4679E-01
  NTH       -0.59478E-02     0.22965       -0.026   0.97934  0.6016E-01
  YH        -0.43799         0.22072       -1.984   0.04721  0.9358E-01
  SW        -0.27732         0.18502       -1.499   0.13392  0.8690E-01
  EA         0.32514         0.24726        1.315   0.18853  0.4011E-01
  WM        -0.26154         0.17344       -1.508   0.13156  0.1083
  EM        -0.25079         0.21472       -1.168   0.24281  0.6684E-01
  NW        -0.20590         0.18561       -1.109   0.26729  0.9492E-01
          Index    equation for FT
  Constant    1.2496         0.22425        5.572   0.00000
  EX         0.24195E-02     0.38667E-02    0.626   0.53150   22.38
  HHSIZE    -0.28594         0.34396E-01   -8.313   0.00000   3.450
  LOCUNEMP  -0.91138E-02     0.16519E-01   -0.552   0.58114   6.917
  Y89       -0.53074E-01     0.10380       -0.511   0.60915  0.2697
  Y90       -0.22244E-01     0.10967       -0.203   0.83927  0.2604
  Y91       -0.46056E-01     0.10719       -0.430   0.66743  0.2095
          Disturbance correlation and std. dev.
  RHO(1,2)   0.57762         0.29959        1.928   0.05385

ÉÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍ»
Knowledge workers
ÉÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍ»
º Bivariate Probit Sample Selection Model                               º
º Two stage   least squares regression    Weighting variable = ONE      º
º Dependent variable is SALARY    Mean =    2.30590, S.D. =      0.4948 º
º Model size: Observations =     124, Parameters =  10, Deg.Fr. =   114 º
º Residuals:  Sum of squares=    22.8332     Std.Dev. =         0.44754 º
º Fit:        R-squared = 0.17530, Adjusted R-squared =         0.11020 º
º             Note:  Not using OLS.  R-squared is not bounded in [0,1]  º
º Model test: F[  9,    114] =    2.69,    Prob value =         0.00708 º
º Diagnostic: Log-L =    -71.0402, Restricted(á=0) Log-L =     -88.2033 º
º             Amemiya Pr. Crt.=    0.216, Akaike Info. Crt.=      1.307 º
º Full sample contains        1279 observations.                        º
º Selected sample contains     124 observations.                        º
º Estimated disturbance standard deviation =        0.810984            º
º Estimated correlation with selection equation A = -1.11323            º
º Estimated correlation with selection equation B =-0.454970            º
ÈÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍͼ

  Variable  Coefficient   Standard Error  z=b/s.e. P[³Z³òz]   Mean of X
  ÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄ
  Constant    2.1082          3.5632        0.592   0.55409
  SCH1       0.17487E-01     0.21301        0.082   0.93457
  EX         0.48468E-01     0.44652E-01    1.085   0.27771
  EXSQ      -0.12749E-02     0.10326E-02   -1.235   0.21693
  PT         0.27624         0.29713        0.930   0.35253
  Y89        0.80801E-02     0.23760        0.034   0.97287
  Y90        0.28605         0.22165        1.291   0.19686
  Y91        0.19782         0.22626        0.874   0.38197
  Lambda-A  -0.90281         0.94990       -0.950   0.34190
  Lambda-B  -0.36897         0.56400       -0.654   0.51298

ÉÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍ»
Other workers
ÉÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍ»
º Bivariate Probit Sample Selection Model                               º
º Two stage   least squares regression    Weighting variable = ONE      º
º Dependent variable is SALARY    Mean =    1.81907, S.D. =      0.5080 º
º Model size: Observations =     624, Parameters =  10, Deg.Fr. =   614 º
º Residuals:  Sum of squares=    135.600     Std.Dev. =         0.46994 º
º Fit:        R-squared = 0.14291, Adjusted R-squared =         0.13035 º
º             Note:  Not using OLS.  R-squared is not bounded in [0,1]  º
º Model test: F[  9,    614] =   11.38,    Prob value =         0.00000 º
º Diagnostic: Log-L =   -409.1672, Restricted(á=0) Log-L =    -462.3224 º
º             Amemiya Pr. Crt.=    0.224, Akaike Info. Crt.=      1.343 º
º Full sample contains        1279 observations.                        º
º Selected sample contains     624 observations.                        º
º Estimated disturbance standard deviation =        0.918030E-01        º
º Estimated correlation with selection equation A = -6.21481            º
º Estimated correlation with selection equation B = -4.01909            º
ÈÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍͼ

  Variable  Coefficient   Standard Error  z=b/s.e. P[³Z³òz]   Mean of X
  ÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄ
  Constant   0.79013         0.60185        1.313   0.18924
  SCH1       0.64658E-01     0.59434E-01    1.088   0.27664
  EX         0.24836E-01     0.12140E-01    2.046   0.04078
  EXSQ      -0.49819E-03     0.24455E-03   -2.037   0.04163
  PT         0.28575E-01     0.96115E-01    0.297   0.76624
  Y89        0.82678E-01     0.79180E-01    1.044   0.29640
  Y90        0.75594E-01     0.85465E-01    0.885   0.37642
  Y91        0.18285         0.89942E-01    2.033   0.04206
  Lambda-A  -0.57054         0.36230       -1.575   0.11531
  Lambda-B  -0.36897         0.22198       -1.662   0.09648

ÉÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍ»
Structural model
              ÉÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍ»
              º FIML Estimates of Bivariate Probit Model    º
              º Maximum Likelihood Estimates                º
              º Dependent variable                OCCFT     º
              º Number of observations             1279     º
              º Iterations completed                  2     º
              º Log likelihood function       -1151.235     º
              º Selection model based on FT                 º
              º Means for vars. 1- 12 are after selection.  º
              ÈÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍͼ

  Variable  Coefficient   Standard Error  z=b/s.e. P[³Z³òz]   Mean of X
  ÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄ
          Index    equation for OCCUPN
  Constant    1.5097         0.63032        2.395   0.01661
  EXRELW    -0.90399         0.34820       -2.596   0.00943   1.090
  EX        -0.60515E-01     0.93264E-02   -6.489   0.00000   22.75
  SCOT      -0.35844E-01     0.20548       -0.174   0.86152  0.9225E-01
  WALES     -0.15010         0.31335       -0.479   0.63194  0.4679E-01
  NTH        0.39348E-01     0.22420        0.176   0.86068  0.6016E-01
  YH        -0.19680         0.31943       -0.616   0.53782  0.9358E-01
  SW        -0.17404         0.24447       -0.712   0.47652  0.8690E-01
  EA         0.13963         0.24634        0.567   0.57085  0.4011E-01
  WM        -0.18998         0.23159       -0.820   0.41202  0.1083
  EM        -0.16545         0.32781       -0.505   0.61375  0.6684E-01
  NW        -0.15035         0.26376       -0.570   0.56866  0.9492E-01
          Index    equation for FT
  Constant   0.96233         0.21833        4.408   0.00001
  EX        -0.91605E-02     0.38151E-02   -2.401   0.01635   22.38
  HHSIZE    -0.11264         0.33420E-01   -3.370   0.00075   3.450
  LOCUNEMP  -0.10248E-01     0.16315E-01   -0.628   0.52993   6.917
  Y89       -0.23611E-01     0.10285       -0.230   0.81843  0.2697
  Y90       -0.12776E-01     0.10786       -0.118   0.90571  0.2604
  Y91       -0.18384E-01     0.10500       -0.175   0.86101  0.2095
          Disturbance correlation and std. dev.
  RHO(1,2)  -0.10609E-01     0.98914       -0.011   0.99144

ÉÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍ»
    MGLFX=R5*R6*BBB1(2)
    MGLFX   = -0.659863E-01

USA
Men
Reduced form model
              ÉÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍ»
              º FIML Estimates of Bivariate Probit Model    º
              º Maximum Likelihood Estimates                º
              º Dependent variable                OCCFT     º
              º Number of observations             1391     º
              º Iterations completed                  2     º
              º Log likelihood function       -1249.202     º
              º Selection model based on FT                 º
              º Means for vars. 1- 16 are after selection.  º
              ÈÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍͼ

  Variable  Coefficient   Standard Error  z=b/s.e. P[³Z³òz]   Mean of X
  ÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄ
          Index    equation for OCCUPN
  Constant  -0.93445          1.8347       -0.509   0.61053
  SCH1       0.80756E-01     0.16201E-01    4.985   0.00000   13.77
  EX        -0.65628E-02     0.20738E-01   -0.316   0.75165   16.46
  EXSQ      -0.13091E-02     0.48165E-03   -2.718   0.00657   353.9
  PT        -0.94759E-02     0.11832       -0.080   0.93617  0.1492
  Y89        0.35115E-01     0.11165        0.315   0.75314  0.2770
  Y90        0.20728E-01     0.11618        0.178   0.85841  0.2218
  Y91        0.65616E-01     0.12999        0.505   0.61371  0.2366
  NEWENG     0.31412E-01     0.21583        0.146   0.88428  0.4534E-01
  MIDAT      0.18394E-01     0.15617        0.118   0.90624  0.1204
  ENC        0.26156E-01     0.13835        0.189   0.85005  0.1880
  WNC        0.81178E-01     0.16932        0.479   0.63162  0.8821E-01
  SAT        0.74759E-02     0.14328        0.052   0.95839  0.1764
  ESC        0.33934E-01     0.19074        0.178   0.85879  0.6760E-01
  WSC        0.67851E-01     0.16994        0.399   0.68971  0.1031
  MOUNT      0.35072E-01     0.16659        0.211   0.83326  0.7007E-01
          Index    equation for FT
  Constant   0.72927         0.30141        2.420   0.01554
  EX         0.40717E-02     0.39388E-02    1.034   0.30125   16.27
  HHSIZE     0.11900E-01     0.31077E-01    0.383   0.70179   2.922
  LOCUNEMP   0.15917E-01     0.49073E-01    0.324   0.74567   6.170
  Y89        0.12846E-01     0.12717        0.101   0.91954  0.2732
  Y90       -0.14977E-02     0.14436       -0.010   0.99172  0.2185
  Y91       -0.53696E-01     0.15720       -0.342   0.73267  0.2444
          Disturbance correlation and std. dev.
  RHO(1,2)  -0.49049E-04      5.2819        0.000   0.99999

ÉÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍ»
Knowledge workers
ÉÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍ»
º Bivariate Probit Sample Selection Model                               º
º Two stage   least squares regression    Weighting variable = ONE      º
º Dependent variable is SALARY    Mean =    2.31037, S.D. =      0.6816 º
º Model size: Observations =     393, Parameters =  10, Deg.Fr. =   383 º
º Residuals:  Sum of squares=    125.103     Std.Dev. =         0.57152 º
º Fit:        R-squared = 0.29508, Adjusted R-squared =         0.27851 º
º             Note:  Not using OLS.  R-squared is not bounded in [0,1]  º
º Model test: F[  9,    383] =   17.81,    Prob value =         0.00000 º
º Diagnostic: Log-L =   -332.7151, Restricted(á=0) Log-L =    -406.4900 º
º             Amemiya Pr. Crt.=    0.335, Akaike Info. Crt.=      1.744 º
º Full sample contains        1391 observations.                        º
º Selected sample contains     393 observations.                        º
º Estimated disturbance standard deviation =         3.24787            º
º Estimated correlation with selection equation A =  1.28000            º
º Estimated correlation with selection equation B = 0.193421            º
ÈÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍͼ

  Variable  Coefficient   Standard Error  z=b/s.e. P[³Z³òz]   Mean of X
  ÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄ
  Constant   -5.4955          18.542       -0.296   0.76694
  SCH1       0.28214         0.64090        0.440   0.65978
  EX         0.11770         0.12332        0.954   0.33986
  EXSQ      -0.81409E-02     0.11656E-01   -0.698   0.48490
  PT         0.72093E-01     0.64289        0.112   0.91071
  Y89        0.45775E-01     0.65788        0.070   0.94453
  Y90        0.20307         0.65386        0.311   0.75613
  Y91        0.33302         0.81075        0.411   0.68125
  Lambda-A    4.1573          13.944        0.298   0.76560
  Lambda-B   0.62821          12.053        0.052   0.95843

ÉÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍ»
Other workers
ÉÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍ»
º Bivariate Probit Sample Selection Model                               º
º Two stage   least squares regression    Weighting variable = ONE      º
º Dependent variable is SALARY    Mean =    2.33446, S.D. =      0.6262 º
º Model size: Observations =     820, Parameters =  10, Deg.Fr. =   810 º
º Residuals:  Sum of squares=    264.420     Std.Dev. =         0.57135 º
º Fit:        R-squared = 0.16660, Adjusted R-squared =         0.15734 º
º             Note:  Not using OLS.  R-squared is not bounded in [0,1]  º
º Model test: F[  9,    810] =   17.99,    Prob value =         0.00000 º
º Diagnostic: Log-L =   -699.5052, Restricted(á=0) Log-L =    -779.2543 º
º             Amemiya Pr. Crt.=    0.330, Akaike Info. Crt.=      1.731 º
º Full sample contains        1391 observations.                        º
º Selected sample contains     820 observations.                        º
º Estimated disturbance standard deviation =         1.50963            º
º Estimated correlation with selection equation A = 0.198776            º
º Estimated correlation with selection equation B = -1.32951            º
ÈÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍͼ

  Variable  Coefficient   Standard Error  z=b/s.e. P[³Z³òz]   Mean of X
  ÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄ
  Constant    1.4953          1.7026        0.878   0.37983
  SCH1       0.93963E-01     0.23429E-01    4.011   0.00006
  EX         0.31608E-01     0.19859E-01    1.592   0.11147
  EXSQ      -0.68305E-03     0.32590E-03   -2.096   0.03610
  PT         0.20874         0.60972E-01    3.423   0.00062
  Y89        0.86770E-01     0.12313        0.705   0.48099
  Y90        0.36917E-01     0.12726        0.290   0.77174
  Y91        0.16050         0.13188        1.217   0.22358
  Lambda-A   0.30008         0.57361        0.523   0.60088
  Lambda-B   -2.0071          4.8863       -0.411   0.68125

ÉÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍ»
Structural model
              ÉÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍ»
              º FIML Estimates of Bivariate Probit Model    º
              º Maximum Likelihood Estimates                º
              º Dependent variable                OCCFT     º
              º Number of observations             1391     º
              º Iterations completed                  3     º
              º Log likelihood function       -999.7044     º
              º Selection model based on FT                 º
              º Means for vars. 1- 11 are after selection.  º
              ÈÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍͼ

  Variable  Coefficient   Standard Error  z=b/s.e. P[³Z³òz]   Mean of X
  ÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄ
          Index    equation for OCCUPN
  Constant   0.87719          1.7885        0.490   0.62381
  EXRELW     0.12179E-01     0.96152E-01    0.127   0.89921  -5.575
  EX        -0.98367E-01     0.27920E-01   -3.523   0.00043   16.46
  NEWENG    -0.92680E-02     0.24912       -0.037   0.97032  0.4534E-01
  MIDAT     -0.61779E-01     0.19404       -0.318   0.75020  0.1204
  ENC       -0.24814E-01     0.16166       -0.153   0.87801  0.1880
  WNC        0.24925E-01     0.20803        0.120   0.90463  0.8821E-01
  SAT       -0.68821E-01     0.16882       -0.408   0.68353  0.1764
  ESC       -0.10728         0.27175       -0.395   0.69300  0.6760E-01
  WSC       -0.92512E-01     0.20523       -0.451   0.65215  0.1031
  MOUNT     -0.35483E-01     0.22889       -0.155   0.87681  0.7007E-01
          Index    equation for FT
  Constant   0.72995         0.28758        2.538   0.01114
  EX         0.13002E-01     0.36483E-02    3.564   0.00037   16.27
  HHSIZE     0.14412E-01     0.29653E-01    0.486   0.62695   2.922
  LOCUNEMP   0.20122E-01     0.47127E-01    0.427   0.66939   6.170
  Y89        0.13100E-01     0.12214        0.107   0.91459  0.2732
  Y90       -0.12859E-02     0.13870       -0.009   0.99260  0.2185
  Y91       -0.53692E-01     0.15264       -0.352   0.72501  0.2444
          Disturbance correlation and std. dev.
  RHO(1,2)   0.31502E-02      5.1987        0.001   0.99952

ÉÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍ»

    MGLFX=R5*R6*BBB1(2)
    MGLFX   =  0.542969E-03

Single women
Reduced form model
              ÉÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍ»
              º FIML Estimates of Bivariate Probit Model    º
              º Maximum Likelihood Estimates                º
              º Dependent variable                OCCFT     º
              º Number of observations              784     º
              º Iterations completed                 17     º
              º Log likelihood function       -726.1587     º
              º Selection model based on FT                 º
              º Means for vars. 1- 16 are after selection.  º
              ÈÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍͼ

  Variable  Coefficient   Standard Error  z=b/s.e. P[³Z³òz]   Mean of X
  ÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄ
          Index    equation for OCCUPN
  Constant   -4.1471         0.92877       -4.465   0.00001
  SCH1       0.27557         0.48364E-01    5.698   0.00000   13.89
  EX        -0.13926E-01     0.26089E-01   -0.534   0.59350   15.46
  EXSQ       0.50912E-03     0.77806E-03    0.654   0.51288   328.2
  PT        -0.31607         0.13704       -2.306   0.02109  0.2891
  Y89        0.19221         0.17792        1.080   0.28001  0.2527
  Y90       -0.15304         0.18665       -0.820   0.41223  0.2145
  Y91        0.30461         0.17485        1.742   0.08148  0.2509
  NEWENG     0.33577         0.29890        1.123   0.26129  0.6364E-01
  MIDAT     -0.90312E-01     0.20989       -0.430   0.66700  0.1473
  ENC        0.28360E-01     0.20152        0.141   0.88808  0.1818
  WNC       -0.11604         0.23732       -0.489   0.62487  0.8364E-01
  SAT        0.12064         0.21865        0.552   0.58111  0.1618
  ESC       -0.13795E-01     0.26289       -0.052   0.95815  0.7273E-01
  WSC        0.35617E-01     0.21255        0.168   0.86692  0.9091E-01
  MOUNT      0.39378E-01     0.43707        0.090   0.92821  0.3636E-01
          Index    equation for FT
  Constant    1.2054         0.31035        3.884   0.00010
  EX        -0.85469E-03     0.45449E-02   -0.188   0.85083   15.38
  HHSIZE    -0.18383         0.32822E-01   -5.601   0.00000   2.546
  LOCUNEMP  -0.31244E-01     0.50589E-01   -0.618   0.53684   6.130
  Y89       -0.12074         0.13051       -0.925   0.35488  0.2679
  Y90        0.12722         0.15696        0.811   0.41761  0.2003
  Y91        0.22340E-01     0.16901        0.132   0.89484  0.2538
          Disturbance correlation and std. dev.
  RHO(1,2)  -0.45925         0.41860       -1.097   0.27259

ÉÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍ»
Knowledge workers
ÉÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍ»
º Bivariate Probit Sample Selection Model                               º
º Two stage   least squares regression    Weighting variable = ONE      º
º Dependent variable is SALARY    Mean =    2.34215, S.D. =      0.5514 º
º Model size: Observations =     147, Parameters =  10, Deg.Fr. =   137 º
º Residuals:  Sum of squares=    33.2394     Std.Dev. =         0.49257 º
º Fit:        R-squared = 0.19641, Adjusted R-squared =         0.14362 º
º             Note:  Not using OLS.  R-squared is not bounded in [0,1]  º
º Model test: F[  9,    137] =    3.72,    Prob value =         0.00033 º
º Diagnostic: Log-L =    -99.3118, Restricted(á=0) Log-L =    -120.5620 º
º             Amemiya Pr. Crt.=    0.259, Akaike Info. Crt.=      1.487 º
º Full sample contains         784 observations.                        º
º Selected sample contains     147 observations.                        º
º Estimated disturbance standard deviation =        0.966873            º
º Estimated correlation with selection equation A = -1.14217            º
º Estimated correlation with selection equation B = 0.112865            º
ÈÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍͼ

  Variable  Coefficient   Standard Error  z=b/s.e. P[³Z³òz]   Mean of X
  ÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄ
  Constant    4.0632          3.0802        1.319   0.18712
  SCH1      -0.99145E-01     0.15288       -0.649   0.51665
  EX         0.64784E-01     0.27721E-01    2.337   0.01944
  EXSQ      -0.13441E-02     0.89530E-03   -1.501   0.13329
  PT         0.20499E-01     0.28879        0.071   0.94341
  Y89       -0.75568E-02     0.23764       -0.032   0.97463
  Y90        0.14718         0.27950        0.527   0.59850
  Y91        0.98920E-01     0.26143        0.378   0.70515
  Lambda-A   -1.1043          1.0781       -1.024   0.30567
  Lambda-B   0.10913         0.44615        0.245   0.80677

ÉÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍ»
Other workers
ÉÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍ»
º Bivariate Probit Sample Selection Model                               º
º Two stage   least squares regression    Weighting variable = ONE      º
º Dependent variable is SALARY    Mean =    1.94666, S.D. =      0.6139 º
º Model size: Observations =     403, Parameters =  10, Deg.Fr. =   393 º
º Residuals:  Sum of squares=    130.396     Std.Dev. =         0.57602 º
º Fit:        R-squared = 0.11752, Adjusted R-squared =         0.09731 º
º             Note:  Not using OLS.  R-squared is not bounded in [0,1]  º
º Model test: F[  9,    393] =    5.82,    Prob value =         0.00000 º
º Diagnostic: Log-L =   -344.4677, Restricted(á=0) Log-L =    -374.7219 º
º             Amemiya Pr. Crt.=    0.340, Akaike Info. Crt.=      1.759 º
º Full sample contains         784 observations.                        º
º Selected sample contains     403 observations.                        º
º Estimated disturbance standard deviation =        0.569569            º
º Estimated correlation with selection equation A =-0.535971            º
º Estimated correlation with selection equation B =-0.585540            º
ÈÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍͼ

  Variable  Coefficient   Standard Error  z=b/s.e. P[³Z³òz]   Mean of X
  ÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄ
  Constant   0.94847         0.47317        2.005   0.04501
  SCH1       0.40314E-01     0.46800E-01    0.861   0.38901
  EX         0.57985E-01     0.12095E-01    4.794   0.00000
  EXSQ      -0.13833E-02     0.32929E-03   -4.201   0.00003
  PT         0.12951         0.80761E-01    1.604   0.10881
  Y89        0.37599E-01     0.99375E-01    0.378   0.70517
  Y90        0.39342E-01     0.93503E-01    0.421   0.67393
  Y91       -0.38778E-01     0.10522       -0.369   0.71247
  Lambda-A  -0.30527         0.42401       -0.720   0.47154
  Lambda-B  -0.33351         0.27347       -1.220   0.22264

ÉÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍ»
Structural model
              ÉÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍ»
              º FIML Estimates of Bivariate Probit Model    º
              º Maximum Likelihood Estimates                º
              º Dependent variable                OCCFT     º
              º Number of observations              784     º
              º Iterations completed                  3     º
              º Log likelihood function       -712.9696     º
              º Selection model based on FT                 º
              º Means for vars. 1- 11 are after selection.  º
              ÈÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍͼ

  Variable  Coefficient   Standard Error  z=b/s.e. P[³Z³òz]   Mean of X
  ÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄ
          Index    equation for OCCUPN
  Constant    1.3157         0.69480        1.894   0.05827
  EXRELW    -0.93045         0.28588       -3.255   0.00114   1.310
  EX        -0.45070E-01     0.85501E-02   -5.271   0.00000   15.46
  NEWENG     0.31692E-02     0.29689        0.011   0.99148  0.6364E-01
  MIDAT     -0.70198E-01     0.26974       -0.260   0.79467  0.1473
  ENC        0.35344E-01     0.22680        0.156   0.87616  0.1818
  WNC       -0.38483E-01     0.27660       -0.139   0.88935  0.8364E-01
  SAT        0.10046         0.23942        0.420   0.67477  0.1618
  ESC        0.35126E-01     0.37137        0.095   0.92464  0.7273E-01
  WSC        0.82203E-01     0.30214        0.272   0.78557  0.9091E-01
  MOUNT      0.11463         0.38370        0.299   0.76513  0.3636E-01
          Index    equation for FT
  Constant   0.92671         0.31286        2.962   0.00306
  EX        -0.58030E-02     0.46125E-02   -1.258   0.20835   15.38
  HHSIZE    -0.68016E-01     0.32303E-01   -2.106   0.03525   2.546
  LOCUNEMP  -0.12801E-01     0.50866E-01   -0.252   0.80130   6.130
  Y89       -0.48832E-01     0.12875       -0.379   0.70447  0.2679
  Y90        0.44397E-01     0.15686        0.283   0.77714  0.2003
  Y91       -0.17488E-02     0.16638       -0.011   0.99161  0.2538
          Disturbance correlation and std. dev.
  RHO(1,2)  -0.29991E-02      1.5265       -0.002   0.99843

ÉÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍ»
    MGLFX=R5*R6*BBB1(2)
    MGLFX   = -0.390364E-13

Married women
Reduced form model
              ÉÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍ»
              º FIML Estimates of Bivariate Probit Model    º
              º Maximum Likelihood Estimates                º
              º Dependent variable                OCCFT     º
              º Number of observations              912     º
              º Iterations completed                 18     º
              º Log likelihood function       -877.8689     º
              º Selection model based on FT                 º
              º Means for vars. 1- 16 are after selection.  º
              ÈÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍͼ

  Variable  Coefficient   Standard Error  z=b/s.e. P[³Z³òz]   Mean of X
  ÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄ
          Index    equation for OCCUPN
  Constant   -4.5289          1.1989       -3.778   0.00016
  SCH1       0.29520         0.53985E-01    5.468   0.00000   13.49
  EX         0.18250E-01     0.28240E-01    0.646   0.51812   17.96
  EXSQ      -0.37911E-03     0.71379E-03   -0.531   0.59534   396.8
  PT        -0.22532         0.13635       -1.652   0.09844  0.4198
  Y89        0.24013         0.17416        1.379   0.16796  0.2628
  Y90       -0.86948E-02     0.17577       -0.049   0.96055  0.2287
  Y91        0.20684E-01     0.17578        0.118   0.90633  0.2696
  NEWENG    -0.12905E-01     0.28018       -0.046   0.96326  0.6997E-01
  MIDAT      0.34640         0.22352        1.550   0.12119  0.1365
  ENC       -0.80398E-01     0.21131       -0.380   0.70359  0.1826
  WNC       -0.19583         0.28058       -0.698   0.48520  0.9215E-01
  SAT       -0.41997E-01     0.24362       -0.172   0.86313  0.1485
  ESC        0.81292E-01     0.28875        0.282   0.77831  0.6655E-01
  WSC       -0.21571E-01     0.25700       -0.084   0.93311  0.9044E-01
  MOUNT      0.10933E-01     0.27912        0.039   0.96876  0.7509E-01
          Index    equation for FT
  Constant   0.57705         0.31353        1.841   0.06569
  EX        -0.62185E-02     0.46287E-02   -1.343   0.17912   18.21
  HHSIZE    -0.89191E-01     0.32413E-01   -2.752   0.00593   3.563
  LOCUNEMP   0.34458E-01     0.46522E-01    0.741   0.45889   6.196
  Y89        0.46003E-01     0.12174        0.378   0.70551  0.2654
  Y90        0.72304E-02     0.13579        0.053   0.95754  0.2259
  Y91       -0.17258E-01     0.14702       -0.117   0.90655  0.2675
          Disturbance correlation and std. dev.
  RHO(1,2)  -0.27630         0.89759       -0.308   0.75822

ÉÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍ»
Knowledge workers
ÉÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍ»
º Bivariate Probit Sample Selection Model                               º
º Two stage   least squares regression    Weighting variable = ONE      º
º Dependent variable is SALARY    Mean =    2.38377, S.D. =      0.6443 º
º Model size: Observations =     137, Parameters =  10, Deg.Fr. =   127 º
º Residuals:  Sum of squares=    46.4549     Std.Dev. =         0.60480 º
º Fit:        R-squared = 0.11229, Adjusted R-squared =         0.04938 º
º             Note:  Not using OLS.  R-squared is not bounded in [0,1]  º
º Model test: F[  9,    127] =    1.78,    Prob value =         0.07722 º
º Diagnostic: Log-L =   -120.3120, Restricted(á=0) Log-L =    -133.6627 º
º             Amemiya Pr. Crt.=    0.392, Akaike Info. Crt.=      1.902 º
º Full sample contains         912 observations.                        º
º Selected sample contains     137 observations.                        º
º Estimated disturbance standard deviation =        0.650300            º
º Estimated correlation with selection equation A =-0.502075            º
º Estimated correlation with selection equation B = 0.148567            º
ÈÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍͼ

  Variable  Coefficient   Standard Error  z=b/s.e. P[³Z³òz]   Mean of X
  ÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄ
  Constant   0.96471          2.0751        0.465   0.64200
  SCH1       0.56873E-01     0.98109E-01    0.580   0.56212
  EX         0.63312E-01     0.29439E-01    2.151   0.03151
  EXSQ      -0.17424E-02     0.82213E-03   -2.119   0.03406
  PT        -0.46179E-01     0.13416       -0.344   0.73068
  Y89        0.15696         0.17356        0.904   0.36579
  Y90        0.27753         0.17809        1.558   0.11916
  Y91        0.18782         0.16389        1.146   0.25179
  Lambda-A  -0.32650         0.68360       -0.478   0.63292
  Lambda-B   0.96613E-01     0.85769        0.113   0.91031

ÉÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍ»
Other workers
ÉÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍ»
º Bivariate Probit Sample Selection Model                               º
º Two stage   least squares regression    Weighting variable = ONE      º
º Dependent variable is SALARY    Mean =    1.91032, S.D. =      0.6550 º
º Model size: Observations =     449, Parameters =  10, Deg.Fr. =   439 º
º Residuals:  Sum of squares=    169.140     Std.Dev. =         0.62071 º
º Fit:        R-squared = 0.10006, Adjusted R-squared =         0.08161 º
º             Note:  Not using OLS.  R-squared is not bounded in [0,1]  º
º Model test: F[  9,    439] =    5.42,    Prob value =         0.00000 º
º Diagnostic: Log-L =   -417.9248, Restricted(á=0) Log-L =    -446.6507 º
º             Amemiya Pr. Crt.=    0.394, Akaike Info. Crt.=      1.906 º
º Full sample contains         912 observations.                        º
º Selected sample contains     449 observations.                        º
º Estimated disturbance standard deviation =        0.620163            º
º Estimated correlation with selection equation A =-0.214027            º
º Estimated correlation with selection equation B =-0.333796E-01        º
ÈÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍͼ

  Variable  Coefficient   Standard Error  z=b/s.e. P[³Z³òz]   Mean of X
  ÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄ
  Constant   0.32489         0.49136        0.661   0.50848
  SCH1       0.83236E-01     0.42996E-01    1.936   0.05288
  EX         0.47434E-01     0.15744E-01    3.013   0.00259
  EXSQ      -0.96349E-03     0.38400E-03   -2.509   0.01210
  PT         0.10131         0.66322E-01    1.528   0.12664
  Y89        0.30430E-01     0.90024E-01    0.338   0.73534
  Y90       -0.12939         0.85713E-01   -1.510   0.13116
  Y91        0.35008E-01     0.83356E-01    0.420   0.67450
  Lambda-A  -0.13273         0.37278       -0.356   0.72180
  Lambda-B  -0.20701E-01     0.49684       -0.042   0.96677

ÉÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍ»
Structural model
              ÉÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍ»
              º FIML Estimates of Bivariate Probit Model    º
              º Maximum Likelihood Estimates                º
              º Dependent variable                OCCFT     º
              º Number of observations              912     º
              º Iterations completed                  4     º
              º Log likelihood function       -911.5689     º
              º Selection model based on FT                 º
              º Means for vars. 1- 11 are after selection.  º
              ÈÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍͼ

  Variable  Coefficient   Standard Error  z=b/s.e. P[³Z³òz]   Mean of X
  ÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄ
          Index    equation for OCCUPN
  Constant   0.63708          1.0271        0.620   0.53508
  EXRELW    -0.58280         0.35766       -1.629   0.10321  0.3659
  EX        -0.59180E-01     0.19578E-01   -3.023   0.00250   17.96
  NEWENG     0.76092E-01     0.24482        0.311   0.75594  0.6997E-01
  MIDAT      0.15525         0.20234        0.767   0.44292  0.1365
  ENC       -0.67087E-01     0.21550       -0.311   0.75557  0.1826
  WNC       -0.69073E-01     0.25306       -0.273   0.78489  0.9215E-01
  SAT        0.61585E-01     0.20576        0.299   0.76470  0.1485
  ESC       -0.27444E-01     0.26612       -0.103   0.91786  0.6655E-01
  WSC        0.28633E-01     0.25330        0.113   0.91000  0.9044E-01
  MOUNT     -0.23138E-02     0.26126       -0.009   0.99293  0.7509E-01
          Index    equation for FT
  Constant   0.72080         0.30880        2.334   0.01958
  EX        -0.14121E-01     0.46712E-02   -3.023   0.00250   18.21
  HHSIZE    -0.40463E-01     0.32246E-01   -1.255   0.20954   3.563
  LOCUNEMP   0.88771E-02     0.46092E-01    0.193   0.84728   6.196
  Y89        0.19741E-01     0.12153        0.162   0.87096  0.2654
  Y90        0.31753E-02     0.13536        0.023   0.98128  0.2259
  Y91       -0.65361E-02     0.14574       -0.045   0.96423  0.2675
          Disturbance correlation and std. dev.
  RHO(1,2)  -0.51093E-02      2.3141       -0.002   0.99824

ÉÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍ»
    MGLFX=R5*R6*BBB1(2)
    MGLFX   = -0.170170E-13


Information about occupation definitions in the 3 countries:

The sources of occupation codes are as follows:

Germany: International Standard Classification of Occupations (ISCO), 1969.
Great Britain: The 1988-90 data refer to the Classification of Occupations Summary - the Key Occupations for Statistical Purposes (KOS), 1980; the 1991 data refer to the Standard Occupational Classification (SOC) of the Office for Population Censuses and Surveys.
USA: The 1988-90 data refer to the Census of Population, Classified index of industries and occupations, 1970; the 1991 data refer to the 1980 Classified index of industries and occupations.

Using the above sources, the following occupations are supposed to comprise knowledge workers:
Germany: occupation codes less than or equal to 299
Great Britain: 
1988-90:occupation codes less than or equal to 39 
1991: occupation codes less than or equal to 399
USA: 
1988-90: occupation codes less than or equal to 235 except 150-151, 153-156, 162, 165, 171, 173, 220, 225-233
1991: occupation codes less than or equal to 242 except 205, 208-216, 223-225, 235, 243.