. cd /Users/robertyaffee/Documents/data/research/chwk/phase3/Htests/H6
/Users/robertyaffee/Documents/data/research/chwk/phase3/Htests/H6

. use chwide22jul2012, replace
(Chornobyl sequelae for both sexes   women=362  men = 340 wide format)

. save chwide22jul2012, replace
file chwide22jul2012.dta saved

.
. cap gen BFbsidep1 =     max(0, BSIposymp - 50);

. cap gen BFbsidep2 = max(0, WHPsociso - 42.14);

. cap gen BFbsidep3 = max(0, 42.14 - WHPsociso);

. cap gen BFbsidep4 = max(0, BSIhos - 4);

. cap gen BFbsidep5 = max(0, BSIsoma - 6);

. cap gen BFbsidep6w3 = max(0, polprw3 - 20);

. cap gen BFbsidep7 = max(0, BSIpar - 4);

. cap gen BFbsidep11 = max(0, 65.06 - WHPsleep);

. cap gen BFbsidep12 = max(0, BSIoc - 5);

. cap gen BFbsidep13 = max(0, BSIanx - 5);

. cap gen BFbsidep14 = max(0, BSIips - 4);

. cap gen BFbsidep15 = max(0, BSIphanx - 5);

. cap gen BFbsidep16 = max(0, BSIpsyc - 4);

. cap gen BFbsidep19w2 = max(0, neiw2 - 90);

. cap gen BFbsidep21 = max(0, BSItotal - 108);

. cap gen BFbsidep23 = max(0, 108 - BSItotal);

. cap gen BFbsidep25 = CSprbslv;

. cap gen BFbsidep27  = max(0, CSprbslv - 11)* BFbsidep25;

. cap gen BFbsidep29 = max(0, 126.757 - havmil);

. cap gen BFbsidep30 = max(0, aborw2 + 1.40268e-008);

. cap gen BFps1 = max(0, bf4m - 21);

. cap gen BFps2 = max(0, 21-bf4m);

. cap gen BFps3 = max(0, WHPer - 33.7);

. cap gen BFps4 = max(0, 33.7 - WHPer);

. cap gen BFps6w3 = max(0, 20-fdferw3);

. cap gen BFps8 = max(0, 20-CSsocspt);

. cap gen BFps8w2 = max(0, 5- depagw2);

. cap gen BFdepx1 = max(0, WHPer + 3.14933e-007);

. cap gen BFdepx2 = max(0, bf4m - 9);

. cap gen BFptsd1a = max(0, BSIsoma - 6);

. cap gen BFptsd2a = max(0, fdferw2 - 7.17707e-007);

. cap gen BFptsd4a = max(0, 40- depagw3);

. cap gen BFptsd5a = max(0, airw2 - 10);

. cap gen BFptsd7a = max(0, polprw3 - 10);

. cap gen BFptsd13a = max(0, WHPer + 6.5771e-008);

. cap gen BFbsipstd1c = max(0, BSItotal-53);

. cap gen BFbsiptsd2w1c = max(0, 36- shjobw1);

. cap gen BFbsiptsdw2d = max(0, 70-fdferw2);

. cap gen BFbsiptsdd = max(0, 23.48 - WHPpa);

.
. // does it explain predict it at each wave?
. // does it explain predict it at final wave
. // how  well does it do so
. // by assumption fulfillment
. // by encompassing
. // over time
.
. des crhtw*

              storage  display     value
variable name   type   format      label      variable label
-------------


crhtw1            float    %9.0g                         Chornobyl related health threat: wave 1 alpha = .835


crhtw2            float    %9.0g                        Chornobyl related health threat in wave 2 alpha=.822


crhtw3            float    %9.0g                        Chornobyl rleated health threat in wave 3 alpha=0.833


crhtw1Xc1         float    %9.0g


. label var crhtw1 " Chornobyl related health threat: wave 1 alpha = .796"

. label var crhtw1 " Chornobyl related health threat: wave 1 alpha = .822"

. label var crhtw1 " Chornobyl related health threat: wave 1 alpha = .835"

. // zero-order analysis
. cap drop ptsd1g1pred

. cap drop ptsd1g2pred

. cap drop ptsd2g1pred

. cap drop ptsd2g2pred

. cap drop ptsd3g1pred

. cap drop ptsd3g2pred

. cap drop ptsd1g1res

. cap drop ptsd1g2res

. cap drop ptsd2g1res

. cap drop ptsd2g2res

. cap drop ptsd3g1res

. cap drop ptsd3g2res

.
. forvalues k=1/2 {
  2. forvalues i=1/3 {
  3. regress PTSDw`i' crhtw`i' if gender==`k'
  4. cap predict srptsd`i'g`k'pred, xb
  5. cap predict srptsd`i'g`k'res, residual
  6. scatter PTSDw`i' crhtw`i' if gender==`k', color(sand) ///
>   || lowess PTSDw`i' crhtw`i' if gender==`k', ///
>   title("PTSDw`i'and Chornobyl related health risk" "in wave `i' for gender==`k'") ///
>   ytitle(PTSD self-report)
  7. }
  8. }

      Source |       SS       df       MS              Number of obs =     339
-------------+------------------------------           F(  1,   337) =   70.60
       Model |  88945.9816     1  88945.9816           Prob > F      =  0.0000
    Residual |  424566.113   337   1259.8401           R-squared     =  0.1732
-------------+------------------------------           Adj R-squared =  0.1708
       Total |  513512.094   338  1519.26655           Root MSE      =  35.494

-------------+----------------------------------------------------------------
       PTSDw1 |      Coef.    Std. Err.       t    P>|t|      [95% Con f. Interval]
-------------+----------------------------------------------------------------
      crhtw1 |    17.44697    2.076417      8.40    0.000       13.3626     21.53134
       _cons |    26.78632    1.950237     13.73    0.000      22.95015     30.62249
-------------+----------------------------------------------------------------

      Source |       SS       df       MS              Number of obs =     339
-------------+------------------------------           F(  1,   337) =   51.84
       Model |  7220.64144     1  7220.64144           Prob > F      =  0.0000
    Residual |  46936.7285   337  139.278126           R-squared     =  0.1333
-------------+------------------------------           Adj R-squared =  0.1308
       Total |    54157.37   338  160.228905           Root MSE      =  11.802

-------------+----------------------------------------------------------------
       PTSDw2 |      Coef.    Std. Err.       t    P>|t|      [95% Con f. Interval]
-------------+----------------------------------------------------------------
      crhtw2 |    5.014186    .6963923      7.20    0.000      3.644362     6.384009
       _cons |    6.022352     .653659      9.21    0.000      4.736586     7.308118
-------------+----------------------------------------------------------------

      Source |       SS       df       MS              Number of obs =     339
-------------+------------------------------           F(  1,   337) =   53.02
       Model |  6616.91939     1  6616.91939           Prob > F      =  0.0000
    Residual |  42060.1163   337  124.807467           R-squared     =  0.1359
-------------+------------------------------           Adj R-squared =  0.1334
       Total |  48677.0357   338  144.014898           Root MSE      =  11.172

-------------+----------------------------------------------------------------
       PTSDw3 |      Coef.    Std. Err.       t    P>|t|      [95% Con f. Interval]
-------------+----------------------------------------------------------------
      crhtw3 |    4.804477    .6598397      7.28    0.000      3.506553       6.1024
       _cons |    6.115832    .6188266      9.88    0.000      4.898582     7.333081
-------------+----------------------------------------------------------------

      Source |       SS       df       MS              Number of obs =     362
-------------+------------------------------           F(  1,   360) =    0.59
       Model |  653.189301     1  653.189301           Prob > F      =  0.4431
    Residual |  398819.808   360   1107.8328           R-squared     =  0.0016
-------------+------------------------------           Adj R-squared = -0.0011
       Total |  399472.997   361   1106.5734           Root MSE      =  33.284

-------------+----------------------------------------------------------------
       PTSDw1 |      Coef.    Std. Err.       t    P>|t|      [95% Con f. Interval]
-------------+----------------------------------------------------------------
      crhtw1 |    1.423881    1.854348      0.77    0.443     -2.222834     5.070596
       _cons |    17.80774    1.766697     10.08    0.000      14.33339     21.28208
-------------+----------------------------------------------------------------

      Source |       SS       df       MS              Number of obs =     363
-------------+------------------------------           F(  1,   361) =    5.71
       Model |    384.3236     1    384.3236           Prob > F      =  0.0174
    Residual |  24311.2565   361  67.3442008           R-squared     =  0.0156
-------------+------------------------------           Adj R-squared =  0.0128
       Total |  24695.5801   362  68.2198345           Root MSE      =  8.2064

-------------+----------------------------------------------------------------
       PTSDw2 |      Coef.    Std. Err.       t    P>|t|      [95% Con f. Interval]
-------------+----------------------------------------------------------------
      crhtw2 |    1.182857    .4951466      2.39    0.017      .2091229     2.156591
       _cons |     3.30197    .4390455      7.52    0.000      2.438562     4.165378
-------------+----------------------------------------------------------------

      Source |       SS       df       MS              Number of obs =     363
-------------+------------------------------           F(  1,   361) =    7.63
       Model |  542.295421     1  542.295421           Prob > F      =  0.0060
    Residual |  25643.2934   361  71.0340538           R-squared     =  0.0207
-------------+------------------------------           Adj R-squared =  0.0180
       Total |  26185.5888   362  72.3358808           Root MSE      =  8.4282

-------------+----------------------------------------------------------------
       PTSDw3 |      Coef.    Std. Err.       t    P>|t|      [95% Con f. Interval]
-------------+----------------------------------------------------------------
      crhtw3 |    1.373599    .4971362      2.76    0.006      .3959526     2.351246
       _cons |    4.535616    .4505602     10.07    0.000      3.649563     5.421668
-------------+----------------------------------------------------------------

.
. // summary of findings:  self-report r^2 are very low in all but wave 1 to say
. // that there is a prediction here.
.
.
.
. // wave three retrospective analysis
.
. forvalues k=1/2 {
  2. forvalues i=1/3 {
  3. set more off
  4. regress MiPTSD crhtw`i' if gender==`k'
  5. cap predict ptsd`i'g`k'pred, xb
  6. cap predict ptsd`i'g`k'res, residual
  7. scatter MiPTSD crhtw`i' if gender==`k', color(sand) ///
>   || lowess MiPTSD crhtw`i' if gender==`k', ///
>   title("MiPTSD and Chornobyl related health risk" "in wave `i' for gender==`k'") ///
>   ytitle(MiPTSD)
  8. }
  9. }

      Source |       SS       df       MS              Number of obs =     339
-------------+------------------------------           F(  1,   337) =   62.27
       Model |  7481.97182     1  7481.97182           Prob > F      =  0.0000
    Residual |  40488.6889   337  120.144478           R-squared     =  0.1560
-------------+------------------------------           Adj R-squared =  0.1535
       Total |  47970.6608   338  141.925032           Root MSE      =  10.961

-------------+----------------------------------------------------------------
       MiPTSD |      Coef.    Std. Err.       t    P>|t|      [95% Con f. Interval]
-------------+----------------------------------------------------------------
      crhtw1 |    5.060167    .6412225      7.89    0.000      3.798865      6.32147
       _cons |    47.90203    .6022567     79.54    0.000      46.71738     49.08669
-------------+----------------------------------------------------------------

      Source |       SS       df       MS              Number of obs =     339
-------------+------------------------------           F(  1,   337) =  120.40
       Model |  12627.3854     1  12627.3854           Prob > F      =  0.0000
    Residual |  35343.2753   337  104.876188           R-squared     =  0.2632
-------------+------------------------------           Adj R-squared =  0.2610
       Total |  47970.6608   338  141.925032           Root MSE      =  10.241

-------------+----------------------------------------------------------------
       MiPTSD |      Coef.    Std. Err.       t    P>|t|      [95% Con f. Interval]
-------------+----------------------------------------------------------------
      crhtw2 |     6.63085    .6042977     10.97    0.000       5.44218     7.819521
       _cons |    48.40302    .5672157     85.33    0.000      47.28729     49.51875
-------------+----------------------------------------------------------------

      Source |       SS       df       MS              Number of obs =     339
-------------+------------------------------           F(  1,   337) =  132.86
       Model |  13564.0698     1  13564.0698           Prob > F      =  0.0000
    Residual |   34406.591   337  102.096709           R-squared     =  0.2828
-------------+------------------------------           Adj R-squared =  0.2806
       Total |  47970.6608   338  141.925032           Root MSE      =  10.104

-------------+----------------------------------------------------------------
       MiPTSD |      Coef.    Std. Err.       t    P>|t|      [95% Con f. Interval]
-------------+----------------------------------------------------------------
      crhtw3 |    6.878808    .5967934     11.53    0.000      5.704898     8.052717
       _cons |     48.4504     .559699     86.57    0.000      47.34946     49.55135
-------------+----------------------------------------------------------------

      Source |       SS       df       MS              Number of obs =     362
-------------+------------------------------           F(  1,   360) =    0.02
       Model |  3.61764538     1  3.61764538           Prob > F      =  0.8748
    Residual |  52388.6807   360  145.524113           R-squared     =  0.0001
-------------+------------------------------           Adj R-squared = -0.0027
       Total |  52392.2983   361   145.13102           Root MSE      =  12.063

-------------+----------------------------------------------------------------
       MiPTSD |      Coef.    Std. Err.       t    P>|t|      [95% Con f. Interval]
-------------+----------------------------------------------------------------
      crhtw1 |   -.1059661    .6720809     -0.16    0.875     -1.427664     1.215732
       _cons |    49.62184    .6403132     77.50    0.000      48.36261     50.88106
-------------+----------------------------------------------------------------

      Source |       SS       df       MS              Number of obs =     363
-------------+------------------------------           F(  1,   361) =   29.86
       Model |   4018.3944     1   4018.3944           Prob > F      =  0.0000
    Residual |   48586.702   361  134.589202           R-squared     =  0.0764
-------------+------------------------------           Adj R-squared =  0.0738
       Total |  52605.0964   362  145.317946           Root MSE      =  11.601

-------------+----------------------------------------------------------------
       MiPTSD |      Coef.    Std. Err.       t    P>|t|      [95% Con f. Interval]
-------------+----------------------------------------------------------------
      crhtw2 |    3.824811    .6999852      5.46    0.000       2.44825     5.201372
       _cons |    48.91023    .6206754     78.80    0.000      47.68964     50.13083
-------------+----------------------------------------------------------------

      Source |       SS       df       MS              Number of obs =     363
-------------+------------------------------           F(  1,   361) =   38.71
       Model |   5094.8758     1   5094.8758           Prob > F      =  0.0000
    Residual |  47510.2206   361  131.607259           R-squared     =  0.0969
-------------+------------------------------           Adj R-squared =  0.0943
       Total |  52605.0964   362  145.317946           Root MSE      =  11.472

-------------+----------------------------------------------------------------
       MiPTSD |      Coef.    Std. Err.       t    P>|t|      [95% Con f. Interval]
-------------+----------------------------------------------------------------
      crhtw3 |     4.21026     .676678      6.22    0.000      2.879534     5.540986
       _cons |    48.84299    .6132809     79.64    0.000      47.63694     50.04904
-------------+----------------------------------------------------------------

. // R^2 are very small except in wave 1 to say that there is a prediction in waves
. // 2 and 3
.
. // Therefore we will examine wave 1 in more detail but MiPTSD is not a valid scale there.
. cap drop r3

. regress MiPTSD crhtw1 if gender==1, vce(cluster id)    // males

Linear regression                                      Number of obs =     339
                                                       F(  1,   338) =   47.55
                                                       Prob > F      =  0.0000
                                                       R-squared     =  0.1560
                                                       Root MSE      =  10.961

(Std. Err. adjusted for 339 clusters in id)
-------------+----------------------------------------------------------------
             |               Robust
       MiPTSD |      Coef.    Std. Err.       t    P>|t|      [95% Con f. Interval]
-------------+----------------------------------------------------------------
      crhtw1 |    5.060167    .7338265      6.90    0.000      3.616725     6.503609
       _cons |    47.90203    .6420479     74.61    0.000      46.63912     49.16495
-------------+----------------------------------------------------------------

. rdiagrbst r3     // only 50% of assumptions fulfilled


*******************************************************************************
*******************************************************************************
*****                                                                     *****
*****                                                                     *****
*****                   Robust Residual Diagnostics :                     *****
*****                                                                     *****
*****                                                                     *****
*****                                            26 Jul 2012    11:21:02  *****
*******************************************************************************
*******************************************************************************


(option xb assumed; fitted values)
(1 missing value generated)
yvar =  MiPTSD
crhtw1 _cons
12
6
xvar = ,crhtw1

   | Problems
-------------+-----------
POuts   |         0
NOuts   |         0
Badouts   |         0
Collin   |         0
SW   |         0
KS   |         0
Spec   |         0
Strbrk   |         0
Coltot   |         0
Pcttot   |         0
-------------
Outlier dignostics
-------------

                        Fitted values
-------------------------------------------------------------
      Percentiles      Smallest
 1%     40.47459       40.30397
 5%     40.73051       40.30397
10%     41.19661       40.30397       Obs                 701
25%     43.73628       40.30397       Sum of Wgt.         701

50%     47.43007                      Mean           47.90203
                        Largest       Std. Dev.      4.791489
75%     53.78838       54.21492
90%     54.21492       54.21492       Variance       22.95837
95%     54.21492       54.21492       Skewness       .0578457
99%     54.21492       54.21492       Kurtosis       1.577906
(175 real changes made)
(175 real changes made)
(175 real changes made)
(175 real changes made)
(176 real changes made)
(176 real changes made)
(176 real changes made)
(176 real changes made)

    Variable |       Obs        Mean    Std. Dev.       Min        Max
-------------+--------------------------------------------------------
    __000001 |       702    .2492877    .4329092          0          1

    Variable |       Obs        Mean    Std. Dev.       Min        Max
-------------+--------------------------------------------------------
    __000003 |       702    .2492877    .4329092          0          1

    Variable |       Obs        Mean    Std. Dev.       Min        Max
-------------+--------------------------------------------------------
    __000005 |       702    .2507123    .4337322          0          1

    Variable |       Obs        Mean    Std. Dev.       Min        Max
-------------+--------------------------------------------------------
    __000007 |       702    .2507123    .4337322          0          1
.18741034
.1881236
Split sample analysis
There are 4 subsamples of equal sample size.

Variance 1 =  .18741034
Variance 2 =  .18741034
Variance 3 =  .1881236
Variance 4 =  .1881236

Minimum variance =  .18741034

Maximum variance =  .1881236

Computing Box's F test as indication of stuctural break:

for segment 1  Box's F =  1
for segment 2  Box's F =  1
for segment 3  Box's F =  1.0038059
for segment 4  Box's F =  1.0038059

There appears to be no structural break.

  +------------------+
  | obs:       r3std |
  +------------------+
  |  32.   -1.585742 |
  |  52.   -1.585742 |
  |  85.   -1.585742 |
  | 366.   -1.585742 |
  | 552.   -1.585742 |
  +------------------+

  +-----------------+
  | 679.   1.317521 |
  | 680.   1.317521 |
  | 684.   1.317521 |
  | 686.   1.317521 |
  | 690.   1.317521 |
  +-----------------+

note: 6 values of -1.585742
note: 129 values of 1.317521
There are 0 positive outliers
There are 0 negative outliers:
(702 real changes made, 702 to missing)
    0
    0
(1 real change made)
    1

                           outlier
-------------------------------------------------------------
      Percentiles      Smallest
 1%            0              0
 5%            0              0
10%            0              0       Obs                 702
25%            0              0       Sum of Wgt.         702

50%            0                      Mean           .0014245
                        Largest       Std. Dev.      .0377426
75%            0              0
90%            0              0       Variance       .0014245
95%            0              0       Skewness       26.43864
99%            0              1       Kurtosis       700.0014
(351 real changes made)

   __00000N |      Freq.     Percent        Cum.
------------+-----------------------------------
          0 |        351       50.00       50.00
          1 |        351       50.00      100.00
------------+-----------------------------------
      Total |        702      100.00

Stem-and-leaf plot for r3std

r3std rounded to nearest multiple of .01
plot in units of .01

 -15* | 9999997533333331111000000000000000000
 -14* | 8888865444444443333311111111111110
 -13* | 86666643311
 -12* | 8777776433222222111
 -11* | 9987755444444444440
 -10* | 9999999774443100000000
  -9* | 99998877776555555333222221100000
  -8* | 8776666666665433322100
  -7* | 9985433222111100
  -6* | 99988888888888887775433110000
  -5* | 999999998665555533320000
  -4* | 66666655332211111111100
  -3* | 987777665543222222222222221
  -2* | 9877777222221
  -1* | 98877755533333333333220000000
  -0* | 888888888887666665555555555555555555511
   0* | 0116666
   1* | 4556777799999
   2* | 0012333556788
   3* | 6999
   4* | 011334555566888
   5* | 00000005688
   6* | 0222888889
   7* | 0124444445777777
   8* | 27
   9* | 011156677999
  10* | 112334444446
  11* | 04446888889
  12* | 13333333333333335566666666668888888888888888
  13* | 0000012222222222222222222222222222222222222222222222222222222 ... (135)


Computing the leverage h

                            r3std
-------------------------------------------------------------
      Percentiles      Smallest
 1%    -1.550133      -1.585742
 5%     -1.49672      -1.585742
10%    -1.399444      -1.585742       Obs                 701
25%    -.8694065      -1.585742       Sum of Wgt.         701

50%    -.0985015                      Mean          -1.32e-09
                        Largest       Std. Dev.             1
75%       1.2285       1.317521
90%     1.317521       1.317521       Variance              1
95%     1.317521       1.317521       Skewness       .0578457
99%     1.317521       1.317521       Kurtosis       1.577906

standard error of snewness=  .09231881
(702 real changes made, 702 to missing)

t test for skewness =  .
badout=  1
Positive outliers are :
    0
Negative outliers are:
    0

-------------
multicollinearity tests:

    Variable |       VIF       1/VIF
-------------+----------------------
      crhtw1 |      1.00    1.000000
-------------+----------------------
    Mean VIF |      1.00

Condition number using scaled variables =         1.16

Condition Indexes and Variance-Decomposition Proportions

condition
   index  _cons crhtw1
1  1.00   0.42   0.42
2  1.16   0.58   0.58

-------------
Residual normality tests:

                   Shapiro-Wilk W test for normal data

    Variable |    Obs       W           V         z       Prob>z
-------------+--------------------------------------------------
       r3std |    701    0.94229     26.350     7.982    0.00000

                    Skewness/Kurtosis tests for Normality
                                                         ------- joint ------
    Variable |    Obs   Pr(Skewness)   Pr(Kurtosis)  adj chi2(2)    Prob>chi2
-------------+---------------------------------------------------------------
       r3std |    701      0.5279             .             .              .
Jarque-Bera normality test:  59.46 Chi(2)  1.2e-13
Jarque-Bera test for Ho: normality:
(file qnr3std.gph saved)
(file bxr3std.gph saved)
(bin=26, start=-1.5857415, width=.11166396)
(file histr3std.gph saved)
(file kdr3std.gph saved)
(file combi4.gph saved)
-------------
Specification error tests:

Ramsey RESET test using powers of the fitted values of MiPTSD
       Ho:  model has no omitted variables
                 F(3, 334) =      7.65
                  Prob > F =      0.0001
-------------

Assessment of fulfillment of cross-sectional statistical regression assumptions:
Caveat: This score is a heuristicc only as these scores are not strictly additive


Regression model assumption test result matrix


                  +--------------------+--------------+
                  |    | Problems=1   |
                  +--------------------+--------------+
                  | Positive-outliers   |            0   |
                  +--------------------+--------------+
                  | Negative-outliers   |            0   |
                  +--------------------+--------------+
                  | Bad-outliers   |            1   |
                  +--------------------+--------------+
                  | collinearity   |            0   |
                  +--------------------+--------------+
                  | Shapiro-Wilk   |            1   |
                  +--------------------+--------------+
                  | Kolmogorov-Smirnov   |            0   |
                  +--------------------+--------------+
                  | Specification   |            1   |
                  +--------------------+--------------+
                  | StructuralBreak   |            0   |
                  +--------------------+--------------+
                  | TotalProbs   |          2.5   |
                  +--------------------+--------------+
                  | Percent-Invalid   |           50   |
                  +--------------------+--------------+
                  | Pct-Pseudo-Valid   |           50   |
                  +--------------------+--------------+

Problem type total score = Bad-Outliers+collinearity + (SW + KS)/2 + Spec + strbrk
Percent_Valid = 100 - (percent of violated assumptions)

-------------
Legend:
PosOutliers = positive outlier issue
NegOutliers = negative outlier issue
Badoutliers = Number of outliers with bad leverage
collinearity =  Collinearity problem when condition number > 20
SW = Significant Shapiro-Wilk test
KS = Significant Kolmogorov-Smirnov test
Spec = Significant Ramset23 reset test
strbrk = Existence of structural break in residual variance between 4 sample segments
Column total = the total number of counted problems in the column vector
Percent invalid: Percent of 6 regression assumptions violated
Percent Pseudo-Validity =  percent of regression model assumptions not violated
-------------

.
. cap drop r4

. regress MiPTSD crhtw1 if gender==2, vce(cluster id)  // only 50% of assumptions fufilled

Linear regression                                      Number of obs =     362
                                                       F(  1,   361) =    0.02
                                                       Prob > F      =  0.8874
                                                       R-squared     =  0.0001
                                                       Root MSE      =  12.063

(Std. Err. adjusted for 362 clusters in id)
-------------+----------------------------------------------------------------
             |               Robust
       MiPTSD |      Coef.    Std. Err.       t    P>|t|      [95% Con f. Interval]
-------------+----------------------------------------------------------------
      crhtw1 |   -.1059661    .7481201     -0.14    0.887     -1.577187     1.365255
       _cons |    49.62184     .642563     77.22    0.000       48.3582     50.88547
-------------+----------------------------------------------------------------

. rdiagrbst r4


*******************************************************************************
*******************************************************************************
*****                                                                     *****
*****                                                                     *****
*****                   Robust Residual Diagnostics :                     *****
*****                                                                     *****
*****                                                                     *****
*****                                            26 Jul 2012    11:21:14  *****
*******************************************************************************
*******************************************************************************


(option xb assumed; fitted values)
(1 missing value generated)
yvar =  MiPTSD
crhtw1 _cons
12
6
xvar = ,crhtw1

   | Problems
-------------+-----------
POuts   |         0
NOuts   |         0
Badouts   |         0
Collin   |         0
SW   |         0
KS   |         0
Spec   |         0
Strbrk   |         0
Coltot   |         0
Pcttot   |         0
-------------
Outlier dignostics
-------------

                        Fitted values
-------------------------------------------------------------
      Percentiles      Smallest
 1%     49.48964       49.48964
 5%     49.48964       49.48964
10%     49.48964       49.48964       Obs                 701
25%     49.49857       49.48964       Sum of Wgt.         701

50%     49.63172                      Mean           49.62184
                        Largest       Std. Dev.      .1003394
75%     49.70907       49.78095
90%     49.76226       49.78095       Variance        .010068
95%     49.77202       49.78095       Skewness       -.057841
99%     49.77738       49.78095       Kurtosis       1.577905
(175 real changes made)
(175 real changes made)
(175 real changes made)
(175 real changes made)
(176 real changes made)
(176 real changes made)
(176 real changes made)
(176 real changes made)

    Variable |       Obs        Mean    Std. Dev.       Min        Max
-------------+--------------------------------------------------------
    __000001 |       702    .2492877    .4329092          0          1

    Variable |       Obs        Mean    Std. Dev.       Min        Max
-------------+--------------------------------------------------------
    __000003 |       702    .2492877    .4329092          0          1

    Variable |       Obs        Mean    Std. Dev.       Min        Max
-------------+--------------------------------------------------------
    __000005 |       702    .2507123    .4337322          0          1

    Variable |       Obs        Mean    Std. Dev.       Min        Max
-------------+--------------------------------------------------------
    __000007 |       702    .2507123    .4337322          0          1
.18741034
.1881236
Split sample analysis
There are 4 subsamples of equal sample size.

Variance 1 =  .18741034
Variance 2 =  .18741034
Variance 3 =  .1881236
Variance 4 =  .1881236

Minimum variance =  .18741034

Maximum variance =  .1881236

Computing Box's F test as indication of stuctural break:

for segment 1  Box's F =  1
for segment 2  Box's F =  1
for segment 3  Box's F =  1.0038059
for segment 4  Box's F =  1.0038059

There appears to be no structural break.

  +------------------+
  | obs:       r4std |
  +------------------+
  |  10.   -1.317515 |
  |  11.   -1.317515 |
  |  17.   -1.317515 |
  |  18.   -1.317515 |
  |  21.   -1.317515 |
  +------------------+

  +-----------------+
  |  52.   1.585726 |
  |  85.   1.585726 |
  | 366.   1.585726 |
  | 552.   1.585726 |
  | 649.   1.585726 |
  +-----------------+

note: 129 values of -1.317515
note: 6 values of 1.585726
There are 0 positive outliers
There are 0 negative outliers:
(702 real changes made, 702 to missing)
    0
    0
(1 real change made)
    1

                           outlier
-------------------------------------------------------------
      Percentiles      Smallest
 1%            0              0
 5%            0              0
10%            0              0       Obs                 702
25%            0              0       Sum of Wgt.         702

50%            0                      Mean           .0014245
                        Largest       Std. Dev.      .0377426
75%            0              0
90%            0              0       Variance       .0014245
95%            0              0       Skewness       26.43864
99%            0              1       Kurtosis       700.0014
(351 real changes made)

   __00000N |      Freq.     Percent        Cum.
------------+-----------------------------------
          0 |        351       50.00       50.00
          1 |        351       50.00      100.00
------------+-----------------------------------
      Total |        702      100.00

Stem-and-leaf plot for r4std

r4std rounded to nearest multiple of .01
plot in units of .01

 -13* | 2222222222222222222222222222222222222222222222222222222222222 ... (135)
 -12* | 88888888888888886666666666553333333333333331
 -11* | 98888864440
 -10* | 655555433211
  -9* | 999776651110
  -8* | 72
  -7* | 7777775444444210
  -6* | 9888882220
  -5* | 88650000000
  -4* | 888665555433110
  -3* | 9996
  -2* | 8876553332100
  -1* | 9999977776554
  -0* | 666611
   0* | 0115555555555555555555566666788888888888
   1* | 00000002233333333333555777889
   2* | 1222227777789
   3* | 122222222222222345566777789
   4* | 00111111111223355666666
   5* | 000023335555566899999999
   6* | 00001133457778888888888888999
   7* | 0011112223345899
   8* | 0012233345666666666778
   9* | 00000112222233355555567777889999
  10* | 0000000013444779999999
  11* | 0444444444445577899
  12* | 1112222223346777778
  13* | 11334666668
  14* | 0111111111111133333444444445688888
  15* | 0000000000000000001111333333357999999


Computing the leverage h

                            r4std
-------------------------------------------------------------
      Percentiles      Smallest
 1%    -1.317515      -1.317515
 5%    -1.317515      -1.317515
10%    -1.317515      -1.317515       Obs                 701
25%    -1.228515      -1.317515       Sum of Wgt.         701

50%     .0985018                      Mean           1.06e-08
                        Largest       Std. Dev.             1
75%     .8693919       1.585726
90%     1.399438       1.585726       Variance              1
95%     1.496726       1.585726       Skewness      -.0578409
99%     1.550141       1.585726       Kurtosis       1.577905

standard error of snewness=  .09231881
(702 real changes made, 702 to missing)

t test for skewness =  .
badout=  1
Positive outliers are :
    0
Negative outliers are:
    0

-------------
multicollinearity tests:

    Variable |       VIF       1/VIF
-------------+----------------------
      crhtw1 |      1.00    1.000000
-------------+----------------------
    Mean VIF |      1.00

Condition number using scaled variables =         1.15

Condition Indexes and Variance-Decomposition Proportions

condition
   index  _cons crhtw1
1  1.00   0.43   0.43
2  1.15   0.57   0.57

-------------
Residual normality tests:

                   Shapiro-Wilk W test for normal data

    Variable |    Obs       W           V         z       Prob>z
-------------+--------------------------------------------------
       r4std |    701    0.93196     31.067     8.384    0.00000

                    Skewness/Kurtosis tests for Normality
                                                         ------- joint ------
    Variable |    Obs   Pr(Skewness)   Pr(Kurtosis)  adj chi2(2)    Prob>chi2
-------------+---------------------------------------------------------------
       r4std |    701      0.5279             .             .              .
Jarque-Bera normality test:  59.46 Chi(2)  1.2e-13
Jarque-Bera test for Ho: normality:
(file qnr4std.gph saved)
(file bxr4std.gph saved)
(bin=26, start=-1.3175149, width=.11166312)
(file histr4std.gph saved)
(file kdr4std.gph saved)
(file combi4.gph saved)
-------------
Specification error tests:

Ramsey RESET test using powers of the fitted values of MiPTSD
       Ho:  model has no omitted variables
                 F(3, 357) =      5.75
                  Prob > F =      0.0008
-------------

Assessment of fulfillment of cross-sectional statistical regression assumptions:
Caveat: This score is a heuristicc only as these scores are not strictly additive


Regression model assumption test result matrix


                  +--------------------+--------------+
                  |    | Problems=1   |
                  +--------------------+--------------+
                  | Positive-outliers   |            0   |
                  +--------------------+--------------+
                  | Negative-outliers   |            0   |
                  +--------------------+--------------+
                  | Bad-outliers   |            1   |
                  +--------------------+--------------+
                  | collinearity   |            0   |
                  +--------------------+--------------+
                  | Shapiro-Wilk   |            1   |
                  +--------------------+--------------+
                  | Kolmogorov-Smirnov   |            0   |
                  +--------------------+--------------+
                  | Specification   |            1   |
                  +--------------------+--------------+
                  | StructuralBreak   |            0   |
                  +--------------------+--------------+
                  | TotalProbs   |          2.5   |
                  +--------------------+--------------+
                  | Percent-Invalid   |           50   |
                  +--------------------+--------------+
                  | Pct-Pseudo-Valid   |           50   |
                  +--------------------+--------------+

Problem type total score = Bad-Outliers+collinearity + (SW + KS)/2 + Spec + strbrk
Percent_Valid = 100 - (percent of violated assumptions)

-------------
Legend:
PosOutliers = positive outlier issue
NegOutliers = negative outlier issue
Badoutliers = Number of outliers with bad leverage
collinearity =  Collinearity problem when condition number > 20
SW = Significant Shapiro-Wilk test
KS = Significant Kolmogorov-Smirnov test
Spec = Significant Ramset23 reset test
strbrk = Existence of structural break in residual variance between 4 sample segments
Column total = the total number of counted problems in the column vector
Percent invalid: Percent of 6 regression assumptions violated
Percent Pseudo-Validity =  percent of regression model assumptions not violated
-------------

.
.
. fracpoly regress MiPTSD crhtw1 if gender==1, vce(cluster id) powers(-3 3)  // no significant results
...
-> gen double Icrht__1 = X^3-2.515413384 if e(sample)
-> gen double Icrht__2 = X^3*ln(X)-.7734369237 if e(sample)
   (where: X = (crhtw1+1.502110719680786))

Linear regression                                      Number of obs =     339
                                                       F(  2,   338) =   25.46
                                                       Prob > F      =  0.0000
                                                       R-squared     =  0.1924
                                                       Root MSE      =  10.738

(Std. Err. adjusted for 339 clusters in id)
-------------+----------------------------------------------------------------
             |               Robust
       MiPTSD |      Coef.    Std. Err.       t    P>|t|      [95% Con f. Interval]
-------------+----------------------------------------------------------------
    Icrht__1 |   -.4806786    .6918798     -0.69    0.488     -1.841611      .880254
    Icrht__2 |    1.093705     .689085      1.59    0.113     -.2617305      2.44914
       _cons |    43.82821     .823832     53.20    0.000      42.20773      45.4487
-------------+----------------------------------------------------------------
Deviance:  2568.46. Best powers of crhtw1 among 9 models fit: 3 3.

. fracpoly regress MiPTSD crhtw1 if gender==2, vce(cluster id) powers(-3 3)  // no significant results
...
-> gen double Icrht__1 = X-1.634755478 if e(sample)
-> gen double Icrht__2 = X^3-4.368762183 if e(sample)
   (where: X = (crhtw1+1.501666713505983))

Linear regression                                      Number of obs =     362
                                                       F(  2,   361) =    3.61
                                                       Prob > F      =  0.0282
                                                       R-squared     =  0.0186
                                                       Root MSE      =  11.968

(Std. Err. adjusted for 362 clusters in id)
-------------+----------------------------------------------------------------
             |               Robust
       MiPTSD |      Coef.    Std. Err.       t    P>|t|      [95% Con f. Interval]
-------------+----------------------------------------------------------------
    Icrht__1 |   -5.085849    2.052607     -2.48    0.014     -9.122417     -1.04928
    Icrht__2 |    .5686021    .2118178      2.68    0.008      .1520504     .9851539
       _cons |    47.19258    .8935598     52.81    0.000      45.43534     48.94981
-------------+----------------------------------------------------------------
Deviance:  2821.41. Best powers of crhtw1 among 9 models fit: 1 3.

.
.
.
. title "self-report analysis"


*******************************************************************************
*******************************************************************************
*****                                                                     *****
*****                                                                     *****
*****                        self-report analysis                         *****
*****                                                                     *****
*****                                                                     *****
*****                                            26 Jul 2012    11:21:27  *****
*******************************************************************************
*******************************************************************************


.
. fracpoly regress PTSDw1 crhtw1 if gender==1, vce(cluster id) powers(-3 3) // no tranformation nec.
...
-> gen double Icrht__1 = X-1.35999232 if e(sample)
-> gen double Icrht__2 = X*ln(X)-.4181691498 if e(sample)
   (where: X = (crhtw1+1.502110719680786))

Linear regression                                      Number of obs =     339
                                                       F(  2,   338) =   54.31
                                                       Prob > F      =  0.0000
                                                       R-squared     =  0.1736
                                                       Root MSE      =  35.538

(Std. Err. adjusted for 339 clusters in id)
-------------+----------------------------------------------------------------
             |               Robust
       PTSDw1 |      Coef.    Std. Err.       t    P>|t|      [95% Con f. Interval]
-------------+----------------------------------------------------------------
    Icrht__1 |    20.44299    5.252379      3.89    0.000      10.11152     30.77446
    Icrht__2 |   -2.427416    5.047644     -0.48    0.631     -12.35617     7.501335
       _cons |    25.17053    2.800766      8.99    0.000      19.66141     30.67966
-------------+----------------------------------------------------------------
Deviance:  3379.89. Best powers of crhtw1 among 9 models fit: 1 1.

.
. fracpoly regress PTSDw1 crhtw1 if gender==2, vce(cluster id) powers(-3 3) // no transformation nec.
...
-> gen double Icrht__1 = X-1.634755478 if e(sample)
-> gen double Icrht__2 = X*ln(X)-.8034712642 if e(sample)
   (where: X = (crhtw1+1.501666713505983))

Linear regression                                      Number of obs =     362
                                                       F(  2,   361) =    2.91
                                                       Prob > F      =  0.0560
                                                       R-squared     =  0.0084
                                                       Root MSE      =  33.217

(Std. Err. adjusted for 362 clusters in id)
-------------+----------------------------------------------------------------
             |               Robust
       PTSDw1 |      Coef.    Std. Err.       t    P>|t|      [95% Con f. Interval]
-------------+----------------------------------------------------------------
    Icrht__1 |    13.95545     6.17506      2.26    0.024      1.811847     26.09906
    Icrht__2 |   -9.525836    4.912623     -1.94    0.053     -19.18679     .1351189
       _cons |    21.06585    2.500396      8.43    0.000      16.14868     25.98302
-------------+----------------------------------------------------------------
Deviance:  3560.51. Best powers of crhtw1 among 9 models fit: 1 1.

.
.
. // Full male multivarate regression model
.
. local demogs age i.educ i.marrw1 i.marrw2 i.marrw3 childw1 childw2 childw3 i.emplw1 i.emplw2 ///
>   i.emplw3 inc1w1-inc4w3 occ1w1-occ8w3

. mvreg ptsdw1 ptsdw2 ptsdw3 = `demogs' airw1 radchw1 injothr BSIdep BSIanx healthef crhtw1 if gender==1
note: 3.emplw2 omitted because of collinearity
note: 5.emplw2 omitted because of collinearity
note: occ6w2 omitted because of collinearity
note: occ7w3 omitted because of collinearity
note: occ8w3 omitted because of collinearity

Equation          Obs  Parms        RMSE    "R-sq"          F        P
----------------------------------------------------------------------
ptsdw1            338     76    32.59707    0.4572   2.942886   0.0000
ptsdw2            338     76     10.2595    0.4887   3.338778   0.0000
ptsdw3            338     76    10.00996    0.4604   2.980376   0.0000

-------------+----------------------------------------------------------------
             |      Coef.    Std. Err.       t    P>|t|      [95% Con f. Interval]
-------------+----------------------------------------------------------------
ptsdw1       |
         age |    .0652943    .4077007      0.16    0.873     -.7374927     .8680813
             |
        educ |
          3  |     8.45106    12.34727      0.68    0.494     -15.86145     32.76357
          4  |    9.633239    15.23269      0.63    0.528     -20.36083     39.62731
          5  |    19.20168    13.24657      1.45    0.148     -6.881596     45.28496
          6  |   -.7028323    12.20564     -0.06    0.954     -24.73647      23.3308
          7  |   -1.306986    21.83665     -0.06    0.952     -44.30465     41.69068
          8  |   -27.90554    27.54834     -1.01    0.312     -82.14988     26.33879
             |
      marrw1 |
          1  |   -10.54586    36.21996     -0.29    0.771     -81.86512     60.77341
          2  |    8.831418    39.43265      0.22    0.823     -68.81382     86.47665
          3  |   -2.628605    37.40538     -0.07    0.944     -76.28203     71.02482
          5  |   -38.28629    43.10464     -0.89    0.375     -123.1619     46.58931
             |
      marrw2 |
          1  |    16.03065    54.40489      0.29    0.768     -91.09582     123.1571
          2  |    16.45164    56.09913      0.29    0.770     -94.01089     126.9142
          3  |    20.97518    54.58135      0.38    0.701     -86.49877     128.4491
          4  |    28.60108    58.90431      0.49    0.628     -87.38501     144.5872
          5  |    15.95701    58.27178      0.27    0.784     -98.78361     130.6976
          6  |    6.088073    58.57911      0.10    0.917     -109.2577     121.4338
             |
      marrw3 |
          2  |    4.518975    10.77546      0.42    0.675     -16.69855      25.7365
          3  |    3.616603    8.423116      0.43    0.668     -12.96902     20.20222
          5  |    -7.22613    12.27658     -0.59    0.557     -31.39944     16.94718
          6  |    9.502402    15.80969      0.60    0.548     -21.62782     40.63263
             |
     childw1 |   -3.183885    5.217167     -0.61    0.542      -13.4568     7.089028
     childw2 |   -.0382252    5.639166     -0.01    0.995     -11.14208     11.06563
     childw3 |   -2.566409    3.872762     -0.66    0.508     -10.19211     5.059291
             |
      emplw1 |
          1  |    37.68569     56.2413      0.67    0.503      -73.0568     148.4282
          2  |    57.20176    56.56735      1.01    0.313     -54.18273     168.5863
          3  |    35.09998    66.28334      0.53    0.597     -95.41587     165.6158
          4  |   -12.36829     70.4104     -0.18    0.861     -151.0106      126.274
          5  |     41.1689    56.13701      0.73    0.464     -69.36821      151.706
             |
      emplw2 |
          1  |    2.131723    13.04791      0.16    0.870      -23.5604     27.82384
          2  |   -6.196227     13.2892     -0.47    0.641     -32.36346     19.97101
          3  |           0   (omitted)
          4  |    9.176578    23.21594      0.40    0.693       -36.537     54.89016
          5  |           0   (omitted)
             |
      emplw3 |
          2  |   -2.730603    7.323707     -0.37    0.710     -17.15142     11.69021
          4  |   -4.464715    13.43722     -0.33    0.740     -30.92341     21.99398
          5  |   -21.36273    17.31818     -1.23    0.218     -55.46325      12.7378
             |
      inc1w1 |    .6244982    30.74675      0.02    0.984     -59.91768     61.16668
      inc2w1 |    2.291404    30.68861      0.07    0.941     -58.13631     62.71912
      inc3w1 |    3.363511    30.89121      0.11    0.913     -57.46313     64.19015
      inc4w1 |    7.436507    31.50587      0.24    0.814     -54.60043     69.47345
      inc1w2 |    28.78466    36.53848      0.79    0.432      -43.1618     100.7311
      inc2w2 |    25.86967    35.89286      0.72    0.472     -44.80551     96.54485
      inc3w2 |    33.63758    35.34181      0.95    0.342     -35.95255     103.2277
      inc4w2 |    3.728556     37.5998      0.10    0.921      -70.3077     77.76481
      inc1w3 |     .528931    29.15226      0.02    0.986      -56.8736     57.93146
      inc2w3 |   -2.745049     28.9518     -0.09    0.925     -59.75288     54.26278
      inc3w3 |   -13.75836    28.83279     -0.48    0.634     -70.53184     43.01511
      inc4w3 |    7.111589    30.65358      0.23    0.817     -53.24713     67.47031
      occ1w1 |   -36.72455    26.21504     -1.40    0.162     -88.34353     14.89443
      occ2w1 |   -21.72353    25.45922     -0.85    0.394     -71.85425     28.40719
      occ3w1 |   -33.88273    27.91677     -1.21    0.226     -88.85252     21.08706
      occ4w1 |   -36.60089    27.40753     -1.34    0.183     -90.56795     17.36617
      occ5w1 |   -16.70638    28.93287     -0.58    0.564     -73.67693     40.26417
      occ6w1 |   -51.76305    29.35242     -1.76    0.079     -109.5597     6.033621
      occ7w1 |     3.60118    32.92653      0.11    0.913     -61.23313     68.43549
      occ8w1 |   -31.11893    25.55705     -1.22    0.224     -81.44229     19.20443
      occ1w2 |     3.41252    33.39163      0.10    0.919      -62.3376     69.16264
      occ2w2 |   -9.298653    32.15786     -0.29    0.773      -72.6194     54.02209
      occ3w2 |   -4.910453     34.7961     -0.14    0.888     -73.42604     63.60514
      occ4w2 |   -7.219287    34.47309     -0.21    0.834     -75.09887     60.66029
      occ5w2 |   -5.600562     34.5399     -0.16    0.871     -73.61169     62.41056
      occ6w2 |           0   (omitted)
      occ7w2 |   -32.24874    37.98224     -0.85    0.397      -107.038     42.54056
      occ8w2 |   -3.355426    32.55715     -0.10    0.918     -67.46239     60.75154
      occ1w3 |    8.877136    14.55095      0.61    0.542     -19.77455     37.52882
      occ2w3 |    7.470262    14.55091      0.51    0.608     -21.18134     36.12186
      occ3w3 |   -2.413332    15.79755     -0.15    0.879     -33.51965     28.69298
      occ4w3 |    2.148485    15.67975      0.14    0.891     -28.72587     33.02284
      occ5w3 |   -8.157986    16.64796     -0.49    0.625     -40.93882     24.62285
      occ6w3 |    24.39979    27.52485      0.89    0.376     -29.79829     78.59787
      occ7w3 |           0   (omitted)
      occ8w3 |           0   (omitted)
       airw1 |    .1676251    .0672196      2.49    0.013      .0352657     .2999845
     radchw1 |   -.1472231    .0567279     -2.60    0.010     -.2589236    -.0355225
     injothr |    9.976798    5.152818      1.94    0.054     -.1694088       20.123
      BSIdep |    .1936098    .8516298      0.23    0.820       -1.4833      1.87052
      BSIanx |    .7232534    .9238086      0.78    0.434     -1.095781     2.542288
    healthef |    .1443113    .0763212      1.89    0.060     -.0059697     .2945924
      crhtw1 |    12.91632    2.599956      4.97    0.000      7.796855     18.03579
       _cons |   -52.87153    43.65574     -1.21    0.227     -138.8323     33.08922
-------------+----------------------------------------------------------------
ptsdw2       |
         age |    .2310586    .1283184      1.80    0.073     -.0216079     .4837252
             |
        educ |
          3  |    3.165499    3.886139      0.81    0.416      -4.48654     10.81754
          4  |   -2.721828    4.794286     -0.57    0.571     -12.16206     6.718408
          5  |    2.521853    4.169181      0.60    0.546     -5.687513     10.73122
          6  |    1.326255    3.841563      0.35    0.730     -6.238012     8.890523
          7  |   -3.602441    6.872795     -0.52    0.601     -17.13539     9.930504
          8  |   -9.410264    8.670475     -1.09    0.279     -26.48295      7.66242
             |
      marrw1 |
          1  |    -6.82498    11.39975     -0.60    0.550     -29.27177     15.62181
          2  |   -2.901938     12.4109     -0.23    0.815     -27.33974     21.53587
          3  |   -1.106922    11.77285     -0.09    0.925     -24.28836     22.07451
          5  |   -4.848846    13.56661     -0.36    0.721     -31.56232     21.86462
             |
      marrw2 |
          1  |   -.1126359    17.12322     -0.01    0.995     -33.82927       33.604
          2  |    1.584356    17.65646      0.09    0.929     -33.18226     36.35097
          3  |    .3527078    17.17876      0.02    0.984     -33.47329      34.1787
          4  |   -6.879523    18.53935     -0.37    0.711     -43.38461     29.62556
          5  |   -4.229833    18.34027     -0.23    0.818     -40.34292     31.88325
          6  |   -.3580644      18.437     -0.02    0.985     -36.66161     35.94548
             |
      marrw3 |
          2  |   -1.242089    3.391433     -0.37    0.714     -7.920023     5.435845
          3  |   -1.193206    2.651064     -0.45    0.653     -6.413309     4.026898
          5  |   -4.044118     3.86389     -1.05    0.296     -11.65235     3.564112
          6  |   -.1114104     4.97589     -0.02    0.982     -9.909235     9.686414
             |
     childw1 |   -1.822153    1.642034     -1.11    0.268     -5.055416      1.41111
     childw2 |    -2.74153    1.774853     -1.54    0.124     -6.236321     .7532612
     childw3 |    1.529954      1.2189      1.26    0.211      -.870133     3.930042
             |
      emplw1 |
          1  |   -1.154357     17.7012     -0.07    0.948     -36.00908     33.70037
          2  |    3.622412    17.80382      0.20    0.839     -31.43438      38.6792
          3  |   -6.827845     20.8618     -0.33    0.744     -47.90598     34.25029
          4  |    5.114071    22.16074      0.23    0.818     -38.52174     48.74989
          5  |   -1.352957    17.66838     -0.08    0.939     -36.14305     33.43713
             |
      emplw2 |
          1  |    2.504313    4.106658      0.61    0.543     -5.581941     10.59057
          2  |   -1.102691      4.1826     -0.26    0.792     -9.338481     7.133098
          3  |           0   (omitted)
          4  |   -6.491405     7.30691     -0.89    0.375     -20.87915     7.896337
          5  |           0   (omitted)
             |
      emplw3 |
          2  |    2.267587     2.30504      0.98    0.326     -2.271174     6.806348
          4  |     14.8827    4.229187      3.52    0.001      6.555181     23.21023
          5  |     12.0808    5.450666      2.22    0.028      1.348114     22.81349
             |
      inc1w1 |    1.942021     9.67713      0.20    0.841     -17.11283     20.99687
      inc2w1 |   -.6601056    9.658833     -0.07    0.946     -19.67893     18.35871
      inc3w1 |    1.438328    9.722599      0.15    0.883     -17.70605     20.58271
      inc4w1 |    7.591937    9.916054      0.77    0.445     -11.93337     27.11724
      inc1w2 |   -3.495585        11.5     -0.30    0.761     -26.13978     19.14861
      inc2w2 |    -.411374     11.2968     -0.04    0.971     -22.65545      21.8327
      inc3w2 |    -.530566    11.12336     -0.05    0.962     -22.43313       21.372
      inc4w2 |    -10.7932    11.83404     -0.91    0.363     -34.09512     12.50873
      inc1w3 |    3.128384    9.175285      0.34    0.733      -14.9383     21.19507
      inc2w3 |    1.022405    9.112196      0.11    0.911     -16.92005     18.96486
      inc3w3 |   -1.670245    9.074736     -0.18    0.854     -19.53894     16.19845
      inc4w3 |    5.174347    9.647806      0.54    0.592     -13.82276     24.17145
      occ1w1 |   -45.70994    8.250836     -5.54    0.000     -61.95633    -29.46355
      occ2w1 |     -51.295     8.01295     -6.40    0.000     -67.07298    -35.51703
      occ3w1 |   -46.66025    8.786434     -5.31    0.000     -63.96126    -29.35924
      occ4w1 |   -47.35592    8.626156     -5.49    0.000     -64.34134    -30.37051
      occ5w1 |   -47.50986    9.106237     -5.22    0.000     -65.44058    -29.57913
      occ6w1 |   -9.282865    9.238285     -1.00    0.316      -27.4736     8.907869
      occ7w1 |   -51.30649    10.36319     -4.95    0.000     -71.71223    -30.90075
      occ8w1 |   -45.59506    8.043742     -5.67    0.000     -61.43367    -29.75645
      occ1w2 |    45.63725    10.50957      4.34    0.000      24.94327     66.33122
      occ2w2 |     45.8119    10.12126      4.53    0.000      25.88254     65.74126
      occ3w2 |    43.42398    10.95161      3.97    0.000      21.85961     64.98836
      occ4w2 |    50.28913    10.84995      4.63    0.000      28.92493     71.65333
      occ5w2 |    46.37659    10.87097      4.27    0.000      24.97099     67.78219
      occ6w2 |           0   (omitted)
      occ7w2 |    47.16876    11.95441      3.95    0.000      23.62982      70.7077
      occ8w2 |     47.3328    10.24693      4.62    0.000      27.15599     67.50962
      occ1w3 |    13.91456    4.579718      3.04    0.003      4.896821      22.9323
      occ2w3 |    15.12116    4.579704      3.30    0.001       6.10345     24.13887
      occ3w3 |    13.14072    4.972069      2.64    0.009      3.350423     22.93102
      occ4w3 |     8.27188    4.934992      1.68    0.095     -1.445414     17.98917
      occ5w3 |    10.80197    5.239725      2.06    0.040      .4846328      21.1193
      occ6w3 |    15.60603    8.663082      1.80    0.073       -1.4521     32.66415
      occ7w3 |           0   (omitted)
      occ8w3 |           0   (omitted)
       airw1 |   -.0011224    .0211565     -0.05    0.958     -.0427808      .040536
     radchw1 |   -.0530823    .0178543     -2.97    0.003     -.0882386    -.0179261
     injothr |    .0766283    1.621781      0.05    0.962     -3.116755     3.270012
      BSIdep |    .5254682    .2680392      1.96    0.051      -.002317     1.053253
      BSIanx |    .0483793    .2907565      0.17    0.868     -.5241376     .6208962
    healthef |    .0869073    .0240211      3.62    0.000      .0396083     .1342063
      crhtw1 |     2.73004    .8183016      3.34    0.001      1.118755     4.341324
       _cons |   -22.83294    13.74006     -1.66    0.098     -49.88795     4.222068
-------------+----------------------------------------------------------------
ptsdw3       |
         age |    .2090005    .1251974      1.67    0.096     -.0375206     .4555216
             |
        educ |
          3  |    4.742954    3.791619      1.25    0.212      -2.72297     12.20888
          4  |   -.9345637    4.677678     -0.20    0.842     -10.14519     8.276064
          5  |    2.645224    4.067777      0.65    0.516     -5.364471     10.65492
          6  |    4.778665    3.748128      1.27    0.203     -2.601622     12.15895
          7  |   -2.029033    6.705633     -0.30    0.762     -15.23283     11.17476
          8  |    .8419554     8.45959      0.10    0.921     -15.81548     17.49939
             |
      marrw1 |
          1  |   -3.378602    11.12248     -0.30    0.762     -25.27944     18.52223
          2  |   -.3675445    12.10904     -0.03    0.976     -24.21097     23.47588
          3  |   -1.527753     11.4865     -0.13    0.894     -24.14536     21.08986
          5  |   -5.119576    13.23664     -0.39    0.699     -31.18331     20.94416
             |
      marrw2 |
          1  |   -24.61126    16.70674     -1.47    0.142     -57.50782     8.285311
          2  |   -26.14135    17.22701     -1.52    0.130     -60.06236     7.779665
          3  |   -22.68701    16.76093     -1.35    0.177     -55.69028     10.31626
          4  |   -22.40469    18.08843     -1.24    0.217     -58.02188     13.21251
          5  |   -22.58323    17.89419     -1.26    0.208     -57.81796     12.65151
          6  |   -22.83657    17.98857     -1.27    0.205     -58.25714     12.58399
             |
      marrw3 |
          2  |   -1.044106    3.308946     -0.32    0.753     -7.559617     5.471405
          3  |     -.19582    2.586584     -0.08    0.940     -5.288958     4.897318
          5  |   -3.778226    3.769911     -1.00    0.317     -11.20141     3.644954
          6  |   -6.258289    4.854865     -1.29    0.199     -15.81781      3.30123
             |
     childw1 |   -.2752986    1.602096     -0.17    0.864     -3.429921     2.879324
     childw2 |   -2.558599    1.731684     -1.48    0.141     -5.968389     .8511908
     childw3 |    1.150207    1.189254      0.97    0.334     -1.191505     3.491919
             |
      emplw1 |
          1  |    7.151848    17.27067      0.41    0.679     -26.85513     41.15883
          2  |    11.72014    17.37079      0.67    0.500     -22.48399     45.92426
          3  |   -2.206656    20.35439     -0.11    0.914     -42.28567     37.87236
          4  |    9.791588    21.62174      0.45    0.651      -32.7829     52.36608
          5  |     8.47269    17.23864      0.49    0.623     -25.47123     42.41661
             |
      emplw2 |
          1  |   -1.390877    4.006774     -0.35    0.729     -9.280455     6.498701
          2  |   -1.431189     4.08087     -0.35    0.726     -9.466666     6.604287
          3  |           0   (omitted)
          4  |   -11.38989    7.129189     -1.60    0.111     -25.42769     2.647909
          5  |           0   (omitted)
             |
      emplw3 |
          2  |    .5537141    2.248976      0.25    0.806     -3.874654     4.982082
          4  |    16.43099    4.126324      3.98    0.000      8.306012     24.55597
          5  |    9.282131    5.318094      1.75    0.082     -1.189513     19.75377
             |
      inc1w1 |   -7.814599     9.44176     -0.83    0.409     -26.40599     10.77679
      inc2w1 |   -8.136713    9.423909     -0.86    0.389     -26.69295     10.41953
      inc3w1 |   -7.264892    9.486123     -0.77    0.444     -25.94364     11.41385
      inc4w1 |   -1.614201    9.674873     -0.17    0.868      -20.6646      17.4362
      inc1w2 |    7.283338     11.2203      0.65    0.517      -14.8101     29.37677
      inc2w2 |    8.624572    11.02204      0.78    0.435     -13.07848     30.32762
      inc3w2 |    6.402235    10.85282      0.59    0.556     -14.96761     27.77208
      inc4w2 |    6.341181    11.54621      0.55    0.583     -16.39399     29.07635
      inc1w3 |    2.431876    8.952121      0.27    0.786     -15.19539     20.05914
      inc2w3 |   -2.052201    8.890567     -0.23    0.818     -19.55826     15.45386
      inc3w3 |   -1.839143    8.854018     -0.21    0.836     -19.27323     15.59495
      inc4w3 |   -1.448364     9.41315     -0.15    0.878     -19.98342     17.08669
      occ1w1 |   -48.62188    8.050157     -6.04    0.000     -64.47312    -32.77064
      occ2w1 |    -50.7422    7.818057     -6.49    0.000     -66.13642    -35.34798
      occ3w1 |    -45.8444    8.572728     -5.35    0.000     -62.72461    -28.96418
      occ4w1 |   -47.90807    8.416348     -5.69    0.000     -64.48036    -31.33577
      occ5w1 |   -44.38041    8.884752     -5.00    0.000     -61.87502    -26.88581
      occ6w1 |   -10.93728    9.013588     -1.21    0.226     -28.68557     6.811014
      occ7w1 |   -48.00161    10.11113     -4.75    0.000     -67.91104    -28.09219
      occ8w1 |   -46.90386      7.8481     -5.98    0.000     -62.35724    -31.45048
      occ1w2 |    44.81688    10.25396      4.37    0.000      24.62623     65.00753
      occ2w2 |    40.99928    9.875087      4.15    0.000      21.55464     60.44392
      occ3w2 |    38.42027    10.68524      3.60    0.000      17.38039     59.46015
      occ4w2 |    45.54421    10.58605      4.30    0.000      24.69964     66.38878
      occ5w2 |    43.99805    10.60657      4.15    0.000      23.11309     64.88302
      occ6w2 |           0   (omitted)
      occ7w2 |    46.07323    11.66365      3.95    0.000      23.10681     69.03965
      occ8w2 |    42.72728      9.9977      4.27    0.000      23.04121     62.41335
      occ1w3 |    14.57818    4.468329      3.26    0.001      5.779775     23.37659
      occ2w3 |    15.19153    4.468315      3.40    0.001      6.393154     23.98991
      occ3w3 |    14.52958    4.851137      3.00    0.003      4.977398     24.08175
      occ4w3 |    10.16896    4.814962      2.11    0.036      .6880154     19.64991
      occ5w3 |    12.40044    5.112283      2.43    0.016      2.334045     22.46683
      occ6w3 |     17.6468    8.452376      2.09    0.038      1.003565     34.29003
      occ7w3 |           0   (omitted)
      occ8w3 |           0   (omitted)
       airw1 |   -.0126733    .0206419     -0.61    0.540     -.0533184     .0279718
     radchw1 |   -.0505459    .0174201     -2.90    0.004     -.0848471    -.0162447
     injothr |    1.996836    1.582336      1.26    0.208     -1.118878     5.112549
      BSIdep |    .5644152    .2615199      2.16    0.032      .0494669     1.079363
      BSIanx |    .2087656    .2836846      0.74    0.462     -.3498264     .7673576
    healthef |    .0473322    .0234368      2.02    0.044      .0011837     .0934808
      crhtw1 |     2.64555    .7983986      3.31    0.001      1.073456     4.217645
       _cons |   -3.521169    13.40587     -0.26    0.793     -29.91814      22.8758
-------------+----------------------------------------------------------------

.
.
. // trimmed male multivarate regression model
. mvreg ptsdw1 ptsdw2 ptsdw3 = i.emplw3  inc4w1 airw1 radchw1 injothr BSIanx BSIdep healthef crhtw1 if gender==1

Equation          Obs  Parms        RMSE    "R-sq"          F        P
----------------------------------------------------------------------
ptsdw1            339     12    33.84412    0.2706   11.02873   0.0000
ptsdw2            339     12    11.11676    0.2538   10.11173   0.0000
ptsdw3            339     12    10.65981    0.2367   9.216002   0.0000

-------------+----------------------------------------------------------------
             |      Coef.    Std. Err.       t    P>|t|      [95% Con f. Interval]
-------------+----------------------------------------------------------------
ptsdw1       |
      emplw3 |
          2  |   -.0113905    6.406339     -0.00    0.999     -12.61423     12.59145
          4  |   -6.659418    4.819431     -1.38    0.168     -16.14042     2.821585
          5  |   -13.83359     12.2725     -1.13    0.260     -37.97662     10.30943
             |
      inc4w1 |   -1.821212    6.572141     -0.28    0.782     -14.75022      11.1078
       airw1 |    .2350042    .0587181      4.00    0.000      .1194912     .3505172
     radchw1 |   -.1901568    .0526722     -3.61    0.000     -.2937759    -.0865377
     injothr |     5.54087    4.429199      1.25    0.212     -3.172451     14.25419
      BSIanx |    1.144396    .8549514      1.34    0.182     -.5375029     2.826295
      BSIdep |    .0121899    .7858275      0.02    0.988     -1.533725     1.558105
    healthef |    .0991087     .069477      1.43    0.155     -.0375696      .235787
      crhtw1 |    12.00601    2.357078      5.09    0.000       7.36906     16.64296
       _cons |    4.089466    8.827384      0.46    0.643     -13.27616     21.45509
-------------+----------------------------------------------------------------
ptsdw2       |
      emplw3 |
          2  |    1.630785    2.104287      0.77    0.439     -2.508862     5.770433
          4  |    4.669949    1.583036      2.95    0.003      1.555729     7.784169
          5  |    3.821873    4.031143      0.95    0.344     -4.108373     11.75212
             |
      inc4w1 |    3.955105    2.158748      1.83    0.068     -.2916809     8.201891
       airw1 |     .023443    .0192871      1.22    0.225     -.0144995     .0613854
     radchw1 |   -.0764601    .0173012     -4.42    0.000     -.1104959    -.0424244
     injothr |    1.101504    1.454857      0.76    0.450     -1.760556     3.963564
      BSIanx |    .1672143    .2808254      0.60    0.552     -.3852381     .7196667
      BSIdep |     .533448    .2581204      2.07    0.040       .025662     1.041234
    healthef |    .0680506    .0228211      2.98    0.003      .0231559     .1129452
      crhtw1 |    2.296706    .7742281      2.97    0.003      .7736095     3.819802
       _cons |   -4.507289    2.899526     -1.55    0.121     -10.21137      1.19679
-------------+----------------------------------------------------------------
ptsdw3       |
      emplw3 |
          2  |    .4506435    2.017791      0.22    0.823     -3.518846     4.420133
          4  |    5.170559    1.517966      3.41    0.001      2.184348      8.15677
          5  |    2.628076    3.865445      0.68    0.497     -4.976202     10.23235
             |
      inc4w1 |    3.784561    2.070014      1.83    0.068      -.287663     7.856785
       airw1 |    .0052049    .0184943      0.28    0.779      -.031178     .0415878
     radchw1 |     -.06379      .01659     -3.85    0.000     -.0964267    -.0311533
     injothr |    2.517006    1.395056      1.80    0.072     -.2274103     5.261422
      BSIanx |    .2469762    .2692822      0.92    0.360      -.282768     .7767203
      BSIdep |    .5836647    .2475104      2.36    0.019       .096751     1.070578
    healthef |    .0338223     .021883      1.55    0.123      -.009227     .0768716
      crhtw1 |    2.289914    .7424038      3.08    0.002      .8294242     3.750405
       _cons |   -3.679208    2.780343     -1.32    0.187     -9.148823     1.790408
-------------+----------------------------------------------------------------

.
. // Full male multivariate regression model
. local demogs age i.educ i.marrw1 i.marrw2 i.marrw3 childw1 childw2 childw3 i.emplw1 i.emplw2 ///
>   i.emplw3 inc1w1-inc4w3 occ1w1-occ8w3

. mvreg ptsdw1 ptsdw2 ptsdw3 = `demogs' airw1 radchw1 injothr BSIdep BSIanx healthef crhtw1 if gender==1
note: 3.emplw2 omitted because of collinearity
note: 5.emplw2 omitted because of collinearity
note: occ6w2 omitted because of collinearity
note: occ7w3 omitted because of collinearity
note: occ8w3 omitted because of collinearity

Equation          Obs  Parms        RMSE    "R-sq"          F        P
----------------------------------------------------------------------
ptsdw1            338     76    32.59707    0.4572   2.942886   0.0000
ptsdw2            338     76     10.2595    0.4887   3.338778   0.0000
ptsdw3            338     76    10.00996    0.4604   2.980376   0.0000

-------------+----------------------------------------------------------------
             |      Coef.    Std. Err.       t    P>|t|      [95% Con f. Interval]
-------------+----------------------------------------------------------------
ptsdw1       |
         age |    .0652943    .4077007      0.16    0.873     -.7374927     .8680813
             |
        educ |
          3  |     8.45106    12.34727      0.68    0.494     -15.86145     32.76357
          4  |    9.633239    15.23269      0.63    0.528     -20.36083     39.62731
          5  |    19.20168    13.24657      1.45    0.148     -6.881596     45.28496
          6  |   -.7028323    12.20564     -0.06    0.954     -24.73647      23.3308
          7  |   -1.306986    21.83665     -0.06    0.952     -44.30465     41.69068
          8  |   -27.90554    27.54834     -1.01    0.312     -82.14988     26.33879
             |
      marrw1 |
          1  |   -10.54586    36.21996     -0.29    0.771     -81.86512     60.77341
          2  |    8.831418    39.43265      0.22    0.823     -68.81382     86.47665
          3  |   -2.628605    37.40538     -0.07    0.944     -76.28203     71.02482
          5  |   -38.28629    43.10464     -0.89    0.375     -123.1619     46.58931
             |
      marrw2 |
          1  |    16.03065    54.40489      0.29    0.768     -91.09582     123.1571
          2  |    16.45164    56.09913      0.29    0.770     -94.01089     126.9142
          3  |    20.97518    54.58135      0.38    0.701     -86.49877     128.4491
          4  |    28.60108    58.90431      0.49    0.628     -87.38501     144.5872
          5  |    15.95701    58.27178      0.27    0.784     -98.78361     130.6976
          6  |    6.088073    58.57911      0.10    0.917     -109.2577     121.4338
             |
      marrw3 |
          2  |    4.518975    10.77546      0.42    0.675     -16.69855      25.7365
          3  |    3.616603    8.423116      0.43    0.668     -12.96902     20.20222
          5  |    -7.22613    12.27658     -0.59    0.557     -31.39944     16.94718
          6  |    9.502402    15.80969      0.60    0.548     -21.62782     40.63263
             |
     childw1 |   -3.183885    5.217167     -0.61    0.542      -13.4568     7.089028
     childw2 |   -.0382252    5.639166     -0.01    0.995     -11.14208     11.06563
     childw3 |   -2.566409    3.872762     -0.66    0.508     -10.19211     5.059291
             |
      emplw1 |
          1  |    37.68569     56.2413      0.67    0.503      -73.0568     148.4282
          2  |    57.20176    56.56735      1.01    0.313     -54.18273     168.5863
          3  |    35.09998    66.28334      0.53    0.597     -95.41587     165.6158
          4  |   -12.36829     70.4104     -0.18    0.861     -151.0106      126.274
          5  |     41.1689    56.13701      0.73    0.464     -69.36821      151.706
             |
      emplw2 |
          1  |    2.131723    13.04791      0.16    0.870      -23.5604     27.82384
          2  |   -6.196227     13.2892     -0.47    0.641     -32.36346     19.97101
          3  |           0   (omitted)
          4  |    9.176578    23.21594      0.40    0.693       -36.537     54.89016
          5  |           0   (omitted)
             |
      emplw3 |
          2  |   -2.730603    7.323707     -0.37    0.710     -17.15142     11.69021
          4  |   -4.464715    13.43722     -0.33    0.740     -30.92341     21.99398
          5  |   -21.36273    17.31818     -1.23    0.218     -55.46325      12.7378
             |
      inc1w1 |    .6244982    30.74675      0.02    0.984     -59.91768     61.16668
      inc2w1 |    2.291404    30.68861      0.07    0.941     -58.13631     62.71912
      inc3w1 |    3.363511    30.89121      0.11    0.913     -57.46313     64.19015
      inc4w1 |    7.436507    31.50587      0.24    0.814     -54.60043     69.47345
      inc1w2 |    28.78466    36.53848      0.79    0.432      -43.1618     100.7311
      inc2w2 |    25.86967    35.89286      0.72    0.472     -44.80551     96.54485
      inc3w2 |    33.63758    35.34181      0.95    0.342     -35.95255     103.2277
      inc4w2 |    3.728556     37.5998      0.10    0.921      -70.3077     77.76481
      inc1w3 |     .528931    29.15226      0.02    0.986      -56.8736     57.93146
      inc2w3 |   -2.745049     28.9518     -0.09    0.925     -59.75288     54.26278
      inc3w3 |   -13.75836    28.83279     -0.48    0.634     -70.53184     43.01511
      inc4w3 |    7.111589    30.65358      0.23    0.817     -53.24713     67.47031
      occ1w1 |   -36.72455    26.21504     -1.40    0.162     -88.34353     14.89443
      occ2w1 |   -21.72353    25.45922     -0.85    0.394     -71.85425     28.40719
      occ3w1 |   -33.88273    27.91677     -1.21    0.226     -88.85252     21.08706
      occ4w1 |   -36.60089    27.40753     -1.34    0.183     -90.56795     17.36617
      occ5w1 |   -16.70638    28.93287     -0.58    0.564     -73.67693     40.26417
      occ6w1 |   -51.76305    29.35242     -1.76    0.079     -109.5597     6.033621
      occ7w1 |     3.60118    32.92653      0.11    0.913     -61.23313     68.43549
      occ8w1 |   -31.11893    25.55705     -1.22    0.224     -81.44229     19.20443
      occ1w2 |     3.41252    33.39163      0.10    0.919      -62.3376     69.16264
      occ2w2 |   -9.298653    32.15786     -0.29    0.773      -72.6194     54.02209
      occ3w2 |   -4.910453     34.7961     -0.14    0.888     -73.42604     63.60514
      occ4w2 |   -7.219287    34.47309     -0.21    0.834     -75.09887     60.66029
      occ5w2 |   -5.600562     34.5399     -0.16    0.871     -73.61169     62.41056
      occ6w2 |           0   (omitted)
      occ7w2 |   -32.24874    37.98224     -0.85    0.397      -107.038     42.54056
      occ8w2 |   -3.355426    32.55715     -0.10    0.918     -67.46239     60.75154
      occ1w3 |    8.877136    14.55095      0.61    0.542     -19.77455     37.52882
      occ2w3 |    7.470262    14.55091      0.51    0.608     -21.18134     36.12186
      occ3w3 |   -2.413332    15.79755     -0.15    0.879     -33.51965     28.69298
      occ4w3 |    2.148485    15.67975      0.14    0.891     -28.72587     33.02284
      occ5w3 |   -8.157986    16.64796     -0.49    0.625     -40.93882     24.62285
      occ6w3 |    24.39979    27.52485      0.89    0.376     -29.79829     78.59787
      occ7w3 |           0   (omitted)
      occ8w3 |           0   (omitted)
       airw1 |    .1676251    .0672196      2.49    0.013      .0352657     .2999845
     radchw1 |   -.1472231    .0567279     -2.60    0.010     -.2589236    -.0355225
     injothr |    9.976798    5.152818      1.94    0.054     -.1694088       20.123
      BSIdep |    .1936098    .8516298      0.23    0.820       -1.4833      1.87052
      BSIanx |    .7232534    .9238086      0.78    0.434     -1.095781     2.542288
    healthef |    .1443113    .0763212      1.89    0.060     -.0059697     .2945924
      crhtw1 |    12.91632    2.599956      4.97    0.000      7.796855     18.03579
       _cons |   -52.87153    43.65574     -1.21    0.227     -138.8323     33.08922
-------------+----------------------------------------------------------------
ptsdw2       |
         age |    .2310586    .1283184      1.80    0.073     -.0216079     .4837252
             |
        educ |
          3  |    3.165499    3.886139      0.81    0.416      -4.48654     10.81754
          4  |   -2.721828    4.794286     -0.57    0.571     -12.16206     6.718408
          5  |    2.521853    4.169181      0.60    0.546     -5.687513     10.73122
          6  |    1.326255    3.841563      0.35    0.730     -6.238012     8.890523
          7  |   -3.602441    6.872795     -0.52    0.601     -17.13539     9.930504
          8  |   -9.410264    8.670475     -1.09    0.279     -26.48295      7.66242
             |
      marrw1 |
          1  |    -6.82498    11.39975     -0.60    0.550     -29.27177     15.62181
          2  |   -2.901938     12.4109     -0.23    0.815     -27.33974     21.53587
          3  |   -1.106922    11.77285     -0.09    0.925     -24.28836     22.07451
          5  |   -4.848846    13.56661     -0.36    0.721     -31.56232     21.86462
             |
      marrw2 |
          1  |   -.1126359    17.12322     -0.01    0.995     -33.82927       33.604
          2  |    1.584356    17.65646      0.09    0.929     -33.18226     36.35097
          3  |    .3527078    17.17876      0.02    0.984     -33.47329      34.1787
          4  |   -6.879523    18.53935     -0.37    0.711     -43.38461     29.62556
          5  |   -4.229833    18.34027     -0.23    0.818     -40.34292     31.88325
          6  |   -.3580644      18.437     -0.02    0.985     -36.66161     35.94548
             |
      marrw3 |
          2  |   -1.242089    3.391433     -0.37    0.714     -7.920023     5.435845
          3  |   -1.193206    2.651064     -0.45    0.653     -6.413309     4.026898
          5  |   -4.044118     3.86389     -1.05    0.296     -11.65235     3.564112
          6  |   -.1114104     4.97589     -0.02    0.982     -9.909235     9.686414
             |
     childw1 |   -1.822153    1.642034     -1.11    0.268     -5.055416      1.41111
     childw2 |    -2.74153    1.774853     -1.54    0.124     -6.236321     .7532612
     childw3 |    1.529954      1.2189      1.26    0.211      -.870133     3.930042
             |
      emplw1 |
          1  |   -1.154357     17.7012     -0.07    0.948     -36.00908     33.70037
          2  |    3.622412    17.80382      0.20    0.839     -31.43438      38.6792
          3  |   -6.827845     20.8618     -0.33    0.744     -47.90598     34.25029
          4  |    5.114071    22.16074      0.23    0.818     -38.52174     48.74989
          5  |   -1.352957    17.66838     -0.08    0.939     -36.14305     33.43713
             |
      emplw2 |
          1  |    2.504313    4.106658      0.61    0.543     -5.581941     10.59057
          2  |   -1.102691      4.1826     -0.26    0.792     -9.338481     7.133098
          3  |           0   (omitted)
          4  |   -6.491405     7.30691     -0.89    0.375     -20.87915     7.896337
          5  |           0   (omitted)
             |
      emplw3 |
          2  |    2.267587     2.30504      0.98    0.326     -2.271174     6.806348
          4  |     14.8827    4.229187      3.52    0.001      6.555181     23.21023
          5  |     12.0808    5.450666      2.22    0.028      1.348114     22.81349
             |
      inc1w1 |    1.942021     9.67713      0.20    0.841     -17.11283     20.99687
      inc2w1 |   -.6601056    9.658833     -0.07    0.946     -19.67893     18.35871
      inc3w1 |    1.438328    9.722599      0.15    0.883     -17.70605     20.58271
      inc4w1 |    7.591937    9.916054      0.77    0.445     -11.93337     27.11724
      inc1w2 |   -3.495585        11.5     -0.30    0.761     -26.13978     19.14861
      inc2w2 |    -.411374     11.2968     -0.04    0.971     -22.65545      21.8327
      inc3w2 |    -.530566    11.12336     -0.05    0.962     -22.43313       21.372
      inc4w2 |    -10.7932    11.83404     -0.91    0.363     -34.09512     12.50873
      inc1w3 |    3.128384    9.175285      0.34    0.733      -14.9383     21.19507
      inc2w3 |    1.022405    9.112196      0.11    0.911     -16.92005     18.96486
      inc3w3 |   -1.670245    9.074736     -0.18    0.854     -19.53894     16.19845
      inc4w3 |    5.174347    9.647806      0.54    0.592     -13.82276     24.17145
      occ1w1 |   -45.70994    8.250836     -5.54    0.000     -61.95633    -29.46355
      occ2w1 |     -51.295     8.01295     -6.40    0.000     -67.07298    -35.51703
      occ3w1 |   -46.66025    8.786434     -5.31    0.000     -63.96126    -29.35924
      occ4w1 |   -47.35592    8.626156     -5.49    0.000     -64.34134    -30.37051
      occ5w1 |   -47.50986    9.106237     -5.22    0.000     -65.44058    -29.57913
      occ6w1 |   -9.282865    9.238285     -1.00    0.316      -27.4736     8.907869
      occ7w1 |   -51.30649    10.36319     -4.95    0.000     -71.71223    -30.90075
      occ8w1 |   -45.59506    8.043742     -5.67    0.000     -61.43367    -29.75645
      occ1w2 |    45.63725    10.50957      4.34    0.000      24.94327     66.33122
      occ2w2 |     45.8119    10.12126      4.53    0.000      25.88254     65.74126
      occ3w2 |    43.42398    10.95161      3.97    0.000      21.85961     64.98836
      occ4w2 |    50.28913    10.84995      4.63    0.000      28.92493     71.65333
      occ5w2 |    46.37659    10.87097      4.27    0.000      24.97099     67.78219
      occ6w2 |           0   (omitted)
      occ7w2 |    47.16876    11.95441      3.95    0.000      23.62982      70.7077
      occ8w2 |     47.3328    10.24693      4.62    0.000      27.15599     67.50962
      occ1w3 |    13.91456    4.579718      3.04    0.003      4.896821      22.9323
      occ2w3 |    15.12116    4.579704      3.30    0.001       6.10345     24.13887
      occ3w3 |    13.14072    4.972069      2.64    0.009      3.350423     22.93102
      occ4w3 |     8.27188    4.934992      1.68    0.095     -1.445414     17.98917
      occ5w3 |    10.80197    5.239725      2.06    0.040      .4846328      21.1193
      occ6w3 |    15.60603    8.663082      1.80    0.073       -1.4521     32.66415
      occ7w3 |           0   (omitted)
      occ8w3 |           0   (omitted)
       airw1 |   -.0011224    .0211565     -0.05    0.958     -.0427808      .040536
     radchw1 |   -.0530823    .0178543     -2.97    0.003     -.0882386    -.0179261
     injothr |    .0766283    1.621781      0.05    0.962     -3.116755     3.270012
      BSIdep |    .5254682    .2680392      1.96    0.051      -.002317     1.053253
      BSIanx |    .0483793    .2907565      0.17    0.868     -.5241376     .6208962
    healthef |    .0869073    .0240211      3.62    0.000      .0396083     .1342063
      crhtw1 |     2.73004    .8183016      3.34    0.001      1.118755     4.341324
       _cons |   -22.83294    13.74006     -1.66    0.098     -49.88795     4.222068
-------------+----------------------------------------------------------------
ptsdw3       |
         age |    .2090005    .1251974      1.67    0.096     -.0375206     .4555216
             |
        educ |
          3  |    4.742954    3.791619      1.25    0.212      -2.72297     12.20888
          4  |   -.9345637    4.677678     -0.20    0.842     -10.14519     8.276064
          5  |    2.645224    4.067777      0.65    0.516     -5.364471     10.65492
          6  |    4.778665    3.748128      1.27    0.203     -2.601622     12.15895
          7  |   -2.029033    6.705633     -0.30    0.762     -15.23283     11.17476
          8  |    .8419554     8.45959      0.10    0.921     -15.81548     17.49939
             |
      marrw1 |
          1  |   -3.378602    11.12248     -0.30    0.762     -25.27944     18.52223
          2  |   -.3675445    12.10904     -0.03    0.976     -24.21097     23.47588
          3  |   -1.527753     11.4865     -0.13    0.894     -24.14536     21.08986
          5  |   -5.119576    13.23664     -0.39    0.699     -31.18331     20.94416
             |
      marrw2 |
          1  |   -24.61126    16.70674     -1.47    0.142     -57.50782     8.285311
          2  |   -26.14135    17.22701     -1.52    0.130     -60.06236     7.779665
          3  |   -22.68701    16.76093     -1.35    0.177     -55.69028     10.31626
          4  |   -22.40469    18.08843     -1.24    0.217     -58.02188     13.21251
          5  |   -22.58323    17.89419     -1.26    0.208     -57.81796     12.65151
          6  |   -22.83657    17.98857     -1.27    0.205     -58.25714     12.58399
             |
      marrw3 |
          2  |   -1.044106    3.308946     -0.32    0.753     -7.559617     5.471405
          3  |     -.19582    2.586584     -0.08    0.940     -5.288958     4.897318
          5  |   -3.778226    3.769911     -1.00    0.317     -11.20141     3.644954
          6  |   -6.258289    4.854865     -1.29    0.199     -15.81781      3.30123
             |
     childw1 |   -.2752986    1.602096     -0.17    0.864     -3.429921     2.879324
     childw2 |   -2.558599    1.731684     -1.48    0.141     -5.968389     .8511908
     childw3 |    1.150207    1.189254      0.97    0.334     -1.191505     3.491919
             |
      emplw1 |
          1  |    7.151848    17.27067      0.41    0.679     -26.85513     41.15883
          2  |    11.72014    17.37079      0.67    0.500     -22.48399     45.92426
          3  |   -2.206656    20.35439     -0.11    0.914     -42.28567     37.87236
          4  |    9.791588    21.62174      0.45    0.651      -32.7829     52.36608
          5  |     8.47269    17.23864      0.49    0.623     -25.47123     42.41661
             |
      emplw2 |
          1  |   -1.390877    4.006774     -0.35    0.729     -9.280455     6.498701
          2  |   -1.431189     4.08087     -0.35    0.726     -9.466666     6.604287
          3  |           0   (omitted)
          4  |   -11.38989    7.129189     -1.60    0.111     -25.42769     2.647909
          5  |           0   (omitted)
             |
      emplw3 |
          2  |    .5537141    2.248976      0.25    0.806     -3.874654     4.982082
          4  |    16.43099    4.126324      3.98    0.000      8.306012     24.55597
          5  |    9.282131    5.318094      1.75    0.082     -1.189513     19.75377
             |
      inc1w1 |   -7.814599     9.44176     -0.83    0.409     -26.40599     10.77679
      inc2w1 |   -8.136713    9.423909     -0.86    0.389     -26.69295     10.41953
      inc3w1 |   -7.264892    9.486123     -0.77    0.444     -25.94364     11.41385
      inc4w1 |   -1.614201    9.674873     -0.17    0.868      -20.6646      17.4362
      inc1w2 |    7.283338     11.2203      0.65    0.517      -14.8101     29.37677
      inc2w2 |    8.624572    11.02204      0.78    0.435     -13.07848     30.32762
      inc3w2 |    6.402235    10.85282      0.59    0.556     -14.96761     27.77208
      inc4w2 |    6.341181    11.54621      0.55    0.583     -16.39399     29.07635
      inc1w3 |    2.431876    8.952121      0.27    0.786     -15.19539     20.05914
      inc2w3 |   -2.052201    8.890567     -0.23    0.818     -19.55826     15.45386
      inc3w3 |   -1.839143    8.854018     -0.21    0.836     -19.27323     15.59495
      inc4w3 |   -1.448364     9.41315     -0.15    0.878     -19.98342     17.08669
      occ1w1 |   -48.62188    8.050157     -6.04    0.000     -64.47312    -32.77064
      occ2w1 |    -50.7422    7.818057     -6.49    0.000     -66.13642    -35.34798
      occ3w1 |    -45.8444    8.572728     -5.35    0.000     -62.72461    -28.96418
      occ4w1 |   -47.90807    8.416348     -5.69    0.000     -64.48036    -31.33577
      occ5w1 |   -44.38041    8.884752     -5.00    0.000     -61.87502    -26.88581
      occ6w1 |   -10.93728    9.013588     -1.21    0.226     -28.68557     6.811014
      occ7w1 |   -48.00161    10.11113     -4.75    0.000     -67.91104    -28.09219
      occ8w1 |   -46.90386      7.8481     -5.98    0.000     -62.35724    -31.45048
      occ1w2 |    44.81688    10.25396      4.37    0.000      24.62623     65.00753
      occ2w2 |    40.99928    9.875087      4.15    0.000      21.55464     60.44392
      occ3w2 |    38.42027    10.68524      3.60    0.000      17.38039     59.46015
      occ4w2 |    45.54421    10.58605      4.30    0.000      24.69964     66.38878
      occ5w2 |    43.99805    10.60657      4.15    0.000      23.11309     64.88302
      occ6w2 |           0   (omitted)
      occ7w2 |    46.07323    11.66365      3.95    0.000      23.10681     69.03965
      occ8w2 |    42.72728      9.9977      4.27    0.000      23.04121     62.41335
      occ1w3 |    14.57818    4.468329      3.26    0.001      5.779775     23.37659
      occ2w3 |    15.19153    4.468315      3.40    0.001      6.393154     23.98991
      occ3w3 |    14.52958    4.851137      3.00    0.003      4.977398     24.08175
      occ4w3 |    10.16896    4.814962      2.11    0.036      .6880154     19.64991
      occ5w3 |    12.40044    5.112283      2.43    0.016      2.334045     22.46683
      occ6w3 |     17.6468    8.452376      2.09    0.038      1.003565     34.29003
      occ7w3 |           0   (omitted)
      occ8w3 |           0   (omitted)
       airw1 |   -.0126733    .0206419     -0.61    0.540     -.0533184     .0279718
     radchw1 |   -.0505459    .0174201     -2.90    0.004     -.0848471    -.0162447
     injothr |    1.996836    1.582336      1.26    0.208     -1.118878     5.112549
      BSIdep |    .5644152    .2615199      2.16    0.032      .0494669     1.079363
      BSIanx |    .2087656    .2836846      0.74    0.462     -.3498264     .7673576
    healthef |    .0473322    .0234368      2.02    0.044      .0011837     .0934808
      crhtw1 |     2.64555    .7983986      3.31    0.001      1.073456     4.217645
       _cons |   -3.521169    13.40587     -0.26    0.793     -29.91814      22.8758
-------------+----------------------------------------------------------------

.
.
. //  trimmed female multivariate regression model
.
. mvreg ptsdw1 ptsdw2 ptsdw3 =   inc4w1  airw1 injothr BSIanx  crhtw1 if gender==2

Equation          Obs  Parms        RMSE    "R-sq"          F        P
----------------------------------------------------------------------
ptsdw1            361      6    31.81516    0.0998   7.867107   0.0000
ptsdw2            361      6    7.921598    0.0970   7.630115   0.0000
ptsdw3            361      6    7.974865    0.1370    11.2747   0.0000

-------------+----------------------------------------------------------------
             |      Coef.    Std. Err.       t    P>|t|      [95% Con f. Interval]
-------------+----------------------------------------------------------------
ptsdw1       |
      inc4w1 |    .8753648    7.747675      0.11    0.910     -14.36175     16.11248
       airw1 |    .1195544    .0520178      2.30    0.022      .0172526     .2218562
     injothr |    -2.93559    4.834838     -0.61    0.544     -12.44411     6.572935
      BSIanx |    2.813658    .4711959      5.97    0.000      1.886972     3.740345
      crhtw1 |    1.555315    1.923495      0.81    0.419     -2.227562     5.338193
       _cons |    -12.9657    6.545948     -1.98    0.048     -25.83941    -.0919906
-------------+----------------------------------------------------------------
ptsdw2       |
      inc4w1 |    4.961738    1.929079      2.57    0.011      1.167878     8.755598
       airw1 |    .0180924    .0129518      1.40    0.163     -.0073795     .0435643
     injothr |    1.365859    1.203817      1.13    0.257     -1.001651     3.733369
      BSIanx |    .5736461    .1173222      4.89    0.000      .3429122       .80438
      crhtw1 |    .7742437    .4789274      1.62    0.107     -.1676479     1.716135
       _cons |    -4.33193    1.629863     -2.66    0.008     -7.537332    -1.126529
-------------+----------------------------------------------------------------
ptsdw3       |
      inc4w1 |    6.921718    1.942051      3.56    0.000      3.102347     10.74109
       airw1 |    .0115745    .0130389      0.89    0.375     -.0140687     .0372177
     injothr |    2.573177    1.211912      2.12    0.034      .1897471     4.956606
      BSIanx |    .6921086    .1181111      5.86    0.000      .4598232      .924394
      crhtw1 |   -.2735022    .4821478     -0.57    0.571     -1.221727     .6747229
       _cons |   -4.674755    1.640823     -2.85    0.005     -7.901711    -1.447799
-------------+----------------------------------------------------------------

.
. // Full male model  sig main effect
. set showbaselevels

. regress MiPTSD age i.marrw2 i.marrw3 childw1 inc4w1 airw1 airw2 airw3 radchw1 ///
>    radchw2 radchw3 injothr BSIanx BSIdep relselfr ///
>    healthef i.sector crhtw1 if gender==1, vce(cluster id)

Linear regression                                      Number of obs =     338
                                                       F( 29,   337) =       .
                                                       Prob > F      =       .
                                                       R-squared     =  0.5472
                                                       Root MSE      =  8.4109

(Std. Err. adjusted for 338 clusters in id)
-------------+----------------------------------------------------------------
             |               Robust
       MiPTSD |      Coef.    Std. Err.       t    P>|t|      [95% Con f. Interval]
-------------+----------------------------------------------------------------
         age |    .1091791    .0801459      1.36    0.174     -.0484701     .2668283
             |
      marrw2 |
          1  |      6.2149    3.115333      1.99    0.047      .0869518     12.34285
          2  |    5.537184    3.966739      1.40    0.164     -2.265505     13.33987
          3  |    6.539504    3.251598      2.01    0.045      .1435187     12.93549
          4  |     14.8022    4.008059      3.69    0.000       6.91823     22.68616
          5  |    5.410028    4.674814      1.16    0.248     -3.785464     14.60552
          6  |    11.21927    6.874484      1.63    0.104     -2.303038     24.74157
             |
      marrw3 |
          2  |   -1.683829    2.241364     -0.75    0.453     -6.092655     2.724998
          3  |   -.0270861    1.578352     -0.02    0.986     -3.131749     3.077577
          5  |   -3.235415    2.417974     -1.34    0.182     -7.991638     1.520807
          6  |   -2.371841    3.970662     -0.60    0.551     -10.18225     5.438564
             |
     childw1 |   -.7231683    .8932318     -0.81    0.419     -2.480181     1.033844
      inc4w1 |    3.885624    1.786077      2.18    0.030      .3723589     7.398888
       airw1 |    .0476542    .0168174      2.83    0.005      .0145739     .0807346
       airw2 |   -.0155555    .0267894     -0.58    0.562     -.0682511     .0371401
       airw3 |   -.0261086    .0229065     -1.14    0.255     -.0711663     .0189492
     radchw1 |   -.0283403    .0192272     -1.47    0.141     -.0661607     .0094801
     radchw2 |   -.0133153    .0458134     -0.29    0.772     -.1034316     .0768009
     radchw3 |   -.0361871    .0406907     -0.89    0.374     -.1162269     .0438527
     injothr |    4.225543    1.055368      4.00    0.000      2.149605     6.301481
      BSIanx |    1.100249    .2369373      4.64    0.000      .6341871     1.566312
      BSIdep |     .823542    .2326279      3.54    0.000      .3659564     1.281128
    relselfr |    4.964519    15.30107      0.32    0.746     -25.13313     35.06216
    healthef |    .0500498    .0215382      2.32    0.021      .0076836     .0924161
             |
      sector |
          2  |    .0955233    4.550332      0.02    0.983     -8.855109     9.046155
          3  |    .1872184    3.309873      0.06    0.955     -6.323396     6.697833
          4  |   -2.862056    3.453724     -0.83    0.408     -9.655629     3.931517
          5  |    5.710966    3.716135      1.54    0.125     -1.598777     13.02071
          6  |    4.908595    4.050679      1.21    0.226     -3.059205      12.8764
             |
      crhtw1 |    2.711253    .7578778      3.58    0.000      1.220486      4.20202
       _cons |    16.36812    6.992514      2.34    0.020      2.613647     30.12259
-------------+----------------------------------------------------------------

.
. // Trimmed male model
. set showbaselevels

. regress MiPTSD  marrw11 i.marrw2 occ1w1 occ2w1 occ1w2 occ2w2-occ4w2 ///
>      occ8w2 airw1  radchw1 injothr anxagw1 anxagw2 depagw1 depagw2 PTSDw1 PTSDw2  ///
>     healthef crhtw1 if gender==1, vce(cluster id)

Linear regression                                      Number of obs =     339
                                                       F( 24,   338) =       .
                                                       Prob > F      =       .
                                                       R-squared     =  0.5102
                                                       Root MSE      =  8.6642

(Std. Err. adjusted for 339 clusters in id)
-------------+----------------------------------------------------------------
             |               Robust
       MiPTSD |      Coef.    Std. Err.       t    P>|t|      [95% Con f. Interval]
-------------+----------------------------------------------------------------
     marrw11 |   -2.188698     1.33479     -1.64    0.102     -4.814241     .4368439
             |
      marrw2 |
          1  |     6.24635    2.723444      2.29    0.022      .8893172     11.60338
          2  |    4.877014    3.270689      1.49    0.137     -1.556456     11.31048
          3  |    6.188568    2.394776      2.58    0.010      1.478026     10.89911
          4  |    10.27922    3.832278      2.68    0.008      2.741104     17.81734
          5  |    10.01666    4.662052      2.15    0.032      .8463671     19.18695
          6  |    12.77898    9.185847      1.39    0.165     -5.289643     30.84761
             |
      occ1w1 |     -1.4465    1.904388     -0.76    0.448     -5.192445     2.299445
      occ2w1 |   -1.964074    1.974584     -0.99    0.321     -5.848094     1.919947
      occ1w2 |    1.975713    2.145161      0.92    0.358     -2.243834      6.19526
      occ2w2 |    1.035533    2.171528      0.48    0.634     -3.235879     5.306944
      occ3w2 |    1.106289    2.249143      0.49    0.623     -3.317791     5.530369
      occ4w2 |    1.287651    2.034089      0.63    0.527     -2.713418     5.288719
      occ8w2 |    .1775306    2.000134      0.09    0.929     -3.756747     4.111808
       airw1 |    .0047002    .0169325      0.28    0.781     -.0286061     .0380066
     radchw1 |    -.030431    .0136216     -2.23    0.026     -.0572247    -.0036373
     injothr |    3.500812    1.014296      3.45    0.001      1.505684     5.495941
     anxagw1 |    .0266297    .0266968      1.00    0.319      -.025883     .0791424
     anxagw2 |    .0326808     .076395      0.43    0.669     -.1175886     .1829503
     depagw1 |     .040214    .0285685      1.41    0.160     -.0159806     .0964085
     depagw2 |    .0713408    .0639669      1.12    0.266     -.0544825     .1971641
      PTSDw1 |    .0399857    .0217061      1.84    0.066     -.0027103     .0826817
      PTSDw2 |    .1674753    .0638744      2.62    0.009      .0418338     .2931167
    healthef |    .0437509    .0199579      2.19    0.029      .0044936     .0830083
      crhtw1 |    1.329444    .6786417      1.96    0.051     -.0054496     2.664337
       _cons |    33.53686    3.072292     10.92    0.000      27.49364     39.58008
-------------+----------------------------------------------------------------

.
. // Full female model sig main effect
. regress MiPTSD  marrw11 i.marrw2 occ1w1 occ2w1 occ1w2 occ2w2-occ4w2 ///
>      occ8w2 airw1  radchw1 injothr anxagw1 anxagw2 depagw1 depagw2 PTSDw1 PTSDw2  ///
>     healthef crhtw1 if gender==1, vce(cluster id)

Linear regression                                      Number of obs =     339
                                                       F( 24,   338) =       .
                                                       Prob > F      =       .
                                                       R-squared     =  0.5102
                                                       Root MSE      =  8.6642

(Std. Err. adjusted for 339 clusters in id)
-------------+----------------------------------------------------------------
             |               Robust
       MiPTSD |      Coef.    Std. Err.       t    P>|t|      [95% Con f. Interval]
-------------+----------------------------------------------------------------
     marrw11 |   -2.188698     1.33479     -1.64    0.102     -4.814241     .4368439
             |
      marrw2 |
          1  |     6.24635    2.723444      2.29    0.022      .8893172     11.60338
          2  |    4.877014    3.270689      1.49    0.137     -1.556456     11.31048
          3  |    6.188568    2.394776      2.58    0.010      1.478026     10.89911
          4  |    10.27922    3.832278      2.68    0.008      2.741104     17.81734
          5  |    10.01666    4.662052      2.15    0.032      .8463671     19.18695
          6  |    12.77898    9.185847      1.39    0.165     -5.289643     30.84761
             |
      occ1w1 |     -1.4465    1.904388     -0.76    0.448     -5.192445     2.299445
      occ2w1 |   -1.964074    1.974584     -0.99    0.321     -5.848094     1.919947
      occ1w2 |    1.975713    2.145161      0.92    0.358     -2.243834      6.19526
      occ2w2 |    1.035533    2.171528      0.48    0.634     -3.235879     5.306944
      occ3w2 |    1.106289    2.249143      0.49    0.623     -3.317791     5.530369
      occ4w2 |    1.287651    2.034089      0.63    0.527     -2.713418     5.288719
      occ8w2 |    .1775306    2.000134      0.09    0.929     -3.756747     4.111808
       airw1 |    .0047002    .0169325      0.28    0.781     -.0286061     .0380066
     radchw1 |    -.030431    .0136216     -2.23    0.026     -.0572247    -.0036373
     injothr |    3.500812    1.014296      3.45    0.001      1.505684     5.495941
     anxagw1 |    .0266297    .0266968      1.00    0.319      -.025883     .0791424
     anxagw2 |    .0326808     .076395      0.43    0.669     -.1175886     .1829503
     depagw1 |     .040214    .0285685      1.41    0.160     -.0159806     .0964085
     depagw2 |    .0713408    .0639669      1.12    0.266     -.0544825     .1971641
      PTSDw1 |    .0399857    .0217061      1.84    0.066     -.0027103     .0826817
      PTSDw2 |    .1674753    .0638744      2.62    0.009      .0418338     .2931167
    healthef |    .0437509    .0199579      2.19    0.029      .0044936     .0830083
      crhtw1 |    1.329444    .6786417      1.96    0.051     -.0054496     2.664337
       _cons |    33.53686    3.072292     10.92    0.000      27.49364     39.58008
-------------+----------------------------------------------------------------

.
. local demogs age i.educ   ///
>   i.emplw3 inc1w3-inc4w3 occ1w2-occ8w3

. set more off

.
. regress MiPTSD `demogs'  injothr BSIanx BSIdep relselfr ///
>    healthef crhtw1 if gender==2, vce(cluster id)

Linear regression                                      Number of obs =     362
                                                       F( 33,   361) =       .
                                                       Prob > F      =       .
                                                       R-squared     =  0.4728
                                                       Root MSE      =  9.2327

(Std. Err. adjusted for 362 clusters in id)
-------------+----------------------------------------------------------------
             |               Robust
       MiPTSD |      Coef.    Std. Err.       t    P>|t|      [95% Con f. Interval]
-------------+----------------------------------------------------------------
         age |    .0502488    .0642135      0.78    0.434     -.0760308     .1765284
             |
        educ |
          2  |    7.824926    2.953542      2.65    0.008      2.016616     13.63324
          3  |    11.26955    3.270005      3.45    0.001      4.838895      17.7002
          4  |    14.42489    3.406223      4.23    0.000      7.726355     21.12342
          5  |    13.52343     3.00366      4.50    0.000      7.616563      19.4303
          6  |    10.62282    2.868779      3.70    0.000      4.981203     16.26444
          7  |    4.144434    4.538329      0.91    0.362     -4.780449     13.06932
             |
      emplw3 |
          2  |    2.432312    2.035357      1.20    0.233     -1.570333     6.434957
          3  |    5.854351    3.475821      1.68    0.093     -.9810491     12.68975
          4  |    4.153586    3.132757      1.33    0.186      -2.00716     10.31433
          5  |    2.572525    2.863655      0.90    0.370     -3.059017     8.204066
             |
      inc1w3 |    4.964655    2.424058      2.05    0.041      .1976072     9.731703
      inc2w3 |    4.723632    1.865166      2.53    0.012      1.055677     8.391586
      inc3w3 |    3.898666    2.021646      1.93    0.055     -.0770171     7.874349
      inc4w3 |     6.01322    2.914847      2.06    0.040      .2810062     11.74543
      occ1w2 |    7.075529    1.580779      4.48    0.000      3.966836     10.18422
      occ2w2 |     10.4231    3.512467      2.97    0.003      3.515627     17.33056
      occ3w2 |    4.726024    2.620019      1.80    0.072     -.4263938     9.878441
      occ4w2 |    9.565371    3.930337      2.43    0.015      1.836139      17.2946
      occ5w2 |    9.072122     4.02432      2.25    0.025      1.158067     16.98618
      occ6w2 |     10.8325    4.132123      2.62    0.009       2.70644     18.95855
      occ7w2 |    7.456106    2.373359      3.14    0.002       2.78876     12.12345
      occ8w2 |    5.658227    2.647741      2.14    0.033      .4512936     10.86516
      occ1w3 |   -5.830224    1.811973     -3.22    0.001     -9.393572    -2.266876
      occ2w3 |    -9.61282    3.451889     -2.78    0.006     -16.40116    -2.824483
      occ3w3 |   -4.994549    3.121217     -1.60    0.110      -11.1326     1.143502
      occ4w3 |   -7.383398    4.187963     -1.76    0.079     -15.61927     .8524694
      occ5w3 |   -4.542986     5.20887     -0.87    0.384     -14.78653     5.700554
      occ6w3 |   -6.762455    5.760536     -1.17    0.241     -18.09088     4.565968
      occ7w3 |   -7.757961    3.443716     -2.25    0.025     -14.53023    -.9856963
      occ8w3 |   -17.66271    4.956917     -3.56    0.000     -27.41078    -7.914653
     injothr |     4.64749    1.297788      3.58    0.000      2.095317     7.199663
      BSIanx |    1.307435    .2071739      6.31    0.000      .9000161     1.714855
      BSIdep |    .4902927    .1559024      3.14    0.002      .1837018     .7968836
    relselfr |    5.821054    2.651975      2.19    0.029      .6057937     11.03631
    healthef |    .0364296    .0173664      2.10    0.037      .0022777     .0705816
      crhtw1 |    -.932488    .6512383     -1.43    0.153     -2.213185     .3482093
       _cons |    7.010912    4.481811      1.56    0.119     -1.802825     15.82465
-------------+----------------------------------------------------------------

.
. // trimmed female model  no significant main effect
.
. title "potential moderators"


*******************************************************************************
*******************************************************************************
*****                                                                     *****
*****                                                                     *****
*****                        potential moderators                         *****
*****                                                                     *****
*****                                                                     *****
*****                                            26 Jul 2012    11:21:28  *****
*******************************************************************************
*******************************************************************************


.
.
. // Trimmed male model
. set showbaselevels

. regress MiPTSD  marrw11 i.marrw2 occ1w1 occ2w1 occ1w2 occ2w2-occ4w2 ///
>      occ8w2 airw1  radchw1 injothr BSIanx BSIdep  ///
>     healthef crhtw1 if gender==1, vce(cluster id)

Linear regression                                      Number of obs =     339
                                                       F( 20,   338) =       .
                                                       Prob > F      =       .
                                                       R-squared     =  0.5126
                                                       Root MSE      =   8.588

(Std. Err. adjusted for 339 clusters in id)
-------------+----------------------------------------------------------------
             |               Robust
       MiPTSD |      Coef.    Std. Err.       t    P>|t|      [95% Con f. Interval]
-------------+----------------------------------------------------------------
     marrw11 |    -3.69917    1.355988     -2.73    0.007     -6.366408    -1.031932
             |
      marrw2 |
          1  |    14.93768    2.500149      5.97    0.000      10.01987     19.85549
          2  |    14.76366    3.245953      4.55    0.000      8.378851     21.14848
          3  |    15.92636    2.005236      7.94    0.000      11.98204     19.87067
          4  |    22.32133    2.871905      7.77    0.000      16.67227     27.97039
          5  |    10.48137    3.960545      2.65    0.009      2.690947     18.27179
          6  |    19.02093    6.379827      2.98    0.003      6.471761     31.57009
             |
      occ1w1 |   -3.490328    1.669202     -2.09    0.037      -6.77366    -.2069956
      occ2w1 |   -3.333493    1.714247     -1.94    0.053     -6.705429     .0384442
      occ1w2 |     4.13909    1.858169      2.23    0.027      .4840587     7.794121
      occ2w2 |    4.207078    1.931395      2.18    0.030      .4080093     8.006147
      occ3w2 |    2.052336    2.366203      0.87    0.386     -2.602003     6.706675
      occ4w2 |    2.065244    2.086134      0.99    0.323     -2.038197     6.168685
      occ8w2 |    3.821917    2.041862      1.87    0.062     -.1944404     7.838275
       airw1 |    .0405391    .0162265      2.50    0.013      .0086214     .0724569
     radchw1 |   -.0578935    .0136334     -4.25    0.000     -.0847104    -.0310765
     injothr |     3.98173    .9593866      4.15    0.000      2.094609      5.86885
      BSIanx |    1.246141    .2418172      5.15    0.000      .7704853     1.721798
      BSIdep |    .8599627    .2421029      3.55    0.000      .3837445     1.336181
    healthef |    .0514048    .0193715      2.65    0.008      .0133009     .0895087
      crhtw1 |    2.225844    .6683531      3.33    0.001      .9111889       3.5405
       _cons |    9.717619    3.343422      2.91    0.004      3.141084     16.29415
-------------+----------------------------------------------------------------

. foreach xvar in marrw11 i.marrw2 occ1w1 occ2w1 occ1w2 occ2w2-occ4w2 ///
>      occ8w2 airw1  radchw1 injothr BSIanx BSIdep  ///
>     healthef crhtw1 {
  2.  foreach wvar in  marrw11 i.marrw2 occ1w1 occ2w1 occ1w2 occ2w2-occ4w2 ///
>      occ8w2 airw1  radchw1 injothr BSIanx BSIdep  ///
>     healthef crhtw1  {
  3.    cap gen `xvar'Xc1 = `xvar'*crhtw1
  4.    }
  5.    }

.
. cap gen radchw1Xhealthef = radchw1*healthef

. // Trimmed male moderator model
. regress MiPTSD i.marrw2 ///
>      occ8w2 airw1  radchw1 injothr BSIanx BSIdep  ///
>     healthef crhtw1 radchw1Xhealthef ///
>     injothrXc1    ///
>     if gender==1, vce(cluster id)

Linear regression                                      Number of obs =     339
                                                       F( 15,   338) =       .
                                                       Prob > F      =       .
                                                       R-squared     =  0.5546
                                                       Root MSE      =  8.1463

(Std. Err. adjusted for 339 clusters in id)
-----------------+----------------------------------------------------------------
                 |               Robust
           MiPTSD |      Coef.    Std. Err.       t    P>|t|      [95% Con f. Interval]
-----------------+----------------------------------------------------------------
          marrw2 |
              1  |    15.19538    2.045778      7.43    0.000      11.17132     19.21944
              2  |     15.3422    3.042512      5.04    0.000      9.357553     21.32684
              3  |    17.04739    1.689697     10.09    0.000      13.72374     20.37103
              4  |    25.03557    2.055594     12.18    0.000       20.9922     29.07894
              5  |    14.03628    3.409788      4.12    0.000      7.329206     20.74336
              6  |    23.21055    6.529182      3.55    0.000       10.3676      36.0535
                 |
          occ8w2 |    .5905334    1.293692      0.46    0.648     -1.954168     3.135235
           airw1 |    .0340698    .0156069      2.18    0.030      .0033709     .0647688
         radchw1 |     .117905    .0357208      3.30    0.001       .047642      .188168
         injothr |    5.570789    1.071611      5.20    0.000      3.462922     7.678656
          BSIanx |    1.204177    .2235038      5.39    0.000      .7645432     1.643811
          BSIdep |    .9183206    .2202302      4.17    0.000      .4851261     1.351515
        healthef |     .189287    .0285942      6.62    0.000      .1330419      .245532
          crhtw1 |   -1.397745     .954407     -1.46    0.144      -3.27507     .4795808
radchw1Xhealthef |   -.0024874    .0004566     -5.45    0.000     -.0033855    -.0015893
      injothrXc1 |    4.233159    1.081137      3.92    0.000      2.106555     6.359763
           _cons |    -2.82663    3.648554     -0.77    0.439     -10.00336     4.350102
-----------------+----------------------------------------------------------------

.
.
.  //trimmed female moderator model
.
.  regress MiPTSD i.educ occ1w2-occ8w3 inc1w3-inc4w3 injothr healthef ///
>     relselfr BSIdep BSIanx  ///
>      if gender==2, vce(cluster id)

Linear regression                                      Number of obs =     363
                                                       F( 28,   362) =       .
                                                       Prob > F      =       .
                                                       R-squared     =  0.4617
                                                       Root MSE      =  9.2491

(Std. Err. adjusted for 363 clusters in id)
-------------+----------------------------------------------------------------
             |               Robust
       MiPTSD |      Coef.    Std. Err.       t    P>|t|      [95% Con f. Interval]
-------------+----------------------------------------------------------------
        educ |
          2  |    5.820597    2.197737      2.65    0.008      1.498662     10.14253
          3  |    8.557859    2.217783      3.86    0.000      4.196502     12.91922
          4  |    11.49647    2.709472      4.24    0.000      6.168188     16.82475
          5  |    10.60789       2.125      4.99    0.000      6.428997     14.78679
          6  |    7.777783    1.860158      4.18    0.000      4.119711     11.43586
          7  |    1.051078    4.626026      0.23    0.820     -8.046182     10.14834
             |
      occ1w2 |    8.731455    1.389735      6.28    0.000      5.998488     11.46442
      occ2w2 |    12.16897     3.32078      3.66    0.000      5.638522     18.69941
      occ3w2 |    6.136189    2.529987      2.43    0.016      1.160872     11.11151
      occ4w2 |     11.3704    3.877347      2.93    0.004      3.745448     18.99536
      occ5w2 |    9.690954      3.6034      2.69    0.007      2.604727     16.77718
      occ6w2 |     12.0527    4.179716      2.88    0.004       3.83313     20.27228
      occ7w2 |    8.595732    2.286393      3.76    0.000      4.099451     13.09201
      occ8w2 |    6.619405     2.39341      2.77    0.006      1.912671     11.32614
      occ1w3 |   -8.428897    1.599397     -5.27    0.000     -11.57417    -5.283622
      occ2w3 |   -12.30266    3.201279     -3.84    0.000      -18.5981    -6.007221
      occ3w3 |   -7.176093    2.980536     -2.41    0.017     -13.03743    -1.314753
      occ4w3 |   -10.52378    4.019201     -2.62    0.009      -18.4277    -2.619867
      occ5w3 |    -6.18308    4.886825     -1.27    0.207     -15.79321     3.427052
      occ6w3 |   -7.000901     6.15707     -1.14    0.256     -19.10902     5.107215
      occ7w3 |   -6.645748    2.342885     -2.84    0.005     -11.25312    -2.038374
      occ8w3 |   -19.71297    4.625454     -4.26    0.000      -28.8091    -10.61683
      inc1w3 |    4.828286    2.204215      2.19    0.029      .4936128      9.16296
      inc2w3 |    4.698021    1.460079      3.22    0.001      1.826719     7.569323
      inc3w3 |    3.730748    1.556939      2.40    0.017      .6689661     6.792529
      inc4w3 |    5.574113    2.476172      2.25    0.025       .704624      10.4436
     injothr |    4.198758    1.230167      3.41    0.001      1.779587     6.617929
    healthef |    .0349093    .0164606      2.12    0.035      .0025389     .0672797
    relselfr |    5.867825    2.654574      2.21    0.028      .6475019     11.08815
      BSIdep |       .5038    .1557731      3.23    0.001      .1974662     .8101338
      BSIanx |    1.368257    .2023052      6.76    0.000      .9704159     1.766098
       _cons |     13.4466    2.623317      5.13    0.000      8.287744     18.60545
-------------+----------------------------------------------------------------

.
.
.   foreach xvar in  i.educ occ1w2-occ8w3 inc1w3-inc4w3 injothr healthef ///
>     relselfr BSIdep BSIanx {
  2.      cap gen `xvar'Xc1 = `xvar'*crhtw1
  3.      }

.
.  // Significant moderatore effects for females of air and relocate self
.    regress MiPTSD i.educ occ1w2-occ8w3 inc1w3-inc4w3 injothr healthef ///
>     relselfr BSIdep BSIanx  ///
>     airw1Xc1 relselfrXc1  crhtw1 if gender==2, vce(cluster id)

Linear regression                                      Number of obs =     361
                                                       F( 31,   360) =       .
                                                       Prob > F      =       .
                                                       R-squared     =  0.4791
                                                       Root MSE      =  9.1493

(Std. Err. adjusted for 361 clusters in id)
-------------+----------------------------------------------------------------
             |               Robust
       MiPTSD |      Coef.    Std. Err.       t    P>|t|      [95% Con f. Interval]
-------------+----------------------------------------------------------------
        educ |
          2  |    1.503984    2.305584      0.65    0.515      -3.03012     6.038089
          3  |    3.925554    2.334682      1.68    0.094     -.6657741     8.516882
          4  |    6.379501    2.951497      2.16    0.031      .5751603     12.18384
          5  |    5.682827    2.436853      2.33    0.020      .8905708     10.47508
          6  |    3.324189    2.112783      1.57    0.117     -.8307582     7.479136
          7  |    -1.06893    3.423387     -0.31    0.755     -7.801279      5.66342
             |
      occ1w2 |    9.641632    1.561659      6.17    0.000      6.570511     12.71275
      occ2w2 |    13.07145    3.585007      3.65    0.000       6.02126     20.12163
      occ3w2 |    6.974412    2.614122      2.67    0.008      1.833543     12.11528
      occ4w2 |    12.07693      4.0078      3.01    0.003      4.195286     19.95857
      occ5w2 |    9.883227    3.962389      2.49    0.013      2.090889     17.67557
      occ6w2 |    12.61973    4.149506      3.04    0.003      4.459415     20.78005
      occ7w2 |    9.675388    2.300889      4.21    0.000      5.150515     14.20026
      occ8w2 |     6.74932    2.656016      2.54    0.011      1.526064     11.97258
      occ1w3 |   -8.385778    1.581519     -5.30    0.000     -11.49595    -5.275602
      occ2w3 |   -12.58695     3.39746     -3.70    0.000     -19.26831    -5.905584
      occ3w3 |    -7.37086    2.980232     -2.47    0.014     -13.23171    -1.510009
      occ4w3 |   -11.04816    4.021294     -2.75    0.006     -18.95633    -3.139978
      occ5w3 |   -7.545585    5.726028     -1.32    0.188     -18.80625     3.715081
      occ6w3 |   -7.393468     6.45954     -1.14    0.253     -20.09664     5.309705
      occ7w3 |   -6.724345    2.323231     -2.89    0.004     -11.29315    -2.155536
      occ8w3 |   -18.40974    4.572921     -4.03    0.000     -27.40273    -9.416745
      inc1w3 |      4.0668    2.239695      1.82    0.070     -.3377298      8.47133
      inc2w3 |     3.70331    1.575915      2.35    0.019      .6041539     6.802467
      inc3w3 |    2.985094    1.638948      1.82    0.069     -.2380221      6.20821
      inc4w3 |    4.265061    2.581923      1.65    0.099      -.812485     9.342607
     injothr |    4.458633    1.226623      3.63    0.000      2.046386      6.87088
    healthef |    .0270717    .0162242      1.67    0.096     -.0048345     .0589779
    relselfr |     3.81397    2.294743      1.66    0.097     -.6988157     8.326756
      BSIdep |    .5503419     .158567      3.47    0.001      .2385079     .8621758
      BSIanx |    1.353913    .2092183      6.47    0.000      .9424698     1.765357
    airw1Xc1 |    .0435927    .0164779      2.65    0.009      .0111877     .0759978
 relselfrXc1 |    4.783147    2.075833      2.30    0.022      .7008638     8.865429
      crhtw1 |   -3.471504    1.174961     -2.95    0.003     -5.782153    -1.160854
       _cons |    17.82546    2.836786      6.28    0.000      12.24671     23.40422
-------------+----------------------------------------------------------------

.
.
. title "Mediator analysis begins here"


*******************************************************************************
*******************************************************************************
*****                                                                     *****
*****                                                                     *****
*****                   Mediator analysis begins here                     *****
*****                                                                     *****
*****                                                                     *****
*****                                            26 Jul 2012    11:21:28  *****
*******************************************************************************
*******************************************************************************


.
.  // two wave path models for males and females
. // male model 1
. title "H6: Common perceptions of Chornobyl" "related health risk and PTSD for males"


*******************************************************************************
*******************************************************************************
*****                                                                     *****
*****                                                                     *****
*****                H6: Common perceptions of Chornobyl                  *****
*****               related health risk and PTSD for males                *****
*****                                                                     *****
*****                                                                     *****
*****                                            26 Jul 2012    11:21:28  *****
*******************************************************************************
*******************************************************************************


. set more off

. sem(age avgcumdosew1 illw1 havmil icdxcnt radchw1->avgcumdosew2) ///
>    (age avgcumdosew1 airw1 radw1 radchw1 illw1 havmil  ->crhtw1)  ///
>    (avgcumdosew1 avgcumdosew2 crhtw1 radchw1 radchw2 illw2 radw2 icdxcnt->crhtw2) ///
>    (age havmil  icdxcnt airw1 radchw1 crhtw1->PTSDw1) ///
>    (crhtw2 illw2 radw2->PTSDw2)(age PTSDw1 radw2 radchw2 avgcumdosew1 icdxcnt->PTSDw2) ///
>    (avgcumdosew2 havmil age PTSDw1 PTSDw2 crhtw2 radchw2 radw2 illw1 illw2 ->MiPTSD) ///
>    if gender==1, nocapslatent  cov(e.crhtw1*e.PTSDw1)
(8 observations with missing values excluded;
 specify option 'method(mlmv)' to use all observations)

Endogenous variables


Observed:  avgcumdosew2 crhtw1 crhtw2 PTSDw1 PTSDw2 MiPTSD


Exogenous variables


Observed:  age avgcumdosew1 illw1 havmil icdxcnt radchw1 airw1 radw1 radchw2 illw2 radw2


Fitting target model:

Iteration 0:   log likelihood = -17736.211  (not concave)
Iteration 1:   log likelihood = -17734.628
Iteration 2:   log likelihood = -17732.988
Iteration 3:   log likelihood = -17732.741
Iteration 4:   log likelihood = -17732.484  (not concave)
Iteration 5:   log likelihood =  -17732.46
Iteration 6:   log likelihood = -17732.407
Iteration 7:   log likelihood = -17732.364
Iteration 8:   log likelihood = -17732.345
Iteration 9:   log likelihood = -17732.344
Iteration 10:  log likelihood = -17732.343

Structural equation model                       Number of obs      =       331
 Estimation method = ml
 Log likelihood    = -17732.343

------------------+----------------------------------------------------------------
                  |                 OIM
                  |      Coef.    Std. Err.       z    P>|z|      [95% Con f. Interval]
------------------+----------------------------------------------------------------
 Structural        |
   avgcumdosew2 <- |
              age |    .0042741    .0030971      1.38    0.168     -.0017961     .0103442
     avgcumdosew1 |    1.150662    .0219955     52.31    0.000      1.107552     1.193772
            illw1 |    .1267012    .0977114      1.30    0.195     -.0648096     .3182121
           havmil |   -.0004914    .0003429     -1.43    0.152     -.0011634     .0001806
          icdxcnt |    .0337108    .0229542      1.47    0.142     -.0112785     .0787002
          radchw1 |    .0022718    .0009887      2.30    0.022       .000334     .0042097
            _cons |    .0410652    .1602539      0.26    0.798     -.2730267      .355157
  ----------------+----------------------------------------------------------------
   crhtw1 <-       |
              age |    .0162549    .0037448      4.34    0.000      .0089152     .0235945
     avgcumdosew1 |    .0264232    .0166899      1.58    0.113     -.0062883     .0591348
            illw1 |    .0471139    .0587031      0.80    0.422     -.0679421     .1621699
           havmil |   -.0009598      .00042     -2.29    0.022     -.0017831    -.0001365
          radchw1 |   -.0007053    .0012619     -0.56    0.576     -.0031786      .001768
            airw1 |    .0098957    .0012683      7.80    0.000        .00741     .0123815
            radw1 |    .0030811    .0011504      2.68    0.007      .0008263     .0053359
            _cons |   -1.608145    .2009874     -8.00    0.000     -2.002073    -1.214217
  ----------------+----------------------------------------------------------------
   crhtw2 <-       |
     avgcumdosew2 |     .065308    .0352059      1.86    0.064     -.0036943     .1343104
           crhtw1 |    .7972845    .0253508     31.45    0.000      .7475978     .8469712
     avgcumdosew1 |   -.1103735    .0431472     -2.56    0.011     -.1949404    -.0258065
          icdxcnt |    .0598666    .0150157      3.99    0.000      .0304364     .0892969
          radchw1 |   -.0077508    .0008868     -8.74    0.000      -.009489    -.0060127
          radchw2 |    .0088279    .0009278      9.51    0.000      .0070094     .0106463
            illw2 |    .1897646    .0424434      4.47    0.000      .1065771     .2729521
            radw2 |     .002624    .0007346      3.57    0.000      .0011842     .0040638
            _cons |   -.3867569    .0597536     -6.47    0.000     -.5038717     -.269642
  ----------------+----------------------------------------------------------------
   PTSDw1 <-       |
           crhtw1 |    88.53432    33.19479      2.67    0.008      23.47373     153.5949
              age |   -1.366398     .635724     -2.15    0.032     -2.612394    -.1204018
           havmil |    .0318959    .0514347      0.62    0.535     -.0689143      .132706
          icdxcnt |    2.894117    1.166754      2.48    0.013      .6073206     5.180914
          radchw1 |   -.2361654    .1068714     -2.21    0.027     -.4456294    -.0267013
            airw1 |   -.5533222    .3594229     -1.54    0.124     -1.257778     .1511338
            _cons |    142.6942    53.01585      2.69    0.007      38.78504     246.6034
  ----------------+----------------------------------------------------------------
   PTSDw2 <-       |
           crhtw2 |    2.834522    .7464686      3.80    0.000      1.371471     4.297574
           PTSDw1 |    .0921938    .0167231      5.51    0.000       .059417     .1249706
              age |    .1234019    .0523709      2.36    0.018      .0207569     .2260469
     avgcumdosew1 |   -.6086565    .3591329     -1.69    0.090     -1.312544     .0952311
          icdxcnt |   -1.037163     .393548     -2.64    0.008     -1.808503    -.2658233
          radchw2 |   -.0685893    .0177563     -3.86    0.000     -.1033911    -.0337875
            illw2 |    3.023728    1.059671      2.85    0.004      .9468102     5.100645
            radw2 |    .0909383    .0189519      4.80    0.000      .0537933     .1280832
            _cons |   -1.499636    2.928032     -0.51    0.609     -7.238475     4.239202
  ----------------+----------------------------------------------------------------
   MiPTSD <-       |
     avgcumdosew2 |   -.5350549    .2171563     -2.46    0.014     -.9606735    -.1094364
           crhtw2 |    3.015964    .5486406      5.50    0.000      1.940648      4.09128
           PTSDw1 |    .0847532    .0127838      6.63    0.000      .0596974     .1098089
           PTSDw2 |    .1565365    .0398818      3.93    0.000      .0783697     .2347033
              age |    .0880756    .0378223      2.33    0.020      .0139451      .162206
            illw1 |    3.325304     1.15483      2.88    0.004      1.061879     5.588729
           havmil |    .0089795    .0040689      2.21    0.027      .0010047     .0169544
          radchw2 |   -.0418445    .0133722     -3.13    0.002     -.0680534    -.0156355
            illw2 |    3.172976    .7514112      4.22    0.000      1.700237     4.645715
            radw2 |    .0835874    .0143676      5.82    0.000      .0554274     .1117475
            _cons |    37.53437    2.178869     17.23    0.000      33.26387     41.80488
------------------+----------------------------------------------------------------
 Variance          |
   e.avgcumdosew2 |    .4060684    .0315646                    .3486851     .4728953
         e.crhtw1 |    .6150598    .0478353                    .5281004     .7163385
         e.crhtw2 |    .1677849    .0130423                    .1440745     .1953973
         e.PTSDw1 |    4735.537    3191.375                    1263.942     17742.36
         e.PTSDw2 |    105.3381    8.188165                    90.45234     122.6737
         e.MiPTSD |    56.61918    4.401134                    48.61808     65.93703
------------------+----------------------------------------------------------------
 Covariance        |
   e.crhtw1        |
         e.PTSDw1 |   -47.75493    20.80817     -2.30    0.022      -88.5382    -6.971659
------------------+----------------------------------------------------------------
LR test of model vs. saturated: chi2(39)  =     50.60, Prob > chi2 = 0.1010

. estat gof

---------------------+------------------------------------------------------
Fit statistic  |      Value   Description
---------------------+------------------------------------------------------
Likelihood ratio  |
chi2_ms(39)   |     50.603     model vs. saturated
p > chi2   |      0.101
chi2_bs(81)   |   2023.828     baseline vs. saturated
p > chi2   |      0.000
---------------------+------------------------------------------------------

. estat stable

Stability analysis of simultaneous equation systems

   Eigenvalue stability condition
  +--------------------------+-------------+
  |        Eigenvalue        |   Modulus   |
  +--------------------------+-------------+
  |  .00002939 + .00002941i  |   .000042   |
  |  .00002939 - .00002941i  |   .000042   |
  | -.00002939 + .00002937i  |   .000042   |
  | -.00002939 - .00002937i  |   .000042   |
  | -1.374e-08               |   1.4e-08   |
  |  1.374e-08               |   1.4e-08   |
  +--------------------------+-------------+
   stability index =  .0000416
   All the eigenvalues lie inside the unit circle.
   SEM satisfies stability condition.

. estat mindices

Modification indices

-------------+-----------------------------------------------
             |                                       Standard
             |         MI      df    P>MI         EPC         EPC
-------------+-----------------------------------------------
 Structural   |
   crhtw1 <-  |
       illw2 |       4.523       1    0.03   -.0785421   -.0501654
  -----------+-----------------------------------------------
   crhtw2 <-  |
       radw1 |      11.783       1    0.00   -.0034118   -.1378039
  -----------+-----------------------------------------------
   PTSDw1 <-  |
     radchw2 |       4.455       1    0.03   -.1556654   -.1360587
-------------+-----------------------------------------------
EPC = expected parameter change

. estat teffects


Direct effects
------------------+----------------------------------------------------------------
                  |                 OIM
                  |      Coef.    Std. Err.       z    P>|z|      [95% Con f. Interval]
------------------+----------------------------------------------------------------
 Structural        |
   avgcumdosew2 <- |
              age |    .0042741    .0030971      1.38    0.168     -.0017961     .0103442
     avgcumdosew1 |    1.150662    .0219955     52.31    0.000      1.107552     1.193772
            illw1 |    .1267012    .0977114      1.30    0.195     -.0648096     .3182121
           havmil |   -.0004914    .0003429     -1.43    0.152     -.0011634     .0001806
          icdxcnt |    .0337108    .0229542      1.47    0.142     -.0112785     .0787002
          radchw1 |    .0022718    .0009887      2.30    0.022       .000334     .0042097
  ----------------+----------------------------------------------------------------
   crhtw1 <-       |
              age |    .0162549    .0037448      4.34    0.000      .0089152     .0235945
     avgcumdosew1 |    .0264232    .0166899      1.58    0.113     -.0062883     .0591348
            illw1 |    .0471139    .0587031      0.80    0.422     -.0679421     .1621699
           havmil |   -.0009598      .00042     -2.29    0.022     -.0017831    -.0001365
          radchw1 |   -.0007053    .0012619     -0.56    0.576     -.0031786      .001768
            airw1 |    .0098957    .0012683      7.80    0.000        .00741     .0123815
            radw1 |    .0030811    .0011504      2.68    0.007      .0008263     .0053359
  ----------------+----------------------------------------------------------------
   crhtw2 <-       |
     avgcumdosew2 |     .065308    .0352059      1.86    0.064     -.0036943     .1343104
           crhtw1 |    .7972845    .0253508     31.45    0.000      .7475978     .8469712
              age |           0   (no path)
     avgcumdosew1 |   -.1103735    .0431472     -2.56    0.011     -.1949404    -.0258065
            illw1 |           0   (no path)
           havmil |           0   (no path)
          icdxcnt |    .0598666    .0150157      3.99    0.000      .0304364     .0892969
          radchw1 |   -.0077508    .0008868     -8.74    0.000      -.009489    -.0060127
            airw1 |           0   (no path)
            radw1 |           0   (no path)
          radchw2 |    .0088279    .0009278      9.51    0.000      .0070094     .0106463
            illw2 |    .1897646    .0424434      4.47    0.000      .1065771     .2729521
            radw2 |     .002624    .0007346      3.57    0.000      .0011842     .0040638
  ----------------+----------------------------------------------------------------
   PTSDw1 <-       |
           crhtw1 |    88.53432    33.19479      2.67    0.008      23.47373     153.5949
              age |   -1.366398     .635724     -2.15    0.032     -2.612394    -.1204018
     avgcumdosew1 |           0   (no path)
            illw1 |           0   (no path)
           havmil |    .0318959    .0514347      0.62    0.535     -.0689143      .132706
          icdxcnt |    2.894117    1.166754      2.48    0.013      .6073206     5.180914
          radchw1 |   -.2361654    .1068714     -2.21    0.027     -.4456294    -.0267013
            airw1 |   -.5533222    .3594229     -1.54    0.124     -1.257778     .1511338
            radw1 |           0   (no path)
  ----------------+----------------------------------------------------------------
   PTSDw2 <-       |
     avgcumdosew2 |           0   (no path)
           crhtw1 |           0   (no path)
           crhtw2 |    2.834522    .7464686      3.80    0.000      1.371471     4.297574
           PTSDw1 |    .0921938    .0167231      5.51    0.000       .059417     .1249706
              age |    .1234019    .0523709      2.36    0.018      .0207569     .2260469
     avgcumdosew1 |   -.6086565    .3591329     -1.69    0.090     -1.312544     .0952311
            illw1 |           0   (no path)
           havmil |           0   (no path)
          icdxcnt |   -1.037163     .393548     -2.64    0.008     -1.808503    -.2658233
          radchw1 |           0   (no path)
            airw1 |           0   (no path)
            radw1 |           0   (no path)
          radchw2 |   -.0685893    .0177563     -3.86    0.000     -.1033911    -.0337875
            illw2 |    3.023728    1.059671      2.85    0.004      .9468102     5.100645
            radw2 |    .0909383    .0189519      4.80    0.000      .0537933     .1280832
  ----------------+----------------------------------------------------------------
   MiPTSD <-       |
     avgcumdosew2 |   -.5350549    .2171563     -2.46    0.014     -.9606735    -.1094364
           crhtw1 |           0   (no path)
           crhtw2 |    3.015964    .5486406      5.50    0.000      1.940648      4.09128
           PTSDw1 |    .0847532    .0127838      6.63    0.000      .0596974     .1098089
           PTSDw2 |    .1565365    .0398818      3.93    0.000      .0783697     .2347033
              age |    .0880756    .0378223      2.33    0.020      .0139451      .162206
     avgcumdosew1 |           0   (no path)
            illw1 |    3.325304     1.15483      2.88    0.004      1.061879     5.588729
           havmil |    .0089795    .0040689      2.21    0.027      .0010047     .0169544
          icdxcnt |           0   (no path)
          radchw1 |           0   (no path)
            airw1 |           0   (no path)
            radw1 |           0   (no path)
          radchw2 |   -.0418445    .0133722     -3.13    0.002     -.0680534    -.0156355
            illw2 |    3.172976    .7514112      4.22    0.000      1.700237     4.645715
            radw2 |    .0835874    .0143676      5.82    0.000      .0554274     .1117475
------------------+----------------------------------------------------------------


Indirect effects
------------------+----------------------------------------------------------------
                  |                 OIM
                  |      Coef.    Std. Err.       z    P>|z|      [95% Con f. Interval]
------------------+----------------------------------------------------------------
 Structural        |
   avgcumdosew2 <- |
              age |           0   (no path)
     avgcumdosew1 |           0   (no path)
            illw1 |           0   (no path)
           havmil |           0   (no path)
          icdxcnt |           0   (no path)
          radchw1 |           0   (no path)
  ----------------+----------------------------------------------------------------
   crhtw1 <-       |
              age |           0   (no path)
     avgcumdosew1 |           0   (no path)
            illw1 |           0   (no path)
           havmil |           0   (no path)
          radchw1 |           0   (no path)
            airw1 |           0   (no path)
            radw1 |           0   (no path)
  ----------------+----------------------------------------------------------------
   crhtw2 <-       |
     avgcumdosew2 |           0   (no path)
           crhtw1 |           0   (no path)
              age |    .0132389     .003023      4.38    0.000      .0073139     .0191639
     avgcumdosew1 |    .0962143    .0426232      2.26    0.024      .0126743     .1797543
            illw1 |    .0458378    .0474532      0.97    0.334     -.0471688     .1388444
           havmil |   -.0007973    .0003369     -2.37    0.018     -.0014576    -.0001371
          icdxcnt |    .0022016     .001912      1.15    0.250     -.0015459     .0059491
          radchw1 |    -.000414    .0010116     -0.41    0.682     -.0023967     .0015687
            airw1 |    .0078897    .0010418      7.57    0.000      .0058478     .0099316
            radw1 |    .0024565    .0009205      2.67    0.008      .0006523     .0042607
          radchw2 |           0   (no path)
            illw2 |           0   (no path)
            radw2 |           0   (no path)
  ----------------+----------------------------------------------------------------
   PTSDw1 <-       |
           crhtw1 |           0   (no path)
              age |    1.439113    .6382586      2.25    0.024      .1881494     2.690077
     avgcumdosew1 |    2.339364    1.122108      2.08    0.037       .140073     4.538654
            illw1 |    4.171198    5.005803      0.83    0.405     -5.639996     13.98239
           havmil |   -.0849753    .0513746     -1.65    0.098     -.1856677     .0157172
          icdxcnt |           0   (no path)
          radchw1 |   -.0624437    .1086569     -0.57    0.566     -.2754073     .1505199
            airw1 |    .8761098      .35776      2.45    0.014      .1749131     1.577306
            radw1 |    .2727822    .0530673      5.14    0.000      .1687721     .3767922
  ----------------+----------------------------------------------------------------
   PTSDw2 <-       |
     avgcumdosew2 |    .1851171     .099792      1.86    0.064     -.0104717     .3807058
           crhtw1 |    10.42224    3.061197      3.40    0.001      4.422399     16.42207
           crhtw2 |           0   (no path)
           PTSDw1 |           0   (no path)
              age |    .0442298    .0216272      2.05    0.041      .0018414     .0866182
     avgcumdosew1 |    .1755403    .1482977      1.18    0.237     -.1151178     .4661985
            illw1 |    .5144869    .5936437      0.87    0.386     -.6490333     1.678007
           havmil |   -.0071536    .0022963     -3.12    0.002     -.0116543     -.002653
          icdxcnt |    .4427534    .1283494      3.45    0.001      .1911932     .6943135
          radchw1 |   -.0506732    .0094877     -5.34    0.000     -.0692687    -.0320777
            airw1 |    .0521225    .0092774      5.62    0.000      .0339391      .070306
            radw1 |    .0321118     .007493      4.29    0.000      .0174258     .0467979
          radchw2 |    .0250228    .0070951      3.53    0.000      .0111167      .038929
            illw2 |    .5378921    .1858477      2.89    0.004      .1736373     .9021469
            radw2 |    .0074377    .0028587      2.60    0.009      .0018347     .0130407
  ----------------+----------------------------------------------------------------
   MiPTSD <-       |
     avgcumdosew2 |    .2259443    .1218009      1.86    0.064     -.0127812     .4646697
           crhtw1 |    11.53961     3.29359      3.50    0.000      5.084291     17.99493
           crhtw2 |    .4437063    .1168496      3.80    0.000      .2146853     .6727273
           PTSDw1 |    .0144317    .0026178      5.51    0.000      .0093009     .0195625
           PTSDw2 |           0   (no path)
              age |    .0700445    .0246427      2.84    0.004      .0217457     .1183432
     avgcumdosew1 |   -.5279009    .3006099     -1.76    0.079     -1.117086     .0612838
            illw1 |    .5045114    .6558992      0.77    0.442     -.7810275      1.79005
           havmil |   -.0077602    .0024565     -3.16    0.002     -.0125749    -.0029456
          icdxcnt |    .3213972    .1512048      2.13    0.034      .0250413     .6177531
          radchw1 |   -.0590806    .0091661     -6.45    0.000     -.0770459    -.0411153
            airw1 |    .0593114    .0087387      6.79    0.000      .0421839     .0764389
            radw1 |    .0355546    .0075358      4.72    0.000      .0207846     .0503245
          radchw2 |    .0198048    .0068949      2.87    0.004      .0062911     .0333185
            illw2 |    1.129847    .2735332      4.13    0.000      .5937316     1.665962
            radw2 |    .0233133    .0055199      4.22    0.000      .0124945      .034132
------------------+----------------------------------------------------------------


Total effects
------------------+----------------------------------------------------------------
                  |                 OIM
                  |      Coef.    Std. Err.       z    P>|z|      [95% Con f. Interval]
------------------+----------------------------------------------------------------
 Structural        |
   avgcumdosew2 <- |
              age |    .0042741    .0030971      1.38    0.168     -.0017961     .0103442
     avgcumdosew1 |    1.150662    .0219955     52.31    0.000      1.107552     1.193772
            illw1 |    .1267012    .0977114      1.30    0.195     -.0648096     .3182121
           havmil |   -.0004914    .0003429     -1.43    0.152     -.0011634     .0001806
          icdxcnt |    .0337108    .0229542      1.47    0.142     -.0112785     .0787002
          radchw1 |    .0022718    .0009887      2.30    0.022       .000334     .0042097
  ----------------+----------------------------------------------------------------
   crhtw1 <-       |
              age |    .0162549    .0037448      4.34    0.000      .0089152     .0235945
     avgcumdosew1 |    .0264232    .0166899      1.58    0.113     -.0062883     .0591348
            illw1 |    .0471139    .0587031      0.80    0.422     -.0679421     .1621699
           havmil |   -.0009598      .00042     -2.29    0.022     -.0017831    -.0001365
          radchw1 |   -.0007053    .0012619     -0.56    0.576     -.0031786      .001768
            airw1 |    .0098957    .0012683      7.80    0.000        .00741     .0123815
            radw1 |    .0030811    .0011504      2.68    0.007      .0008263     .0053359
  ----------------+----------------------------------------------------------------
   crhtw2 <-       |
     avgcumdosew2 |     .065308    .0352059      1.86    0.064     -.0036943     .1343104
           crhtw1 |    .7972845    .0253508     31.45    0.000      .7475978     .8469712
              age |    .0132389     .003023      4.38    0.000      .0073139     .0191639
     avgcumdosew1 |   -.0141592     .019643     -0.72    0.471     -.0526588     .0243404
            illw1 |    .0458378    .0474532      0.97    0.334     -.0471688     .1388444
           havmil |   -.0007973    .0003369     -2.37    0.018     -.0014576    -.0001371
          icdxcnt |    .0620682    .0149948      4.14    0.000      .0326789     .0914575
          radchw1 |   -.0081648     .001342     -6.08    0.000      -.010795    -.0055346
            airw1 |    .0078897    .0010418      7.57    0.000      .0058478     .0099316
            radw1 |    .0024565    .0009205      2.67    0.008      .0006523     .0042607
          radchw2 |    .0088279    .0009278      9.51    0.000      .0070094     .0106463
            illw2 |    .1897646    .0424434      4.47    0.000      .1065771     .2729521
            radw2 |     .002624    .0007346      3.57    0.000      .0011842     .0040638
  ----------------+----------------------------------------------------------------
   PTSDw1 <-       |
           crhtw1 |    88.53432    33.19479      2.67    0.008      23.47373     153.5949
              age |    .0727153     .166997      0.44    0.663     -.2545928     .4000235
     avgcumdosew1 |    2.339364    1.122108      2.08    0.037       .140073     4.538654
            illw1 |    4.171198    5.005803      0.83    0.405     -5.639996     13.98239
           havmil |   -.0530794    .0179578     -2.96    0.003      -.088276    -.0178828
          icdxcnt |    2.894117    1.166754      2.48    0.013      .6073206     5.180914
          radchw1 |   -.2986091      .05428     -5.50    0.000     -.4049959    -.1922223
            airw1 |    .3227876    .0548931      5.88    0.000       .215199     .4303761
            radw1 |    .2727822    .0530673      5.14    0.000      .1687721     .3767922
  ----------------+----------------------------------------------------------------
   PTSDw2 <-       |
     avgcumdosew2 |    .1851171     .099792      1.86    0.064     -.0104717     .3807058
           crhtw1 |    10.42224    3.061197      3.40    0.001      4.422399     16.42207
           crhtw2 |    2.834522    .7464686      3.80    0.000      1.371471     4.297574
           PTSDw1 |    .0921938    .0167231      5.51    0.000       .059417     .1249706
              age |    .1676317    .0540941      3.10    0.002      .0616092     .2736542
     avgcumdosew1 |   -.4331162    .3845588     -1.13    0.260     -1.186838     .3206053
            illw1 |    .5144869    .5936437      0.87    0.386     -.6490333     1.678007
           havmil |   -.0071536    .0022963     -3.12    0.002     -.0116543     -.002653
          icdxcnt |   -.5944098    .4005187     -1.48    0.138     -1.379412     .1905925
          radchw1 |   -.0506732    .0094877     -5.34    0.000     -.0692687    -.0320777
            airw1 |    .0521225    .0092774      5.62    0.000      .0339391      .070306
            radw1 |    .0321118     .007493      4.29    0.000      .0174258     .0467979
          radchw2 |   -.0435665    .0178722     -2.44    0.015     -.0785954    -.0085375
            illw2 |     3.56162    1.069887      3.33    0.001       1.46468      5.65856
            radw2 |     .098376    .0188778      5.21    0.000      .0613763     .1353757
  ----------------+----------------------------------------------------------------
   MiPTSD <-       |
     avgcumdosew2 |   -.3091107    .2489826     -1.24    0.214     -.7971076     .1788862
           crhtw1 |    11.53961     3.29359      3.50    0.000      5.084291     17.99493
           crhtw2 |     3.45967     .560946      6.17    0.000      2.360237     4.559104
           PTSDw1 |    .0991849    .0130491      7.60    0.000      .0736092     .1247606
           PTSDw2 |    .1565365    .0398818      3.93    0.000      .0783697     .2347033
              age |      .15812    .0426605      3.71    0.000      .0745071      .241733
     avgcumdosew1 |   -.5279009    .3006099     -1.76    0.079     -1.117086     .0612838
            illw1 |    3.829815    1.325412      2.89    0.004      1.232056     6.427574
           havmil |    .0012193    .0046767      0.26    0.794     -.0079469     .0103855
          icdxcnt |    .3213972    .1512048      2.13    0.034      .0250413     .6177531
          radchw1 |   -.0590806    .0091661     -6.45    0.000     -.0770459    -.0411153
            airw1 |    .0593114    .0087387      6.79    0.000      .0421839     .0764389
            radw1 |    .0355546    .0075358      4.72    0.000      .0207846     .0503245
          radchw2 |   -.0220397    .0137525     -1.60    0.109     -.0489941     .0049147
            illw2 |    4.302823    .7770512      5.54    0.000      2.779831     5.825815
            radw2 |    .1069007    .0142987      7.48    0.000      .0788758     .1349256
------------------+----------------------------------------------------------------

.
. // regular male path model
. sem (crhtw1 -> crhtw2) (crhtw1 -> crhtw3) (crhtw1 -> PTSDw1) (crhtw1 -> PTSDw2) (crhtw1 -> MiPTSD) (crhtw2 -> crhtw3) ///
>  (crhtw2 -> PTSDw2) (crhtw2 -> MiPTSD) (crhtw3 -> MiPTSD) (PTSDw1 -> PTSDw2) (PTSDw1 -> MiPTSD) ///
>  (PTSDw2 -> MiPTSD) if female==0,  wrap(1) nocapslatent

Endogenous variables


Observed:  crhtw2 crhtw3 PTSDw1 PTSDw2 MiPTSD


Exogenous variables


Observed:  crhtw1


Fitting target model:

Iteration 0:   log likelihood = -4935.7279
Iteration 1:   log likelihood = -4935.7279

Structural equation model                       Number of obs      =       339
 Estimation method = ml
 Log likelihood    = -4935.7279

-------------+----------------------------------------------------------------
             |                 OIM
             |      Coef.    Std. Err.       z    P>|z|      [95% Con f. Interval]
-------------+----------------------------------------------------------------
 Structural   |
   crhtw2 <-  |
      crhtw1 |     .808673    .0311482     25.96    0.000      .7476236     .8697223
       _cons |   -.0690801    .0292554     -2.36    0.018     -.1264196    -.0117406
  -----------+----------------------------------------------------------------
   crhtw3 <-  |
      crhtw2 |    1.055212    .0258798     40.77    0.000      1.004489     1.105936
      crhtw1 |   -.1183935    .0256569     -4.61    0.000     -.1686801    -.0681068
       _cons |   -.0069225    .0140543     -0.49    0.622     -.0344684     .0206233
  -----------+----------------------------------------------------------------
   PTSDw1 <-  |
      crhtw1 |    17.44697    2.070283      8.43    0.000      13.38929     21.50465
       _cons |    26.78632    1.944476     13.78    0.000      22.97522     30.59742
  -----------+----------------------------------------------------------------
   PTSDw2 <-  |
      crhtw2 |    4.185639    1.124376      3.72    0.000      1.981903     6.389375
      PTSDw1 |    .1174034    .0169167      6.94    0.000      .0842474     .1505595
      crhtw1 |   -1.207601    1.140348     -1.06    0.290     -3.442642      1.02744
       _cons |    2.844572     .762135      3.73    0.000      1.350815     4.338329
  -----------+----------------------------------------------------------------
   MiPTSD <-  |
      crhtw2 |   -.2129708    2.145157     -0.10    0.921     -4.417402     3.991461
      crhtw3 |    5.742549    1.857988      3.09    0.002       2.10096     9.384138
      PTSDw1 |    .0899427    .0142199      6.33    0.000      .0620723     .1178131
      PTSDw2 |    .2641445    .0428121      6.17    0.000      .1802343     .3480547
      crhtw1 |   -1.673352    .9217186     -1.82    0.069     -3.479888     .1331828
       _cons |    44.43075     .611012     72.72    0.000      43.23319     45.62832
-------------+----------------------------------------------------------------
 Variance     |
    e.crhtw2 |    .2834994    .0217755                    .2438775     .3295585
    e.crhtw3 |    .0643684    .0049441                    .0553723     .0748261
    e.PTSDw1 |    1252.407    96.19683                    1077.371     1455.881
    e.PTSDw2 |    121.1661    9.306713                    104.2319     140.8515
    e.MiPTSD |    74.76247    5.742471                    64.31367     86.90886
-------------+----------------------------------------------------------------
LR test of model vs. saturated: chi2(3)   =      3.49, Prob > chi2 = 0.3223

. estat gof

---------------------+------------------------------------------------------
Fit statistic  |      Value   Description
---------------------+------------------------------------------------------
Likelihood ratio  |
chi2_ms(3)   |      3.488     model vs. saturated
p > chi2   |      0.322
chi2_bs(15)   |   1622.160     baseline vs. saturated
p > chi2   |      0.000
---------------------+------------------------------------------------------

. estat stable

Stability analysis of simultaneous equation systems

   Eigenvalue stability condition
  +--------------------------+-------------+
  |        Eigenvalue        |   Modulus   |
  +--------------------------+-------------+
  |  7.392e-06               |   7.4e-06   |
  | -3.696e-06 + 6.401e-06i  |   7.4e-06   |
  | -3.696e-06 - 6.401e-06i  |   7.4e-06   |
  |  1.039e-17 + 9.511e-10i  |   9.5e-10   |
  |  1.039e-17 - 9.511e-10i  |   9.5e-10   |
  +--------------------------+-------------+
   stability index =  7.39e-06
   All the eigenvalues lie inside the unit circle.
   SEM satisfies stability condition.

. estat teffects


Direct effects
-------------+----------------------------------------------------------------
             |                 OIM
             |      Coef.    Std. Err.       z    P>|z|      [95% Con f. Interval]
-------------+----------------------------------------------------------------
 Structural   |
   crhtw2 <-  |
      crhtw1 |     .808673    .0311482     25.96    0.000      .7476236     .8697223
  -----------+----------------------------------------------------------------
   crhtw3 <-  |
      crhtw2 |    1.055212    .0258798     40.77    0.000      1.004489     1.105936
      crhtw1 |   -.1183935    .0256569     -4.61    0.000     -.1686801    -.0681068
  -----------+----------------------------------------------------------------
   PTSDw1 <-  |
      crhtw1 |    17.44697    2.070283      8.43    0.000      13.38929     21.50465
  -----------+----------------------------------------------------------------
   PTSDw2 <-  |
      crhtw2 |    4.185639    1.124376      3.72    0.000      1.981903     6.389375
      PTSDw1 |    .1174034    .0169167      6.94    0.000      .0842474     .1505595
      crhtw1 |   -1.207601    1.140348     -1.06    0.290     -3.442642      1.02744
  -----------+----------------------------------------------------------------
   MiPTSD <-  |
      crhtw2 |   -.2129708    2.145157     -0.10    0.921     -4.417402     3.991461
      crhtw3 |    5.742549    1.857988      3.09    0.002       2.10096     9.384138
      PTSDw1 |    .0899427    .0142199      6.33    0.000      .0620723     .1178131
      PTSDw2 |    .2641445    .0428121      6.17    0.000      .1802343     .3480547
      crhtw1 |   -1.673352    .9217186     -1.82    0.069     -3.479888     .1331828
-------------+----------------------------------------------------------------


Indirect effects
-------------+----------------------------------------------------------------
             |                 OIM
             |      Coef.    Std. Err.       z    P>|z|      [95% Con f. Interval]
-------------+----------------------------------------------------------------
 Structural   |
   crhtw2 <-  |
      crhtw1 |           0   (no path)
  -----------+----------------------------------------------------------------
   crhtw3 <-  |
      crhtw2 |           0   (no path)
      crhtw1 |    .8533216    .0389653     21.90    0.000       .776951     .9296922
  -----------+----------------------------------------------------------------
   PTSDw1 <-  |
      crhtw1 |           0   (no path)
  -----------+----------------------------------------------------------------
   PTSDw2 <-  |
      crhtw2 |           0   (no path)
      PTSDw1 |           0   (no path)
      crhtw1 |    5.433147     .980717      5.54    0.000      3.510977     7.355317
  -----------+----------------------------------------------------------------
   MiPTSD <-  |
      crhtw2 |    7.165222    .3321059     21.58    0.000      6.514306     7.816138
      crhtw3 |           0   (no path)
      PTSDw1 |    .0310115    .0044684      6.94    0.000      .0222535     .0397695
      PTSDw2 |           0   (no path)
      crhtw1 |     6.73352    .8599438      7.83    0.000      5.048061     8.418979
-------------+----------------------------------------------------------------


Total effects
-------------+----------------------------------------------------------------
             |                 OIM
             |      Coef.    Std. Err.       z    P>|z|      [95% Con f. Interval]
-------------+----------------------------------------------------------------
 Structural   |
   crhtw2 <-  |
      crhtw1 |     .808673    .0311482     25.96    0.000      .7476236     .8697223
  -----------+----------------------------------------------------------------
   crhtw3 <-  |
      crhtw2 |    1.055212    .0258798     40.77    0.000      1.004489     1.105936
      crhtw1 |    .7349282    .0360637     20.38    0.000      .6642447     .8056117
  -----------+----------------------------------------------------------------
   PTSDw1 <-  |
      crhtw1 |    17.44697    2.070283      8.43    0.000      13.38929     21.50465
  -----------+----------------------------------------------------------------
   PTSDw2 <-  |
      crhtw2 |    4.185639    1.124376      3.72    0.000      1.981903     6.389375
      PTSDw1 |    .1174034    .0169167      6.94    0.000      .0842474     .1505595
      crhtw1 |    4.225546    .7005266      6.03    0.000      2.852539     5.598553
  -----------+----------------------------------------------------------------
   MiPTSD <-  |
      crhtw2 |    6.952251    2.170713      3.20    0.001      2.697732     11.20677
      crhtw3 |    5.742549    1.857988      3.09    0.002       2.10096     9.384138
      PTSDw1 |    .1209542    .0149054      8.11    0.000      .0917401     .1501683
      PTSDw2 |    .2641445    .0428121      6.17    0.000      .1802343     .3480547
      crhtw1 |    5.060167    .6337603      7.98    0.000       3.81802     6.302315
-------------+----------------------------------------------------------------

.
.
. // Robust path model
. sem (crhtw1 -> crhtw2) (crhtw1 -> crhtw3) (crhtw1 -> PTSDw1) (crhtw1 -> PTSDw2) (crhtw1 -> MiPTSD) (crhtw2 -> crhtw3) ///
>  (crhtw2 -> PTSDw2) (crhtw2 -> MiPTSD) (crhtw3 -> MiPTSD) (PTSDw1 -> PTSDw2) (PTSDw1 -> MiPTSD) ///
>  (PTSDw2 -> MiPTSD) if female==0, vce(cluster id) wrap(1) nocapslatent

Endogenous variables


Observed:  crhtw2 crhtw3 PTSDw1 PTSDw2 MiPTSD


Exogenous variables


Observed:  crhtw1


Fitting target model:

Iteration 0:   log pseudolikelihood = -4935.7279
Iteration 1:   log pseudolikelihood = -4935.7279

Structural equation model                       Number of obs      =       339
 Estimation method = ml
 Log pseudolikelihood = -4935.7279

(Std. Err. adjusted for 339 clusters in id)
-------------+----------------------------------------------------------------
             |               Robust
             |      Coef.    Std. Err.       z    P>|z|      [95% Con f. Interval]
-------------+----------------------------------------------------------------
 Structural   |
   crhtw2 <-  |
      crhtw1 |     .808673    .0331388     24.40    0.000      .7437222     .8736237
       _cons |   -.0690801    .0275271     -2.51    0.012     -.1230322    -.0151279
  -----------+----------------------------------------------------------------
   crhtw3 <-  |
      crhtw2 |    1.055212    .0329604     32.01    0.000       .990611     1.119814
      crhtw1 |   -.1183935    .0379136     -3.12    0.002     -.1927027    -.0440842
       _cons |   -.0069225    .0143437     -0.48    0.629     -.0350358     .0211907
  -----------+----------------------------------------------------------------
   PTSDw1 <-  |
      crhtw1 |    17.44697    2.143377      8.14    0.000      13.24603     21.64791
       _cons |    26.78632     2.10841     12.70    0.000      22.65391     30.91873
  -----------+----------------------------------------------------------------
   PTSDw2 <-  |
      crhtw2 |    4.185639    .8883731      4.71    0.000       2.44446     5.926818
      PTSDw1 |    .1174034    .0225907      5.20    0.000      .0731264     .1616804
      crhtw1 |   -1.207601    .8573454     -1.41    0.159     -2.887967     .4727654
       _cons |    2.844572    .4519448      6.29    0.000      1.958776     3.730367
  -----------+----------------------------------------------------------------
   MiPTSD <-  |
      crhtw2 |   -.2129708    2.284107     -0.09    0.926     -4.689737     4.263796
      crhtw3 |    5.742549    1.906455      3.01    0.003      2.005967     9.479132
      PTSDw1 |    .0899427     .016917      5.32    0.000       .056786     .1230994
      PTSDw2 |    .2641445    .0679108      3.89    0.000      .1310418     .3972472
      crhtw1 |   -1.673352    1.044293     -1.60    0.109     -3.720128     .3734236
       _cons |    44.43075    .6539724     67.94    0.000      43.14899     45.71252
-------------+----------------------------------------------------------------
 Variance     |
    e.crhtw2 |    .2834994    .0390594                    .2164098     .3713876
    e.crhtw3 |    .0643684    .0182932                    .0368779     .1123517
    e.PTSDw1 |    1252.407    86.18473                    1094.385     1433.248
    e.PTSDw2 |    121.1661     30.4715                    74.01444     198.3562
    e.MiPTSD |    74.76247    7.173899                    61.94494     90.23218
-------------+----------------------------------------------------------------

. estat gof


Note: model was fit with vce(cluster); only stats(residuals) valid.


. estat stable

Stability analysis of simultaneous equation systems

   Eigenvalue stability condition
  +--------------------------+-------------+
  |        Eigenvalue        |   Modulus   |
  +--------------------------+-------------+
  |  7.392e-06               |   7.4e-06   |
  | -3.696e-06 + 6.401e-06i  |   7.4e-06   |
  | -3.696e-06 - 6.401e-06i  |   7.4e-06   |
  |  1.039e-17 + 9.511e-10i  |   9.5e-10   |
  |  1.039e-17 - 9.511e-10i  |   9.5e-10   |
  +--------------------------+-------------+
   stability index =  7.39e-06
   All the eigenvalues lie inside the unit circle.
   SEM satisfies stability condition.

. estat teffects


Direct effects
(Std. Err. adjusted for 339 clusters in id)
-------------+----------------------------------------------------------------
             |               Robust
             |      Coef.    Std. Err.       z    P>|z|      [95% Con f. Interval]
-------------+----------------------------------------------------------------
 Structural   |
   crhtw2 <-  |
      crhtw1 |     .808673    .0331388     24.40    0.000      .7437222     .8736237
  -----------+----------------------------------------------------------------
   crhtw3 <-  |
      crhtw2 |    1.055212    .0329604     32.01    0.000       .990611     1.119814
      crhtw1 |   -.1183935    .0379136     -3.12    0.002     -.1927027    -.0440842
  -----------+----------------------------------------------------------------
   PTSDw1 <-  |
      crhtw1 |    17.44697    2.143377      8.14    0.000      13.24603     21.64791
  -----------+----------------------------------------------------------------
   PTSDw2 <-  |
      crhtw2 |    4.185639    .8883731      4.71    0.000       2.44446     5.926818
      PTSDw1 |    .1174034    .0225907      5.20    0.000      .0731264     .1616804
      crhtw1 |   -1.207601    .8573454     -1.41    0.159     -2.887967     .4727654
  -----------+----------------------------------------------------------------
   MiPTSD <-  |
      crhtw2 |   -.2129708    2.284107     -0.09    0.926     -4.689737     4.263796
      crhtw3 |    5.742549    1.906455      3.01    0.003      2.005967     9.479132
      PTSDw1 |    .0899427     .016917      5.32    0.000       .056786     .1230994
      PTSDw2 |    .2641445    .0679108      3.89    0.000      .1310418     .3972472
      crhtw1 |   -1.673352    1.044293     -1.60    0.109     -3.720128     .3734236
-------------+----------------------------------------------------------------


Indirect effects
(Std. Err. adjusted for 339 clusters in id)
-------------+----------------------------------------------------------------
             |               Robust
             |      Coef.    Std. Err.       z    P>|z|      [95% Con f. Interval]
-------------+----------------------------------------------------------------
 Structural   |
   crhtw2 <-  |
      crhtw1 |           0   (no path)
  -----------+----------------------------------------------------------------
   crhtw3 <-  |
      crhtw2 |           0   (no path)
      crhtw1 |    .8533216     .050095     17.03    0.000      .7551372     .9515061
  -----------+----------------------------------------------------------------
   PTSDw1 <-  |
      crhtw1 |           0   (no path)
  -----------+----------------------------------------------------------------
   PTSDw2 <-  |
      crhtw2 |           0   (no path)
      PTSDw1 |           0   (no path)
      crhtw1 |    5.433147    1.028607      5.28    0.000      3.417115     7.449179
  -----------+----------------------------------------------------------------
   MiPTSD <-  |
      crhtw2 |    7.165222    .3311412     21.64    0.000      6.516197     7.814247
      crhtw3 |           0   (no path)
      PTSDw1 |    .0310115    .0059672      5.20    0.000      .0193159      .042707
      PTSDw2 |           0   (no path)
      crhtw1 |     6.73352     1.01958      6.60    0.000      4.735179      8.73186
-------------+----------------------------------------------------------------


Total effects
(Std. Err. adjusted for 339 clusters in id)
-------------+----------------------------------------------------------------
             |               Robust
             |      Coef.    Std. Err.       z    P>|z|      [95% Con f. Interval]
-------------+----------------------------------------------------------------
 Structural   |
   crhtw2 <-  |
      crhtw1 |     .808673    .0331388     24.40    0.000      .7437222     .8736237
  -----------+----------------------------------------------------------------
   crhtw3 <-  |
      crhtw2 |    1.055212    .0329604     32.01    0.000       .990611     1.119814
      crhtw1 |    .7349282    .0384391     19.12    0.000      .6595889     .8102674
  -----------+----------------------------------------------------------------
   PTSDw1 <-  |
      crhtw1 |    17.44697    2.143377      8.14    0.000      13.24603     21.64791
  -----------+----------------------------------------------------------------
   PTSDw2 <-  |
      crhtw2 |    4.185639    .8883731      4.71    0.000       2.44446     5.926818
      PTSDw1 |    .1174034    .0225907      5.20    0.000      .0731264     .1616804
      crhtw1 |    4.225546    .9614717      4.39    0.000      2.341096     6.109996
  -----------+----------------------------------------------------------------
   MiPTSD <-  |
      crhtw2 |    6.952251    2.267127      3.07    0.002      2.508764     11.39574
      crhtw3 |    5.742549    1.906455      3.01    0.003      2.005967     9.479132
      PTSDw1 |    .1209542    .0189164      6.39    0.000      .0838787     .1580296
      PTSDw2 |    .2641445    .0679108      3.89    0.000      .1310418     .3972472
      crhtw1 |    5.060167    .7327402      6.91    0.000      3.624023     6.496312
-------------+----------------------------------------------------------------

.
.
. title "H6  Common perceptions of Chornobyl health" "related risk and PTSD for females"


*******************************************************************************
*******************************************************************************
*****                                                                     *****
*****                                                                     *****
*****             H6  Common perceptions of Chornobyl health              *****
*****                 related risk and PTSD for females                   *****
*****                                                                     *****
*****                                                                     *****
*****                                            26 Jul 2012    11:21:30  *****
*******************************************************************************
*******************************************************************************


. set more off

. sem( kmacc kmwork age->avgcumdosew1)(airw1 kmacc kmwork radhlw1 ->radchw1) ///
>     (radchw1->radchw2)(airw1->airw2)(age ->crhtw1) (PTSDw1->PTSDw2)  ///
>     (avgcumdosew1 age airw2 crhtw1 radchw1 radchw2->crhtw2) (havmil->kmwork) ///
>     ( radchw1 radhlw1->crhtw1)(kmacc->radchw2) (radchw1 kmacc  kmwork ->radchw2) ///
>     (age->radchw1) (avgcumdosew1->radchw2) ///
>      (age ->airw2) (crhtw1->airw2)(age avgcumdosew1 PTSDw1 PTSDw2 crhtw1 crhtw2->MiPTSD) ///
>       if gender==2, nocapslatent iterate(100)
(2 observations with missing values excluded;
 specify option 'method(mlmv)' to use all observations)

Endogenous variables


Observed:  kmwork avgcumdosew1 radchw1 radchw2 airw2 crhtw1 PTSDw2 crhtw2 MiPTSD


Exogenous variables


Observed:  kmacc age airw1 radhlw1 PTSDw1 havmil


Fitting target model:

Iteration 0:   log likelihood =  -22130.55
Iteration 1:   log likelihood =  -22130.55

Structural equation model                       Number of obs      =       361
 Estimation method = ml
 Log likelihood    =  -22130.55

------------------+----------------------------------------------------------------
                  |                 OIM
                  |      Coef.    Std. Err.       z    P>|z|      [95% Con f. Interval]
------------------+----------------------------------------------------------------
 Structural        |
   kmwork <-       |
           havmil |    .7678895    .0343703     22.34    0.000      .7005248     .8352541
            _cons |    115.5287    8.498615     13.59    0.000      98.87176     132.1857
  ----------------+----------------------------------------------------------------
   avgcumdosew1 <- |
           kmwork |   -.0003884    .0010885     -0.36    0.721     -.0025219     .0017451
            kmacc |    .0001138    .0010881      0.10    0.917     -.0020187     .0022464
              age |    .0083599    .0024082      3.47    0.001      .0036399       .01308
            _cons |   -.0307946    .1274701     -0.24    0.809     -.2806314     .2190422
  ----------------+----------------------------------------------------------------
   radchw1 <-      |
           kmwork |    .0429111    .0663036      0.65    0.518     -.0870416     .1728638
            kmacc |   -.0528118    .0662331     -0.80    0.425     -.1826262     .0770026
              age |    .0673804    .1508917      0.45    0.655     -.2283619     .3631227
            airw1 |    .1670343    .0530069      3.15    0.002      .0631426      .270926
          radhlw1 |    .3888016    .0505154      7.70    0.000      .2897933     .4878099
            _cons |    26.64206    8.263142      3.22    0.001       10.4466     42.83752
  ----------------+----------------------------------------------------------------
   radchw2 <-      |
           kmwork |    .0079469    .0465706      0.17    0.865     -.0833298     .0992235
     avgcumdosew1 |    4.420976    2.229713      1.98    0.047      .0508189     8.791134
          radchw1 |    .6597837    .0330808     19.94    0.000      .5949466     .7246208
            kmacc |    .0251429    .0465594      0.54    0.589     -.0661119     .1163977
            _cons |    14.13323    2.813062      5.02    0.000      8.619724     19.64673
  ----------------+----------------------------------------------------------------
   airw2 <-        |
           crhtw1 |    4.771939    2.084931      2.29    0.022      .6855498     8.858328
              age |    .3126984    .1592591      1.96    0.050      .0005563     .6248405
            airw1 |    .3413792    .0558761      6.11    0.000      .2318641     .4508942
            _cons |   -3.519867    8.850243     -0.40    0.691     -20.86602     13.82629
  ----------------+----------------------------------------------------------------
   crhtw1 <-       |
          radchw1 |    .0008685    .0003081      2.82    0.005      .0002646     .0014723
              age |    .0017182    .0008924      1.93    0.054      -.000031     .0034673
          radhlw1 |    .0242765    .0003149     77.09    0.000      .0236593     .0248937
            _cons |   -1.382818    .0459768    -30.08    0.000     -1.472931    -1.292705
  ----------------+----------------------------------------------------------------
   PTSDw2 <-       |
           PTSDw1 |    .0674057    .0125948      5.35    0.000      .0427203     .0920911
            _cons |    2.308176    .4765058      4.84    0.000      1.374242      3.24211
  ----------------+----------------------------------------------------------------
   crhtw2 <-       |
     avgcumdosew1 |    .1152761    .0477203      2.42    0.016      .0217462     .2088061
          radchw1 |   -.0123678    .0010597    -11.67    0.000     -.0144448    -.0102909
          radchw2 |       .0146    .0010744     13.59    0.000      .0124942     .0167058
            airw2 |    .0025974    .0007155      3.63    0.000      .0011951     .0039998
           crhtw1 |    .6674235    .0325239     20.52    0.000      .6036778     .7311692
              age |    .0061971     .002302      2.69    0.007      .0016853     .0107089
            _cons |    -.510907    .1255166     -4.07    0.000      -.756915    -.2648989
  ----------------+----------------------------------------------------------------
   MiPTSD <-       |
     avgcumdosew1 |    1.551965    1.044057      1.49    0.137     -.4943493      3.59828
           crhtw1 |   -4.871468     .798249     -6.10    0.000     -6.436007    -3.306928
           PTSDw2 |    .2589542     .069413      3.73    0.000      .1229072     .3950013
           crhtw2 |    6.261799    .8793748      7.12    0.000      4.538256     7.985342
              age |    .1595218    .0493551      3.23    0.001      .0627876      .256256
           PTSDw1 |    .0473405    .0174816      2.71    0.007      .0130772     .0816038
            _cons |    38.86229    2.471961     15.72    0.000      34.01734     43.70725
------------------+----------------------------------------------------------------
 Variance          |
         e.kmwork |    21773.32    1620.638                    18817.75      25193.1
   e.avgcumdosew1 |    .2883424     .021462                     .249202     .3336302
        e.radchw1 |    1064.138    79.20622                     919.689     1231.274
        e.radchw2 |    534.0213    39.74843                    461.5318     617.8961
          e.airw2 |    1193.041    88.80078                    1031.095     1380.423
         e.crhtw1 |    .0377576    .0028104                    .0326323     .0436879
         e.PTSDw2 |     63.3167    4.712807                    54.72192     73.26138
         e.crhtw2 |    .2389773    .0177876                    .2065379     .2765117
         e.MiPTSD |     108.173    8.051563                     93.4893     125.1629
------------------+----------------------------------------------------------------
LR test of model vs. saturated: chi2(60)  =   1328.32, Prob > chi2 = 0.0000

. estat gof

---------------------+------------------------------------------------------
Fit statistic  |      Value   Description
---------------------+------------------------------------------------------
Likelihood ratio  |
chi2_ms(60)   |   1328.323     model vs. saturated
p > chi2   |      0.000
chi2_bs(90)   |   3781.969     baseline vs. saturated
p > chi2   |      0.000
---------------------+------------------------------------------------------

. estat stable

Stability analysis of simultaneous equation systems

   Eigenvalue stability condition
  +--------------------------+-------------+
  |        Eigenvalue        |   Modulus   |
  +--------------------------+-------------+
  |  .00020949               |   .000209   |
  |  .00009893 + .00018164i  |   .000207   |
  |  .00009893 - .00018164i  |   .000207   |
  | -.00010516 +  .0001709i  |   .000201   |
  | -.00010516 -  .0001709i  |   .000201   |
  | -.00019703               |   .000197   |
  | -1.042e-16 + 1.201e-09i  |   1.2e-09   |
  | -1.042e-16 - 1.201e-09i  |   1.2e-09   |
  |  9.920e-19               |   9.9e-19   |
  +--------------------------+-------------+
   stability index =  .0002095
   All the eigenvalues lie inside the unit circle.
   SEM satisfies stability condition.

. estat mindices

Modification indices

------------------+-----------------------------------------------
                  |                                       Standard
                  |         MI      df    P>MI         EPC         EPC
------------------+-----------------------------------------------
 Structural        |
   kmwork <-       |
          radchw1 |      29.162       1    0.00   -1.502379   -.2565479
           crhtw1 |       9.168       1    0.00   -25.04916   -.1038661
            kmacc |     349.712       1    0.00    .9854817    .9864471
          radhlw1 |       9.569       1    0.00   -.6490071   -.1059972
  ----------------+-----------------------------------------------
   avgcumdosew1 <- |
           MiPTSD |       7.928       1    0.00    .0224692     .484807
           PTSDw1 |      15.402       1    0.00    .0033577    .2022973
  ----------------+-----------------------------------------------
   radchw1 <-      |
           havmil |       4.159       1    0.04   -.0302535   -.1757583
  ----------------+-----------------------------------------------
   radchw2 <-      |
            airw2 |       6.262       1    0.01    .0820942    .0899958
           PTSDw1 |       4.324       1    0.04   -.0761589    -.073878
  ----------------+-----------------------------------------------
   airw2 <-        |
           crhtw2 |       9.285       1    0.00    13.63817    .3188338
  ----------------+-----------------------------------------------
   PTSDw2 <-       |
           crhtw2 |       4.538       1    0.03    1.016017    .1079692
------------------+-----------------------------------------------
 Covariance        |
   e.kmwork        |
        e.radchw1 |       4.159       1    0.04    857.8305    .1782133
  ----------------+-----------------------------------------------
   e.radchw2       |
          e.airw2 |       9.565       1    0.00    129.9239    .1627729
------------------+-----------------------------------------------
EPC = expected parameter change

. estat teffects


Direct effects
------------------+----------------------------------------------------------------
                  |                 OIM
                  |      Coef.    Std. Err.       z    P>|z|      [95% Con f. Interval]
------------------+----------------------------------------------------------------
 Structural        |
   kmwork <-       |
           havmil |    .7678895    .0343703     22.34    0.000      .7005248     .8352541
  ----------------+----------------------------------------------------------------
   avgcumdosew1 <- |
           kmwork |   -.0003884    .0010885     -0.36    0.721     -.0025219     .0017451
            kmacc |    .0001138    .0010881      0.10    0.917     -.0020187     .0022464
              age |    .0083599    .0024082      3.47    0.001      .0036399       .01308
           havmil |           0   (no path)
  ----------------+----------------------------------------------------------------
   radchw1 <-      |
           kmwork |    .0429111    .0663036      0.65    0.518     -.0870416     .1728638
            kmacc |   -.0528118    .0662331     -0.80    0.425     -.1826262     .0770026
              age |    .0673804    .1508917      0.45    0.655     -.2283619     .3631227
            airw1 |    .1670343    .0530069      3.15    0.002      .0631426      .270926
          radhlw1 |    .3888016    .0505154      7.70    0.000      .2897933     .4878099
           havmil |           0   (no path)
  ----------------+----------------------------------------------------------------
   radchw2 <-      |
           kmwork |    .0079469    .0465706      0.17    0.865     -.0833298     .0992235
     avgcumdosew1 |    4.420976    2.229713      1.98    0.047      .0508189     8.791134
          radchw1 |    .6597837    .0330808     19.94    0.000      .5949466     .7246208
            kmacc |    .0251429    .0465594      0.54    0.589     -.0661119     .1163977
              age |           0   (no path)
            airw1 |           0   (no path)
          radhlw1 |           0   (no path)
           havmil |           0   (no path)
  ----------------+----------------------------------------------------------------
   airw2 <-        |
           kmwork |           0   (no path)
          radchw1 |           0   (no path)
           crhtw1 |    4.771939    2.084931      2.29    0.022      .6855498     8.858328
            kmacc |           0   (no path)
              age |    .3126984    .1592591      1.96    0.050      .0005563     .6248405
            airw1 |    .3413792    .0558761      6.11    0.000      .2318641     .4508942
          radhlw1 |           0   (no path)
           havmil |           0   (no path)
  ----------------+----------------------------------------------------------------
   crhtw1 <-       |
           kmwork |           0   (no path)
          radchw1 |    .0008685    .0003081      2.82    0.005      .0002646     .0014723
            kmacc |           0   (no path)
              age |    .0017182    .0008924      1.93    0.054      -.000031     .0034673
            airw1 |           0   (no path)
          radhlw1 |    .0242765    .0003149     77.09    0.000      .0236593     .0248937
           havmil |           0   (no path)
  ----------------+----------------------------------------------------------------
   PTSDw2 <-       |
           PTSDw1 |    .0674057    .0125948      5.35    0.000      .0427203     .0920911
  ----------------+----------------------------------------------------------------
   crhtw2 <-       |
           kmwork |           0   (no path)
     avgcumdosew1 |    .1152761    .0477203      2.42    0.016      .0217462     .2088061
          radchw1 |   -.0123678    .0010597    -11.67    0.000     -.0144448    -.0102909
          radchw2 |       .0146    .0010744     13.59    0.000      .0124942     .0167058
            airw2 |    .0025974    .0007155      3.63    0.000      .0011951     .0039998
           crhtw1 |    .6674235    .0325239     20.52    0.000      .6036778     .7311692
            kmacc |           0   (no path)
              age |    .0061971     .002302      2.69    0.007      .0016853     .0107089
            airw1 |           0   (no path)
          radhlw1 |           0   (no path)
           havmil |           0   (no path)
  ----------------+----------------------------------------------------------------
   MiPTSD <-       |
           kmwork |           0   (no path)
     avgcumdosew1 |    1.551965    1.044057      1.49    0.137     -.4943493      3.59828
          radchw1 |           0   (no path)
          radchw2 |           0   (no path)
            airw2 |           0   (no path)
           crhtw1 |   -4.871468     .798249     -6.10    0.000     -6.436007    -3.306928
           PTSDw2 |    .2589542     .069413      3.73    0.000      .1229072     .3950013
           crhtw2 |    6.261799    .8793748      7.12    0.000      4.538256     7.985342
            kmacc |           0   (no path)
              age |    .1595218    .0493551      3.23    0.001      .0627876      .256256
            airw1 |           0   (no path)
          radhlw1 |           0   (no path)
           PTSDw1 |    .0473405    .0174816      2.71    0.007      .0130772     .0816038
           havmil |           0   (no path)
------------------+----------------------------------------------------------------


Indirect effects
------------------+----------------------------------------------------------------
                  |                 OIM
                  |      Coef.    Std. Err.       z    P>|z|      [95% Con f. Interval]
------------------+----------------------------------------------------------------
 Structural        |
   kmwork <-       |
           havmil |           0   (no path)
  ----------------+----------------------------------------------------------------
   avgcumdosew1 <- |
           kmwork |           0   (no path)
            kmacc |           0   (no path)
              age |           0   (no path)
           havmil |   -.0002983     .000836     -0.36    0.721     -.0019368     .0013402
  ----------------+----------------------------------------------------------------
   radchw1 <-      |
           kmwork |           0   (no path)
            kmacc |           0   (no path)
              age |           0   (no path)
            airw1 |           0   (no path)
          radhlw1 |           0   (no path)
           havmil |     .032951    .0509352      0.65    0.518     -.0668802     .1327821
  ----------------+----------------------------------------------------------------
   radchw2 <-      |
           kmwork |    .0265949    .0440099      0.60    0.546      -.059663     .1128528
     avgcumdosew1 |           0   (no path)
          radchw1 |           0   (no path)
            kmacc |   -.0343412    .0439992     -0.78    0.435     -.1205782     .0518957
              age |    .0814155     .101853      0.80    0.424     -.1182126     .2810436
            airw1 |    .1102065    .0354069      3.11    0.002      .0408101     .1796028
          radhlw1 |     .256525    .0357248      7.18    0.000      .1865056     .3265444
           havmil |    .0265242    .0492117      0.54    0.590     -.0699289     .1229774
  ----------------+----------------------------------------------------------------
   airw2 <-        |
           kmwork |    .0001778    .0002748      0.65    0.518     -.0003607     .0007164
          radchw1 |    .0041443    .0014702      2.82    0.005      .0012628     .0070257
           crhtw1 |           0   (no path)
            kmacc |   -.0002189    .0003009     -0.73    0.467     -.0008085     .0003708
              age |    .0084783    .0056774      1.49    0.135     -.0026492     .0196058
            airw1 |    .0006922    .0004473      1.55    0.122     -.0001844     .0015688
          radhlw1 |    .1174574    .0513372      2.29    0.022      .0168383     .2180764
           havmil |    .0001366    .0002246      0.61    0.543     -.0003037     .0005769
  ----------------+----------------------------------------------------------------
   crhtw1 <-       |
           kmwork |    .0000373    .0000576      0.65    0.518     -.0000756     .0001501
          radchw1 |           0   (no path)
            kmacc |   -.0000459    .0000598     -0.77    0.443      -.000163     .0000713
              age |    .0000585    .0001327      0.44    0.659     -.0002015     .0003186
            airw1 |    .0001451     .000069      2.10    0.036      9.74e-06     .0002804
          radhlw1 |    .0003377    .0001276      2.65    0.008      .0000876     .0005877
           havmil |    .0000286    .0000454      0.63    0.528     -.0000603     .0001176
  ----------------+----------------------------------------------------------------
   PTSDw2 <-       |
           PTSDw1 |           0   (no path)
  ----------------+----------------------------------------------------------------
   crhtw2 <-       |
           kmwork |   -.0000458    .0007217     -0.06    0.949     -.0014603     .0013686
     avgcumdosew1 |    .0645464    .0325539      1.98    0.047       .000742     .1283508
          radchw1 |    .0102233    .0005264     19.42    0.000      .0091915     .0112551
          radchw2 |           0   (no path)
            airw2 |           0   (no path)
           crhtw1 |    .0123948    .0054155      2.29    0.022      .0017807      .023009
            kmacc |    .0005008    .0007227      0.69    0.488     -.0009156     .0019172
              age |    .0033391    .0010555      3.16    0.002      .0012703     .0054078
            airw1 |    .0005285    .0003446      1.53    0.125     -.0001469     .0012039
          radhlw1 |    .0156698    .0007737     20.25    0.000      .0141534     .0171863
           havmil |   -.0000352    .0005549     -0.06    0.949     -.0011229     .0010525
  ----------------+----------------------------------------------------------------
   MiPTSD <-       |
           kmwork |   -.0010714    .0052911     -0.20    0.840     -.0114418      .009299
     avgcumdosew1 |    1.126013    .3617227      3.11    0.002      .4170491     1.834976
          radchw1 |   -.0176595    .0072947     -2.42    0.015     -.0319568    -.0033621
          radchw2 |    .0914224    .0067277     13.59    0.000      .0782364     .1046085
            airw2 |    .0162647    .0044802      3.63    0.000      .0074836     .0250458
           crhtw1 |    4.256886    .2064622     20.62    0.000      3.852227     4.661544
           PTSDw2 |           0   (no path)
           crhtw2 |           0   (no path)
            kmacc |     .003536    .0053157      0.67    0.506     -.0068826     .0139546
              age |    .0640327    .0192951      3.32    0.001       .026215     .1018504
            airw1 |    .0026027    .0022665      1.15    0.251     -.0018395     .0070448
          radhlw1 |   -.0217859    .0155778     -1.40    0.162     -.0523178     .0087459
           PTSDw1 |     .017455    .0057034      3.06    0.002      .0062765     .0286335
           havmil |   -.0008227    .0040787     -0.20    0.840     -.0088169     .0071714
------------------+----------------------------------------------------------------


Total effects
------------------+----------------------------------------------------------------
                  |                 OIM
                  |      Coef.    Std. Err.       z    P>|z|      [95% Con f. Interval]
------------------+----------------------------------------------------------------
 Structural        |
   kmwork <-       |
           havmil |    .7678895    .0343703     22.34    0.000      .7005248     .8352541
  ----------------+----------------------------------------------------------------
   avgcumdosew1 <- |
           kmwork |   -.0003884    .0010885     -0.36    0.721     -.0025219     .0017451
            kmacc |    .0001138    .0010881      0.10    0.917     -.0020187     .0022464
              age |    .0083599    .0024082      3.47    0.001      .0036399       .01308
           havmil |   -.0002983     .000836     -0.36    0.721     -.0019368     .0013402
  ----------------+----------------------------------------------------------------
   radchw1 <-      |
           kmwork |    .0429111    .0663036      0.65    0.518     -.0870416     .1728638
            kmacc |   -.0528118    .0662331     -0.80    0.425     -.1826262     .0770026
              age |    .0673804    .1508917      0.45    0.655     -.2283619     .3631227
            airw1 |    .1670343    .0530069      3.15    0.002      .0631426      .270926
          radhlw1 |    .3888016    .0505154      7.70    0.000      .2897933     .4878099
           havmil |     .032951    .0509352      0.65    0.518     -.0668802     .1327821
  ----------------+----------------------------------------------------------------
   radchw2 <-      |
           kmwork |    .0345418    .0640757      0.54    0.590     -.0910443     .1601278
     avgcumdosew1 |    4.420976    2.229713      1.98    0.047      .0508189     8.791134
          radchw1 |    .6597837    .0330808     19.94    0.000      .5949466     .7246208
            kmacc |   -.0091983    .0640048     -0.14    0.886     -.1346455     .1162489
              age |    .0814155     .101853      0.80    0.424     -.1182126     .2810436
            airw1 |    .1102065    .0354069      3.11    0.002      .0408101     .1796028
          radhlw1 |     .256525    .0357248      7.18    0.000      .1865056     .3265444
           havmil |    .0265242    .0492117      0.54    0.590     -.0699289     .1229774
  ----------------+----------------------------------------------------------------
   airw2 <-        |
           kmwork |    .0001778    .0002748      0.65    0.518     -.0003607     .0007164
          radchw1 |    .0041443    .0014702      2.82    0.005      .0012628     .0070257
           crhtw1 |    4.771939    2.084931      2.29    0.022      .6855498     8.858328
            kmacc |   -.0002189    .0003009     -0.73    0.467     -.0008085     .0003708
              age |    .3211767    .1584558      2.03    0.043       .010609     .6317444
            airw1 |    .3420714      .05579      6.13    0.000      .2327249     .4514179
          radhlw1 |    .1174574    .0513372      2.29    0.022      .0168383     .2180764
           havmil |    .0001366    .0002246      0.61    0.543     -.0003037     .0005769
  ----------------+----------------------------------------------------------------
   crhtw1 <-       |
           kmwork |    .0000373    .0000576      0.65    0.518     -.0000756     .0001501
          radchw1 |    .0008685    .0003081      2.82    0.005      .0002646     .0014723
            kmacc |   -.0000459    .0000598     -0.77    0.443      -.000163     .0000713
              age |    .0017767    .0009016      1.97    0.049      9.57e-06     .0035438
            airw1 |    .0001451     .000069      2.10    0.036      9.74e-06     .0002804
          radhlw1 |    .0246142    .0002874     85.65    0.000      .0240509     .0251775
           havmil |    .0000286    .0000454      0.63    0.528     -.0000603     .0001176
  ----------------+----------------------------------------------------------------
   PTSDw2 <-       |
           PTSDw1 |    .0674057    .0125948      5.35    0.000      .0427203     .0920911
  ----------------+----------------------------------------------------------------
   crhtw2 <-       |
           kmwork |   -.0000458    .0007217     -0.06    0.949     -.0014603     .0013686
     avgcumdosew1 |    .1798225    .0577666      3.11    0.002      .0666021      .293043
          radchw1 |   -.0021446    .0011832     -1.81    0.070     -.0044637     .0001746
          radchw2 |       .0146    .0010744     13.59    0.000      .0124942     .0167058
            airw2 |    .0025974    .0007155      3.63    0.000      .0011951     .0039998
           crhtw1 |    .6798183    .0329717     20.62    0.000       .615195     .7444417
            kmacc |    .0005008    .0007227      0.69    0.488     -.0009156     .0019172
              age |    .0095362    .0024347      3.92    0.000      .0047642     .0143081
            airw1 |    .0005285    .0003446      1.53    0.125     -.0001469     .0012039
          radhlw1 |    .0156698    .0007737     20.25    0.000      .0141534     .0171863
           havmil |   -.0000352    .0005549     -0.06    0.949     -.0011229     .0010525
  ----------------+----------------------------------------------------------------
   MiPTSD <-       |
           kmwork |   -.0010714    .0052911     -0.20    0.840     -.0114418      .009299
     avgcumdosew1 |    2.677978    1.104943      2.42    0.015      .5123295     4.843626
          radchw1 |   -.0176595    .0072947     -2.42    0.015     -.0319568    -.0033621
          radchw2 |    .0914224    .0067277     13.59    0.000      .0782364     .1046085
            airw2 |    .0162647    .0044802      3.63    0.000      .0074836     .0250458
           crhtw1 |   -.6145821    .8245169     -0.75    0.456     -2.230606     1.001441
           PTSDw2 |    .2589542     .069413      3.73    0.000      .1229072     .3950013
           crhtw2 |    6.261799    .8793748      7.12    0.000      4.538256     7.985342
            kmacc |     .003536    .0053157      0.67    0.506     -.0068826     .0139546
              age |    .2235545    .0506688      4.41    0.000      .1242455     .3228635
            airw1 |    .0026027    .0022665      1.15    0.251     -.0018395     .0070448
          radhlw1 |   -.0217859    .0155778     -1.40    0.162     -.0523178     .0087459
           PTSDw1 |    .0647955    .0172434      3.76    0.000      .0309991     .0985919
           havmil |   -.0008227    .0040787     -0.20    0.840     -.0088169     .0071714
------------------+----------------------------------------------------------------

.
. // regular female path
. set more off

. sem (crhtw1 age kmwork goferw1  fdferw2-> crhtw2) (crhtw1  goferw1-> crhtw3) (crhtw1 fdferw1 -> PTSDw1) ///
>  (crhtw1 fdferw1-> PTSDw2) (crhtw1 fdferw2 -> MiPTSD) (crhtw2 -> crhtw3) ///
>  (crhtw2 -> PTSDw2) (crhtw2 -> MiPTSD)  (PTSDw1 -> PTSDw2) (PTSDw1 -> MiPTSD) ///
>  (age havmil kmwork PTSDw2 -> MiPTSD) if female==1,  wrap(1) nocapslatent
(1 observations with missing values excluded;
 specify option 'method(mlmv)' to use all observations)

Endogenous variables


Observed:  crhtw2 crhtw3 PTSDw1 PTSDw2 MiPTSD


Exogenous variables


Observed:  crhtw1 age kmwork goferw1 fdferw2 fdferw1 havmil


Fitting target model:

Iteration 0:   log likelihood = -16555.948
Iteration 1:   log likelihood = -16555.948

Structural equation model                       Number of obs      =       362
 Estimation method = ml
 Log likelihood    = -16555.948

-------------+----------------------------------------------------------------
             |                 OIM
             |      Coef.    Std. Err.       z    P>|z|      [95% Con f. Interval]
-------------+----------------------------------------------------------------
 Structural   |
   crhtw2 <-  |
      crhtw1 |    .5978997    .0352374     16.97    0.000      .5288356     .6669638
         age |    .0067204    .0028199      2.38    0.017      .0011935     .0122473
      kmwork |     .000302    .0001433      2.11    0.035      .0000212     .0005829
     goferw1 |     .003697    .0009446      3.91    0.000      .0018457     .0055483
     fdferw2 |    .0031687    .0013177      2.40    0.016      .0005861     .0057513
       _cons |   -.4607364    .1448095     -3.18    0.001     -.7445579    -.1769149
  -----------+----------------------------------------------------------------
   crhtw3 <-  |
      crhtw2 |    1.035551    .0241144     42.94    0.000      .9882879     1.082815
      crhtw1 |   -.1126471    .0215249     -5.23    0.000     -.1548352     -.070459
     goferw1 |    .0010505    .0004297      2.44    0.015      .0002082     .0018928
       _cons |   -.0216496    .0191719     -1.13    0.259     -.0592257     .0159266
  -----------+----------------------------------------------------------------
   PTSDw1 <-  |
      crhtw1 |    .1802023    1.830742      0.10    0.922     -3.407986     3.768391
     fdferw1 |    .1905237    .0456693      4.17    0.000      .1010135     .2800339
       _cons |    10.84177    2.397857      4.52    0.000       6.14206     15.54149
  -----------+----------------------------------------------------------------
   PTSDw2 <-  |
      crhtw2 |    .8650436    .6682401      1.29    0.195     -.4446829      2.17477
      PTSDw1 |    .0638212    .0127881      4.99    0.000       .038757     .0888854
      crhtw1 |    .0940193    .6013894      0.16    0.876     -1.084682     1.272721
     fdferw1 |    .0093967    .0116693      0.81    0.421     -.0134748     .0322682
       _cons |    1.852971    .6005676      3.09    0.002      .6758803     3.030062
  -----------+----------------------------------------------------------------
   MiPTSD <-  |
      crhtw2 |     5.61472     .863879      6.50    0.000      3.921548     7.307892
      PTSDw1 |    .0337288    .0169678      1.99    0.047      .0004725     .0669851
      PTSDw2 |    .2458808    .0672094      3.66    0.000      .1141528     .3776087
      crhtw1 |   -4.723429    .7902744     -5.98    0.000     -6.272338    -3.174519
         age |    .1975397    .0474226      4.17    0.000      .1045931     .2904863
      kmwork |   -.0067987    .0036698     -1.85    0.064     -.0139914      .000394
     fdferw2 |    .1074096     .022037      4.87    0.000      .0642178     .1506014
      havmil |    .0065438    .0036331      1.80    0.072     -.0005769     .0136646
       _cons |    36.88496     2.46019     14.99    0.000      32.06307     41.70684
-------------+----------------------------------------------------------------
 Variance     |
    e.crhtw2 |    .3595649    .0267262                    .3108193     .4159553
    e.crhtw3 |    .0793505    .0058981                    .0685931      .091795
    e.PTSDw1 |    1051.174    78.13318                    908.6686      1216.03
    e.PTSDw2 |    62.22138    4.624878                    53.78613     71.97952
    e.MiPTSD |     100.934    7.502365                    87.25055     116.7634
-------------+----------------------------------------------------------------
LR test of model vs. saturated: chi2(23)  =     27.08, Prob > chi2 = 0.2527

. estat gof

---------------------+------------------------------------------------------
Fit statistic  |      Value   Description
---------------------+------------------------------------------------------
Likelihood ratio  |
chi2_ms(23)   |     27.077     model vs. saturated
p > chi2   |      0.253
chi2_bs(45)   |   1312.214     baseline vs. saturated
p > chi2   |      0.000
---------------------+------------------------------------------------------

. estat stable

Stability analysis of simultaneous equation systems

   Eigenvalue stability condition
  +--------------------------+-------------+
  |        Eigenvalue        |   Modulus   |
  +--------------------------+-------------+
  | -2.198e-16 + 8.325e-09i  |   8.3e-09   |
  | -2.198e-16 - 8.325e-09i  |   8.3e-09   |
  | -2.356e-10               |   2.4e-10   |
  |  2.356e-10               |   2.4e-10   |
  |          0               |         0   |
  +--------------------------+-------------+
   stability index =  8.32e-09
   All the eigenvalues lie inside the unit circle.
   SEM satisfies stability condition.

. estat teffects


Direct effects
-------------+----------------------------------------------------------------
             |                 OIM
             |      Coef.    Std. Err.       z    P>|z|      [95% Con f. Interval]
-------------+----------------------------------------------------------------
 Structural   |
   crhtw2 <-  |
      crhtw1 |    .5978997    .0352374     16.97    0.000      .5288356     .6669638
         age |    .0067204    .0028199      2.38    0.017      .0011935     .0122473
      kmwork |     .000302    .0001433      2.11    0.035      .0000212     .0005829
     goferw1 |     .003697    .0009446      3.91    0.000      .0018457     .0055483
     fdferw2 |    .0031687    .0013177      2.40    0.016      .0005861     .0057513
  -----------+----------------------------------------------------------------
   crhtw3 <-  |
      crhtw2 |    1.035551    .0241144     42.94    0.000      .9882879     1.082815
      crhtw1 |   -.1126471    .0215249     -5.23    0.000     -.1548352     -.070459
         age |           0   (no path)
      kmwork |           0   (no path)
     goferw1 |    .0010505    .0004297      2.44    0.015      .0002082     .0018928
     fdferw2 |           0   (no path)
  -----------+----------------------------------------------------------------
   PTSDw1 <-  |
      crhtw1 |    .1802023    1.830742      0.10    0.922     -3.407986     3.768391
     fdferw1 |    .1905237    .0456693      4.17    0.000      .1010135     .2800339
  -----------+----------------------------------------------------------------
   PTSDw2 <-  |
      crhtw2 |    .8650436    .6682401      1.29    0.195     -.4446829      2.17477
      PTSDw1 |    .0638212    .0127881      4.99    0.000       .038757     .0888854
      crhtw1 |    .0940193    .6013894      0.16    0.876     -1.084682     1.272721
         age |           0   (no path)
      kmwork |           0   (no path)
     goferw1 |           0   (no path)
     fdferw2 |           0   (no path)
     fdferw1 |    .0093967    .0116693      0.81    0.421     -.0134748     .0322682
  -----------+----------------------------------------------------------------
   MiPTSD <-  |
      crhtw2 |     5.61472     .863879      6.50    0.000      3.921548     7.307892
      PTSDw1 |    .0337288    .0169678      1.99    0.047      .0004725     .0669851
      PTSDw2 |    .2458808    .0672094      3.66    0.000      .1141528     .3776087
      crhtw1 |   -4.723429    .7902744     -5.98    0.000     -6.272338    -3.174519
         age |    .1975397    .0474226      4.17    0.000      .1045931     .2904863
      kmwork |   -.0067987    .0036698     -1.85    0.064     -.0139914      .000394
     goferw1 |           0   (no path)
     fdferw2 |    .1074096     .022037      4.87    0.000      .0642178     .1506014
     fdferw1 |           0   (no path)
      havmil |    .0065438    .0036331      1.80    0.072     -.0005769     .0136646
-------------+----------------------------------------------------------------


Indirect effects
-------------+----------------------------------------------------------------
             |                 OIM
             |      Coef.    Std. Err.       z    P>|z|      [95% Con f. Interval]
-------------+----------------------------------------------------------------
 Structural   |
   crhtw2 <-  |
      crhtw1 |           0   (no path)
         age |           0   (no path)
      kmwork |           0   (no path)
     goferw1 |           0   (no path)
     fdferw2 |           0   (no path)
  -----------+----------------------------------------------------------------
   crhtw3 <-  |
      crhtw2 |           0   (no path)
      crhtw1 |    .6191558    .0392353     15.78    0.000      .5422559     .6960556
         age |    .0069593    .0029246      2.38    0.017      .0012271     .0126915
      kmwork |    .0003128    .0001486      2.11    0.035      .0000216      .000604
     goferw1 |    .0038285    .0009822      3.90    0.000      .0019034     .0057535
     fdferw2 |    .0032814    .0013667      2.40    0.016      .0006028       .00596
  -----------+----------------------------------------------------------------
   PTSDw1 <-  |
      crhtw1 |           0   (no path)
     fdferw1 |           0   (no path)
  -----------+----------------------------------------------------------------
   PTSDw2 <-  |
      crhtw2 |           0   (no path)
      PTSDw1 |           0   (no path)
      crhtw1 |    .5287101    .4174198      1.27    0.205     -.2894178     1.346838
         age |    .0058134    .0051106      1.14    0.255     -.0042031       .01583
      kmwork |    .0002613    .0002369      1.10    0.270      -.000203     .0007255
     goferw1 |    .0031981    .0026021      1.23    0.219      -.001902     .0082981
     fdferw2 |    .0027411    .0024048      1.14    0.254     -.0019722     .0074544
     fdferw1 |    .0121595    .0037989      3.20    0.001      .0047138     .0196051
  -----------+----------------------------------------------------------------
   MiPTSD <-  |
      crhtw2 |    .2126976    .1643074      1.29    0.195      -.109339     .5347342
      PTSDw1 |    .0156924    .0031443      4.99    0.000      .0095296     .0218552
      PTSDw2 |           0   (no path)
      crhtw1 |    3.516235    .5733546      6.13    0.000       2.39248     4.639989
         age |    .0391626    .0174609      2.24    0.025      .0049399     .0733853
      kmwork |    .0017601    .0008761      2.01    0.045      .0000429     .0034773
     goferw1 |    .0215441    .0063909      3.37    0.001      .0090181     .0340701
     fdferw2 |    .0184655    .0081675      2.26    0.024      .0024575     .0344735
     fdferw1 |    .0117264    .0048194      2.43    0.015      .0022805     .0211722
      havmil |           0   (no path)
-------------+----------------------------------------------------------------


Total effects
-------------+----------------------------------------------------------------
             |                 OIM
             |      Coef.    Std. Err.       z    P>|z|      [95% Con f. Interval]
-------------+----------------------------------------------------------------
 Structural   |
   crhtw2 <-  |
      crhtw1 |    .5978997    .0352374     16.97    0.000      .5288356     .6669638
         age |    .0067204    .0028199      2.38    0.017      .0011935     .0122473
      kmwork |     .000302    .0001433      2.11    0.035      .0000212     .0005829
     goferw1 |     .003697    .0009446      3.91    0.000      .0018457     .0055483
     fdferw2 |    .0031687    .0013177      2.40    0.016      .0005861     .0057513
  -----------+----------------------------------------------------------------
   crhtw3 <-  |
      crhtw2 |    1.035551    .0241144     42.94    0.000      .9882879     1.082815
      crhtw1 |    .5065086    .0397895     12.73    0.000      .4285227     .5844946
         age |    .0069593    .0029246      2.38    0.017      .0012271     .0126915
      kmwork |    .0003128    .0001486      2.11    0.035      .0000216      .000604
     goferw1 |    .0048789    .0010626      4.59    0.000      .0027964     .0069615
     fdferw2 |    .0032814    .0013667      2.40    0.016      .0006028       .00596
  -----------+----------------------------------------------------------------
   PTSDw1 <-  |
      crhtw1 |    .1802023    1.830742      0.10    0.922     -3.407986     3.768391
     fdferw1 |    .1905237    .0456693      4.17    0.000      .1010135     .2800339
  -----------+----------------------------------------------------------------
   PTSDw2 <-  |
      crhtw2 |    .8650436    .6682401      1.29    0.195     -.4446829      2.17477
      PTSDw1 |    .0638212    .0127881      4.99    0.000       .038757     .0888854
      crhtw1 |    .6227294    .4615092      1.35    0.177      -.281812     1.527271
         age |    .0058134    .0051106      1.14    0.255     -.0042031       .01583
      kmwork |    .0002613    .0002369      1.10    0.270      -.000203     .0007255
     goferw1 |    .0031981    .0026021      1.23    0.219      -.001902     .0082981
     fdferw2 |    .0027411    .0024048      1.14    0.254     -.0019722     .0074544
     fdferw1 |    .0215562    .0117725      1.83    0.067     -.0015175     .0446299
  -----------+----------------------------------------------------------------
   MiPTSD <-  |
      crhtw2 |    5.827418    .8793656      6.63    0.000      4.103893     7.550943
      PTSDw1 |    .0494212    .0172567      2.86    0.004      .0155987     .0832437
      PTSDw2 |    .2458808    .0672094      3.66    0.000      .1141528     .3776087
      crhtw1 |   -1.207194    .6409172     -1.88    0.060     -2.463368     .0489807
         age |    .2367023     .049625      4.77    0.000       .139439     .3339655
      kmwork |   -.0050386    .0037553     -1.34    0.180     -.0123989     .0023217
     goferw1 |    .0215441    .0063909      3.37    0.001      .0090181     .0340701
     fdferw2 |    .1258751    .0230024      5.47    0.000      .0807912      .170959
     fdferw1 |    .0117264    .0048194      2.43    0.015      .0022805     .0211722
      havmil |    .0065438    .0036331      1.80    0.072     -.0005769     .0136646
-------------+----------------------------------------------------------------

.
. // robust path model
. sem (crhtw1 -> crhtw2) (crhtw1 -> crhtw3) (crhtw1 -> PTSDw1) (crhtw1 -> PTSDw2) (crhtw1 -> MiPTSD) (crhtw2 -> crhtw3) ///
>  (crhtw2 -> PTSDw2) (crhtw2 -> MiPTSD) (crhtw3 -> MiPTSD) (PTSDw1 -> PTSDw2) (PTSDw1 -> MiPTSD) ///
>  (PTSDw2 -> MiPTSD) if female==1, vce(cluster id) wrap(1) nocapslatent
(1 observations with missing values excluded;
 specify option 'method(mlmv)' to use all observations)

Endogenous variables


Observed:  crhtw2 crhtw3 PTSDw1 PTSDw2 MiPTSD


Exogenous variables


Observed:  crhtw1


Fitting target model:

Iteration 0:   log pseudolikelihood = -5310.8935
Iteration 1:   log pseudolikelihood = -5310.8935

Structural equation model                       Number of obs      =       362
 Estimation method = ml
 Log pseudolikelihood = -5310.8935

(Std. Err. adjusted for 362 clusters in id)
-------------+----------------------------------------------------------------
             |               Robust
             |      Coef.    Std. Err.       z    P>|z|      [95% Con f. Interval]
-------------+----------------------------------------------------------------
 Structural   |
   crhtw2 <-  |
      crhtw1 |    .6298151    .0400677     15.72    0.000      .5512839     .7083464
       _cons |    .0885394    .0360796      2.45    0.014      .0178246     .1592541
  -----------+----------------------------------------------------------------
   crhtw3 <-  |
      crhtw2 |    1.051234    .0277728     37.85    0.000      .9967999     1.105667
      crhtw1 |   -.1166559    .0316122     -3.69    0.000     -.1786147     -.054697
       _cons |    .0072487    .0158942      0.46    0.648     -.0239034     .0384008
  -----------+----------------------------------------------------------------
   PTSDw1 <-  |
      crhtw1 |    1.423881    1.804336      0.79    0.430     -2.112553     4.960315
       _cons |    17.80774    1.740305     10.23    0.000       14.3968     21.21867
  -----------+----------------------------------------------------------------
   PTSDw2 <-  |
      crhtw2 |    .9851263    .5144778      1.91    0.056     -.0232317     1.993484
      PTSDw1 |    .0659966    .0151775      4.35    0.000      .0362494     .0957439
      crhtw1 |    .0766308    .5442594      0.14    0.888     -.9900981      1.14336
       _cons |    2.147164    .2859438      7.51    0.000      1.586725     2.707604
  -----------+----------------------------------------------------------------
   MiPTSD <-  |
      crhtw2 |    3.611711     2.45526      1.47    0.141      -1.20051     8.423931
      crhtw3 |    3.099064    2.173102      1.43    0.154     -1.160138     7.358265
      PTSDw1 |    .0583282    .0187505      3.11    0.002      .0215779     .0950785
      PTSDw2 |    .2669377       .0937      2.85    0.004      .0832891     .4505864
      crhtw1 |   -4.365202    .9294167     -4.70    0.000     -6.186826    -2.543579
       _cons |    47.04229    .5888074     79.89    0.000      45.88825     48.19634
-------------+----------------------------------------------------------------
 Variance     |
    e.crhtw2 |    .4056681    .0388891                    .3361793     .4895204
    e.crhtw3 |    .0806603    .0123716                    .0597178     .1089472
    e.PTSDw1 |    1101.712     98.1412                    925.2148     1311.879
    e.PTSDw2 |    62.33283    25.07038                    28.33753     137.1108
    e.MiPTSD |    111.4816    11.44536                    91.16201     136.3304
-------------+----------------------------------------------------------------

. estat gof


Note: model was fit with vce(cluster); only stats(residuals) valid.


. estat stable

Stability analysis of simultaneous equation systems

   Eigenvalue stability condition
  +--------------------------+-------------+
  |        Eigenvalue        |   Modulus   |
  +--------------------------+-------------+
  |  2.597e-08               |   2.6e-08   |
  | -2.597e-08               |   2.6e-08   |
  | -2.796e-17 + 1.313e-09i  |   1.3e-09   |
  | -2.796e-17 - 1.313e-09i  |   1.3e-09   |
  |  8.882e-16               |   8.9e-16   |
  +--------------------------+-------------+
   stability index =  2.60e-08
   All the eigenvalues lie inside the unit circle.
   SEM satisfies stability condition.

. estat teffects


Direct effects
(Std. Err. adjusted for 362 clusters in id)
-------------+----------------------------------------------------------------
             |               Robust
             |      Coef.    Std. Err.       z    P>|z|      [95% Con f. Interval]
-------------+----------------------------------------------------------------
 Structural   |
   crhtw2 <-  |
      crhtw1 |    .6298151    .0400677     15.72    0.000      .5512839     .7083464
  -----------+----------------------------------------------------------------
   crhtw3 <-  |
      crhtw2 |    1.051234    .0277728     37.85    0.000      .9967999     1.105667
      crhtw1 |   -.1166559    .0316122     -3.69    0.000     -.1786147     -.054697
  -----------+----------------------------------------------------------------
   PTSDw1 <-  |
      crhtw1 |    1.423881    1.804336      0.79    0.430     -2.112553     4.960315
  -----------+----------------------------------------------------------------
   PTSDw2 <-  |
      crhtw2 |    .9851263    .5144778      1.91    0.056     -.0232317     1.993484
      PTSDw1 |    .0659966    .0151775      4.35    0.000      .0362494     .0957439
      crhtw1 |    .0766308    .5442594      0.14    0.888     -.9900981      1.14336
  -----------+----------------------------------------------------------------
   MiPTSD <-  |
      crhtw2 |    3.611711     2.45526      1.47    0.141      -1.20051     8.423931
      crhtw3 |    3.099064    2.173102      1.43    0.154     -1.160138     7.358265
      PTSDw1 |    .0583282    .0187505      3.11    0.002      .0215779     .0950785
      PTSDw2 |    .2669377       .0937      2.85    0.004      .0832891     .4505864
      crhtw1 |   -4.365202    .9294167     -4.70    0.000     -6.186826    -2.543579
-------------+----------------------------------------------------------------


Indirect effects
(Std. Err. adjusted for 362 clusters in id)
-------------+----------------------------------------------------------------
             |               Robust
             |      Coef.    Std. Err.       z    P>|z|      [95% Con f. Interval]
-------------+----------------------------------------------------------------
 Structural   |
   crhtw2 <-  |
      crhtw1 |           0   (no path)
  -----------+----------------------------------------------------------------
   crhtw3 <-  |
      crhtw2 |           0   (no path)
      crhtw1 |    .6620828    .0497097     13.32    0.000      .5646536      .759512
  -----------+----------------------------------------------------------------
   PTSDw1 <-  |
      crhtw1 |           0   (no path)
  -----------+----------------------------------------------------------------
   PTSDw2 <-  |
      crhtw2 |           0   (no path)
      PTSDw1 |           0   (no path)
      crhtw1 |    .7144187    .3638013      1.96    0.050      .0013813     1.427456
  -----------+----------------------------------------------------------------
   MiPTSD <-  |
      crhtw2 |    3.520807    .1764301     19.96    0.000       3.17501     3.866604
      crhtw3 |           0   (no path)
      PTSDw1 |     .017617    .0040514      4.35    0.000      .0096763     .0255577
      PTSDw2 |           0   (no path)
      crhtw1 |    4.259236    .7494689      5.68    0.000      2.790304     5.728168
-------------+----------------------------------------------------------------


Total effects
(Std. Err. adjusted for 362 clusters in id)
-------------+----------------------------------------------------------------
             |               Robust
             |      Coef.    Std. Err.       z    P>|z|      [95% Con f. Interval]
-------------+----------------------------------------------------------------
 Structural   |
   crhtw2 <-  |
      crhtw1 |    .6298151    .0400677     15.72    0.000      .5512839     .7083464
  -----------+----------------------------------------------------------------
   crhtw3 <-  |
      crhtw2 |    1.051234    .0277728     37.85    0.000      .9967999     1.105667
      crhtw1 |    .5454269    .0449103     12.14    0.000      .4574044     .6334495
  -----------+----------------------------------------------------------------
   PTSDw1 <-  |
      crhtw1 |    1.423881    1.804336      0.79    0.430     -2.112553     4.960315
  -----------+----------------------------------------------------------------
   PTSDw2 <-  |
      crhtw2 |    .9851263    .5144778      1.91    0.056     -.0232317     1.993484
      PTSDw1 |    .0659966    .0151775      4.35    0.000      .0362494     .0957439
      crhtw1 |    .7910495    .4830154      1.64    0.101     -.1556434     1.737742
  -----------+----------------------------------------------------------------
   MiPTSD <-  |
      crhtw2 |    7.132518    2.483768      2.87    0.004      2.264423     12.00061
      crhtw3 |    3.099064    2.173102      1.43    0.154     -1.160138     7.358265
      PTSDw1 |    .0759452    .0187906      4.04    0.000      .0391163     .1127741
      PTSDw2 |    .2669377       .0937      2.85    0.004      .0832891     .4505864
      crhtw1 |   -.1059661    .7470832     -0.14    0.887     -1.570222      1.35829
-------------+----------------------------------------------------------------

.
.
.  save chwide25jul2012, replace
file chwide25jul2012.dta saved

.  use ch3wMaster22jul2012, clear

.  mi xtset id wave
       panel variable:  id (strongly balanced)
        time variable:   wave, 1 to 3
                delta:  1 unit

.  use chwide25jul2012, clear
(Chornobyl sequelae for both sexes   women=362  men = 340 wide format)

. label var crhtw2 "Chornobyl related health threat in wave 2 alpha=.822"

. label var crhtw3 "Chornobyl rleated health threat in wave 3 alpha=0.833"

.
.  save chwide22jul2012, replace
file chwide22jul2012.dta saved

.
.
.  use ch3w22jul2012, clear

. sjlog using PTSDpanel, replace