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

.
.
. forvalues k=1/2 {
  2. forvalues i=1/3 {
  3. set more off
  4. mvreg  WHPpain WHPsleep WHPsociso WHPpa WHPel = crhtw`i' avgcumdosew`i' icdxcnt  ///
>    airw`i' radchw`i' kzchorn if gender==`k'
  5. }
  6. }

Equation          Obs  Parms        RMSE    "R-sq"          F        P
----------------------------------------------------------------------
WHPpain           339      7    15.30871    0.1510    9.84158   0.0000
WHPsleep          339      7    22.82093    0.1713   11.43585   0.0000
WHPsociso         339      7    14.74458    0.0336   1.926532   0.0759
WHPpa             339      7    13.69334    0.1315   8.378596   0.0000
WHPel             339      7    27.85346    0.1553   10.17627   0.0000

-------------+----------------------------------------------------------------
             |      Coef.    Std. Err.       t    P>|t|      [95% Con f. Interval]
-------------+----------------------------------------------------------------
WHPpain      |
      crhtw1 |    1.712645    1.133175      1.51    0.132     -.5164626     3.941752
avgcumdosew1 |    .9218099    .5025185      1.83    0.067     -.0667118     1.910332
     icdxcnt |    2.124077    .5243019      4.05    0.000      1.092705      3.15545
       airw1 |   -.0504658     .026276     -1.92    0.056     -.1021542     .0012225
     radchw1 |    .0364728    .0236814      1.54    0.124     -.0101117     .0830572
     kzchorn |    .1378862    .0376103      3.67    0.000      .0639016     .2118707
       _cons |    -2.98708    3.522282     -0.85    0.397     -9.915885     3.941725
-------------+----------------------------------------------------------------
WHPsleep     |
      crhtw1 |    5.190844    1.689241      3.07    0.002       1.86788     8.513808
avgcumdosew1 |    1.481022     .749112      1.98    0.049       .007418     2.954627
     icdxcnt |    2.637959    .7815847      3.38    0.001      1.100476     4.175441
       airw1 |   -.0268993      .03917     -0.69    0.493      -.103952     .0501533
     radchw1 |   -.0174828    .0353022     -0.50    0.621     -.0869269     .0519613
     kzchorn |    .1640065    .0560663      2.93    0.004      .0537166     .2742964
       _cons |    2.894135     5.25072      0.55    0.582      -7.43474     13.22301
-------------+----------------------------------------------------------------
WHPsociso    |
      crhtw1 |    1.019617    1.091416      0.93    0.351     -1.127347      3.16658
avgcumdosew1 |    .3525934    .4840003      0.73    0.467     -.5995006     1.304687
     icdxcnt |    1.484703    .5049809      2.94    0.004      .4913377     2.478069
       airw1 |   -.0006631    .0253077     -0.03    0.979     -.0504468     .0491205
     radchw1 |    .0031856    .0228087      0.14    0.889     -.0416822     .0480534
     kzchorn |   -.0253454    .0362244     -0.70    0.485     -.0966036     .0459127
       _cons |    4.801267    3.392484      1.42    0.158     -1.872206     11.47474
-------------+----------------------------------------------------------------
WHPpa        |
      crhtw1 |    2.171147    1.013602      2.14    0.033      .1772546      4.16504
avgcumdosew1 |   -.2369793    .4494928     -0.53    0.598     -1.121192     .6472338
     icdxcnt |    2.078993    .4689776      4.43    0.000       1.15645     3.001535
       airw1 |   -.0092656    .0235033     -0.39    0.694     -.0554998     .0369687
     radchw1 |    .0372417    .0211825      1.76    0.080     -.0044271     .0789106
     kzchorn |    .0650998    .0336417      1.94    0.054      -.001078     .1312775
       _cons |   -.4878267    3.150612     -0.15    0.877     -6.685506     5.709852
-------------+----------------------------------------------------------------
WHPel        |
      crhtw1 |    6.619025    2.061756      3.21    0.001      2.563273     10.67478
avgcumdosew1 |    .9815096     .914308      1.07    0.284     -.8170577     2.780077
     icdxcnt |    5.727662    .9539417      6.00    0.000       3.85113     7.604195
       airw1 |   -.0699187    .0478079     -1.46    0.145     -.1639632     .0241258
     radchw1 |    .0061817    .0430871      0.14    0.886     -.0785764     .0909399
     kzchorn |    -.004252    .0684301     -0.06    0.950     -.1388633     .1303593
       _cons |    15.72922    6.408622      2.45    0.015      3.122598     28.33585
-------------+----------------------------------------------------------------

Equation          Obs  Parms        RMSE    "R-sq"          F        P
----------------------------------------------------------------------
WHPpain           339      7    15.05612    0.1788   12.04679   0.0000
WHPsleep          339      7    22.65333    0.1834   12.42748   0.0000
WHPsociso         339      7    14.61233    0.0509   2.967633   0.0078
WHPpa             339      7    13.41739    0.1662   11.02624   0.0000
WHPel             339      7    27.43416    0.1806   12.19404   0.0000

-------------+----------------------------------------------------------------
             |      Coef.    Std. Err.       t    P>|t|      [95% Con f. Interval]
-------------+----------------------------------------------------------------
WHPpain      |
      crhtw2 |    3.958005    1.250566      3.16    0.002      1.497972     6.418038
avgcumdosew2 |    .5012325    .3308508      1.51    0.131     -.1495957     1.152061
     icdxcnt |    1.586542    .5255613      3.02    0.003      .5526924     2.620393
       airw2 |   -.0011869     .027823     -0.04    0.966     -.0559185     .0535447
     radchw2 |     .057236    .0246006      2.33    0.021      .0088433     .1056286
     kzchorn |     .062269    .0409127      1.52    0.129     -.0182118     .1427498
       _cons |      -.0537    3.576045     -0.02    0.988     -7.088264     6.980864
-------------+----------------------------------------------------------------
WHPsleep     |
      crhtw2 |     7.58543    1.881594      4.03    0.000      3.884082     11.28678
avgcumdosew2 |     .773677    .4977958      1.55    0.121     -.2055545     1.752909
     icdxcnt |    1.950945    .7907559      2.47    0.014      .3954214     3.506469
       airw2 |   -.0106997    .0418623     -0.26    0.798     -.0930485      .071649
     radchw2 |    -.022871    .0370138     -0.62    0.537     -.0956822     .0499402
     kzchorn |    .0956969     .061557      1.55    0.121      -.025394     .2167878
       _cons |    8.483896    5.380493      1.58    0.116      -2.10026     19.06805
-------------+----------------------------------------------------------------
WHPsociso    |
      crhtw2 |    3.063891    1.213706      2.52    0.012      .6763678     5.451413
avgcumdosew2 |    .2070704    .3210988      0.64    0.519     -.4245743     .8387151
     icdxcnt |    1.227748    .5100702      2.41    0.017      .2243708     2.231125
       airw2 |   -.0140995    .0270029     -0.52    0.602     -.0672179     .0390188
     radchw2 |    .0116536    .0238755      0.49    0.626     -.0353127     .0586198
     kzchorn |   -.0678827    .0397068     -1.71    0.088     -.1459913     .0102259
       _cons |    8.592364     3.47064      2.48    0.014      1.765146     15.41958
-------------+----------------------------------------------------------------
WHPpa        |
      crhtw2 |    3.383965    1.114453      3.04    0.003      1.191685     5.576244
avgcumdosew2 |   -.1380143    .2948405     -0.47    0.640     -.7180053     .4419768
     icdxcnt |     1.71787    .4683585      3.67    0.000      .7965454     2.639194
       airw2 |    .0215353    .0247947      0.87    0.386     -.0272392     .0703098
     radchw2 |    .0606458     .021923      2.77    0.006      .0175202     .1037713
     kzchorn |    .0145793    .0364597      0.40    0.690     -.0571419     .0863005
       _cons |    2.128083    3.186823      0.67    0.505     -4.140829     8.396994
-------------+----------------------------------------------------------------
WHPel        |
      crhtw2 |    8.942615    2.278691      3.92    0.000      4.460122     13.42511
avgcumdosew2 |    1.235125    .6028522      2.05    0.041      .0492331     2.421017
     icdxcnt |    4.884435    .9576396      5.10    0.000      3.000628     6.768241
       airw2 |   -.0740653     .050697     -1.46    0.145     -.1737932     .0256626
     radchw2 |    .0326922    .0448254      0.73    0.466     -.0554853     .1208697
     kzchorn |   -.0765414    .0745482     -1.03    0.305     -.2231877     .0701049
       _cons |    18.70985     6.51601      2.87    0.004      5.891981     31.52772
-------------+----------------------------------------------------------------

Equation          Obs  Parms        RMSE    "R-sq"          F        P
----------------------------------------------------------------------
WHPpain           339      7    14.84559    0.2016   13.97143   0.0000
WHPsleep          339      7    22.63905    0.1844   12.51301   0.0000
WHPsociso         339      7    14.60505    0.0518   3.025755   0.0068
WHPpa             339      7     13.3822    0.1705   11.37572   0.0000
WHPel             339      7    27.72325    0.1632    10.7931   0.0000

-------------+----------------------------------------------------------------
             |      Coef.    Std. Err.       t    P>|t|      [95% Con f. Interval]
-------------+----------------------------------------------------------------
WHPpain      |
      crhtw3 |     4.59055    1.250811      3.67    0.000      2.130036     7.051063
avgcumdosew3 |    .4459064    .3057968      1.46    0.146     -.1556372      1.04745
     icdxcnt |    1.495087    .5214908      2.87    0.004      .4692447      2.52093
       airw3 |     .004801    .0251673      0.19    0.849     -.0447065     .0543084
     radchw3 |    .0724206    .0244798      2.96    0.003      .0242655     .1205757
     kzchorn |    .0440838     .040416      1.09    0.276     -.0354199     .1235876
       _cons |    .7401192     3.57644      0.21    0.836     -6.295221     7.775459
-------------+----------------------------------------------------------------
WHPsleep     |
      crhtw3 |    7.236111    1.907447      3.79    0.000      3.483906     10.98832
avgcumdosew3 |     .711297    .4663304      1.53    0.128     -.2060379     1.628632
     icdxcnt |    1.822198    .7952569      2.29    0.023      .2578206     3.386576
       airw3 |    .0237258    .0383793      0.62    0.537     -.0517715     .0992231
     radchw3 |   -.0435792    .0373309     -1.17    0.244     -.1170142     .0298558
     kzchorn |     .094378    .0616332      1.53    0.127     -.0268627     .2156188
       _cons |    8.492664    5.453957      1.56    0.120     -2.236005     19.22133
-------------+----------------------------------------------------------------
WHPsociso    |
      crhtw3 |    3.019644    1.230545      2.45    0.015      .5989969     5.440292
avgcumdosew3 |    .1713677    .3008422      0.57    0.569     -.4204295     .7631649
     icdxcnt |    1.265912    .5130414      2.47    0.014      .2566907     2.275134
       airw3 |   -.0219494    .0247595     -0.89    0.376     -.0706547     .0267559
     radchw3 |    .0177519    .0240832      0.74    0.462     -.0296229     .0651268
     kzchorn |   -.0634542    .0397612     -1.60    0.111     -.1416698     .0147614
       _cons |    8.358894    3.518493      2.38    0.018      1.437543     15.28024
-------------+----------------------------------------------------------------
WHPpa        |
      crhtw3 |    3.937198    1.127513      3.49    0.001      1.719227     6.155168
avgcumdosew3 |   -.1214166    .2756531     -0.44    0.660     -.6636635     .4208304
     icdxcnt |    1.684881    .4700853      3.58    0.000      .7601598     2.609602
       airw3 |    .0010018    .0226864      0.04    0.965     -.0436254     .0456291
     radchw3 |    .0614401    .0220667      2.78    0.006      .0180318     .1048483
     kzchorn |    .0078929     .036432      0.22    0.829     -.0637739     .0795596
       _cons |    3.381091    3.223895      1.05    0.295     -2.960746     9.722928
-------------+----------------------------------------------------------------
WHPel        |
      crhtw3 |    7.086379    2.335814      3.03    0.003      2.491516     11.68124
avgcumdosew3 |    1.197562    .5710573      2.10    0.037      .0742156     2.320909
     icdxcnt |    4.882645    .9738529      5.01    0.000      2.966945     6.798345
       airw3 |   -.0172007    .0469984     -0.37    0.715     -.1096529     .0752515
     radchw3 |    .0026988    .0457146      0.06    0.953      -.087228     .0926255
     kzchorn |   -.0502755    .0754745     -0.67    0.506      -.198744     .0981931
       _cons |    16.59266    6.678788      2.48    0.013      3.454586     29.73074
-------------+----------------------------------------------------------------

Equation          Obs  Parms        RMSE    "R-sq"          F        P
----------------------------------------------------------------------
WHPpain           361      7    21.50405    0.0895   5.801383   0.0000
WHPsleep          361      7    29.65421    0.0986   6.450315   0.0000
WHPsociso         361      7    18.32113    0.0658   4.158485   0.0005
WHPpa             361      7    20.39215    0.1109   7.359559   0.0000
WHPel             361      7    33.57185    0.0665   4.203936   0.0004

-------------+----------------------------------------------------------------
             |      Coef.    Std. Err.       t    P>|t|      [95% Con f. Interval]
-------------+----------------------------------------------------------------
WHPpain      |
      crhtw1 |    2.423568    1.460139      1.66    0.098     -.4480699     5.295206
avgcumdosew1 |    6.790181    2.078671      3.27    0.001      2.702084     10.87828
     icdxcnt |    1.873775    .5096138      3.68    0.000      .8715237     2.876026
       airw1 |   -.0571147    .0353731     -1.61    0.107     -.1266825     .0124531
     radchw1 |   -.0061522    .0355415     -0.17    0.863     -.0760511     .0637468
     kzchorn |    .0352114    .0450506      0.78    0.435      -.053389     .1238118
       _cons |    10.98627     4.77626      2.30    0.022      1.592857     20.37968
-------------+----------------------------------------------------------------
WHPsleep     |
      crhtw1 |    6.014866     2.01354      2.99    0.003      2.054861      9.97487
avgcumdosew1 |    8.255107    2.866498      2.88    0.004        2.6176     13.89261
     icdxcnt |    2.305302    .7027602      3.28    0.001      .9231917     3.687412
       airw1 |   -.0298568    .0487797     -0.61    0.541     -.1257911     .0660776
     radchw1 |   -.0737494    .0490119     -1.50    0.133     -.1701405     .0226416
     kzchorn |    .0824447     .062125      1.33    0.185     -.0397357     .2046252
       _cons |    15.53004    6.586488      2.36    0.019      2.576475     28.48361
-------------+----------------------------------------------------------------
WHPsociso    |
      crhtw1 |    3.334535    1.244016      2.68    0.008      .8879429     5.781127
avgcumdosew1 |    4.373699    1.770996      2.47    0.014      .8907026     7.856695
     icdxcnt |    .5932136    .4341832      1.37    0.173     -.2606893     1.447116
       airw1 |   -.0324346    .0301373     -1.08    0.283     -.0917053     .0268361
     radchw1 |   -.0089379    .0302808     -0.30    0.768     -.0684907      .050615
     kzchorn |    .0612291    .0383824      1.60    0.112     -.0142571     .1367153
       _cons |    4.489256    4.069301      1.10    0.271     -3.513788      12.4923
-------------+----------------------------------------------------------------
WHPpa        |
      crhtw1 |    4.247966     1.38464      3.07    0.002       1.52481     6.971122
avgcumdosew1 |    4.408479     1.97119      2.24    0.026      .5317635     8.285194
     icdxcnt |    2.029621    .4832634      4.20    0.000      1.079193      2.98005
       airw1 |   -.0267979    .0335441     -0.80    0.425     -.0927686     .0391728
     radchw1 |    .0326285    .0337037      0.97    0.334     -.0336562     .0989133
     kzchorn |      .03977    .0427212      0.93    0.353     -.0442492     .1237892
       _cons |    6.769598    4.529296      1.49    0.136     -2.138113     15.67731
-------------+----------------------------------------------------------------
WHPel        |
      crhtw1 |    4.148502     2.27955      1.82    0.070     -.3346612     8.631666
avgcumdosew1 |    5.920427    3.245193      1.82    0.069     -.4618559     12.30271
     icdxcnt |     2.30064    .7956023      2.89    0.004      .7359385     3.865341
       airw1 |    -.081897     .055224     -1.48    0.139     -.1905053     .0267113
     radchw1 |   -.0212054    .0554869     -0.38    0.703     -.1303307       .08792
     kzchorn |    .1229607    .0703324      1.75    0.081     -.0153611     .2612825
       _cons |    19.36588    7.456634      2.60    0.010      4.701009     34.03075
-------------+----------------------------------------------------------------

Equation          Obs  Parms        RMSE    "R-sq"          F        P
----------------------------------------------------------------------
WHPpain           363      7    21.17946    0.1150   7.711937   0.0000
WHPsleep          363      7    29.24663    0.1192   8.027174   0.0000
WHPsociso         363      7    18.13818    0.0808   5.213322   0.0000
WHPpa             363      7    20.21939    0.1246   8.444649   0.0000
WHPel             363      7    33.26445    0.0827   5.351154   0.0000

-------------+----------------------------------------------------------------
             |      Coef.    Std. Err.       t    P>|t|      [95% Con f. Interval]
-------------+----------------------------------------------------------------
WHPpain      |
      crhtw2 |    7.410312    1.695029      4.37    0.000      4.076782     10.74384
avgcumdosew2 |     1.63405     .823119      1.99    0.048       .015263     3.252837
     icdxcnt |    1.771595    .4941783      3.58    0.000      .7997187      2.74347
       airw2 |   -.0645679    .0315139     -2.05    0.041     -.1265447    -.0025911
     radchw2 |    -.046027    .0370114     -1.24    0.214     -.1188154     .0267614
     kzchorn |   -.0405731    .0500955     -0.81    0.419     -.1390935     .0579472
       _cons |    17.96384    4.677356      3.84    0.000      8.765119     27.16256
-------------+----------------------------------------------------------------
WHPsleep     |
      crhtw2 |    10.03493    2.340659      4.29    0.000       5.43167     14.63818
avgcumdosew2 |    2.275549    1.136642      2.00    0.046      .0401725     4.510925
     icdxcnt |     1.96937    .6824088      2.89    0.004      .6273106     3.311429
       airw2 |     .004131    .0435174      0.09    0.924     -.0814525     .0897145
     radchw2 |   -.1276013    .0511088     -2.50    0.013     -.2281145    -.0270881
     kzchorn |    .0028551    .0691767      0.04    0.967     -.1331912     .1389015
       _cons |    23.86073    6.458942      3.69    0.000      11.15826     36.56321
-------------+----------------------------------------------------------------
WHPsociso    |
      crhtw2 |    5.534162     1.45163      3.81    0.000      2.679313     8.389011
avgcumdosew2 |    1.366355    .7049226      1.94    0.053     -.0199807     2.752691
     icdxcnt |    .4147829    .4232164      0.98    0.328     -.4175356     1.247102
       airw2 |   -.0449927    .0269886     -1.67    0.096     -.0980699     .0080845
     radchw2 |   -.0322396    .0316967     -1.02    0.310     -.0945759     .0300967
     kzchorn |    .0199242     .042902      0.46    0.643      -.064449     .1042975
       _cons |    8.855995    4.005708      2.21    0.028      .9781698     16.73382
-------------+----------------------------------------------------------------
WHPpa        |
      crhtw2 |    6.426416    1.618193      3.97    0.000      3.243997     9.608835
avgcumdosew2 |    1.348209    .7858066      1.72    0.087      -.197198     2.893615
     icdxcnt |    1.716006     .471777      3.64    0.000      .7881856     2.643826
       airw2 |    .0208287    .0300854      0.69    0.489     -.0383387      .079996
     radchw2 |    .0079681    .0353336      0.23    0.822     -.0615207      .077457
     kzchorn |   -.0327161    .0478247     -0.68    0.494     -.1267705     .0613383
       _cons |    12.01987    4.465329      2.69    0.007      3.238127     20.80161
-------------+----------------------------------------------------------------
WHPel        |
      crhtw2 |    8.981656    2.662212      3.37    0.001      3.746017      14.2173
avgcumdosew2 |    1.806486     1.29279      1.40    0.163     -.7359799     4.348952
     icdxcnt |    2.063021    .7761562      2.66    0.008      .5365934     3.589449
       airw2 |   -.0265804    .0494957     -0.54    0.592     -.1239212     .0707603
     radchw2 |   -.0686238      .05813     -1.18    0.239     -.1829452     .0456976
     kzchorn |    .0285538      .07868      0.36    0.717     -.1261822     .1832899
       _cons |    25.16568    7.346253      3.43    0.001      10.71817     39.61319
-------------+----------------------------------------------------------------

Equation          Obs  Parms        RMSE    "R-sq"          F        P
----------------------------------------------------------------------
WHPpain           363      7    21.08499    0.1229   8.314122   0.0000
WHPsleep          363      7    29.11717    0.1269    8.62752   0.0000
WHPsociso         363      7    18.02687    0.0920   6.012917   0.0000
WHPpa             363      7    20.20668    0.1257   8.529938   0.0000
WHPel             363      7    33.04538    0.0948   6.211657   0.0000

-------------+----------------------------------------------------------------
             |      Coef.    Std. Err.       t    P>|t|      [95% Con f. Interval]
-------------+----------------------------------------------------------------
WHPpain      |
      crhtw3 |    8.335972    1.701378      4.90    0.000      4.989957     11.68199
avgcumdosew3 |    1.091759    .6420331      1.70    0.090     -.1708951     2.354414
     icdxcnt |    1.809017    .4917099      3.68    0.000      .8419959     2.776038
       airw3 |   -.0407478    .0315584     -1.29    0.197     -.1028122     .0213166
     radchw3 |    -.029961     .036462     -0.82    0.412     -.1016689      .041747
     kzchorn |   -.0797151    .0509318     -1.57    0.118     -.1798802       .02045
       _cons |    19.32865    4.746513      4.07    0.000      9.993922     28.66338
-------------+----------------------------------------------------------------
WHPsleep     |
      crhtw3 |    11.17041    2.349506      4.75    0.000      6.549758     15.79107
avgcumdosew3 |    1.598725    .8866112      1.80    0.072      -.144929     3.342379
     icdxcnt |    2.007304    .6790235      2.96    0.003      .6719023     3.342705
       airw3 |    .0077981    .0435804      0.18    0.858     -.0779093     .0935055
     radchw3 |   -.1110622    .0503519     -2.21    0.028     -.2100868    -.0120376
     kzchorn |   -.0343486     .070334     -0.49    0.626     -.1726709     .1039738
       _cons |    25.19239    6.554665      3.84    0.000      12.30166     38.08312
-------------+----------------------------------------------------------------
WHPsociso    |
      crhtw3 |    6.118552    1.454614      4.21    0.000      3.257835     8.979269
avgcumdosew3 |    .9515525    .5489141      1.73    0.084     -.1279695     2.031074
     icdxcnt |    .4875994    .4203935      1.16    0.247     -.3391675     1.314366
       airw3 |   -.0669775    .0269813     -2.48    0.014     -.1200402    -.0139148
     radchw3 |   -.0132443    .0311736     -0.42    0.671     -.0745519     .0480633
     kzchorn |    .0011118    .0435448      0.03    0.980     -.0845256     .0867492
       _cons |    10.19502     4.05809      2.51    0.012      2.214177     18.17586
-------------+----------------------------------------------------------------
WHPpa        |
      crhtw3 |    6.497179    1.630506      3.98    0.000      3.290544     9.703814
avgcumdosew3 |    .9891908    .6152888      1.61    0.109      -.220867     2.199249
     icdxcnt |    1.746575    .4712275      3.71    0.000      .8198352     2.673314
       airw3 |    .0076068    .0302438      0.25    0.802     -.0518723     .0670858
     radchw3 |    .0247044    .0349431      0.71    0.480     -.0440165     .0934253
     kzchorn |   -.0460092    .0488102     -0.94    0.347     -.1420019     .0499834
       _cons |    12.33897    4.548794      2.71    0.007      3.393084     21.28486
-------------+----------------------------------------------------------------
WHPel        |
      crhtw3 |    11.13369    2.666479      4.18    0.000      5.889655     16.37772
avgcumdosew3 |    .9760696    1.006224      0.97    0.333     -1.002821     2.954961
     icdxcnt |     2.07994    .7706308      2.70    0.007       .564379     3.595501
       airw3 |   -.0339945    .0494598     -0.69    0.492     -.1312647     .0632757
     radchw3 |   -.0772145    .0571449     -1.35    0.177     -.1895985     .0351696
     kzchorn |   -.0143613    .0798228     -0.18    0.857     -.1713447     .1426221
       _cons |    29.50761    7.438957      3.97    0.000      14.87778     44.13743
-------------+----------------------------------------------------------------

.
.
. title "Zero order pain relationships for men and women in different waves"


*******************************************************************************
*******************************************************************************
*****                                                                     *****
*****                                                                     *****
***** Zero order pain relationships for men and women in different waves  *****
*****                                                                     *****
*****                                                                     *****
*****                                            12 Jul 2012    23:56:13  *****
*******************************************************************************
*******************************************************************************


. forvalues k=1/2 {
  2. regress WHPpain crhtw1 crhtw2 crhtw3 if gender==`k'
  3. // zero order relationships in wave 1 and 3 for men
. // zero order relationship in wave three for women
. }

      Source |       SS       df       MS              Number of obs =     339
-------------+------------------------------           F(  3,   335) =   22.62
       Model |  15434.8724     3  5144.95746           Prob > F      =  0.0000
    Residual |   76210.199   335  227.493131           R-squared     =  0.1684
-------------+------------------------------           Adj R-squared =  0.1610
       Total |  91645.0713   338  271.139264           Root MSE      =  15.083

-------------+----------------------------------------------------------------
      WHPpain |      Coef.    Std. Err.       t    P>|t|      [95% Con f. Interval]
-------------+----------------------------------------------------------------
      crhtw1 |   -3.425163    1.572465     -2.18    0.030     -6.518312    -.3320138
      crhtw2 |    1.273679    3.738393      0.34    0.734     -6.080003     8.627362
      crhtw3 |    8.362157     3.22885      2.59    0.010      2.010781     14.71353
       _cons |     11.4981    .8358178     13.76    0.000      9.853991     13.14222
-------------+----------------------------------------------------------------

      Source |       SS       df       MS              Number of obs =     362
-------------+------------------------------           F(  3,   358) =    9.60
       Model |  13415.6503     3  4471.88345           Prob > F      =  0.0000
    Residual |    166705.9   358   465.65894           R-squared     =  0.0745
-------------+------------------------------           Adj R-squared =  0.0667
       Total |  180121.551   361  498.951664           Root MSE      =  21.579

-------------+----------------------------------------------------------------
      WHPpain |      Coef.    Std. Err.       t    P>|t|      [95% Con f. Interval]
-------------+----------------------------------------------------------------
      crhtw1 |    -2.44253    1.708858     -1.43    0.154     -5.803192     .9181314
      crhtw2 |   -1.162584    4.560118     -0.25    0.799     -10.13057     7.805402
      crhtw3 |     8.99153    3.993464      2.25    0.025      1.137933     16.84513
       _cons |    17.03202    1.156567     14.73    0.000       14.7575     19.30654
-------------+----------------------------------------------------------------

.
.
. title "perception of risk => the Male Sleep model "


*******************************************************************************
*******************************************************************************
*****                                                                     *****
*****                                                                     *****
*****            perception of risk => the Male Sleep model               *****
*****                                                                     *****
*****                                                                     *****
*****                                            12 Jul 2012    23:56:13  *****
*******************************************************************************
*******************************************************************************


.
.
. *--------------------------------  Male Sleep Model
. /* candiate illnesses  for male sleep model
>
> foreach i in 7 {
> foreach j in 1 5 6 {
> foreach k in 3 {
> foreach m in 2 {
> foreach n in 4 {
> regress whpsleep icdx`i'nr1-icdx`i'nr28 if gender==1
> regress whpsleep icdx`j'nr1-icdx`j'nr18 if gender==1
> regress whpsleep icdx`k'nr1-icdx`k'nr11 if gender==1
> regress whpsleep icdx`m'nr1-icdx`m'nr8 if gender==1
> regress whpsleep icdx`n'nr1-icdx`n'nr20 if gender==1
> }
> }
> }
> }
> }
>
> */
.
. /*
>
> EQ(62) Modelling WHPsleep by OLS-CS
>        The dataset is: /Users/robertyaffee/Documents/data/research/chwk/phase3/
>        data/ox/workingdatasets/MARS/guys.dta
>        The estimation sample is: 1 - 340
>        Dropped 4 observation(s) with missing values from the sample
>
>                   Coefficient  Std.Error      HACSE  t-HACSE  t-prob Part.R^2
> PTSDw1               0.105433    0.03256    0.04230     2.49  0.0132   0.0191
> emplw22              -5.98474      2.470      2.541    -2.36  0.0191   0.0171
> emplw24               49.4794      19.42      4.842     10.2  0.0000   0.2466
> inc1w2               -25.4472      5.061      4.732    -5.38  0.0000   0.0831
> inc2w2               -23.5018      3.534      3.399    -6.91  0.0000   0.1303
> inc3w2               -21.0581      3.592      3.619    -5.82  0.0000   0.0959
> inc4w2               -19.9443      6.400      6.626    -3.01  0.0028   0.0276
> shrelaw1            0.0772105    0.03973    0.04003     1.93  0.0546   0.0115
> shrelaw2           -0.0973671    0.04382    0.04555    -2.14  0.0333   0.0141
> dafter              0.0590987    0.02348   0.008785     6.73  0.0000   0.1242
> dauthw1            -0.0750523    0.03842    0.03704    -2.03  0.0436   0.0127
> dauthw3             0.0758573    0.03144    0.02970     2.55  0.0111   0.0200
> BSIposymp            0.491856    0.04520    0.04746     10.4  0.0000   0.2519
> avgcumdosew1  U       1.43906     0.6365     0.4173     3.45  0.0006   0.0359
> radhw1       U   -0.00976929    0.05730    0.06965   -0.140  0.8885   0.0001
> radhw2       U     0.0495865     0.1175     0.1131    0.438  0.6613   0.0006
> radhw3       U   0.000243982     0.1004    0.07006  0.00348  0.9972   0.0000
>
> sigma                 18.9984  RSS                115139.816
> log-likelihood       -1457.34
> no. of observations       336  no. of parameters          17
> mean(WHPsleep)        17.5028  se(WHPsleep)          24.9014
> When the log-likelihood constant is NOT included:
> AIC                   5.93798  SC                    6.13111
> HQ                    6.01497  FPE                   379.202
> When the log-likelihood constant is included:
> AIC                   8.77586  SC                    8.96899
> HQ                    8.85284  FPE                   6476.56
>
> Normality test:   Chi^2(2)  =   91.495 [0.0000]**
> Hetero test:      F(27,307) =   1.7727 [0.0119]*
> Hetero-X test:    F(82,252) =   1.0160 [0.4528]
> RESET23 test:     F(2,317)  =   3.4531 [0.0328]*
>
> */
.
. *--------Male Sleep crhtw models --------------------------------------------
. cap drop rmslp1

. cap drop H8MaleSLpred

. cap drop h8MSLsepred

. cap drop MSLrsires

. cap drop upbSLm

. cap drop lpbSLm

.
.
. des icdx1nr8 icdx2nr8 icdx3nr7 icdx4nr7 icdx4nr12

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


icdx1nr8          double   %8.0g                        icdx1nr==466 ac bronchitis/brnchial


icdx2nr8          double   %8.0g                        icdx2nr==732 osteochondropathies


icdx3nr7          double   %8.0g                        icdx3nr==434.91 crbrl art ocl nos w infarc


icdx4nr7          double   %8.0g                        icdx4nr==acute myocardial infarct


icdx4nr12         double   %8.0g                        icdx4nr==gastritis/duodenitis


.
. // Full male sleep regression model
. regress WHPsleep age emplw24 inc1w2-inc4w2 ///
>   shrelaw1 shrelaw2 dafter  ///
>   BSIposymp WHPpain illw1 illw3 icdxcnt ///
>   avgcumdosew1 crhtw1-crhtw3 ///
>   icdx1nr8 icdx2nr8 icdx3nr7 icdx4nr7 icdx4nr12  ///
>   if gender==1, vce(cluster id)

Linear regression                                      Number of obs =     337
                                                       F( 20,   336) =       .
                                                       Prob > F      =       .
                                                       R-squared     =  0.4938
                                                       Root MSE      =  18.343

(Std. Err. adjusted for 337 clusters in id)
-------------+----------------------------------------------------------------
             |               Robust
     WHPsleep |      Coef.    Std. Err.       t    P>|t|      [95% Con f. Interval]
-------------+----------------------------------------------------------------
         age |    .0641877    .1091018      0.59    0.557     -.1504208     .2787963
     emplw24 |    6.851287     11.8757      0.58    0.564      -16.5088     30.21137
      inc1w2 |   -13.86611    6.250046     -2.22    0.027     -26.16026     -1.57196
      inc2w2 |   -12.76203    5.392418     -2.37    0.019     -23.36918    -2.154877
      inc3w2 |   -12.47387    5.389799     -2.31    0.021     -23.07587    -1.871866
      inc4w2 |   -12.94876    6.680462     -1.94    0.053     -26.08956     .1920365
    shrelaw1 |    .1014936    .0418903      2.42    0.016      .0190933     .1838939
    shrelaw2 |   -.0929358    .0455514     -2.04    0.042     -.1825376    -.0033339
      dafter |    .0483449    .0100404      4.82    0.000      .0285949     .0680949
   BSIposymp |     .482837     .069221      6.98    0.000       .346676     .6189981
     WHPpain |    .2951184    .0988011      2.99    0.003      .1007718      .489465
       illw1 |    6.565633    2.837541      2.31    0.021        .98405     12.14722
       illw3 |   -4.196016    1.704056     -2.46    0.014     -7.547978    -.8440535
     icdxcnt |    .5958248    .9136271      0.65    0.515     -1.201325     2.392974
avgcumdosew1 |    1.243387    .3978586      3.13    0.002      .4607799     2.025995
      crhtw1 |    .6335098    2.471172      0.26    0.798     -4.227407     5.494427
      crhtw2 |    .3798887    4.392074      0.09    0.931     -8.259537     9.019314
      crhtw3 |    .6355449    2.968463      0.21    0.831     -5.203568     6.474658
    icdx1nr8 |     11.7858    7.732264      1.52    0.128     -3.423943     26.99555
    icdx2nr8 |    2.232118    6.240594      0.36    0.721     -10.04344     14.50768
    icdx3nr7 |    18.45049    8.659509      2.13    0.034      1.416808     35.48417
    icdx4nr7 |    40.61513    6.346595      6.40    0.000      28.13106     53.09919
   icdx4nr12 |    35.71014    5.858498      6.10    0.000      24.18618     47.23409
       _cons |   -14.79354    8.256445     -1.79    0.074     -31.03438     1.447294
-------------+----------------------------------------------------------------

.
.
. // Trimmed male sleep regression model
. cap drop H8MaleSLpred

. cap drop h8MSLsepred

. cap drop MSLrsires

. cap drop upbSLm

. cap drop lpbSLm

. cap drop rmslp1

.
. regress WHPsleep age emplw24 inc1w2-inc4w2 ///
>   shrelaw1 shrelaw2 dafter  ///
>   BSIposymp WHPpain illw1 illw3  ///
>   avgcumdosew1 crhtw1-crhtw3 ///
>    icdx3nr7 icdx4nr7 icdx4nr12  ///
>   if gender==1, vce(cluster id)

Linear regression                                      Number of obs =     337
                                                       F( 17,   336) =       .
                                                       Prob > F      =       .
                                                       R-squared     =  0.4870
                                                       Root MSE      =  18.377

(Std. Err. adjusted for 337 clusters in id)
-------------+----------------------------------------------------------------
             |               Robust
     WHPsleep |      Coef.    Std. Err.       t    P>|t|      [95% Con f. Interval]
-------------+----------------------------------------------------------------
         age |    .0685693    .1034704      0.66    0.508      -.134962     .2721006
     emplw24 |    18.32908    8.414051      2.18    0.030      1.778224     34.87993
      inc1w2 |   -13.68522    6.176679     -2.22    0.027     -25.83505    -1.535386
      inc2w2 |   -13.09289    5.403288     -2.42    0.016     -23.72143    -2.464356
      inc3w2 |   -12.67035    5.406765     -2.34    0.020     -23.30573    -2.034979
      inc4w2 |   -12.77384    6.679958     -1.91    0.057     -25.91364      .365971
    shrelaw1 |    .1092769    .0419831      2.60    0.010      .0266941     .1918598
    shrelaw2 |   -.0955916    .0449321     -2.13    0.034     -.1839752     -.007208
      dafter |      .05152    .0070831      7.27    0.000      .0375873     .0654527
   BSIposymp |    .4932165    .0723817      6.81    0.000      .3508381      .635595
     WHPpain |    .2811363    .0987587      2.85    0.005       .086873     .4753996
       illw1 |    7.225872    2.783662      2.60    0.010      1.750271     12.70147
       illw3 |   -3.669956    1.542068     -2.38    0.018     -6.703281    -.6366316
avgcumdosew1 |    1.277645    .3456958      3.70    0.000       .597644     1.957646
      crhtw1 |     1.05903    2.451491      0.43    0.666     -3.763175     5.881234
      crhtw2 |   -.5059819    4.337534     -0.12    0.907     -9.038126     8.026162
      crhtw3 |    1.319219    2.981739      0.44    0.658     -4.546009     7.184448
    icdx3nr7 |    18.04353    8.510413      2.12    0.035      1.303132     34.78394
    icdx4nr7 |    41.93536    5.890901      7.12    0.000      30.34766     53.52305
   icdx4nr12 |    36.07907    5.513198      6.54    0.000      25.23433      46.9238
       _cons |    -14.2854    8.200349     -1.74    0.082      -30.4159     1.845086
-------------+----------------------------------------------------------------

.
. *-+ no male sleep relationships between health threat sleep
.
. // prep for residual and prediction interval analysis
.     predict H8MaleSLpred if gender==2, xb
(348 missing values generated)

.      predict h8MSLsepred if gender==2, stdp
(348 missing values generated)

.      predict MSLrsires if gender==2, residual
(348 missing values generated)

.      gen upbSLm = H8MaleSLpred + 1.96*h8MSLsepred
(348 missing values generated)

.      gen lpbSLm = H8MaleSLpred - 1.96*h8MSLsepred
(348 missing values generated)

.  // residual diagnostics
.
. rdiagrbst rmslp1
(option xb assumed; fitted values)
(11 missing values generated)
yvar =  WHPsleep
age emplw24 inc1w2 inc2w2 inc3w2 inc4w2 shrelaw1 shrelaw2 dafter BSIposymp WHPpain illw1 illw3 avgcumdosew1 crhtw1 crhtw2 crhtw3 icdx3nr7 icdx4nr7 icdx4nr12 _cons
162
156
xvar = ,age emplw24 inc1w2 inc2w2 inc3w2 inc4w2 shrelaw1 shrelaw2 dafter BSIposymp WHPpain illw1 illw3 avgcumdosew1 crhtw1 crhtw2 crhtw3 icdx3nr7 icdx4nr7 icdx4nr12

   | 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%    -3.375236       -9.71623
 5%    -1.434807      -6.597977
10%     1.560467      -5.321205       Obs                 691
25%     7.232412      -4.890507       Sum of Wgt.         691

50%     18.11478                      Mean           22.53079
                        Largest       Std. Dev.      20.35484
75%     33.27406        94.9482
90%     50.52992       102.6599       Variance       414.3194
95%     63.19973       108.7799       Skewness       1.207259
99%     85.01624       114.9269       Kurtosis       4.556202
(11 missing values generated)
(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:   rmslp1std |
  +------------------+
  | 615.   -1.584244 |
  |  16.   -1.431049 |
  | 217.   -1.368323 |
  | 460.   -1.347164 |
  | 646.   -1.344078 |
  +------------------+

  +-----------------+
  | 288.   3.268812 |
  | 344.    3.55775 |
  |  11.   3.936613 |
  | 335.   4.237278 |
  | 178.   4.539272 |
  +-----------------+
There are 4 positive outliers

     +----------------+
     |  id   rmslp1~d |
     +----------------+
 11. | 152   3.936613 |
178. | 126   4.539272 |
335. | 133   4.237278 |
344. | 452    3.55775 |
     +----------------+
There are 6 negative outliers:

     +-----------------+
     |  id   rmslp1std |
     +-----------------+
 16. | 122   -1.431049 |
217. | 110   -1.368323 |
460. | 598   -1.347164 |
570. | 506   -1.308082 |
615. |   8   -1.584244 |
     +-----------------+
646. |   6   -1.344078 |
     +-----------------+
(702 real changes made, 702 to missing)
    4
    6
(15 real changes made)
   21

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

50%            0                      Mean           .0213675
                        Largest       Std. Dev.      .1447093
75%            0              1
90%            0              1       Variance       .0209408
95%            0              1       Skewness       6.619806
99%            1              1       Kurtosis       44.82183
(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 rmslp1std

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

 -15* | 8
 -14* | 3
 -13* | 7541
 -12* | 7655443332222222111100000
 -11* | 99887654433331111100000
 -10* | 999987777766433332222221110000000
  -9* | 9999988777777666655554444433333222222211000
  -8* | 999877777666644444332221111110
  -7* | 9888877766665555444433222222111111110000000
  -6* | 998877666666555443333322100000
  -5* | 99999977766555533321100
  -4* | 999988888777766666655444433333222111000000
  -3* | 99877777666555544333322111110000
  -2* | 9988766555544322111000
  -1* | 998888887777766544333211110000
  -0* | 99877665444322211
   0* | 1123344455556666666777788999999
   1* | 000001111233333355567777788899
   2* | 00000011333456666677899
   3* | 00112224444556777799
   4* | 01123333566
   5* | 1122345557788899
   6* | 022333444556677778999
   7* | 013344566888
   8* | 03344445567777889
   9* | 0457889
  10* | 014455778
  11* | 22334579
  12* | 0245566777
  13* | 12233568
  14* | 011266
  15* | 33368
  16* | 02567
  17* | 0224688
  18* | 01357788
  19* | 034
  20* | 045
  21* | 2556
  22* | 4467
  23* | 289
  24* | 57
  25* | 056
  26* | 28
  27* | 111
  28* | 13
  29* | 2
  30* | 37
  31* | 6
  32* | 7
  33* |
  34* |
  35* | 6
  36* |
  37* |
  38* |
  39* | 4
  40* |
  41* |
  42* | 4
  43* |
  44* |
  45* | 4


Computing the leverage h

                          rmslp1std
-------------------------------------------------------------
      Percentiles      Smallest
 1%    -1.272721      -1.584244
 5%    -1.177391      -1.431049
10%    -1.030238      -1.368323       Obs                 691
25%    -.7515843      -1.347164       Sum of Wgt.         691

50%    -.2169512                      Mean           5.01e-10
                        Largest       Std. Dev.             1
75%     .5277992        3.55775
90%     1.375552       3.936613       Variance              1
95%     1.997999       4.237278       Skewness       1.207259
99%     3.069808       4.539272       Kurtosis       4.556202

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

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

     +-----+
     |  id |
     +-----+
 11. | 152 |
178. | 126 |
335. | 133 |
344. | 452 |
     +-----+
    4
Negative outliers are:

     +-----+
     |  id |
     +-----+
 16. | 122 |
217. | 110 |
460. | 598 |
570. | 506 |
615. |   8 |
     +-----+
646. |   6 |
     +-----+
    6

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

    Variable |       VIF       1/VIF
-------------+----------------------
      crhtw2 |     18.48    0.054118
      crhtw3 |     13.89    0.072003
      inc2w2 |      5.62    0.177934
      inc3w2 |      5.29    0.188912
      crhtw1 |      4.17    0.239696
      inc1w2 |      2.28    0.438400
    shrelaw2 |      2.23    0.449395
    shrelaw1 |      2.06    0.484809
      inc4w2 |      1.76    0.569178
   BSIposymp |      1.69    0.591963
       illw3 |      1.68    0.596771
     WHPpain |      1.63    0.612886
         age |      1.39    0.719163
       illw1 |      1.20    0.830233
     emplw24 |      1.17    0.855084
    icdx3nr7 |      1.10    0.908253
      dafter |      1.09    0.918950
   icdx4nr12 |      1.08    0.924193
    icdx4nr7 |      1.07    0.931910
avgcumdosew1 |      1.06    0.941070
-------------+----------------------
    Mean VIF |      3.50

Condition number using scaled variables =        21.50

Condition Indexes and Variance-Decomposition Proportions

condition
   index        _cons          age      emplw24       inc1w2       inc2w2       inc3w2       inc4w2     shrelaw1     shrelaw2       dafter
1   1.00         0.00         0.00         0.00         0.00         0.00         0.00         0.00         0.00         0.00         0.00
2   1.41         0.00         0.00         0.00         0.00         0.00         0.00         0.00         0.00         0.00         0.00
3   2.06         0.00         0.00         0.04         0.01         0.00         0.01         0.02         0.01         0.01         0.11
4   2.19         0.00         0.00         0.04         0.13         0.00         0.00         0.00         0.00         0.00         0.02
5   2.32         0.00         0.00         0.39         0.00         0.01         0.00         0.02         0.00         0.00         0.02
6   2.36         0.00         0.00         0.11         0.00         0.01         0.00         0.02         0.00         0.00         0.21
7   2.37         0.00         0.00         0.00         0.00         0.00         0.03         0.29         0.00         0.00         0.00
8   2.41         0.00         0.00         0.01         0.00         0.00         0.00         0.03         0.00         0.00         0.07
9   2.53         0.00         0.00         0.02         0.02         0.01         0.00         0.07         0.00         0.00         0.02
10  2.59         0.00         0.00         0.00         0.03         0.01         0.02         0.00         0.02         0.01         0.02
11  2.73         0.00         0.00         0.00         0.11         0.00         0.00         0.01         0.02         0.01         0.16
12  2.74         0.00         0.00         0.05         0.08         0.00         0.00         0.01         0.02         0.03         0.22
13  2.84         0.00         0.00         0.09         0.00         0.00         0.01         0.03         0.02         0.03         0.03
14  3.02         0.00         0.00         0.00         0.01         0.01         0.00         0.01         0.07         0.03         0.09
15  4.01         0.00         0.00         0.21         0.00         0.00         0.01         0.00         0.01         0.00         0.00
16  5.75         0.00         0.00         0.00         0.02         0.00         0.00         0.00         0.35         0.32         0.01
17  5.91         0.00         0.00         0.00         0.00         0.00         0.00         0.00         0.46         0.53         0.00
18 10.49         0.00         0.17         0.00         0.30         0.54         0.56         0.33         0.00         0.00         0.01
19 11.55         0.00         0.38         0.01         0.00         0.00         0.00         0.00         0.01         0.00         0.00
20 13.69         0.00         0.00         0.00         0.00         0.00         0.00         0.00         0.00         0.00         0.01
21 21.50         0.99         0.44         0.02         0.28         0.40         0.35         0.13         0.01         0.01         0.00

condition
   index    BSIposymp      WHPpain        illw1        illw3 avgcumdosew1       crhtw1       crhtw2       crhtw3     icdx3nr7     icdx4nr7
1   1.00         0.00         0.01         0.00         0.01         0.00         0.00         0.00         0.00         0.00         0.00
2   1.41         0.00         0.00         0.00         0.00         0.00         0.02         0.01         0.01         0.01         0.00
3   2.06         0.00         0.02         0.08         0.03         0.01         0.01         0.00         0.00         0.05         0.00
4   2.19         0.00         0.00         0.02         0.00         0.01         0.00         0.00         0.00         0.00         0.00
5   2.32         0.00         0.01         0.05         0.01         0.09         0.00         0.00         0.00         0.03         0.02
6   2.36         0.00         0.00         0.01         0.02         0.05         0.00         0.00         0.00         0.06         0.22
7   2.37         0.00         0.00         0.00         0.00         0.04         0.00         0.00         0.00         0.06         0.00
8   2.41         0.00         0.00         0.01         0.00         0.06         0.00         0.00         0.00         0.00         0.60
9   2.53         0.00         0.01         0.14         0.01         0.08         0.00         0.00         0.00         0.35         0.02
10  2.59         0.00         0.02         0.01         0.01         0.40         0.00         0.00         0.00         0.10         0.00
11  2.73         0.00         0.00         0.02         0.04         0.00         0.01         0.00         0.00         0.14         0.00
12  2.74         0.00         0.02         0.02         0.03         0.21         0.00         0.00         0.00         0.05         0.03
13  2.84         0.00         0.03         0.47         0.01         0.00         0.00         0.00         0.00         0.08         0.03
14  3.02         0.00         0.00         0.04         0.23         0.00         0.03         0.00         0.00         0.04         0.03
15  4.01         0.00         0.79         0.00         0.21         0.00         0.00         0.00         0.00         0.00         0.00
16  5.75         0.00         0.00         0.02         0.17         0.00         0.45         0.01         0.05         0.00         0.00
17  5.91         0.00         0.02         0.00         0.17         0.04         0.28         0.01         0.04         0.02         0.00
18 10.49         0.15         0.03         0.03         0.03         0.00         0.00         0.00         0.01         0.00         0.04
19 11.55         0.68         0.02         0.01         0.00         0.00         0.00         0.00         0.01         0.01         0.01
20 13.69         0.01         0.01         0.01         0.00         0.00         0.19         0.96         0.87         0.00         0.00
21 21.50         0.15         0.01         0.04         0.02         0.00         0.01         0.00         0.00         0.00         0.01

condition
   index    icdx4nr12
1   1.00         0.00
2   1.41         0.00
3   2.06         0.03
4   2.19         0.30
5   2.32         0.01
6   2.36         0.00
7   2.37         0.01
8   2.41         0.00
9   2.53         0.05
10  2.59         0.07
11  2.73         0.36
12  2.74         0.03
13  2.84         0.07
14  3.02         0.00
15  4.01         0.03
16  5.75         0.00
17  5.91         0.00
18 10.49         0.00
19 11.55         0.02
20 13.69         0.00
21 21.50         0.00

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

                   Shapiro-Wilk W test for normal data

    Variable |    Obs       W           V         z       Prob>z
-------------+--------------------------------------------------
   rmslp1std |    691    0.90888     41.065     9.059    0.00000

                    Skewness/Kurtosis tests for Normality
                                                         ------- joint ------
    Variable |    Obs   Pr(Skewness)   Pr(Kurtosis)  adj chi2(2)    Prob>chi2
-------------+---------------------------------------------------------------
   rmslp1std |    691      0.0000         0.0000            .         0.0000
Jarque-Bera normality test:  237.6 Chi(2)  2.6e-52
Jarque-Bera test for Ho: normality:
(file qnrmslp1std.gph saved)
(file bxrmslp1std.gph saved)
(bin=26, start=-1.5842435, width=.23551984)
(file histrmslp1std.gph saved)
(file kdrmslp1std.gph saved)
(file combi4.gph saved)
-------------
Specification error tests:

Ramsey RESET test using powers of the fitted values of WHPsleep
       Ho:  model has no omitted variables
                 F(3, 313) =      2.56
                  Prob > F =      0.0548
-------------

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   |            4   |
                  +--------------------+--------------+
                  | Negative-outliers   |            6   |
                  +--------------------+--------------+
                  | Bad-outliers   |            1   |
                  +--------------------+--------------+
                  | collinearity   |            1   |
                  +--------------------+--------------+
                  | Shapiro-Wilk   |            1   |
                  +--------------------+--------------+
                  | Kolmogorov-Smirnov   |            1   |
                  +--------------------+--------------+
                  | Specification   |            0   |
                  +--------------------+--------------+
                  | StructuralBreak   |            0   |
                  +--------------------+--------------+
                  | TotalProbs   |            3   |
                  +--------------------+--------------+
                  | Percent-Invalid   |           60   |
                  +--------------------+--------------+
                  | Pct-Pseudo-Valid   |           40   |
                  +--------------------+--------------+

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
-------------

.
. // prediction interval
. scatter H8MaleSLpred MSLrsires || lowess H8MaleSLpred  MSLrsires  ///
>   || lowess upbSLm MSLrsires || lowess lpbSLm MSLrsires, ///
>   title(Prediction interval of Male Sleep model)   ///
>   ytitle(Predicted Male Sleep )

.
. // saving graph
. gr save MSLPredInt.gph, replace
(file MSLPredInt.gph saved)

. gr export MSLPredInt.eps, replace
(file MSLPredInt.eps written in EPS format)

. gr use  MSLPredInt.gph

.
.
. //  No signinficant main effect for risk sleep relationship so no moderation
. //  tests performed...
.
. title "tests for possible indirect effects for male sleep model"


*******************************************************************************
*******************************************************************************
*****                                                                     *****
*****                                                                     *****
*****      tests for possible indirect effects for male sleep model       *****
*****                                                                     *****
*****                                                                     *****
*****                                            12 Jul 2012    23:56:27  *****
*******************************************************************************
*******************************************************************************


. set more off

. foreach var in age emplw24 inc1w2 inc3w2 inc4w2 ///
>   shrelaw1 shrelaw2 dafter  ///
>   BSIposymp WHPpain illw1 illw3  ///
>   avgcumdosew1 icdx3nr7 icdx4nr7 icdx4nr12  {
  2. sem (crhtw1-> `var')(`var'-> whpsleep) if gender==1, nocapslatent iterate(50)
  3. sem (crhtw2-> `var')(`var'-> whpsleep) if gender==1, nocapslatent iterate(50)
  4. sem (crhtw3-> `var')(`var'-> whpsleep) if gender==1, nocapslatent iterate(50)
  5. }

Endogenous variables


Observed:  age whpsleep


Exogenous variables


Observed:  crhtw1


Fitting target model:

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

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

--------------+----------------------------------------------------------------
              |                 OIM
              |      Coef.    Std. Err.       z    P>|z|      [95% Con f. Interval]
--------------+----------------------------------------------------------------
 Structural    |
   age <-      |
       crhtw1 |    4.266795    .6763452      6.31    0.000      2.941183     5.592408
        _cons |    49.78043     .635245     78.36    0.000      48.53537     51.02549
  ------------+----------------------------------------------------------------
   whpsleep <- |
          age |    .5911349    .1054724      5.60    0.000      .3844128      .797857
        _cons |   -11.72059    5.344283     -2.19    0.028     -22.19519    -1.245984
--------------+----------------------------------------------------------------
 Variance      |
        e.age |    133.6667    10.26687                    114.9854      155.383
   e.whpsleep |    563.2601    43.26375                    484.5389     654.7709
--------------+----------------------------------------------------------------
LR test of model vs. saturated: chi2(1)   =     19.43, Prob > chi2 = 0.0000

Endogenous variables


Observed:  age whpsleep


Exogenous variables


Observed:  crhtw2


Fitting target model:

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

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

--------------+----------------------------------------------------------------
              |                 OIM
              |      Coef.    Std. Err.       z    P>|z|      [95% Con f. Interval]
--------------+----------------------------------------------------------------
 Structural    |
   age <-      |
       crhtw2 |    4.951499    .6691348      7.40    0.000      3.640019     6.262979
        _cons |    50.08515    .6280741     79.74    0.000      48.85415     51.31616
  ------------+----------------------------------------------------------------
   whpsleep <- |
          age |    .5911349    .1054724      5.60    0.000      .3844128      .797857
        _cons |   -11.72059    5.344283     -2.19    0.028     -22.19519    -1.245984
--------------+----------------------------------------------------------------
 Variance      |
        e.age |    128.5885    9.876824                     110.617     149.4798
   e.whpsleep |    563.2601    43.26375                    484.5389     654.7709
--------------+----------------------------------------------------------------
LR test of model vs. saturated: chi2(1)   =     36.29, Prob > chi2 = 0.0000

Endogenous variables


Observed:  age whpsleep


Exogenous variables


Observed:  crhtw3


Fitting target model:

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

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

--------------+----------------------------------------------------------------
              |                 OIM
              |      Coef.    Std. Err.       z    P>|z|      [95% Con f. Interval]
--------------+----------------------------------------------------------------
 Structural    |
   age <-      |
       crhtw3 |    4.882082    .6713623      7.27    0.000      3.566236     6.197928
        _cons |    50.07363     .629633     79.53    0.000      48.83957     51.30769
  ------------+----------------------------------------------------------------
   whpsleep <- |
          age |    .5911349    .1054724      5.60    0.000      .3844128      .797857
        _cons |   -11.72059    5.344283     -2.19    0.028     -22.19519    -1.245984
--------------+----------------------------------------------------------------
 Variance      |
        e.age |    129.2045    9.924137                    111.1469     150.1959
   e.whpsleep |    563.2601    43.26375                    484.5389     654.7709
--------------+----------------------------------------------------------------
LR test of model vs. saturated: chi2(1)   =     36.18, Prob > chi2 = 0.0000

Endogenous variables


Observed:  emplw24 whpsleep


Exogenous variables


Observed:  crhtw1


Fitting target model:

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

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

--------------+----------------------------------------------------------------
              |                 OIM
              |      Coef.    Std. Err.       z    P>|z|      [95% Con f. Interval]
--------------+----------------------------------------------------------------
 Structural    |
   emplw24 <-  |
       crhtw1 |   -.0024245    .0031699     -0.76    0.444     -.0086373     .0037883
        _cons |    .0026053    .0029772      0.88    0.382       -.00323     .0084406
  ------------+----------------------------------------------------------------
   whpsleep <- |
      emplw24 |    60.46045    24.62702      2.46    0.014      12.19238     108.7285
        _cons |    17.16956    1.337556     12.84    0.000        14.548     19.79112
--------------+----------------------------------------------------------------
 Variance      |
    e.emplw24 |    .0029361    .0002255                    .0025257     .0034131
   e.whpsleep |     604.701     46.4468                    520.1879     702.9445
--------------+----------------------------------------------------------------
LR test of model vs. saturated: chi2(1)   =     35.79, Prob > chi2 = 0.0000

Endogenous variables


Observed:  emplw24 whpsleep


Exogenous variables


Observed:  crhtw2


Fitting target model:

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

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

--------------+----------------------------------------------------------------
              |                 OIM
              |      Coef.    Std. Err.       z    P>|z|      [95% Con f. Interval]
--------------+----------------------------------------------------------------
 Structural    |
   emplw24 <-  |
       crhtw2 |   -.0026848    .0031968     -0.84    0.401     -.0089505     .0035809
        _cons |    .0024558    .0030007      0.82    0.413     -.0034254      .008337
  ------------+----------------------------------------------------------------
   whpsleep <- |
      emplw24 |    60.46045    24.62702      2.46    0.014      12.19238     108.7285
        _cons |    17.16956    1.337556     12.84    0.000        14.548     19.79112
--------------+----------------------------------------------------------------
 Variance      |
    e.emplw24 |     .002935    .0002254                    .0025248     .0034119
   e.whpsleep |     604.701     46.4468                    520.1879     702.9445
--------------+----------------------------------------------------------------
LR test of model vs. saturated: chi2(1)   =     59.55, Prob > chi2 = 0.0000

Endogenous variables


Observed:  emplw24 whpsleep


Exogenous variables


Observed:  crhtw3


Fitting target model:

Iteration 0:   log likelihood = -1511.8786
Iteration 1:   log likelihood = -1511.8786  (backed up)

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

--------------+----------------------------------------------------------------
              |                 OIM
              |      Coef.    Std. Err.       z    P>|z|      [95% Con f. Interval]
--------------+----------------------------------------------------------------
 Structural    |
   emplw24 <-  |
       crhtw3 |   -.0026766    .0031998     -0.84    0.403     -.0089482      .003595
        _cons |    .0024567     .003001      0.82    0.413     -.0034251     .0083384
  ------------+----------------------------------------------------------------
   whpsleep <- |
      emplw24 |    60.46045    24.62702      2.46    0.014      12.19238     108.7285
        _cons |    17.16956    1.337556     12.84    0.000        14.548     19.79112
--------------+----------------------------------------------------------------
 Variance      |
    e.emplw24 |    .0029351    .0002254                    .0025249     .0034119
   e.whpsleep |     604.701     46.4468                    520.1879     702.9445
--------------+----------------------------------------------------------------
LR test of model vs. saturated: chi2(1)   =     59.10, Prob > chi2 = 0.0000

Endogenous variables


Observed:  inc1w2 whpsleep


Exogenous variables


Observed:  crhtw1


Fitting target model:

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

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

--------------+----------------------------------------------------------------
              |                 OIM
              |      Coef.    Std. Err.       z    P>|z|      [95% Con f. Interval]
--------------+----------------------------------------------------------------
 Structural    |
   inc1w2 <-   |
       crhtw1 |    .0027897     .014711      0.19    0.850     -.0260432     .0316227
        _cons |    .0682431     .013817      4.94    0.000      .0411622     .0953239
  ------------+----------------------------------------------------------------
   whpsleep <- |
       inc1w2 |    -.039732    5.357837     -0.01    0.994      -10.5409     10.46144
        _cons |     17.3506    1.395576     12.43    0.000      14.61532     20.08588
--------------+----------------------------------------------------------------
 Variance      |
     e.inc1w2 |    .0632367    .0048572                    .0543988     .0735106
   e.whpsleep |    615.4521    47.27259                    529.4365     715.4424
--------------+----------------------------------------------------------------
LR test of model vs. saturated: chi2(1)   =     33.79, Prob > chi2 = 0.0000

Endogenous variables


Observed:  inc1w2 whpsleep


Exogenous variables


Observed:  crhtw2


Fitting target model:

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

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

--------------+----------------------------------------------------------------
              |                 OIM
              |      Coef.    Std. Err.       z    P>|z|      [95% Con f. Interval]
--------------+----------------------------------------------------------------
 Structural    |
   inc1w2 <-   |
       crhtw2 |    .0251521    .0147765      1.70    0.089     -.0038094     .0541135
        _cons |    .0724748    .0138698      5.23    0.000      .0452905      .099659
  ------------+----------------------------------------------------------------
   whpsleep <- |
       inc1w2 |    -.039732    5.357837     -0.01    0.994      -10.5409     10.46144
        _cons |     17.3506    1.395576     12.43    0.000      14.61532     20.08588
--------------+----------------------------------------------------------------
 Variance      |
     e.inc1w2 |    .0627075    .0048165                    .0539435     .0728954
   e.whpsleep |    615.4521    47.27259                    529.4365     715.4424
--------------+----------------------------------------------------------------
LR test of model vs. saturated: chi2(1)   =     56.92, Prob > chi2 = 0.0000

Endogenous variables


Observed:  inc1w2 whpsleep


Exogenous variables


Observed:  crhtw3


Fitting target model:

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

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

--------------+----------------------------------------------------------------
              |                 OIM
              |      Coef.    Std. Err.       z    P>|z|      [95% Con f. Interval]
--------------+----------------------------------------------------------------
 Structural    |
   inc1w2 <-   |
       crhtw3 |    .0256598    .0147879      1.74    0.083     -.0033239     .0546435
        _cons |    .0725748    .0138687      5.23    0.000      .0453926      .099757
  ------------+----------------------------------------------------------------
   whpsleep <- |
       inc1w2 |    -.039732    5.357837     -0.01    0.994      -10.5409     10.46144
        _cons |     17.3506    1.395576     12.43    0.000      14.61532     20.08588
--------------+----------------------------------------------------------------
 Variance      |
     e.inc1w2 |    .0626867    .0048149                    .0539256     .0728712
   e.whpsleep |    615.4521    47.27259                    529.4365     715.4424
--------------+----------------------------------------------------------------
LR test of model vs. saturated: chi2(1)   =     56.51, Prob > chi2 = 0.0000

Endogenous variables


Observed:  inc3w2 whpsleep


Exogenous variables


Observed:  crhtw1


Fitting target model:

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

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

--------------+----------------------------------------------------------------
              |                 OIM
              |      Coef.    Std. Err.       z    P>|z|      [95% Con f. Interval]
--------------+----------------------------------------------------------------
 Structural    |
   inc3w2 <-   |
       crhtw1 |   -.0125769    .0279666     -0.45    0.653     -.0673905     .0422367
        _cons |    .3521949    .0262672     13.41    0.000      .3007122     .4036776
  ------------+----------------------------------------------------------------
   whpsleep <- |
       inc3w2 |   -2.930759    2.813133     -1.04    0.297     -8.444398     2.582879
        _cons |    18.38534    1.673713     10.98    0.000      15.10493     21.66576
--------------+----------------------------------------------------------------
 Variance      |
     e.inc3w2 |    .2285425    .0175542                    .1966014      .265673
   e.whpsleep |     613.488    47.12173                    527.7469     713.1592
--------------+----------------------------------------------------------------
LR test of model vs. saturated: chi2(1)   =     33.61, Prob > chi2 = 0.0000

Endogenous variables


Observed:  inc3w2 whpsleep


Exogenous variables


Observed:  crhtw2


Fitting target model:

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

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

--------------+----------------------------------------------------------------
              |                 OIM
              |      Coef.    Std. Err.       z    P>|z|      [95% Con f. Interval]
--------------+----------------------------------------------------------------
 Structural    |
   inc3w2 <-   |
       crhtw2 |   -.0284892    .0281755     -1.01    0.312     -.0837122     .0267337
        _cons |    .3487401    .0264465     13.19    0.000      .2969058     .4005743
  ------------+----------------------------------------------------------------
   whpsleep <- |
       inc3w2 |   -2.930759    2.813133     -1.04    0.297     -8.444398     2.582879
        _cons |    18.38534    1.673713     10.98    0.000      15.10493     21.66576
--------------+----------------------------------------------------------------
 Variance      |
     e.inc3w2 |    .2279912    .0175119                    .1961272     .2650321
   e.whpsleep |     613.488    47.12173                    527.7469     713.1592
--------------+----------------------------------------------------------------
LR test of model vs. saturated: chi2(1)   =     55.79, Prob > chi2 = 0.0000

Endogenous variables


Observed:  inc3w2 whpsleep


Exogenous variables


Observed:  crhtw3


Fitting target model:

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

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

--------------+----------------------------------------------------------------
              |                 OIM
              |      Coef.    Std. Err.       z    P>|z|      [95% Con f. Interval]
--------------+----------------------------------------------------------------
 Structural    |
   inc3w2 <-   |
       crhtw3 |   -.0279662    .0282034     -0.99    0.321     -.0832439     .0273116
        _cons |    .3488292    .0264504     13.19    0.000      .2969873     .4006711
  ------------+----------------------------------------------------------------
   whpsleep <- |
       inc3w2 |   -2.930759    2.813133     -1.04    0.297     -8.444398     2.582879
        _cons |    18.38534    1.673713     10.98    0.000      15.10493     21.66576
--------------+----------------------------------------------------------------
 Variance      |
     e.inc3w2 |    .2280175    .0175139                    .1961497     .2650627
   e.whpsleep |     613.488    47.12173                    527.7469     713.1592
--------------+----------------------------------------------------------------
LR test of model vs. saturated: chi2(1)   =     55.38, Prob > chi2 = 0.0000

Endogenous variables


Observed:  inc4w2 whpsleep


Exogenous variables


Observed:  crhtw1


Fitting target model:

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

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

--------------+----------------------------------------------------------------
              |                 OIM
              |      Coef.    Std. Err.       z    P>|z|      [95% Con f. Interval]
--------------+----------------------------------------------------------------
 Structural    |
   inc4w2 <-   |
       crhtw1 |    .0055752      .01123      0.50    0.620     -.0164352     .0275856
        _cons |    .0391404    .0105476      3.71    0.000      .0184676     .0598133
  ------------+----------------------------------------------------------------
   whpsleep <- |
       inc4w2 |   -2.773495    7.014819     -0.40    0.693     -16.52229      10.9753
        _cons |    17.45426    1.373689     12.71    0.000      14.76188     20.14664
--------------+----------------------------------------------------------------
 Variance      |
     e.inc4w2 |    .0368507    .0028305                    .0317005     .0428377
   e.whpsleep |    615.1686    47.25081                    529.1926     715.1128
--------------+----------------------------------------------------------------
LR test of model vs. saturated: chi2(1)   =     33.97, Prob > chi2 = 0.0000

Endogenous variables


Observed:  inc4w2 whpsleep


Exogenous variables


Observed:  crhtw2


Fitting target model:

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

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

--------------+----------------------------------------------------------------
              |                 OIM
              |      Coef.    Std. Err.       z    P>|z|      [95% Con f. Interval]
--------------+----------------------------------------------------------------
 Structural    |
   inc4w2 <-   |
       crhtw2 |   -.0072199    .0113249     -0.64    0.524     -.0294162     .0149764
        _cons |    .0370196    .0106299      3.48    0.000      .0161853     .0578538
  ------------+----------------------------------------------------------------
   whpsleep <- |
       inc4w2 |   -2.773495    7.014819     -0.40    0.693     -16.52229      10.9753
        _cons |    17.45426    1.373689     12.71    0.000      14.76188     20.14664
--------------+----------------------------------------------------------------
 Variance      |
     e.inc4w2 |    .0368333    .0028292                    .0316855     .0428175
   e.whpsleep |    615.1686    47.25081                    529.1926     715.1128
--------------+----------------------------------------------------------------
LR test of model vs. saturated: chi2(1)   =     56.25, Prob > chi2 = 0.0000

Endogenous variables


Observed:  inc4w2 whpsleep


Exogenous variables


Observed:  crhtw3


Fitting target model:

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

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

--------------+----------------------------------------------------------------
              |                 OIM
              |      Coef.    Std. Err.       z    P>|z|      [95% Con f. Interval]
--------------+----------------------------------------------------------------
 Structural    |
   inc4w2 <-   |
       crhtw3 |   -.0109769    .0113266     -0.97    0.332     -.0331765     .0112228
        _cons |    .0363254    .0106225      3.42    0.001      .0155056     .0571453
  ------------+----------------------------------------------------------------
   whpsleep <- |
       inc4w2 |   -2.773495    7.014819     -0.40    0.693     -16.52229      10.9753
        _cons |    17.45426    1.373689     12.71    0.000      14.76188     20.14664
--------------+----------------------------------------------------------------
 Variance      |
     e.inc4w2 |    .0367756    .0028247                    .0316359     .0427504
   e.whpsleep |    615.1686    47.25081                    529.1926     715.1128
--------------+----------------------------------------------------------------
LR test of model vs. saturated: chi2(1)   =     55.80, Prob > chi2 = 0.0000

Endogenous variables


Observed:  shrelaw1 whpsleep


Exogenous variables


Observed:  crhtw1


Fitting target model:

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

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

--------------+----------------------------------------------------------------
              |                 OIM
              |      Coef.    Std. Err.       z    P>|z|      [95% Con f. Interval]
--------------+----------------------------------------------------------------
 Structural    |
   shrelaw1 <- |
       crhtw1 |     6.13074     2.17532      2.82    0.005       1.86719     10.39429
        _cons |    30.47306     2.04313     14.91    0.000       26.4686     34.47752
  ------------+----------------------------------------------------------------
   whpsleep <- |
     shrelaw1 |     .145395    .0349367      4.16    0.000      .0769204     .2138697
        _cons |    13.04396    1.672362      7.80    0.000      9.766188     16.32172
--------------+----------------------------------------------------------------
 Variance      |
   e.shrelaw1 |    1382.715    106.2057                    1189.467      1607.36
   e.whpsleep |    585.5372    44.97484                    503.7024     680.6673
--------------+----------------------------------------------------------------
LR test of model vs. saturated: chi2(1)   =     28.69, Prob > chi2 = 0.0000

Endogenous variables


Observed:  shrelaw1 whpsleep


Exogenous variables


Observed:  crhtw2


Fitting target model:

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

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

--------------+----------------------------------------------------------------
              |                 OIM
              |      Coef.    Std. Err.       z    P>|z|      [95% Con f. Interval]
--------------+----------------------------------------------------------------
 Structural    |
   shrelaw1 <- |
       crhtw2 |    7.651572    2.180522      3.51    0.000      3.377827     11.92532
        _cons |    31.00972    2.046717     15.15    0.000      26.99822     35.02121
  ------------+----------------------------------------------------------------
   whpsleep <- |
     shrelaw1 |     .145395    .0349367      4.16    0.000      .0769204     .2138697
        _cons |    13.04396    1.672362      7.80    0.000      9.766188     16.32172
--------------+----------------------------------------------------------------
 Variance      |
   e.shrelaw1 |    1365.513    104.8845                    1174.669     1587.363
   e.whpsleep |    585.5372    44.97484                    503.7024     680.6673
--------------+----------------------------------------------------------------
LR test of model vs. saturated: chi2(1)   =     48.60, Prob > chi2 = 0.0000

Endogenous variables


Observed:  shrelaw1 whpsleep


Exogenous variables


Observed:  crhtw3


Fitting target model:

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

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

--------------+----------------------------------------------------------------
              |                 OIM
              |      Coef.    Std. Err.       z    P>|z|      [95% Con f. Interval]
--------------+----------------------------------------------------------------
 Structural    |
   shrelaw1 <- |
       crhtw3 |    7.476168    2.184426      3.42    0.001      3.194771     11.75756
        _cons |    30.97935    2.048651     15.12    0.000      26.96407     34.99464
  ------------+----------------------------------------------------------------
   whpsleep <- |
     shrelaw1 |     .145395    .0349367      4.16    0.000      .0769204     .2138697
        _cons |    13.04396    1.672362      7.80    0.000      9.766188     16.32172
--------------+----------------------------------------------------------------
 Variance      |
   e.shrelaw1 |     1367.85    105.0639                    1176.679     1590.079
   e.whpsleep |    585.5372    44.97484                    503.7024     680.6673
--------------+----------------------------------------------------------------
LR test of model vs. saturated: chi2(1)   =     48.40, Prob > chi2 = 0.0000
(1 observations with missing values excluded;
 specify option 'method(mlmv)' to use all observations)

Endogenous variables


Observed:  shrelaw2 whpsleep


Exogenous variables


Observed:  crhtw1


Fitting target model:

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

Structural equation model                       Number of obs      =       338
 Estimation method = ml
 Log likelihood    = -3684.8571

--------------+----------------------------------------------------------------
              |                 OIM
              |      Coef.    Std. Err.       z    P>|z|      [95% Con f. Interval]
--------------+----------------------------------------------------------------
 Structural    |
   shrelaw2 <- |
       crhtw1 |     9.99232    1.983305      5.04    0.000      6.105114     13.87953
        _cons |     28.2478    1.859702     15.19    0.000      24.60285     31.89275
  ------------+----------------------------------------------------------------
   whpsleep <- |
     shrelaw2 |    .1105453    .0380427      2.91    0.004       .035983     .1851077
        _cons |    14.42955    1.680356      8.59    0.000      11.13611     17.72299
--------------+----------------------------------------------------------------
 Variance      |
   e.shrelaw2 |    1143.469    87.95912                    983.4384      1329.54
   e.whpsleep |    601.3572    46.25825                    517.1964     699.2131
--------------+----------------------------------------------------------------
LR test of model vs. saturated: chi2(1)   =     27.23, Prob > chi2 = 0.0000
(1 observations with missing values excluded;
 specify option 'method(mlmv)' to use all observations)

Endogenous variables


Observed:  shrelaw2 whpsleep


Exogenous variables


Observed:  crhtw2


Fitting target model:

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

Structural equation model                       Number of obs      =       338
 Estimation method = ml
 Log likelihood    = -3678.5528

--------------+----------------------------------------------------------------
              |                 OIM
              |      Coef.    Std. Err.       z    P>|z|      [95% Con f. Interval]
--------------+----------------------------------------------------------------
 Structural    |
   shrelaw2 <- |
       crhtw2 |    11.52971     1.97306      5.84    0.000       7.66258     15.39683
        _cons |    28.96435    1.852722     15.63    0.000      25.33308     32.59562
  ------------+----------------------------------------------------------------
   whpsleep <- |
     shrelaw2 |    .1105453    .0380427      2.91    0.004       .035983     .1851077
        _cons |    14.42955    1.680356      8.59    0.000      11.13611     17.72299
--------------+----------------------------------------------------------------
 Variance      |
   e.shrelaw2 |    1116.541    85.88778                    960.2796      1298.23
   e.whpsleep |    601.3572    46.25825                    517.1964     699.2131
--------------+----------------------------------------------------------------
LR test of model vs. saturated: chi2(1)   =     48.24, Prob > chi2 = 0.0000
(1 observations with missing values excluded;
 specify option 'method(mlmv)' to use all observations)

Endogenous variables


Observed:  shrelaw2 whpsleep


Exogenous variables


Observed:  crhtw3


Fitting target model:

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

Structural equation model                       Number of obs      =       338
 Estimation method = ml
 Log likelihood    = -3679.1833

--------------+----------------------------------------------------------------
              |                 OIM
              |      Coef.    Std. Err.       z    P>|z|      [95% Con f. Interval]
--------------+----------------------------------------------------------------
 Structural    |
   shrelaw2 <- |
       crhtw3 |    11.28255    1.977936      5.70    0.000       7.40587     15.15924
        _cons |    28.93806    1.857402     15.58    0.000      25.29762      32.5785
  ------------+----------------------------------------------------------------
   whpsleep <- |
     shrelaw2 |    .1105453    .0380427      2.91    0.004       .035983     .1851077
        _cons |    14.42955    1.680356      8.59    0.000      11.13611     17.72299
--------------+----------------------------------------------------------------
 Variance      |
   e.shrelaw2 |    1121.391    86.26081                    964.4503     1303.869
   e.whpsleep |    601.3572    46.25825                    517.1964     699.2131
--------------+----------------------------------------------------------------
LR test of model vs. saturated: chi2(1)   =     48.16, Prob > chi2 = 0.0000
(1 observations with missing values excluded;
 specify option 'method(mlmv)' to use all observations)

Endogenous variables


Observed:  dafter whpsleep


Exogenous variables


Observed:  crhtw1


Fitting target model:

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

Structural equation model                       Number of obs      =       338
 Estimation method = ml
 Log likelihood    = -3789.8816

--------------+----------------------------------------------------------------
              |                 OIM
              |      Coef.    Std. Err.       z    P>|z|      [95% Con f. Interval]
--------------+----------------------------------------------------------------
 Structural    |
   dafter <-   |
       crhtw1 |   -3.009633    2.698124     -1.12    0.265     -8.297859     2.278593
        _cons |    5.311681     2.53165      2.10    0.036      .3497375     10.27362
  ------------+----------------------------------------------------------------
   whpsleep <- |
       dafter |    .0688925    .0290778      2.37    0.018      .0119011     .1258839
        _cons |      17.003    1.349736     12.60    0.000      14.35756     19.64843
--------------+----------------------------------------------------------------
 Variance      |
     e.dafter |    2113.784    162.5988                    1817.957      2457.75
   e.whpsleep |    606.3108    46.63929                    521.4567     704.9727
--------------+----------------------------------------------------------------
LR test of model vs. saturated: chi2(1)   =     37.30, Prob > chi2 = 0.0000
(1 observations with missing values excluded;
 specify option 'method(mlmv)' to use all observations)

Endogenous variables


Observed:  dafter whpsleep


Exogenous variables


Observed:  crhtw2


Fitting target model:

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

Structural equation model                       Number of obs      =       338
 Estimation method = ml
 Log likelihood    = -3787.0422

--------------+----------------------------------------------------------------
              |                 OIM
              |      Coef.    Std. Err.       z    P>|z|      [95% Con f. Interval]
--------------+----------------------------------------------------------------
 Structural    |
   dafter <-   |
       crhtw2 |    3.030983    2.720912      1.11    0.265     -2.301906     8.363873
        _cons |    6.320788     2.55246      2.48    0.013      1.318059     11.32352
  ------------+----------------------------------------------------------------
   whpsleep <- |
       dafter |    .0688925    .0290778      2.37    0.018      .0119011     .1258839
        _cons |      17.003    1.349736     12.60    0.000      14.35756     19.64843
--------------+----------------------------------------------------------------
 Variance      |
     e.dafter |    2113.805    162.6004                    1817.975     2457.774
   e.whpsleep |    606.3108    46.63929                    521.4567     704.9727
--------------+----------------------------------------------------------------
LR test of model vs. saturated: chi2(1)   =     56.38, Prob > chi2 = 0.0000
(1 observations with missing values excluded;
 specify option 'method(mlmv)' to use all observations)

Endogenous variables


Observed:  dafter whpsleep


Exogenous variables


Observed:  crhtw3


Fitting target model:

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

Structural equation model                       Number of obs      =       338
 Estimation method = ml
 Log likelihood    = -3786.7464

--------------+----------------------------------------------------------------
              |                 OIM
              |      Coef.    Std. Err.       z    P>|z|      [95% Con f. Interval]
--------------+----------------------------------------------------------------
 Structural    |
   dafter <-   |
       crhtw3 |    3.091364    2.722908      1.14    0.256     -2.245438     8.428166
        _cons |    6.332646    2.552453      2.48    0.013       1.32993     11.33536
  ------------+----------------------------------------------------------------
   whpsleep <- |
       dafter |    .0688925    .0290778      2.37    0.018      .0119011     .1258839
        _cons |      17.003    1.349736     12.60    0.000      14.35756     19.64843
--------------+----------------------------------------------------------------
 Variance      |
     e.dafter |    2113.506    162.5774                    1817.717     2457.426
   e.whpsleep |    606.3108    46.63929                    521.4567     704.9727
--------------+----------------------------------------------------------------
LR test of model vs. saturated: chi2(1)   =     55.88, Prob > chi2 = 0.0000

Endogenous variables


Observed:  BSIposymp whpsleep


Exogenous variables


Observed:  crhtw1


Fitting target model:

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

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

---------------+----------------------------------------------------------------
               |                 OIM
               |      Coef.    Std. Err.       z    P>|z|      [95% Con f. Interval]
---------------+----------------------------------------------------------------
 Structural     |
   BSIposymp <- |
        crhtw1 |    8.980154    1.259147      7.13    0.000      6.512271     11.44804
         _cons |     76.1612    1.182631     64.40    0.000      73.84329     78.47911
  -------------+----------------------------------------------------------------
   whpsleep <-  |
     BSIposymp |    .6228113    .0475743     13.09    0.000      .5295673     .7160552
         _cons |   -29.29129    3.728001     -7.86    0.000     -36.59803    -21.98454
---------------+----------------------------------------------------------------
 Variance       |
   e.BSIposymp |    463.2751    35.58394                    398.5277     538.5417
    e.whpsleep |    408.7877     31.3988                    351.6555      475.202
---------------+----------------------------------------------------------------
LR test of model vs. saturated: chi2(1)   =      5.77, Prob > chi2 = 0.0163

Endogenous variables


Observed:  BSIposymp whpsleep


Exogenous variables


Observed:  crhtw2


Fitting target model:

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

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

---------------+----------------------------------------------------------------
               |                 OIM
               |      Coef.    Std. Err.       z    P>|z|      [95% Con f. Interval]
---------------+----------------------------------------------------------------
 Structural     |
   BSIposymp <- |
        crhtw2 |    12.51389    1.180342     10.60    0.000      10.20046     14.82732
         _cons |     77.1876    1.107911     69.67    0.000      75.01614     79.35907
  -------------+----------------------------------------------------------------
   whpsleep <-  |
     BSIposymp |    .6228113    .0475743     13.09    0.000      .5295673     .7160552
         _cons |   -29.29129    3.728001     -7.86    0.000     -36.59803    -21.98454
---------------+----------------------------------------------------------------
 Variance       |
   e.BSIposymp |    400.1201    30.73304                    344.1993     465.1261
    e.whpsleep |    408.7877     31.3988                    351.6555      475.202
---------------+----------------------------------------------------------------
LR test of model vs. saturated: chi2(1)   =      7.18, Prob > chi2 = 0.0074

Endogenous variables


Observed:  BSIposymp whpsleep


Exogenous variables


Observed:  crhtw3


Fitting target model:

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

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

---------------+----------------------------------------------------------------
               |                 OIM
               |      Coef.    Std. Err.       z    P>|z|      [95% Con f. Interval]
---------------+----------------------------------------------------------------
 Structural     |
   BSIposymp <- |
        crhtw3 |    12.69567    1.176077     10.79    0.000       10.3906     15.00074
         _cons |     77.2243    1.102977     70.01    0.000      75.06251      79.3861
  -------------+----------------------------------------------------------------
   whpsleep <-  |
     BSIposymp |    .6228113    .0475743     13.09    0.000      .5295673     .7160552
         _cons |   -29.29129    3.728001     -7.86    0.000     -36.59803    -21.98454
---------------+----------------------------------------------------------------
 Variance       |
   e.BSIposymp |    396.4926    30.45442                    341.0788     460.9094
    e.whpsleep |    408.7877     31.3988                    351.6555      475.202
---------------+----------------------------------------------------------------
LR test of model vs. saturated: chi2(1)   =      6.51, Prob > chi2 = 0.0107

Endogenous variables


Observed:  WHPpain whpsleep


Exogenous variables


Observed:  crhtw1


Fitting target model:

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

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

--------------+----------------------------------------------------------------
              |                 OIM
              |      Coef.    Std. Err.       z    P>|z|      [95% Con f. Interval]
--------------+----------------------------------------------------------------
 Structural    |
   WHPpain <-  |
       crhtw1 |    3.750412    .9400436      3.99    0.000      1.907961     5.592864
        _cons |    10.74268    .8829189     12.17    0.000      9.012189     12.47317
  ------------+----------------------------------------------------------------
   whpsleep <- |
      WHPpain |    .7067514    .0724028      9.76    0.000      .5648446     .8486583
        _cons |     10.1322    1.401283      7.23    0.000       7.38574     12.87867
--------------+----------------------------------------------------------------
 Variance      |
    e.WHPpain |    258.2155    19.83341                    222.1273     300.1668
   e.whpsleep |    480.4183    36.90071                     413.275     558.4701
--------------+----------------------------------------------------------------
LR test of model vs. saturated: chi2(1)   =     20.43, Prob > chi2 = 0.0000

Endogenous variables


Observed:  WHPpain whpsleep


Exogenous variables


Observed:  crhtw2


Fitting target model:

Iteration 0:   log likelihood = -3386.9057
Iteration 1:   log likelihood = -3386.9057  (backed up)

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

--------------+----------------------------------------------------------------
              |                 OIM
              |      Coef.    Std. Err.       z    P>|z|      [95% Con f. Interval]
--------------+----------------------------------------------------------------
 Structural    |
   WHPpain <-  |
       crhtw2 |    6.464798    .9044495      7.15    0.000      4.692109     8.237486
        _cons |    11.39925     .848949     13.43    0.000      9.735336     13.06316
  ------------+----------------------------------------------------------------
   whpsleep <- |
      WHPpain |    .7067514    .0724028      9.76    0.000      .5648446     .8486583
        _cons |     10.1322    1.401283      7.23    0.000       7.38574     12.87867
--------------+----------------------------------------------------------------
 Variance      |
    e.WHPpain |    234.9328    18.04508                    202.0986     273.1014
   e.whpsleep |    480.4183    36.90071                     413.275     558.4701
--------------+----------------------------------------------------------------
LR test of model vs. saturated: chi2(1)   =     25.55, Prob > chi2 = 0.0000

Endogenous variables


Observed:  WHPpain whpsleep


Exogenous variables


Observed:  crhtw3


Fitting target model:

Iteration 0:   log likelihood = -3382.0218
Iteration 1:   log likelihood = -3382.0218  (backed up)

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

--------------+----------------------------------------------------------------
              |                 OIM
              |      Coef.    Std. Err.       z    P>|z|      [95% Con f. Interval]
--------------+----------------------------------------------------------------
 Structural    |
   WHPpain <-  |
       crhtw3 |    7.018416    .8931795      7.86    0.000      5.267816     8.769015
        _cons |    11.50291    .8376629     13.73    0.000      9.861124      13.1447
  ------------+----------------------------------------------------------------
   whpsleep <- |
      WHPpain |    .7067514    .0724028      9.76    0.000      .5648446     .8486583
        _cons |     10.1322    1.401283      7.23    0.000       7.38574     12.87867
--------------+----------------------------------------------------------------
 Variance      |
    e.WHPpain |    228.6868    17.56533                    196.7256     265.8408
   e.whpsleep |    480.4183    36.90071                     413.275     558.4701
--------------+----------------------------------------------------------------
LR test of model vs. saturated: chi2(1)   =     22.57, Prob > chi2 = 0.0000

Endogenous variables


Observed:  illw1 whpsleep


Exogenous variables


Observed:  crhtw1


Fitting target model:

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

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

--------------+----------------------------------------------------------------
              |                 OIM
              |      Coef.    Std. Err.       z    P>|z|      [95% Con f. Interval]
--------------+----------------------------------------------------------------
 Structural    |
   illw1 <-    |
       crhtw1 |   -.0014219    .0213089     -0.07    0.947     -.0431866     .0403427
        _cons |    .0941932     .020014      4.71    0.000      .0549665     .1334199
  ------------+----------------------------------------------------------------
   whpsleep <- |
        illw1 |    12.30936    3.638136      3.38    0.001      5.178746     19.43998
        _cons |    16.18596    1.368988     11.82    0.000      13.50279     18.86913
--------------+----------------------------------------------------------------
 Variance      |
      e.illw1 |    .1326807    .0101911                    .1141372     .1542369
   e.whpsleep |    595.3481    45.72842                    512.1422     692.0722
--------------+----------------------------------------------------------------
LR test of model vs. saturated: chi2(1)   =     35.15, Prob > chi2 = 0.0000

Endogenous variables


Observed:  illw1 whpsleep


Exogenous variables


Observed:  crhtw2


Fitting target model:

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

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

--------------+----------------------------------------------------------------
              |                 OIM
              |      Coef.    Std. Err.       z    P>|z|      [95% Con f. Interval]
--------------+----------------------------------------------------------------
 Structural    |
   illw1 <-    |
       crhtw2 |    .0494761    .0213255      2.32    0.020       .007679     .0912733
        _cons |    .1034992    .0200169      5.17    0.000      .0642669     .1427316
  ------------+----------------------------------------------------------------
   whpsleep <- |
        illw1 |    12.30936    3.638136      3.38    0.001      5.178746     19.43998
        _cons |    16.18596    1.368988     11.82    0.000      13.50279     18.86913
--------------+----------------------------------------------------------------
 Variance      |
      e.illw1 |    .1306087     .010032                    .1123548     .1518282
   e.whpsleep |    595.3481    45.72842                    512.1422     692.0722
--------------+----------------------------------------------------------------
LR test of model vs. saturated: chi2(1)   =     52.26, Prob > chi2 = 0.0000

Endogenous variables


Observed:  illw1 whpsleep


Exogenous variables


Observed:  crhtw3


Fitting target model:

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

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

--------------+----------------------------------------------------------------
              |                 OIM
              |      Coef.    Std. Err.       z    P>|z|      [95% Con f. Interval]
--------------+----------------------------------------------------------------
 Structural    |
   illw1 <-    |
       crhtw3 |     .049437     .021346      2.32    0.021      .0075996     .0912743
        _cons |    .1035047    .0200192      5.17    0.000      .0642678     .1427416
  ------------+----------------------------------------------------------------
   whpsleep <- |
        illw1 |    12.30936    3.638136      3.38    0.001      5.178746     19.43998
        _cons |    16.18596    1.368988     11.82    0.000      13.50279     18.86913
--------------+----------------------------------------------------------------
 Variance      |
      e.illw1 |    .1306158    .0100325                    .1123609     .1518365
   e.whpsleep |    595.3481    45.72842                    512.1422     692.0722
--------------+----------------------------------------------------------------
LR test of model vs. saturated: chi2(1)   =     51.85, Prob > chi2 = 0.0000

Endogenous variables


Observed:  illw3 whpsleep


Exogenous variables


Observed:  crhtw1


Fitting target model:

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

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

--------------+----------------------------------------------------------------
              |                 OIM
              |      Coef.    Std. Err.       z    P>|z|      [95% Con f. Interval]
--------------+----------------------------------------------------------------
 Structural    |
   illw3 <-    |
       crhtw1 |   -.0596154    .0551874     -1.08    0.280     -.1677807     .0485498
        _cons |    .4841529    .0518337      9.34    0.000      .3825607     .5857452
  ------------+----------------------------------------------------------------
   whpsleep <- |
        illw3 |    4.812018    1.401674      3.43    0.001      2.064787      7.55925
        _cons |    14.97738    1.493749     10.03    0.000      12.04969     17.90508
--------------+----------------------------------------------------------------
 Variance      |
      e.illw3 |    .8899501    .0683567                    .7655706     1.034537
   e.whpsleep |     594.774    45.68432                    511.6484     691.4048
--------------+----------------------------------------------------------------
LR test of model vs. saturated: chi2(1)   =     37.79, Prob > chi2 = 0.0000

Endogenous variables


Observed:  illw3 whpsleep


Exogenous variables


Observed:  crhtw2


Fitting target model:

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

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

--------------+----------------------------------------------------------------
              |                 OIM
              |      Coef.    Std. Err.       z    P>|z|      [95% Con f. Interval]
--------------+----------------------------------------------------------------
 Structural    |
   illw3 <-    |
       crhtw2 |    .2464559     .054132      4.55    0.000      .1403591     .3525526
        _cons |    .5379751    .0508102     10.59    0.000      .4383888     .6375613
  ------------+----------------------------------------------------------------
   whpsleep <- |
        illw3 |    4.812018    1.401674      3.43    0.001      2.064787      7.55925
        _cons |    14.97738    1.493749     10.03    0.000      12.04969     17.90508
--------------+----------------------------------------------------------------
 Variance      |
      e.illw3 |    .8415554    .0646395                    .7239396     .9782798
   e.whpsleep |     594.774    45.68432                    511.6484     691.4048
--------------+----------------------------------------------------------------
LR test of model vs. saturated: chi2(1)   =     48.20, Prob > chi2 = 0.0000

Endogenous variables


Observed:  illw3 whpsleep


Exogenous variables


Observed:  crhtw3


Fitting target model:

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

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

--------------+----------------------------------------------------------------
              |                 OIM
              |      Coef.    Std. Err.       z    P>|z|      [95% Con f. Interval]
--------------+----------------------------------------------------------------
 Structural    |
   illw3 <-    |
       crhtw3 |    .3093305     .053226      5.81    0.000      .2050096     .4136515
        _cons |    .5496237    .0499176     11.01    0.000      .4517869     .6474605
  ------------+----------------------------------------------------------------
   whpsleep <- |
        illw3 |    4.812018    1.401674      3.43    0.001      2.064787      7.55925
        _cons |    14.97738    1.493749     10.03    0.000      12.04969     17.90508
--------------+----------------------------------------------------------------
 Variance      |
      e.illw3 |    .8121021    .0623772                    .6986027     .9440414
   e.whpsleep |     594.774    45.68432                    511.6484     691.4048
--------------+----------------------------------------------------------------
LR test of model vs. saturated: chi2(1)   =     46.28, Prob > chi2 = 0.0000

Endogenous variables


Observed:  avgcumdosew1 whpsleep


Exogenous variables


Observed:  crhtw1


Fitting target model:

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

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

------------------+----------------------------------------------------------------
                  |                 OIM
                  |      Coef.    Std. Err.       z    P>|z|      [95% Con f. Interval]
------------------+----------------------------------------------------------------
 Structural        |
   avgcumdosew1 <- |
           crhtw1 |    .1559612    .0972536      1.60    0.109     -.0346524     .3465748
            _cons |    .4522623    .0913437      4.95    0.000      .2732319     .6312927
  ----------------+----------------------------------------------------------------
   whpsleep <-     |
     avgcumdosew1 |    1.897497    .8008306      2.37    0.018      .3278975     3.467096
            _cons |     16.5318    1.380055     11.98    0.000      13.82694     19.23666
------------------+----------------------------------------------------------------
 Variance          |
   e.avgcumdosew1 |    2.763744    .2122819                    2.377483     3.212759
       e.whpsleep |    605.4259    46.50248                    520.8115     703.7872
------------------+----------------------------------------------------------------
LR test of model vs. saturated: chi2(1)   =     32.09, Prob > chi2 = 0.0000

Endogenous variables


Observed:  avgcumdosew1 whpsleep


Exogenous variables


Observed:  crhtw2


Fitting target model:

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

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

------------------+----------------------------------------------------------------
                  |                 OIM
                  |      Coef.    Std. Err.       z    P>|z|      [95% Con f. Interval]
------------------+----------------------------------------------------------------
 Structural        |
   avgcumdosew1 <- |
           crhtw2 |    .1298963    .0982166      1.32    0.186     -.0626047     .3223974
            _cons |    .4539992    .0921897      4.92    0.000      .2733108     .6346877
  ----------------+----------------------------------------------------------------
   whpsleep <-     |
     avgcumdosew1 |    1.897497    .8008306      2.37    0.018      .3278975     3.467096
            _cons |     16.5318    1.380055     11.98    0.000      13.82694     19.23666
------------------+----------------------------------------------------------------
 Variance          |
   e.avgcumdosew1 |    2.770416    .2127943                    2.383222     3.220515
       e.whpsleep |    605.4259    46.50248                    520.8115     703.7872
------------------+----------------------------------------------------------------
LR test of model vs. saturated: chi2(1)   =     54.83, Prob > chi2 = 0.0000

Endogenous variables


Observed:  avgcumdosew1 whpsleep


Exogenous variables


Observed:  crhtw3


Fitting target model:

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

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

------------------+----------------------------------------------------------------
                  |                 OIM
                  |      Coef.    Std. Err.       z    P>|z|      [95% Con f. Interval]
------------------+----------------------------------------------------------------
 Structural        |
   avgcumdosew1 <- |
           crhtw3 |    .1283977    .0983147      1.31    0.192     -.0642955     .3210909
            _cons |    .4537564    .0922038      4.92    0.000      .2730402     .6344725
  ----------------+----------------------------------------------------------------
   whpsleep <-     |
     avgcumdosew1 |    1.897497    .8008306      2.37    0.018      .3278975     3.467096
            _cons |     16.5318    1.380055     11.98    0.000      13.82694     19.23666
------------------+----------------------------------------------------------------
 Variance          |
   e.avgcumdosew1 |     2.77077    .2128215                    2.383527     3.220926
       e.whpsleep |    605.4259    46.50248                    520.8115     703.7872
------------------+----------------------------------------------------------------
LR test of model vs. saturated: chi2(1)   =     54.43, Prob > chi2 = 0.0000

Endogenous variables


Observed:  icdx3nr7 whpsleep


Exogenous variables


Observed:  crhtw1


Fitting target model:

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

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

--------------+----------------------------------------------------------------
              |                 OIM
              |      Coef.    Std. Err.       z    P>|z|      [95% Con f. Interval]
--------------+----------------------------------------------------------------
 Structural    |
   icdx3nr7 <- |
       crhtw1 |    .0263074    .0075801      3.47    0.001      .0114507     .0411641
        _cons |    .0214379    .0071194      3.01    0.003       .007484     .0353917
  ------------+----------------------------------------------------------------
   whpsleep <- |
     icdx3nr7 |    29.31936    10.09395      2.90    0.004       9.53558     49.10313
        _cons |    16.82898    1.342879     12.53    0.000      14.19698     19.46097
--------------+----------------------------------------------------------------
 Variance      |
   e.icdx3nr7 |    .0167893    .0012896                    .0144428      .019517
   e.whpsleep |    600.5069    46.12466                      516.58     698.0691
--------------+----------------------------------------------------------------
LR test of model vs. saturated: chi2(1)   =     29.26, Prob > chi2 = 0.0000

Endogenous variables


Observed:  icdx3nr7 whpsleep


Exogenous variables


Observed:  crhtw2


Fitting target model:

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

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

--------------+----------------------------------------------------------------
              |                 OIM
              |      Coef.    Std. Err.       z    P>|z|      [95% Con f. Interval]
--------------+----------------------------------------------------------------
 Structural    |
   icdx3nr7 <- |
       crhtw2 |    .0249057    .0076621      3.25    0.001      .0098883     .0399231
        _cons |    .0222819    .0071919      3.10    0.002      .0081861     .0363778
  ------------+----------------------------------------------------------------
   whpsleep <- |
     icdx3nr7 |    29.31936    10.09395      2.90    0.004       9.53558     49.10313
        _cons |    16.82898    1.342879     12.53    0.000      14.19698     19.46097
--------------+----------------------------------------------------------------
 Variance      |
   e.icdx3nr7 |    .0168604     .001295                     .014504     .0195996
   e.whpsleep |    600.5069    46.12466                      516.58     698.0691
--------------+----------------------------------------------------------------
LR test of model vs. saturated: chi2(1)   =     51.24, Prob > chi2 = 0.0000

Endogenous variables


Observed:  icdx3nr7 whpsleep


Exogenous variables


Observed:  crhtw3


Fitting target model:

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

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

--------------+----------------------------------------------------------------
              |                 OIM
              |      Coef.    Std. Err.       z    P>|z|      [95% Con f. Interval]
--------------+----------------------------------------------------------------
 Structural    |
   icdx3nr7 <- |
       crhtw3 |    .0234479     .007683      3.05    0.002      .0083895     .0385063
        _cons |    .0220197    .0072054      3.06    0.002      .0078973     .0361421
  ------------+----------------------------------------------------------------
   whpsleep <- |
     icdx3nr7 |    29.31936    10.09395      2.90    0.004       9.53558     49.10313
        _cons |    16.82898    1.342879     12.53    0.000      14.19698     19.46097
--------------+----------------------------------------------------------------
 Variance      |
   e.icdx3nr7 |    .0169209    .0012997                    .0145561       .01967
   e.whpsleep |    600.5069    46.12466                      516.58     698.0691
--------------+----------------------------------------------------------------
LR test of model vs. saturated: chi2(1)   =     51.12, Prob > chi2 = 0.0000

Endogenous variables


Observed:  icdx4nr7 whpsleep


Exogenous variables


Observed:  crhtw1


Fitting target model:

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

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

--------------+----------------------------------------------------------------
              |                 OIM
              |      Coef.    Std. Err.       z    P>|z|      [95% Con f. Interval]
--------------+----------------------------------------------------------------
 Structural    |
   icdx4nr7 <- |
       crhtw1 |    .0047559    .0031621      1.50    0.133     -.0014417     .0109534
        _cons |    .0036257    .0029699      1.22    0.222     -.0021952     .0094467
  ------------+----------------------------------------------------------------
   whpsleep <- |
     icdx4nr7 |    60.46045    24.62702      2.46    0.014      12.19238     108.7285
        _cons |    17.16956    1.337556     12.84    0.000        14.548     19.79112
--------------+----------------------------------------------------------------
 Variance      |
   e.icdx4nr7 |    .0029217    .0002244                    .0025133     .0033963
   e.whpsleep |     604.701     46.4468                    520.1879     702.9445
--------------+----------------------------------------------------------------
LR test of model vs. saturated: chi2(1)   =     32.17, Prob > chi2 = 0.0000

Endogenous variables


Observed:  icdx4nr7 whpsleep


Exogenous variables


Observed:  crhtw2


Fitting target model:

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

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

--------------+----------------------------------------------------------------
              |                 OIM
              |      Coef.    Std. Err.       z    P>|z|      [95% Con f. Interval]
--------------+----------------------------------------------------------------
 Structural    |
   icdx4nr7 <- |
       crhtw2 |    .0048025    .0031895      1.51    0.132     -.0014488     .0110539
        _cons |    .0038336    .0029938      1.28    0.200     -.0020342     .0097013
  ------------+----------------------------------------------------------------
   whpsleep <- |
     icdx4nr7 |    60.46045    24.62702      2.46    0.014      12.19238     108.7285
        _cons |    17.16956    1.337556     12.84    0.000        14.548     19.79112
--------------+----------------------------------------------------------------
 Variance      |
   e.icdx4nr7 |    .0029216    .0002244                    .0025133     .0033963
   e.whpsleep |     604.701     46.4468                    520.1879     702.9445
--------------+----------------------------------------------------------------
LR test of model vs. saturated: chi2(1)   =     54.48, Prob > chi2 = 0.0000

Endogenous variables


Observed:  icdx4nr7 whpsleep


Exogenous variables


Observed:  crhtw3


Fitting target model:

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

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

--------------+----------------------------------------------------------------
              |                 OIM
              |      Coef.    Std. Err.       z    P>|z|      [95% Con f. Interval]
--------------+----------------------------------------------------------------
 Structural    |
   icdx4nr7 <- |
       crhtw3 |    .0046896     .003193      1.47    0.142     -.0015685     .0109478
        _cons |     .003814    .0029945      1.27    0.203     -.0020552     .0096832
  ------------+----------------------------------------------------------------
   whpsleep <- |
     icdx4nr7 |    60.46045    24.62702      2.46    0.014      12.19238     108.7285
        _cons |    17.16956    1.337556     12.84    0.000        14.548     19.79112
--------------+----------------------------------------------------------------
 Variance      |
   e.icdx4nr7 |    .0029226    .0002245                    .0025141     .0033974
   e.whpsleep |     604.701     46.4468                    520.1879     702.9445
--------------+----------------------------------------------------------------
LR test of model vs. saturated: chi2(1)   =     54.12, Prob > chi2 = 0.0000

Endogenous variables


Observed:  icdx4nr12 whpsleep


Exogenous variables


Observed:  crhtw1


Fitting target model:

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

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

---------------+----------------------------------------------------------------
               |                 OIM
               |      Coef.    Std. Err.       z    P>|z|      [95% Con f. Interval]
---------------+----------------------------------------------------------------
 Structural     |
   icdx4nr12 <- |
        crhtw1 |   -.0009735    .0031722     -0.31    0.759     -.0071908     .0052438
         _cons |    .0028115    .0029794      0.94    0.345      -.003028      .008651
  -------------+----------------------------------------------------------------
   whpsleep <-  |
     icdx4nr12 |    60.46045    24.62702      2.46    0.014      12.19238     108.7285
         _cons |    17.16956    1.337556     12.84    0.000        14.548     19.79112
---------------+----------------------------------------------------------------
 Variance       |
   e.icdx4nr12 |    .0029403    .0002258                    .0025294      .003418
    e.whpsleep |     604.701     46.4468                    520.1879     702.9445
---------------+----------------------------------------------------------------
LR test of model vs. saturated: chi2(1)   =     34.95, Prob > chi2 = 0.0000

Endogenous variables


Observed:  icdx4nr12 whpsleep


Exogenous variables


Observed:  crhtw2


Fitting target model:

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

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

---------------+----------------------------------------------------------------
               |                 OIM
               |      Coef.    Std. Err.       z    P>|z|      [95% Con f. Interval]
---------------+----------------------------------------------------------------
 Structural     |
   icdx4nr12 <- |
        crhtw2 |   -.0010612    .0031996     -0.33    0.740     -.0073323       .00521
         _cons |    .0027546    .0030033      0.92    0.359     -.0031318     .0086409
  -------------+----------------------------------------------------------------
   whpsleep <-  |
     icdx4nr12 |    60.46045    24.62702      2.46    0.014      12.19238     108.7285
         _cons |    17.16956    1.337556     12.84    0.000        14.548     19.79112
---------------+----------------------------------------------------------------
 Variance       |
   e.icdx4nr12 |    .0029402    .0002258                    .0025293     .0034179
    e.whpsleep |     604.701     46.4468                    520.1879     702.9445
---------------+----------------------------------------------------------------
LR test of model vs. saturated: chi2(1)   =     58.26, Prob > chi2 = 0.0000

Endogenous variables


Observed:  icdx4nr12 whpsleep


Exogenous variables


Observed:  crhtw3


Fitting target model:

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

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

---------------+----------------------------------------------------------------
               |                 OIM
               |      Coef.    Std. Err.       z    P>|z|      [95% Con f. Interval]
---------------+----------------------------------------------------------------
 Structural     |
   icdx4nr12 <- |
        crhtw3 |   -.0010542    .0032026     -0.33    0.742     -.0073312     .0052229
         _cons |    .0027556    .0030036      0.92    0.359     -.0031313     .0086425
  -------------+----------------------------------------------------------------
   whpsleep <-  |
     icdx4nr12 |    60.46045    24.62702      2.46    0.014      12.19238     108.7285
         _cons |    17.16956    1.337556     12.84    0.000        14.548     19.79112
---------------+----------------------------------------------------------------
 Variance       |
   e.icdx4nr12 |    .0029402    .0002258                    .0025293     .0034179
    e.whpsleep |     604.701     46.4468                    520.1879     702.9445
---------------+----------------------------------------------------------------
LR test of model vs. saturated: chi2(1)   =     57.81, Prob > chi2 = 0.0000

.
. * Possible male sleep indirect effects ----------------------------------------
. *-+- Possible indirect effects through
. * wave 1 icdx3nr7 pain BISposymp shrelaw1 shrelaw2 shrelaw1 age
. * wave 2 icdx3nr7 illw3 illw1 pain BSIposymp shrelaw2 shrelaw1 inc1w2 inc1w age
. * wave 3 icdx3nr7 illw3 illw1 pain BSIposymp shrelaw2 shrelaw1  inc1w2 age
. *------------------------------------------------------------------------------
. /*
>
> EQ(71) Modelling WHPsleep by OLS-CS
>        The dataset is: /Users/robertyaffee/Documents/data/research/chwk/phase3/
>        data/ox/workingdatasets/MARS/gals.dta
>        The estimation sample is: 1 - 363
>        Dropped 2 observation(s) with missing values from the sample
>
>                   Coefficient  Std.Error      HACSE  t-HACSE  t-prob Part.R^2
> sepaw2               -26.8551      8.383      5.060    -5.31  0.0000   0.0745
> shhousw2            0.0659946    0.03288    0.03315     1.99  0.0473   0.0112
> sufamw1            -0.0327914    0.03549    0.02975    -1.10  0.2711   0.0035
> pillw2               -1.10178     0.3324     0.2845    -3.87  0.0001   0.0411
> radchw3             -0.128170    0.03728    0.02801    -4.58  0.0000   0.0564
> WHPpain              0.438681    0.07103    0.07325     5.99  0.0000   0.0930
> HP2probsoc            14.3111      3.544      4.216     3.39  0.0008   0.0319
> BSIposymp            0.219874    0.03616    0.03478     6.32  0.0000   0.1025
> radhw1       U     0.0290023    0.05173    0.04457    0.651  0.5157   0.0012
> radhw2       U     -0.128956     0.1199     0.1024    -1.26  0.2087   0.0045
> radhw3       U      0.187958     0.1026    0.09274     2.03  0.0434   0.0116
>
> sigma                  23.044  RSS                 185858.71
> log-likelihood       -1639.25
> no. of observations       361  no. of parameters          11
> mean(WHPsleep)        26.3552  se(WHPsleep)          30.9492
> When the log-likelihood constant is NOT included:
> AIC                   6.30481  SC                    6.42330
> HQ                    6.35192  FPE                   547.206
> When the log-likelihood constant is included:
> AIC                   9.14268  SC                    9.26118
> HQ                    9.18980  FPE                   9345.98
>
> Normality test:   Chi^2(2)  =   72.322 [0.0000]**
> Hetero test:      F(20,340) =   1.7360 [0.0267]*
> Hetero-X test:    F(56,304) =  0.91870 [0.6407]
> RESET23 test:     F(2,348)  =   10.472 [0.0000]**
> */
.
. /* candiate illnesses  for female sleep model
>
> foreach i in 7 {
> foreach j in 1 5 6 {
> foreach k in 3 {
> foreach m in 2 {
> foreach n in 4 {
> regress whpsleep icdx`i'nr1-icdx`i'nr28 if gender==2
> regress whpsleep icdx`j'nr1-icdx`j'nr18 if gender==2
> regress whpsleep icdx`k'nr1-icdx`k'nr11 if gender==2
> regress whpsleep icdx`m'nr1-icdx`m'nr8 if gender==2
> regress whpsleep icdx`n'nr1-icdx`n'nr20 if gender==2
> }
> }
> }
> }
> }
> */
.
. * for females:
. des icdx1nr9 icdx3nr4 icdx3nr5 icdx4nr6 icdx4nr17 icdx5nr4 icdx5nr10 icdx6nr5 ///
>   icdx6nr6 icdx6nr11 icdx6nr12 icdx6nr15

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


icdx1nr9          double   %8.0g                        icdx1nr==454 chronic t & a dis


icdx3nr4          double   %8.0g                        icdx3nr==rheum fev w/o hrt involv


icdx3nr5          double   %8.0g                        icdx3nr==hypertension


icdx4nr6          double   %8.0g                        icdx4nr==hypertension


icdx4nr17         double   %8.0g                        icdx4nr==renal/ureteral calculus


icdx5nr4          double   %8.0g                        icdx5nr==rheum fev w/o hrt involv


icdx5nr10         double   %8.0g                        icdx5nr==ac bronchitis/brnchial


icdx6nr5          double   %8.0g                        icdx6nr==hypertension


icdx6nr6          double   %8.0g                        icdx6nr==acute myocardial infarct


icdx6nr11         double   %8.0g                        icdx6nr==regional enteritis


icdx6nr12         double   %8.0g                        icdx6nr==575.1 cholecystitis


icdx6nr15         double   %8.0g                        icdx6nr==renal/ureteral calculus


.
. *-----------------   Female Sleep crhtw models--------------------------------------
. // H8: Full Female Sleep Model
.
. des bffpain24 bffpain26 BFdepx1 icdx1nr9 icdx3nr4 icdx3nr5 icdx4nr6 icdx4nr17 ///
>    icdx5nr4 icdx5nr10 icdx6nr5 icdx6nr6 icdx6nr11 icdx6nr12 icdx6nr15

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


bffpain24         float    %9.0g                        max(0, 21 - csprbslv) * bffpain21


bffpain26         float    %9.0g                        max(0, 60 - phlthw1)


BFdepx1           float    %9.0g                         max(0, WHPer + 3.14933e-007)


icdx1nr9          double   %8.0g                        icdx1nr==454 chronic t & a dis


icdx3nr4          double   %8.0g                        icdx3nr==rheum fev w/o hrt involv


icdx3nr5          double   %8.0g                        icdx3nr==hypertension


icdx4nr6          double   %8.0g                        icdx4nr==hypertension


icdx4nr17         double   %8.0g                        icdx4nr==renal/ureteral calculus


icdx5nr4          double   %8.0g                        icdx5nr==rheum fev w/o hrt involv


icdx5nr10         double   %8.0g                        icdx5nr==ac bronchitis/brnchial


icdx6nr5          double   %8.0g                        icdx6nr==hypertension


icdx6nr6          double   %8.0g                        icdx6nr==acute myocardial infarct


icdx6nr11         double   %8.0g                        icdx6nr==regional enteritis


icdx6nr12         double   %8.0g                        icdx6nr==575.1 cholecystitis


icdx6nr15         double   %8.0g                        icdx6nr==renal/ureteral calculus


.
. regress WHPsleep age sepaw2  sufamw1  pillw2  radchw3 WHPpain bffpain24 bffpain26   ///
>   BSIanx HP2probsoc HP2pbfhm BSIposymp BSItotal BSIdep crhtw1 crhtw2 crhtw3  ///
>   BFdepx1 ///
>   icdx1nr9 icdx3nr4 icdx3nr5 icdx4nr6 icdx4nr17 icdx5nr4 icdx5nr10 icdx6nr5 ///
>   icdx6nr6 icdx6nr11 icdx6nr12 icdx6nr15 if gender==2, vce(cluster id)
note: icdx6nr15 omitted because of collinearity

Linear regression                                      Number of obs =     361
                                                       F( 27,   360) =       .
                                                       Prob > F      =       .
                                                       R-squared     =  0.6163
                                                       Root MSE      =  19.992

(Std. Err. adjusted for 361 clusters in id)
-------------+----------------------------------------------------------------
             |               Robust
     WHPsleep |      Coef.    Std. Err.       t    P>|t|      [95% Con f. Interval]
-------------+----------------------------------------------------------------
         age |    .2885584     .108848      2.65    0.008      .0745006     .5026161
      sepaw2 |   -18.78591    4.288654     -4.38    0.000     -27.21987    -10.35195
     sufamw1 |   -.0056148    .0305865     -0.18    0.854     -.0657655     .0545359
      pillw2 |   -.7835463    .3081868     -2.54    0.011     -1.389619    -.1774736
     radchw3 |   -.0644864    .0320315     -2.01    0.045     -.1274788     -.001494
     WHPpain |    .2091758    .0760946      2.75    0.006      .0595301     .3588215
   bffpain24 |   -5.310912    1.011869     -5.25    0.000     -7.300828    -3.320995
   bffpain26 |    -.242837    .1108689     -2.19    0.029     -.4608691    -.0248049
      BSIanx |    .2567236    .6359012      0.40    0.687     -.9938241     1.507271
  HP2probsoc |     7.10925    4.185774      1.70    0.090      -1.12239     15.34089
    HP2pbfhm |   -2.568178     4.57417     -0.56    0.575     -11.56363     6.427272
   BSIposymp |    2.567843     .719277      3.57    0.000      1.153331     3.982356
    BSItotal |   -2.231476    .6698732     -3.33    0.001     -3.548833    -.9141203
      BSIdep |   -.7409658    .5678523     -1.30    0.193      -1.85769     .3757587
      crhtw1 |    3.684884    1.716231      2.15    0.032      .3097866     7.059982
      crhtw2 |   -1.956705    3.852706     -0.51    0.612     -9.533342     5.619931
      crhtw3 |    2.381087    3.424844      0.70    0.487     -4.354128     9.116301
     BFdepx1 |    .6455258    .1099483      5.87    0.000      .4293041     .8617475
    icdx1nr9 |   -.5814064    7.835166     -0.07    0.941     -15.98985     14.82704
    icdx3nr4 |    11.16885     11.3931      0.98    0.328     -11.23655     33.57424
    icdx3nr5 |    11.09276    4.120844      2.69    0.007      2.988807     19.19671
    icdx4nr6 |    26.93376     12.5339      2.15    0.032       2.28489     51.58262
   icdx4nr17 |    23.14041    20.69039      1.12    0.264     -17.54879     63.82962
    icdx5nr4 |    53.13371    6.765575      7.85    0.000      39.82869     66.43872
   icdx5nr10 |   -10.45278    12.66422     -0.83    0.410     -35.35792     14.45237
    icdx6nr5 |    64.92306    12.12459      5.35    0.000      41.07912     88.76699
    icdx6nr6 |   -4.550941    7.811188     -0.58    0.561     -19.91223     10.81035
   icdx6nr11 |    24.76754     19.6241      1.26    0.208     -13.82473     63.35981
   icdx6nr12 |    .0993212     5.82067      0.02    0.986     -11.34746     11.54611
   icdx6nr15 |           0   (omitted)
       _cons |   -11.61046    7.229918     -1.61    0.109     -25.82864     2.607715
-------------+----------------------------------------------------------------

.
. // H8: Trimmed Female Sleep Model
.
. des bffpain24 bffpain26 BFdepx1  icdx3nr5 icdx4nr6 ///
>   icdx5nr4 icdx6nr5

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


bffpain24         float    %9.0g                        max(0, 21 - csprbslv) * bffpain21


bffpain26         float    %9.0g                        max(0, 60 - phlthw1)


BFdepx1           float    %9.0g                         max(0, WHPer + 3.14933e-007)


icdx3nr5          double   %8.0g                        icdx3nr==hypertension


icdx4nr6          double   %8.0g                        icdx4nr==hypertension


icdx5nr4          double   %8.0g                        icdx5nr==rheum fev w/o hrt involv


icdx6nr5          double   %8.0g                        icdx6nr==hypertension


.
. cap drop H8FemaleSLpred

. cap drop h8FSLsepred

. cap drop FSLrsires

. cap drop upbSLF

. cap drop lpbSLF

. cap drop rFslp1

.
.
.
.
.
.
. regress WHPsleep age sepaw2   pillw2  radchw3 WHPpain bffpain24 bffpain26   ///
>   BSIanx HP2probsoc BSIposymp BSItotal BSIdep crhtw1 crhtw2 crhtw3  ///
>   BFdepx1 icdx3nr5 icdx4nr6 icdx5nr4 icdx6nr5 ///
>    if gender==2, vce(cluster id)

Linear regression                                      Number of obs =     362
                                                       F( 19,   361) =       .
                                                       Prob > F      =       .
                                                       R-squared     =  0.6054
                                                       Root MSE      =  19.996

(Std. Err. adjusted for 362 clusters in id)
-------------+----------------------------------------------------------------
             |               Robust
     WHPsleep |      Coef.    Std. Err.       t    P>|t|      [95% Con f. Interval]
-------------+----------------------------------------------------------------
         age |    .2748115    .1074232      2.56    0.011      .0635576     .4860654
      sepaw2 |   -18.77988    4.034896     -4.65    0.000     -26.71473    -10.84502
      pillw2 |   -.8019393    .3119536     -2.57    0.011     -1.415414    -.1884648
     radchw3 |   -.0641173    .0310957     -2.06    0.040     -.1252688    -.0029658
     WHPpain |      .21312    .0751605      2.84    0.005      .0653127     .3609274
   bffpain24 |    -5.31569    .9930722     -5.35    0.000     -7.268624    -3.362757
   bffpain26 |   -.2349557    .0961657     -2.44    0.015      -.424071    -.0458403
      BSIanx |    .3524595    .6178581      0.57    0.569     -.8625938     1.567513
  HP2probsoc |    6.639116    3.844953      1.73    0.085     -.9222049     14.20044
   BSIposymp |    2.501887    .7092851      3.53    0.000      1.107037     3.896737
    BSItotal |   -2.168287    .6649841     -3.26    0.001     -3.476016    -.8605573
      BSIdep |   -.7581339     .567056     -1.34    0.182     -1.873282      .357014
      crhtw1 |    3.307012    1.606036      2.06    0.040      .1486505     6.465373
      crhtw2 |   -1.407417    3.559863     -0.40    0.693     -8.408091     5.593257
      crhtw3 |    2.090788    3.264265      0.64    0.522     -4.328575     8.510152
     BFdepx1 |    .6379949    .1054825      6.05    0.000      .4305576     .8454322
    icdx3nr5 |    10.05342    4.054206      2.48    0.014      2.080598     18.02625
    icdx4nr6 |    20.81846    10.55316      1.97    0.049      .0650634     41.57185
    icdx5nr4 |    50.87478     6.54591      7.77    0.000      38.00187     63.74769
    icdx6nr5 |    58.50618    14.84362      3.94    0.000      29.31535     87.69701
       _cons |   -11.47416    6.993127     -1.64    0.102     -25.22654     2.278222
-------------+----------------------------------------------------------------

.    di e(r2_a)
.58224747

.
.  // residual diagnostics
.
.  rdiagrbst rFslp1
(option xb assumed; fitted values)
(1 missing value generated)
yvar =  WHPsleep
age sepaw2 pillw2 radchw3 WHPpain bffpain24 bffpain26 BSIanx HP2probsoc BSIposymp BSItotal BSIdep crhtw1 crhtw2 crhtw3 BFdepx1 icdx3nr5 icdx4nr6 icdx5nr4 icdx6nr5 _cons
168
162
xvar = ,age sepaw2 pillw2 radchw3 WHPpain bffpain24 bffpain26 BSIanx HP2probsoc BSIposymp BSItotal BSIdep crhtw1 crhtw2 crhtw3 BFdepx1 icdx3nr5 icdx4nr6 icdx5nr4 icdx6nr5

   | 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%    -10.76906      -14.52883
 5%    -3.583781      -13.27693
10%    -.9605731      -13.00432       Obs                 701
25%     6.234726      -12.72207       Sum of Wgt.         701

50%     17.61371                      Mean           22.91103
                        Largest       Std. Dev.      22.26136
75%     35.87301       95.69287
90%     55.72131       100.8995       Variance       495.5683
95%      67.5869            101       Skewness       .9727678
99%     83.28068       121.0178       Kurtosis       3.693234
(1 missing value generated)
(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:   rFslp1std |
  +------------------+
  | 483.   -1.681832 |
  | 580.   -1.625595 |
  | 635.   -1.613349 |
  | 596.    -1.60067 |
  | 664.    -1.57186 |
  +------------------+

  +-----------------+
  | 323.   3.114153 |
  | 335.   3.269424 |
  | 333.   3.503309 |
  | 173.   3.507825 |
  |  55.   4.407042 |
  +-----------------+
There are 3 positive outliers

     +----------------+
     |  id   rFslp1~d |
     +----------------+
 55. | 465   4.407042 |
173. | 266   3.507825 |
333. | 535   3.503309 |
     +----------------+
There are 38 negative outliers:

     +-----------------+
     |  id   rFslp1std |
     +-----------------+
 16. | 122   -1.409681 |
113. | 304   -1.442372 |
134. | 341   -1.226688 |
142. | 421   -1.253992 |
185. | 514   -1.279927 |
     +-----------------+
234. |  93   -1.214576 |
240. | 492   -1.201321 |
241. | 269    -1.19017 |
243. | 326   -1.487744 |
452. | 213   -1.265913 |
     +-----------------+
465. | 671   -1.362111 |
475. | 640   -1.255433 |
482. | 728   -1.478607 |
483. | 673   -1.681832 |
484. | 544   -1.268922 |
     +-----------------+
487. | 585   -1.191651 |
490. | 670   -1.275334 |
500. | 120   -1.211728 |
503. | 743   -1.455113 |
504. | 428   -1.512939 |
     +-----------------+
512. | 582   -1.411956 |
525. | 789   -1.442744 |
531. | 693   -1.433851 |
532. | 312   -1.525312 |
565. | 643   -1.552809 |
     +-----------------+
570. | 506   -1.410424 |
580. | 644   -1.625595 |
593. | 663   -1.256306 |
596. | 638    -1.60067 |
608. | 583   -1.217088 |
     +-----------------+
612. | 163   -1.296646 |
635. | 314   -1.613349 |
650. | 218   -1.188731 |
664. | 214    -1.57186 |
669. | 318   -1.383579 |
     +-----------------+
673. | 136   -1.188209 |
686. | 422   -1.417933 |
691. |  90   -1.198765 |
     +-----------------+
(702 real changes made, 702 to missing)
    3
   38
(4 real changes made)
   42

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

50%            0                      Mean            .005698
                        Largest       Std. Dev.      .0753234
75%            0              1
90%            0              1       Variance       .0056736
95%            0              1       Skewness       13.13414
99%            0              1       Kurtosis       173.5057
(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 rFslp1std

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

 -16* | 8310
 -15* | 7531
 -14* | 9864432111
 -13* | 860
 -12* | 8877665321100
 -11* | 99999866665443333332222111000
 -10* | 999987776654333322211000000
  -9* | 999866666555544443333222111111111
  -8* | 998888876666655555444332222221111110000
  -7* | 999888777665554444433222211111100000
  -6* | 9999988887777776666665555554444443332111100
  -5* | 98888777777666544444333322111100
  -4* | 999887776655544432222111111
  -3* | 999987777665554444333332110000
  -2* | 999999888766655554444444321111000
  -1* | 9997665444433222222111110
  -0* | 99987777666555444433333222
   0* | 000123333555666778999
   1* | 00001144555556666788999
   2* | 1122333455556778889999
   3* | 011122233334445667889
   4* | 0112223466799
   5* | 01123444667899
   6* | 00123446778899
   7* | 11112233667789
   8* | 001122234666678
   9* | 000001234567889
  10* | 001334788889
  11* | 012234579
  12* | 123457
  13* | 112333446778
  14* | 145577
  15* | 00012236888
  16* | 145888
  17* | 17
  18* | 0002234689
  19* | 0333
  20* | 012234
  21* | 1344579
  22* | 01556
  23* | 3569
  24* | 46
  25* | 122
  26* | 1
  27* | 117
  28* |
  29* | 6
  30* |
  31* | 1
  32* | 7
  33* |
  34* |
  35* | 01
  36* |
  37* |
  38* |
  39* |
  40* |
  41* |
  42* |
  43* |
  44* | 1


Computing the leverage h

                          rFslp1std
-------------------------------------------------------------
      Percentiles      Smallest
 1%    -1.512939      -1.681832
 5%     -1.19017      -1.625595
10%    -1.072334      -1.613349       Obs                 701
25%    -.7491144       -1.60067       Sum of Wgt.         701

50%    -.2379604                      Mean          -4.00e-10
                        Largest       Std. Dev.             1
75%     .5822632       3.269424
90%     1.473866       3.503309       Variance              1
95%     2.006879       3.507825       Skewness       .9727678
99%     2.711857       4.407042       Kurtosis       3.693234

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

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

     +-----+
     |  id |
     +-----+
 55. | 465 |
173. | 266 |
333. | 535 |
     +-----+
    3
Negative outliers are:

     +-----+
     |  id |
     +-----+
 16. | 122 |
113. | 304 |
134. | 341 |
142. | 421 |
185. | 514 |
     +-----+
234. |  93 |
240. | 492 |
241. | 269 |
243. | 326 |
452. | 213 |
     +-----+
465. | 671 |
475. | 640 |
482. | 728 |
483. | 673 |
484. | 544 |
     +-----+
487. | 585 |
490. | 670 |
500. | 120 |
503. | 743 |
504. | 428 |
     +-----+
512. | 582 |
525. | 789 |
531. | 693 |
532. | 312 |
565. | 643 |
     +-----+
570. | 506 |
580. | 644 |
593. | 663 |
596. | 638 |
608. | 583 |
     +-----+
612. | 163 |
635. | 314 |
650. | 218 |
664. | 214 |
669. | 318 |
     +-----+
673. | 136 |
686. | 422 |
691. |  90 |
     +-----+
   38

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

    Variable |       VIF       1/VIF
-------------+----------------------
   BSIposymp |    346.17    0.002889
    BSItotal |    328.92    0.003040
      crhtw2 |     12.97    0.077103
      crhtw3 |     10.53    0.094985
      BSIanx |      3.74    0.267321
      BSIdep |      3.35    0.298949
      crhtw1 |      2.45    0.408345
     WHPpain |      2.16    0.463739
     BFdepx1 |      2.08    0.480554
  HP2probsoc |      1.72    0.582832
         age |      1.52    0.656437
     radchw3 |      1.34    0.745517
      pillw2 |      1.14    0.880010
    icdx3nr5 |      1.09    0.916738
    icdx5nr4 |      1.09    0.920349
    icdx4nr6 |      1.08    0.928874
   bffpain24 |      1.07    0.934687
      sepaw2 |      1.04    0.962940
   bffpain26 |      1.04    0.964867
    icdx6nr5 |      1.03    0.969825
-------------+----------------------
    Mean VIF |     36.28

Condition number using scaled variables =       249.04

Condition Indexes and Variance-Decomposition Proportions

condition
    index      _cons        age     sepaw2     pillw2    radchw3    WHPpain  bffpain24  bffpain26     BSIanx HP2probsoc  BSIposymp
1    1.00       0.00       0.00       0.00       0.00       0.00       0.00       0.00       0.00       0.00       0.00       0.00
2    1.94       0.00       0.00       0.00       0.00       0.00       0.00       0.00       0.00       0.00       0.00       0.00
3    2.79       0.00       0.00       0.15       0.08       0.00       0.01       0.27       0.09       0.00       0.04       0.00
4    2.81       0.00       0.00       0.09       0.05       0.00       0.00       0.00       0.22       0.00       0.00       0.00
5    2.90       0.00       0.00       0.09       0.01       0.00       0.00       0.08       0.02       0.00       0.00       0.00
6    2.94       0.00       0.00       0.10       0.00       0.00       0.00       0.04       0.00       0.00       0.00       0.00
7    2.94       0.00       0.00       0.11       0.00       0.00       0.00       0.08       0.00       0.00       0.01       0.00
8    3.05       0.00       0.00       0.04       0.02       0.00       0.00       0.17       0.16       0.00       0.00       0.00
9    3.19       0.00       0.00       0.36       0.02       0.00       0.00       0.04       0.41       0.00       0.01       0.00
10   3.38       0.00       0.00       0.01       0.00       0.01       0.06       0.14       0.06       0.00       0.15       0.00
11   3.52       0.00       0.00       0.00       0.73       0.00       0.00       0.12       0.00       0.00       0.01       0.00
12   4.29       0.00       0.00       0.01       0.02       0.00       0.02       0.02       0.01       0.00       0.23       0.00
13   4.94       0.00       0.00       0.01       0.00       0.05       0.05       0.00       0.00       0.00       0.32       0.00
14   6.11       0.00       0.00       0.01       0.02       0.00       0.63       0.02       0.01       0.00       0.03       0.00
15   6.93       0.00       0.00       0.00       0.01       0.61       0.10       0.00       0.01       0.01       0.03       0.00
16  11.78       0.07       0.23       0.00       0.00       0.23       0.00       0.00       0.00       0.04       0.03       0.00
17  13.22       0.00       0.00       0.00       0.00       0.00       0.02       0.01       0.00       0.30       0.05       0.00
18  14.57       0.00       0.00       0.01       0.01       0.00       0.01       0.00       0.00       0.10       0.02       0.00
19  22.08       0.59       0.65       0.00       0.00       0.05       0.03       0.00       0.00       0.19       0.03       0.00
20  24.94       0.33       0.11       0.00       0.01       0.02       0.04       0.01       0.00       0.34       0.02       0.00
21 249.04       0.00       0.01       0.00       0.01       0.01       0.02       0.00       0.00       0.02       0.00       0.99

condition
    index   BSItotal     BSIdep     crhtw1     crhtw2     crhtw3    BFdepx1   icdx3nr5   icdx4nr6   icdx5nr4   icdx6nr5
1    1.00       0.00       0.00       0.00       0.00       0.00       0.00       0.00       0.00       0.00       0.00
2    1.94       0.00       0.00       0.04       0.01       0.01       0.00       0.00       0.00       0.00       0.00
3    2.79       0.00       0.00       0.00       0.00       0.00       0.00       0.02       0.03       0.00       0.00
4    2.81       0.00       0.00       0.01       0.00       0.00       0.00       0.15       0.19       0.09       0.01
5    2.90       0.00       0.00       0.01       0.00       0.00       0.00       0.23       0.00       0.34       0.06
6    2.94       0.00       0.00       0.00       0.00       0.00       0.00       0.12       0.19       0.37       0.06
7    2.94       0.00       0.00       0.00       0.00       0.00       0.00       0.00       0.00       0.00       0.72
8    3.05       0.00       0.00       0.00       0.00       0.00       0.00       0.01       0.45       0.04       0.03
9    3.19       0.00       0.00       0.00       0.00       0.00       0.00       0.15       0.02       0.00       0.03
10   3.38       0.00       0.00       0.01       0.00       0.00       0.05       0.07       0.00       0.01       0.06
11   3.52       0.00       0.00       0.00       0.00       0.00       0.01       0.18       0.01       0.00       0.00
12   4.29       0.00       0.00       0.24       0.01       0.03       0.04       0.01       0.01       0.09       0.01
13   4.94       0.00       0.00       0.37       0.00       0.01       0.12       0.01       0.04       0.00       0.00
14   6.11       0.00       0.00       0.00       0.00       0.00       0.47       0.00       0.02       0.02       0.00
15   6.93       0.00       0.02       0.01       0.00       0.00       0.09       0.02       0.00       0.00       0.00
16  11.78       0.00       0.07       0.08       0.00       0.03       0.17       0.01       0.01       0.00       0.00
17  13.22       0.00       0.31       0.00       0.16       0.17       0.00       0.01       0.00       0.02       0.00
18  14.57       0.00       0.08       0.22       0.80       0.73       0.02       0.00       0.00       0.00       0.00
19  22.08       0.00       0.10       0.00       0.00       0.00       0.00       0.01       0.01       0.00       0.00
20  24.94       0.01       0.41       0.00       0.00       0.00       0.03       0.00       0.01       0.01       0.01
21 249.04       0.99       0.00       0.00       0.00       0.00       0.00       0.01       0.00       0.00       0.00

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

                   Shapiro-Wilk W test for normal data

    Variable |    Obs       W           V         z       Prob>z
-------------+--------------------------------------------------
   rFslp1std |    701    0.93253     30.806     8.363    0.00000

                    Skewness/Kurtosis tests for Normality
                                                         ------- joint ------
    Variable |    Obs   Pr(Skewness)   Pr(Kurtosis)  adj chi2(2)    Prob>chi2
-------------+---------------------------------------------------------------
   rFslp1std |    701      0.0000         0.0026        71.28         0.0000
Jarque-Bera normality test:  124.6 Chi(2)  8.8e-28
Jarque-Bera test for Ho: normality:
(file qnrFslp1std.gph saved)
(file bxrFslp1std.gph saved)
(bin=26, start=-1.6818317, width=.23418745)
(file histrFslp1std.gph saved)
(file kdrFslp1std.gph saved)
(file combi4.gph saved)
-------------
Specification error tests:

Ramsey RESET test using powers of the fitted values of WHPsleep
       Ho:  model has no omitted variables
                 F(3, 338) =      7.41
                  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   |            3   |
                  +--------------------+--------------+
                  | Negative-outliers   |           38   |
                  +--------------------+--------------+
                  | Bad-outliers   |            1   |
                  +--------------------+--------------+
                  | collinearity   |            1   |
                  +--------------------+--------------+
                  | Shapiro-Wilk   |            1   |
                  +--------------------+--------------+
                  | Kolmogorov-Smirnov   |            1   |
                  +--------------------+--------------+
                  | Specification   |            1   |
                  +--------------------+--------------+
                  | StructuralBreak   |            0   |
                  +--------------------+--------------+
                  | TotalProbs   |            4   |
                  +--------------------+--------------+
                  | Percent-Invalid   |           80   |
                  +--------------------+--------------+
                  | Pct-Pseudo-Valid   |           20   |
                  +--------------------+--------------+

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
-------------

.
. regress WHPsleep age sepaw2   pillw2  radchw3 WHPpain bffpain24 bffpain26   ///
>   BSIanx HP2probsoc BSIposymp BSItotal BSIdep crhtw1 crhtw2 crhtw3  ///
>   BFdepx1 icdx3nr5 icdx4nr6 icdx5nr4 icdx6nr5 ///
>    if gender==2, vce(cluster id)

Linear regression                                      Number of obs =     362
                                                       F( 19,   361) =       .
                                                       Prob > F      =       .
                                                       R-squared     =  0.6054
                                                       Root MSE      =  19.996

(Std. Err. adjusted for 362 clusters in id)
-------------+----------------------------------------------------------------
             |               Robust
     WHPsleep |      Coef.    Std. Err.       t    P>|t|      [95% Con f. Interval]
-------------+----------------------------------------------------------------
         age |    .2748115    .1074232      2.56    0.011      .0635576     .4860654
      sepaw2 |   -18.77988    4.034896     -4.65    0.000     -26.71473    -10.84502
      pillw2 |   -.8019393    .3119536     -2.57    0.011     -1.415414    -.1884648
     radchw3 |   -.0641173    .0310957     -2.06    0.040     -.1252688    -.0029658
     WHPpain |      .21312    .0751605      2.84    0.005      .0653127     .3609274
   bffpain24 |    -5.31569    .9930722     -5.35    0.000     -7.268624    -3.362757
   bffpain26 |   -.2349557    .0961657     -2.44    0.015      -.424071    -.0458403
      BSIanx |    .3524595    .6178581      0.57    0.569     -.8625938     1.567513
  HP2probsoc |    6.639116    3.844953      1.73    0.085     -.9222049     14.20044
   BSIposymp |    2.501887    .7092851      3.53    0.000      1.107037     3.896737
    BSItotal |   -2.168287    .6649841     -3.26    0.001     -3.476016    -.8605573
      BSIdep |   -.7581339     .567056     -1.34    0.182     -1.873282      .357014
      crhtw1 |    3.307012    1.606036      2.06    0.040      .1486505     6.465373
      crhtw2 |   -1.407417    3.559863     -0.40    0.693     -8.408091     5.593257
      crhtw3 |    2.090788    3.264265      0.64    0.522     -4.328575     8.510152
     BFdepx1 |    .6379949    .1054825      6.05    0.000      .4305576     .8454322
    icdx3nr5 |    10.05342    4.054206      2.48    0.014      2.080598     18.02625
    icdx4nr6 |    20.81846    10.55316      1.97    0.049      .0650634     41.57185
    icdx5nr4 |    50.87478     6.54591      7.77    0.000      38.00187     63.74769
    icdx6nr5 |    58.50618    14.84362      3.94    0.000      29.31535     87.69701
       _cons |   -11.47416    6.993127     -1.64    0.102     -25.22654     2.278222
-------------+----------------------------------------------------------------

.    di e(r2_a)
.58224747

.
.
. // prep for residual and prediction interval analysis
.     predict H8FemaleSLpred if gender==2, xb
(340 missing values generated)

.       predict h8FSLsepred if gender==2, stdp
(341 missing values generated)

.      predict FSLrsires if gender==2, residual
(340 missing values generated)

.      gen upbSLF = H8FemaleSLpred + 1.96*h8FSLsepred
(341 missing values generated)

.      gen lpbSLF = H8FemaleSLpred - 1.96*h8FSLsepred
(341 missing values generated)

.  // residual diagnostics
.
. // prediction interval
. scatter H8FemaleSLpred FSLrsires || lowess H8FemaleSLpred  FSLrsires  ///
>   || lowess upbSLF FSLrsires || lowess lpbSLF FSLrsires, ///
>   title(Prediction interval of Female Sleep model)   ///
>   ytitle(Predicted Female Sleep )

.
. // saving graph
. gr save FSLPredInt.gph, replace
(file FSLPredInt.gph saved)

. gr export FSLPredInt.eps, replace
(file FSLPredInt.eps written in EPS format)

. gr use  FSLPredInt.gph

.
. title "Moderation effects testing for female sleep owing to significant main eff"


*******************************************************************************
*******************************************************************************
*****                                                                     *****
*****                                                                     *****
*****Moderation effects testing for female sleep owing to significant main eff*****
*****                                                                     *****
*****                                                                     *****
*****                                            12 Jul 2012    23:56:50  *****
*******************************************************************************
*******************************************************************************


.
. foreach varx in age occ5w3 sepaw2   pillw2  radchw3 WHPpain bffpain24 bffpain26   ///
>   BSIanx HP2probsoc BSIposymp BSItotal BSIdep WHPpa defnw2  ///
>   BFdepx1 icdx3nr5 icdx4nr6 icdx4nr10 icdx5nr4 icdx5nr7 icdx6nr5 {
  2.  cap gen `varx'Xcrht1 = `varx'*crhtw1
  3.  cap gen `varx'Xcrht2 = `varx'*crhtw2
  4.  cap gen `varx'Xcrht3 = `varx'*crhtw3
  5.  }

.
.  // Full moderation model
.  set more off

.  regress WHPsleep age sepaw2   pillw2   ///
>   BSIanx HP2probsoc BSIposymp BSItotal BSIdep WHPpain crhtw1 crhtw2 crhtw3 radchw3 ///
>   icdx3nr5 icdx4nr6 icdx5nr4 icdx6nr5   BFdepx1  bffpain24 bffpain26 ///
>    ageXcrht1-BFdepx1Xcrht3 if gender==2, vce(cluster id)

Linear regression                                      Number of obs =     362
                                                       F( 67,   361) =       .
                                                       Prob > F      =       .
                                                       R-squared     =  0.6719
                                                       Root MSE      =  19.669

(Std. Err. adjusted for 362 clusters in id)
-----------------+----------------------------------------------------------------
                 |               Robust
         WHPsleep |      Coef.    Std. Err.       t    P>|t|      [95% Con f. Interval]
-----------------+----------------------------------------------------------------
             age |    .2190285     .111171      1.97    0.050      .0004043     .4376526
          sepaw2 |   -5.709317     5.82363     -0.98    0.328     -17.16182     5.743183
          pillw2 |   -.5678966    .2708668     -2.10    0.037     -1.100572    -.0352216
          BSIanx |   -.0851406    .6268447     -0.14    0.892     -1.317866     1.147585
      HP2probsoc |    1.398479    3.768446      0.37    0.711     -6.012385     8.809344
       BSIposymp |    3.160442     .760845      4.15    0.000      1.664197     4.656687
        BSItotal |   -2.664706    .7094479     -3.76    0.000     -4.059876    -1.269536
          BSIdep |   -1.109565    .5354898     -2.07    0.039     -2.162637     -.056494
         WHPpain |    .2222525    .0808154      2.75    0.006      .0633243     .3811806
          crhtw1 |   -5.099438    10.93025     -0.47    0.641      -26.5944     16.39552
          crhtw2 |    21.87027    27.66483      0.79    0.430     -32.53421     76.27474
          crhtw3 |   -10.85746    23.81684     -0.46    0.649     -57.69462     35.97971
         radchw3 |    -.071969    .0359809     -2.00    0.046     -.1427276    -.0012105
        icdx3nr5 |    8.039609    4.995557      1.61    0.108      -1.78444     17.86366
        icdx4nr6 |    25.53316    11.13154      2.29    0.022      3.642362     47.42396
        icdx5nr4 |    27.36853    15.63218      1.75    0.081     -3.373048     58.11011
        icdx6nr5 |    61.37043    17.37651      3.53    0.000      27.19853     95.54233
         BFdepx1 |    .6307155    .1130835      5.58    0.000      .4083303     .8531008
       bffpain24 |   -2.998509    1.148595     -2.61    0.009     -5.257287    -.7397304
       bffpain26 |   -.3498598    .0866307     -4.04    0.000     -.5202241    -.1794955
       ageXcrht1 |    .3789823    .1747765      2.17    0.031      .0352743     .7226903
       ageXcrht2 |   -.6395828    .5211399     -1.23    0.221     -1.664434     .3852686
       ageXcrht3 |    .4062316    .4534443      0.90    0.371     -.4854925     1.297956
    occ5w3Xcrht1 |   -24.33741    9.764502     -2.49    0.013     -43.53986    -5.134961
    occ5w3Xcrht2 |   -1099.873    288.8765     -3.81    0.000     -1667.965    -531.7808
    occ5w3Xcrht3 |    1118.276    296.4648      3.77    0.000      535.2607     1701.291
    sepaw2Xcrht1 |   -18.78973    16.03454     -1.17    0.242     -50.32257     12.74311
    sepaw2Xcrht2 |    -14.1871    20.48524     -0.69    0.489      -54.4725      26.0983
    sepaw2Xcrht3 |    32.49299    19.24809      1.69    0.092     -5.359473     70.34545
    pillw2Xcrht1 |   -.1430141    .7918204     -0.18    0.857     -1.700174     1.414146
    pillw2Xcrht2 |   -1.128626    1.691112     -0.67    0.505     -4.454293     2.197041
    pillw2Xcrht3 |    .3988188    1.017318      0.39    0.695     -1.601795     2.399432
   radchw3Xcrht1 |   -.0360494    .0613725     -0.59    0.557     -.1567419     .0846432
   radchw3Xcrht2 |   -.0365228    .1193631     -0.31    0.760     -.2712572     .1982116
   radchw3Xcrht3 |    .0473858    .1083439      0.44    0.662     -.1656786     .2604502
   WHPpainXcrht1 |    .0581803    .1544886      0.38    0.707     -.2456304      .361991
   WHPpainXcrht2 |    .2116283    .4800741      0.44    0.660     -.7324648     1.155721
   WHPpainXcrht3 |   -.2994523    .4151106     -0.72    0.471     -1.115791     .5168864
 bffpain24Xcrht1 |   -14.78661      4.9237     -3.00    0.003     -24.46934    -5.103871
 bffpain24Xcrht2 |     15.1307    8.986859      1.68    0.093     -2.542474     32.80387
 bffpain24Xcrht3 |   -2.166511    7.050661     -0.31    0.759     -16.03204     11.69902
 bffpain26Xcrht1 |    .0000238    .1107514      0.00    1.000      -.217775     .2178227
 bffpain26Xcrht2 |   -.2129464    .2902876     -0.73    0.464     -.7838134     .3579207
 bffpain26Xcrht3 |   -.2304548    .2137754     -1.08    0.282     -.6508563     .1899466
    BSIanxXcrht1 |   -.3551488    .6906359     -0.51    0.607     -1.713324     1.003026
    BSIanxXcrht2 |   -.8744168    2.026639     -0.43    0.666     -4.859918     3.111084
    BSIanxXcrht3 |    2.481644    1.854216      1.34    0.182     -1.164778     6.128065
HP2probsocXcrht1 |   -6.556147    6.896054     -0.95    0.342     -20.11763     7.005336
HP2probsocXcrht2 |    4.356438    17.68665      0.25    0.806     -30.42536     39.13824
HP2probsocXcrht3 |    8.692857    14.40531      0.60    0.547       -19.636     37.02172
 BSIposympXcrht1 |    1.023235    .9538997      1.07    0.284     -.8526627     2.899134
 BSIposympXcrht2 |   -2.210221    3.150975     -0.70    0.483     -8.406793     3.986352
 BSIposympXcrht3 |    .4204572    2.986033      0.14    0.888     -5.451748     6.292662
  BSItotalXcrht1 |   -.9182353    .9196814     -1.00    0.319     -2.726841     .8903707
  BSItotalXcrht2 |    1.632888    2.924043      0.56    0.577      -4.11741     7.383187
  BSItotalXcrht3 |   -.2279941    2.774567     -0.08    0.935     -5.684338      5.22835
    BSIdepXcrht1 |   -1.123317    .7545822     -1.49    0.137     -2.607246      .360612
    BSIdepXcrht2 |    4.946507    2.587706      1.91    0.057     -.1423632     10.03538
    BSIdepXcrht3 |   -3.711385    2.280149     -1.63    0.104     -8.195429     .7726576
     WHPpaXcrht1 |   -.2186439    .1613125     -1.36    0.176     -.5358741     .0985863
     WHPpaXcrht2 |   -.0838338    .5094715     -0.16    0.869     -1.085739      .918071
     WHPpaXcrht3 |    .3427543    .4309854      0.80    0.427     -.5048031     1.190312
    defnw2Xcrht1 |     .040829    .0543161      0.75    0.453     -.0659866     .1476446
    defnw2Xcrht2 |    .0418167    .1315831      0.32    0.751      -.216949     .3005824
    defnw2Xcrht3 |   -.1078665    .1143788     -0.94    0.346      -.332799     .1170661
   BFdepx1Xcrht1 |    .2956623    .1355211      2.18    0.030      .0291525     .5621722
   BFdepx1Xcrht2 |    .3709545    .3133792      1.18    0.237     -.2453235     .9872325
   BFdepx1Xcrht3 |   -.4191597    .3033073     -1.38    0.168     -1.015631     .1773114
           _cons |   -11.08795    8.074243     -1.37    0.171     -26.96641     4.790511
-----------------+----------------------------------------------------------------

.    di e(r2_a)
.59579018

.
. // Trimmed moderation model
. set more off

.  regress WHPsleep age sepaw2   pillw2  radchw3 WHPpain bffpain24 bffpain26   ///
>   BSIanx HP2probsoc BSIposymp BSItotal BSIdep crhtw1 crhtw2 crhtw3  ///
>   BFdepx1 icdx3nr5 icdx4nr6 icdx5nr4 icdx6nr5 ///
>   occ5w3Xcrht1-occ5w3Xcrht3 pillw2Xcrht2  ///
>   WHPpaXcrht3 icdx4nr10Xcrht1 icdx4nr10Xcrht3 ///
>   icdx5nr7Xcrht1  if gender==2, vce(cluster id)

Linear regression                                      Number of obs =     362
                                                       F( 24,   361) =       .
                                                       Prob > F      =       .
                                                       R-squared     =  0.6357
                                                       Root MSE      =  19.443

(Std. Err. adjusted for 362 clusters in id)
----------------+----------------------------------------------------------------
                |               Robust
        WHPsleep |      Coef.    Std. Err.       t    P>|t|      [95% Con f. Interval]
----------------+----------------------------------------------------------------
            age |    .2219826    .1051277      2.11    0.035      .0152429     .4287223
         sepaw2 |   -17.55504    3.656245     -4.80    0.000     -24.74526    -10.36483
         pillw2 |   -.7379301    .2336624     -3.16    0.002     -1.197441    -.2784196
        radchw3 |    -.062173    .0306266     -2.03    0.043     -.1224019    -.0019442
        WHPpain |    .1937617    .0765877      2.53    0.012      .0431477     .3443757
      bffpain24 |   -4.256477    1.018012     -4.18    0.000     -6.258456    -2.254499
      bffpain26 |   -.1903348    .0920855     -2.07    0.039     -.3714262    -.0092434
         BSIanx |    .3657967    .6135035      0.60    0.551     -.8406929     1.572286
     HP2probsoc |    6.143407     3.83199      1.60    0.110     -1.392421     13.67923
      BSIposymp |    2.312992    .7089218      3.26    0.001      .9188564     3.707127
       BSItotal |   -1.961028    .6650981     -2.95    0.003     -3.268981    -.6530747
         BSIdep |   -.9832378    .5841257     -1.68    0.093     -2.131954     .1654789
         crhtw1 |    3.518823    1.594712      2.21    0.028      .3827302     6.654916
         crhtw2 |   -1.309479    3.472693     -0.38    0.706     -8.138727      5.51977
         crhtw3 |    .5895341    3.174691      0.19    0.853     -5.653677     6.832745
        BFdepx1 |    .6391151    .1023493      6.24    0.000      .4378393     .8403909
       icdx3nr5 |    9.160697    4.197352      2.18    0.030      .9063645     17.41503
       icdx4nr6 |    26.47271    9.599979      2.76    0.006      7.593799     45.35161
       icdx5nr4 |    41.73646    7.238337      5.77    0.000      27.50186     55.97107
       icdx6nr5 |     60.5948     17.3984      3.48    0.001      26.37986     94.80974
   occ5w3Xcrht1 |   -27.70124    8.854168     -3.13    0.002     -45.11347    -10.28901
   occ5w3Xcrht2 |   -1012.671    298.6315     -3.39    0.001     -1599.947    -425.3953
   occ5w3Xcrht3 |     1033.58    307.8105      3.36    0.001      428.2537     1638.907
   pillw2Xcrht2 |    -.920573    .2508151     -3.67    0.000     -1.413815    -.4273307
    WHPpaXcrht3 |    .1660406    .0551336      3.01    0.003      .0576173     .2744639
icdx4nr10Xcrht1 |   -150.9534     16.4147     -9.20    0.000     -183.2339     -118.673
icdx4nr10Xcrht3 |   -202.0248     15.0301    -13.44    0.000     -231.5823    -172.4672
 icdx5nr7Xcrht1 |    64.20974    14.03834      4.57    0.000      36.60254     91.81694
          _cons |   -10.15555    6.943583     -1.46    0.144      -23.8105     3.499408
----------------+----------------------------------------------------------------

.    di e(r2_a)
.60503813

. *-+- possible female moderators of sleep include occ5w3 pillw2 pillw3 defnw2 WHPpa
. *-+- and icdx4nr10 and icdx5nr7
.
.
. title "tests for possible indirect effects for female sleep model"


*******************************************************************************
*******************************************************************************
*****                                                                     *****
*****                                                                     *****
*****     tests for possible indirect effects for female sleep model      *****
*****                                                                     *****
*****                                                                     *****
*****                                            12 Jul 2012    23:56:50  *****
*******************************************************************************
*******************************************************************************


. set more off

. foreach var in  age sepaw2   pillw2  radchw3 WHPpain bffpain24 bffpain26   ///
>   BSIanx HP2probsoc BSIposymp BSItotal BSIdep  ///
>   BFdepx1 icdx3nr5 icdx4nr6 icdx5nr4 icdx6nr5 {
  2. sem (crhtw1-> `var')(`var'-> whpsleep) if gender==2, nocapslatent iterate(50)
  3. sem (crhtw2-> `var')(`var'-> whpsleep) if gender==2, nocapslatent iterate(50)
  4. sem (crhtw3-> `var')(`var'-> whpsleep) if gender==2, nocapslatent iterate(50)
  5. }
(1 observations with missing values excluded;
 specify option 'method(mlmv)' to use all observations)

Endogenous variables


Observed:  age whpsleep


Exogenous variables


Observed:  crhtw1


Fitting target model:

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

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

--------------+----------------------------------------------------------------
              |                 OIM
              |      Coef.    Std. Err.       z    P>|z|      [95% Con f. Interval]
--------------+----------------------------------------------------------------
 Structural    |
   age <-      |
       crhtw1 |    3.218873    .6364183      5.06    0.000      1.971516      4.46623
        _cons |    49.84232    .6063362     82.20    0.000      48.65393     51.03072
  ------------+----------------------------------------------------------------
   whpsleep <- |
          age |    1.000626    .1269109      7.88    0.000      .7518852     1.249367
        _cons |   -24.01976    6.553885     -3.66    0.000     -36.86514    -11.17438
--------------+----------------------------------------------------------------
 Variance      |
        e.age |      130.49    9.699243                    112.7997     150.9546
   e.whpsleep |    814.5879    60.54784                    704.1556     942.3392
--------------+----------------------------------------------------------------
LR test of model vs. saturated: chi2(1)   =      0.93, Prob > chi2 = 0.3346

Endogenous variables


Observed:  age whpsleep


Exogenous variables


Observed:  crhtw2


Fitting target model:

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

Structural equation model                       Number of obs      =       363
 Estimation method = ml
 Log likelihood    = -3592.6881

--------------+----------------------------------------------------------------
              |                 OIM
              |      Coef.    Std. Err.       z    P>|z|      [95% Con f. Interval]
--------------+----------------------------------------------------------------
 Structural    |
   age <-      |
       crhtw2 |    3.897443    .6851925      5.69    0.000      2.554491     5.240396
        _cons |    49.54513    .6075589     81.55    0.000      48.35434     50.73593
  ------------+----------------------------------------------------------------
   whpsleep <- |
          age |    .9980208    .1262348      7.91    0.000      .7506051     1.245436
        _cons |   -23.87083    6.513018     -3.67    0.000     -36.63611    -11.10555
--------------+----------------------------------------------------------------
 Variance      |
        e.age |    128.9608     9.57237                    111.5002     149.1557
   e.whpsleep |    812.4615    60.30655                    702.4585     939.6906
--------------+----------------------------------------------------------------
LR test of model vs. saturated: chi2(1)   =      9.45, Prob > chi2 = 0.0021

Endogenous variables


Observed:  age whpsleep


Exogenous variables


Observed:  crhtw3


Fitting target model:

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

Structural equation model                       Number of obs      =       363
 Estimation method = ml
 Log likelihood    = -3602.2653

--------------+----------------------------------------------------------------
              |                 OIM
              |      Coef.    Std. Err.       z    P>|z|      [95% Con f. Interval]
--------------+----------------------------------------------------------------
 Structural    |
   age <-      |
       crhtw3 |    3.646918    .6723377      5.42    0.000       2.32916     4.964675
        _cons |    49.58731    .6093473     81.38    0.000      48.39301     50.78161
  ------------+----------------------------------------------------------------
   whpsleep <- |
          age |    .9980208    .1262348      7.91    0.000      .7506051     1.245436
        _cons |   -23.87083    6.513018     -3.67    0.000     -36.63611    -11.10555
--------------+----------------------------------------------------------------
 Variance      |
        e.age |    129.9244    9.643895                    112.3333     150.2702
   e.whpsleep |    812.4615    60.30655                    702.4585     939.6906
--------------+----------------------------------------------------------------
LR test of model vs. saturated: chi2(1)   =     13.57, Prob > chi2 = 0.0002
(1 observations with missing values excluded;
 specify option 'method(mlmv)' to use all observations)

Endogenous variables


Observed:  sepaw2 whpsleep


Exogenous variables


Observed:  crhtw1


Fitting target model:

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

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

--------------+----------------------------------------------------------------
              |                 OIM
              |      Coef.    Std. Err.       z    P>|z|      [95% Con f. Interval]
--------------+----------------------------------------------------------------
 Structural    |
   sepaw2 <-   |
       crhtw1 |    .0142735    .0081557      1.75    0.080     -.0017115     .0302584
        _cons |    .0201998    .0077702      2.60    0.009      .0049705     .0354292
  ------------+----------------------------------------------------------------
   whpsleep <- |
       sepaw2 |   -23.66287    10.97537     -2.16    0.031      -45.1742    -2.151538
        _cons |    26.80537    1.631586     16.43    0.000      23.60752     30.00322
--------------+----------------------------------------------------------------
 Variance      |
     e.sepaw2 |    .0214297    .0015929                    .0185245     .0247906
   e.whpsleep |    942.3734    70.04606                    814.6174     1090.165
--------------+----------------------------------------------------------------
LR test of model vs. saturated: chi2(1)   =      8.87, Prob > chi2 = 0.0029

Endogenous variables


Observed:  sepaw2 whpsleep


Exogenous variables


Observed:  crhtw2


Fitting target model:

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

Structural equation model                       Number of obs      =       363
 Estimation method = ml
 Log likelihood    = -2040.4802

--------------+----------------------------------------------------------------
              |                 OIM
              |      Coef.    Std. Err.       z    P>|z|      [95% Con f. Interval]
--------------+----------------------------------------------------------------
 Structural    |
   sepaw2 <-   |
       crhtw2 |    .0065809    .0088513      0.74    0.457     -.0107673     .0239291
        _cons |    .0209077    .0078484      2.66    0.008      .0055251     .0362903
  ------------+----------------------------------------------------------------
   whpsleep <- |
       sepaw2 |   -23.62277    10.96315     -2.15    0.031     -45.11014    -2.135396
        _cons |    26.76527    1.627522     16.45    0.000      23.57538     29.95515
--------------+----------------------------------------------------------------
 Variance      |
     e.sepaw2 |    .0215201    .0015974                    .0186064     .0248901
   e.whpsleep |     940.334    69.79814                    813.0178     1087.588
--------------+----------------------------------------------------------------
LR test of model vs. saturated: chi2(1)   =     24.96, Prob > chi2 = 0.0000

Endogenous variables


Observed:  sepaw2 whpsleep


Exogenous variables


Observed:  crhtw3


Fitting target model:

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

Structural equation model                       Number of obs      =       363
 Estimation method = ml
 Log likelihood    = -2048.5056

--------------+----------------------------------------------------------------
              |                 OIM
              |      Coef.    Std. Err.       z    P>|z|      [95% Con f. Interval]
--------------+----------------------------------------------------------------
 Structural    |
   sepaw2 <-   |
       crhtw3 |    .0084506    .0086482      0.98    0.328     -.0084995     .0254007
        _cons |    .0205844    .0078379      2.63    0.009      .0052223     .0359464
  ------------+----------------------------------------------------------------
   whpsleep <- |
       sepaw2 |   -23.62277    10.96315     -2.15    0.031     -45.11014    -2.135396
        _cons |    26.76527    1.627522     16.45    0.000      23.57538     29.95515
--------------+----------------------------------------------------------------
 Variance      |
     e.sepaw2 |    .0214963    .0015956                    .0185858     .0248626
   e.whpsleep |     940.334    69.79814                    813.0178     1087.588
--------------+----------------------------------------------------------------
LR test of model vs. saturated: chi2(1)   =     30.30, Prob > chi2 = 0.0000
(1 observations with missing values excluded;
 specify option 'method(mlmv)' to use all observations)

Endogenous variables


Observed:  pillw2 whpsleep


Exogenous variables


Observed:  crhtw1


Fitting target model:

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

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

--------------+----------------------------------------------------------------
              |                 OIM
              |      Coef.    Std. Err.       z    P>|z|      [95% Con f. Interval]
--------------+----------------------------------------------------------------
 Structural    |
   pillw2 <-   |
       crhtw1 |    .2025624    .2097325      0.97    0.334     -.2085058     .6136306
        _cons |    1.591826     .199819      7.97    0.000      1.200188     1.983464
  ------------+----------------------------------------------------------------
   whpsleep <- |
       pillw2 |    .0227273    .4307801      0.05    0.958     -.8215861     .8670407
        _cons |    26.24564    1.767182     14.85    0.000      22.78203     29.70925
--------------+----------------------------------------------------------------
 Variance      |
     e.pillw2 |    14.17175    1.053378                    12.25051      16.3943
   e.whpsleep |    954.4668    70.94496                    825.0713     1104.155
--------------+----------------------------------------------------------------
LR test of model vs. saturated: chi2(1)   =      7.55, Prob > chi2 = 0.0060

Endogenous variables


Observed:  pillw2 whpsleep


Exogenous variables


Observed:  crhtw2


Fitting target model:

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

Structural equation model                       Number of obs      =       363
 Estimation method = ml
 Log likelihood    = -3219.7581

--------------+----------------------------------------------------------------
              |                 OIM
              |      Coef.    Std. Err.       z    P>|z|      [95% Con f. Interval]
--------------+----------------------------------------------------------------
 Structural    |
   pillw2 <-   |
       crhtw2 |    .3246167    .2265378      1.43    0.152     -.1193892     .7686226
        _cons |    1.558542    .2008706      7.76    0.000      1.164843     1.952242
  ------------+----------------------------------------------------------------
   whpsleep <- |
       pillw2 |    .0270174     .430193      0.06    0.950     -.8161454     .8701803
        _cons |    26.20104    1.762341     14.87    0.000      22.74692     29.65517
--------------+----------------------------------------------------------------
 Variance      |
     e.pillw2 |    14.09659    1.046347                    12.18798     16.30407
   e.whpsleep |    952.3509    70.69012                    823.4077     1101.486
--------------+----------------------------------------------------------------
LR test of model vs. saturated: chi2(1)   =     23.82, Prob > chi2 = 0.0000

Endogenous variables


Observed:  pillw2 whpsleep


Exogenous variables


Observed:  crhtw3


Fitting target model:

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

Structural equation model                       Number of obs      =       363
 Estimation method = ml
 Log likelihood    = -3228.4508

--------------+----------------------------------------------------------------
              |                 OIM
              |      Coef.    Std. Err.       z    P>|z|      [95% Con f. Interval]
--------------+----------------------------------------------------------------
 Structural    |
   pillw2 <-   |
       crhtw3 |    .2342381    .2217468      1.06    0.291     -.2003777      .668854
        _cons |    1.574017    .2009717      7.83    0.000       1.18012     1.967914
  ------------+----------------------------------------------------------------
   whpsleep <- |
       pillw2 |    .0270174     .430193      0.06    0.950     -.8161454     .8701803
        _cons |    26.20104    1.762341     14.87    0.000      22.74692     29.65517
--------------+----------------------------------------------------------------
 Variance      |
     e.pillw2 |    14.13288    1.049041                    12.21936     16.34605
   e.whpsleep |    952.3509    70.69012                    823.4077     1101.486
--------------+----------------------------------------------------------------
LR test of model vs. saturated: chi2(1)   =     28.62, Prob > chi2 = 0.0000
(1 observations with missing values excluded;
 specify option 'method(mlmv)' to use all observations)

Endogenous variables


Observed:  radchw3 whpsleep


Exogenous variables


Observed:  crhtw1


Fitting target model:

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

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

--------------+----------------------------------------------------------------
              |                 OIM
              |      Coef.    Std. Err.       z    P>|z|      [95% Con f. Interval]
--------------+----------------------------------------------------------------
 Structural    |
   radchw3 <-  |
       crhtw1 |    8.629612    1.874593      4.60    0.000      4.955477     12.30375
        _cons |    60.22442    1.785986     33.72    0.000      56.72396     63.72489
  ------------+----------------------------------------------------------------
   whpsleep <- |
      radchw3 |    .0047655    .0469045      0.10    0.919     -.0871657     .0966967
        _cons |    25.98996    3.305045      7.86    0.000      19.51219     32.46773
--------------+----------------------------------------------------------------
 Variance      |
    e.radchw3 |    1132.155    84.15244                    978.6709     1309.711
   e.whpsleep |    954.4469    70.94348                    825.0542     1104.132
--------------+----------------------------------------------------------------
LR test of model vs. saturated: chi2(1)   =      7.85, Prob > chi2 = 0.0051

Endogenous variables


Observed:  radchw3 whpsleep


Exogenous variables


Observed:  crhtw2


Fitting target model:

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

Structural equation model                       Number of obs      =       363
 Estimation method = ml
 Log likelihood    = -3986.9499

--------------+----------------------------------------------------------------
              |                 OIM
              |      Coef.    Std. Err.       z    P>|z|      [95% Con f. Interval]
--------------+----------------------------------------------------------------
 Structural    |
   radchw3 <-  |
       crhtw2 |    17.56613    1.875049      9.37    0.000       13.8911     21.24116
        _cons |    58.46075    1.662602     35.16    0.000      55.20211     61.71939
  ------------+----------------------------------------------------------------
   whpsleep <- |
      radchw3 |    .0035358    .0467727      0.08    0.940     -.0881371     .0952087
        _cons |    26.02728     3.30036      7.89    0.000      19.55869     32.49586
--------------+----------------------------------------------------------------
 Variance      |
    e.radchw3 |    965.7349    71.68357                    834.9796     1116.966
   e.whpsleep |    952.3463    70.68977                    823.4036     1101.481
--------------+----------------------------------------------------------------
LR test of model vs. saturated: chi2(1)   =     29.28, Prob > chi2 = 0.0000

Endogenous variables


Observed:  radchw3 whpsleep


Exogenous variables


Observed:  crhtw3


Fitting target model:

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

Structural equation model                       Number of obs      =       363
 Estimation method = ml
 Log likelihood    = -3992.8598

--------------+----------------------------------------------------------------
              |                 OIM
              |      Coef.    Std. Err.       z    P>|z|      [95% Con f. Interval]
--------------+----------------------------------------------------------------
 Structural    |
   radchw3 <-  |
       crhtw3 |    17.61708    1.821378      9.67    0.000      14.04725     21.18692
        _cons |    58.44778    1.650736     35.41    0.000       55.2124     61.68316
  ------------+----------------------------------------------------------------
   whpsleep <- |
      radchw3 |    .0035358    .0467727      0.08    0.940     -.0881371     .0952087
        _cons |    26.02728     3.30036      7.89    0.000      19.55869     32.49586
--------------+----------------------------------------------------------------
 Variance      |
    e.radchw3 |    953.4896    70.77464                    824.3922     1102.803
   e.whpsleep |    952.3463    70.68977                    823.4036     1101.481
--------------+----------------------------------------------------------------
LR test of model vs. saturated: chi2(1)   =     35.81, Prob > chi2 = 0.0000
(1 observations with missing values excluded;
 specify option 'method(mlmv)' to use all observations)

Endogenous variables


Observed:  WHPpain whpsleep


Exogenous variables


Observed:  crhtw1


Fitting target model:

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

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

--------------+----------------------------------------------------------------
              |                 OIM
              |      Coef.    Std. Err.       z    P>|z|      [95% Con f. Interval]
--------------+----------------------------------------------------------------
 Structural    |
   WHPpain <-  |
       crhtw1 |    1.729479    1.239418      1.40    0.163     -.6997351     4.158694
        _cons |    17.83115    1.180833     15.10    0.000      15.51676     20.14554
  ------------+----------------------------------------------------------------
   whpsleep <- |
      WHPpain |    .7680828    .0605752     12.68    0.000      .6493576      .886808
        _cons |    12.40984    1.738609      7.14    0.000      9.002228     15.81745
--------------+----------------------------------------------------------------
 Variance      |
    e.WHPpain |    494.9113    36.78647                     427.817     572.5279
   e.whpsleep |    660.9301    49.12655                     571.329     764.5834
--------------+----------------------------------------------------------------
LR test of model vs. saturated: chi2(1)   =      5.62, Prob > chi2 = 0.0177

Endogenous variables


Observed:  WHPpain whpsleep


Exogenous variables


Observed:  crhtw2


Fitting target model:

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

Structural equation model                       Number of obs      =       363
 Estimation method = ml
 Log likelihood    = -3790.3382

--------------+----------------------------------------------------------------
              |                 OIM
              |      Coef.    Std. Err.       z    P>|z|      [95% Con f. Interval]
--------------+----------------------------------------------------------------
 Structural    |
   WHPpain <-  |
       crhtw2 |    5.721746    1.311307      4.36    0.000      3.151631     8.291861
        _cons |    17.02834    1.162733     14.65    0.000      14.74942     19.30725
  ------------+----------------------------------------------------------------
   whpsleep <- |
      WHPpain |    .7680668    .0604372     12.71    0.000      .6496121     .8865215
        _cons |    12.41057    1.732256      7.16    0.000      9.015409     15.80573
--------------+----------------------------------------------------------------
 Variance      |
    e.WHPpain |    472.3255    35.05929                    408.3752     546.2904
   e.whpsleep |    659.1094    48.92369                    569.8695     762.3241
--------------+----------------------------------------------------------------
LR test of model vs. saturated: chi2(1)   =      9.11, Prob > chi2 = 0.0025

Endogenous variables


Observed:  WHPpain whpsleep


Exogenous variables


Observed:  crhtw3


Fitting target model:

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

Structural equation model                       Number of obs      =       363
 Estimation method = ml
 Log likelihood    = -3795.4294

--------------+----------------------------------------------------------------
              |                 OIM
              |      Coef.    Std. Err.       z    P>|z|      [95% Con f. Interval]
--------------+----------------------------------------------------------------
 Structural    |
   WHPpain <-  |
       crhtw3 |    6.442221    1.270903      5.07    0.000      3.951297     8.933145
        _cons |    16.90299    1.151834     14.67    0.000      14.64544     19.16054
  ------------+----------------------------------------------------------------
   whpsleep <- |
      WHPpain |    .7680668    .0604372     12.71    0.000      .6496121     .8865215
        _cons |    12.41057    1.732256      7.16    0.000      9.015409     15.80573
--------------+----------------------------------------------------------------
 Variance      |
    e.WHPpain |    464.2378    34.45896                    401.3825     536.9361
   e.whpsleep |    659.1094    48.92369                    569.8695     762.3241
--------------+----------------------------------------------------------------
LR test of model vs. saturated: chi2(1)   =      9.98, Prob > chi2 = 0.0016
(1 observations with missing values excluded;
 specify option 'method(mlmv)' to use all observations)

Endogenous variables


Observed:  bffpain24 whpsleep


Exogenous variables


Observed:  crhtw1


Fitting target model:

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

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

---------------+----------------------------------------------------------------
               |                 OIM
               |      Coef.    Std. Err.       z    P>|z|      [95% Con f. Interval]
---------------+----------------------------------------------------------------
 Structural     |
   bffpain24 <- |
        crhtw1 |   -.0163264    .0295801     -0.55    0.581     -.0743023     .0416495
         _cons |    .0739961    .0281819      2.63    0.009      .0187606     .1292316
  -------------+----------------------------------------------------------------
   whpsleep <-  |
     bffpain24 |   -5.072163    3.045385     -1.67    0.096     -11.04101     .8966822
         _cons |    26.64673    1.632318     16.32    0.000      23.44745     29.84601
---------------+----------------------------------------------------------------
 Variance       |
   e.bffpain24 |     .281897    .0209532                    .2436807     .3261068
    e.whpsleep |    947.2157    70.40599                    818.8033     1095.767
---------------+----------------------------------------------------------------
LR test of model vs. saturated: chi2(1)   =      7.34, Prob > chi2 = 0.0067

Endogenous variables


Observed:  bffpain24 whpsleep


Exogenous variables


Observed:  crhtw2


Fitting target model:

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

Structural equation model                       Number of obs      =       363
 Estimation method = ml
 Log likelihood    = -2507.1158

---------------+----------------------------------------------------------------
               |                 OIM
               |      Coef.    Std. Err.       z    P>|z|      [95% Con f. Interval]
---------------+----------------------------------------------------------------
 Structural     |
   bffpain24 <- |
        crhtw2 |   -.0422022    .0319287     -1.32    0.186     -.1047813     .0203769
         _cons |    .0788774    .0283111      2.79    0.005      .0233887     .1343662
  -------------+----------------------------------------------------------------
   whpsleep <-  |
     bffpain24 |   -5.062292    3.041987     -1.66    0.096     -11.02448     .8998935
         _cons |    26.60724    1.628249     16.34    0.000      23.41594     29.79855
---------------+----------------------------------------------------------------
 Variance       |
   e.bffpain24 |    .2800235    .0207853                    .2421098     .3238743
    e.whpsleep |    945.1506    70.15566                    817.1822     1093.158
---------------+----------------------------------------------------------------
LR test of model vs. saturated: chi2(1)   =     22.86, Prob > chi2 = 0.0000

Endogenous variables


Observed:  bffpain24 whpsleep


Exogenous variables


Observed:  crhtw3


Fitting target model:

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

Structural equation model                       Number of obs      =       363
 Estimation method = ml
 Log likelihood    = -2514.9789

---------------+----------------------------------------------------------------
               |                 OIM
               |      Coef.    Std. Err.       z    P>|z|      [95% Con f. Interval]
---------------+----------------------------------------------------------------
 Structural     |
   bffpain24 <- |
        crhtw3 |   -.0490766    .0311821     -1.57    0.116     -.1101923     .0120392
         _cons |    .0800705    .0282607      2.83    0.005      .0246806     .1354604
  -------------+----------------------------------------------------------------
   whpsleep <-  |
     bffpain24 |   -5.062292    3.041987     -1.66    0.096     -11.02448     .8998935
         _cons |    26.60724    1.628249     16.34    0.000      23.41594     29.79855
---------------+----------------------------------------------------------------
 Variance       |
   e.bffpain24 |    .2794642    .0207438                    .2416262     .3232274
    e.whpsleep |    945.1506    70.15566                    817.1822     1093.158
---------------+----------------------------------------------------------------
LR test of model vs. saturated: chi2(1)   =     27.45, Prob > chi2 = 0.0000
(1 observations with missing values excluded;
 specify option 'method(mlmv)' to use all observations)

Endogenous variables


Observed:  bffpain26 whpsleep


Exogenous variables


Observed:  crhtw1


Fitting target model:

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

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

---------------+----------------------------------------------------------------
               |                 OIM
               |      Coef.    Std. Err.       z    P>|z|      [95% Con f. Interval]
---------------+----------------------------------------------------------------
 Structural     |
   bffpain26 <- |
        crhtw1 |   -.1352462    .4308242     -0.31    0.754     -.9796461     .7091538
         _cons |    1.454464    .4104601      3.54    0.000      .6499767     2.258951
  -------------+----------------------------------------------------------------
   whpsleep <-  |
     bffpain26 |   -.2864247    .2094129     -1.37    0.171     -.6968663      .124017
         _cons |    26.69387    1.647301     16.20    0.000      23.46522     29.92252
---------------+----------------------------------------------------------------
 Variance       |
   e.bffpain26 |    59.79874    4.444806                    51.69193     69.17695
    e.whpsleep |    949.5669    70.58075                    820.8358     1098.487
---------------+----------------------------------------------------------------
LR test of model vs. saturated: chi2(1)   =      7.46, Prob > chi2 = 0.0063

Endogenous variables


Observed:  bffpain26 whpsleep


Exogenous variables


Observed:  crhtw2


Fitting target model:

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

Structural equation model                       Number of obs      =       363
 Estimation method = ml
 Log likelihood    = -3480.6721

---------------+----------------------------------------------------------------
               |                 OIM
               |      Coef.    Std. Err.       z    P>|z|      [95% Con f. Interval]
---------------+----------------------------------------------------------------
 Structural     |
   bffpain26 <- |
        crhtw2 |   -.0283059     .466024     -0.06    0.952     -.9416962     .8850843
         _cons |    1.437371    .4132225      3.48    0.001      .6274697     2.247272
  -------------+----------------------------------------------------------------
   whpsleep <-  |
     bffpain26 |   -.2854904    .2091748     -1.36    0.172     -.6954654     .1244847
         _cons |    26.65362    1.643161     16.22    0.000      23.43309     29.87416
---------------+----------------------------------------------------------------
 Variance       |
   e.bffpain26 |    59.65531    4.428032                    51.57829     68.99716
    e.whpsleep |     947.499    70.32997                    819.2127     1095.875
---------------+----------------------------------------------------------------
LR test of model vs. saturated: chi2(1)   =     23.81, Prob > chi2 = 0.0000

Endogenous variables


Observed:  bffpain26 whpsleep


Exogenous variables


Observed:  crhtw3


Fitting target model:

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

Structural equation model                       Number of obs      =       363
 Estimation method = ml
 Log likelihood    = -3488.8833

---------------+----------------------------------------------------------------
               |                 OIM
               |      Coef.    Std. Err.       z    P>|z|      [95% Con f. Interval]
---------------+----------------------------------------------------------------
 Structural     |
   bffpain26 <- |
        crhtw3 |   -.0830469    .4555636     -0.18    0.855     -.9759351     .8098413
         _cons |    1.446798    .4128824      3.50    0.000       .637563     2.256032
  -------------+----------------------------------------------------------------
   whpsleep <-  |
     bffpain26 |   -.2854904    .2091748     -1.36    0.172     -.6954654     .1244847
         _cons |    26.65362    1.643161     16.22    0.000      23.43309     29.87416
---------------+----------------------------------------------------------------
 Variance       |
   e.bffpain26 |    59.65045    4.427672                     51.5741     68.99154
    e.whpsleep |     947.499    70.32997                    819.2127     1095.875
---------------+----------------------------------------------------------------
LR test of model vs. saturated: chi2(1)   =     28.57, Prob > chi2 = 0.0000
(1 observations with missing values excluded;
 specify option 'method(mlmv)' to use all observations)

Endogenous variables


Observed:  BSIanx whpsleep


Exogenous variables


Observed:  crhtw1


Fitting target model:

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

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

--------------+----------------------------------------------------------------
              |                 OIM
              |      Coef.    Std. Err.       z    P>|z|      [95% Con f. Interval]
--------------+----------------------------------------------------------------
 Structural    |
   BSIanx <-   |
       crhtw1 |   -.4122699    .2021768     -2.04    0.041     -.8085292    -.0160107
        _cons |    9.043819    .1926204     46.95    0.000       8.66629     9.421348
  ------------+----------------------------------------------------------------
   whpsleep <- |
       BSIanx |    4.067341    .3901842     10.42    0.000      3.302594     4.832088
        _cons |    -10.2787    3.785421     -2.72    0.007     -17.69799    -2.859408
--------------+----------------------------------------------------------------
 Variance      |
     e.BSIanx |    13.16905    .9788479                    11.38375     15.23435
   e.whpsleep |    734.1125    54.56615                    634.5901     849.2429
--------------+----------------------------------------------------------------
LR test of model vs. saturated: chi2(1)   =     18.53, Prob > chi2 = 0.0000

Endogenous variables


Observed:  BSIanx whpsleep


Exogenous variables


Observed:  crhtw2


Fitting target model:

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

Structural equation model                       Number of obs      =       363
 Estimation method = ml
 Log likelihood    = -3154.7839

--------------+----------------------------------------------------------------
              |                 OIM
              |      Coef.    Std. Err.       z    P>|z|      [95% Con f. Interval]
--------------+----------------------------------------------------------------
 Structural    |
   BSIanx <-   |
       crhtw2 |    .8193537    .2160312      3.79    0.000      .3959404     1.242767
        _cons |    8.837162    .1915544     46.13    0.000      8.461723     9.212602
  ------------+----------------------------------------------------------------
   whpsleep <- |
       BSIanx |    4.065276    .3890115     10.45    0.000      3.302827     4.827724
        _cons |   -10.25323    3.770224     -2.72    0.007     -17.64274     -2.86373
--------------+----------------------------------------------------------------
 Variance      |
     e.BSIanx |    12.81933    .9515401                    11.08366      14.8268
   e.whpsleep |    732.1075    54.34211                     632.984     846.7534
--------------+----------------------------------------------------------------
LR test of model vs. saturated: chi2(1)   =     12.41, Prob > chi2 = 0.0004

Endogenous variables


Observed:  BSIanx whpsleep


Exogenous variables


Observed:  crhtw3


Fitting target model:

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

Structural equation model                       Number of obs      =       363
 Estimation method = ml
 Log likelihood    = -3158.9796

--------------+----------------------------------------------------------------
              |                 OIM
              |      Coef.    Std. Err.       z    P>|z|      [95% Con f. Interval]
--------------+----------------------------------------------------------------
 Structural    |
   BSIanx <-   |
       crhtw3 |    .9985733    .2088589      4.78    0.000      .5892175     1.407929
        _cons |    8.806126    .1892912     46.52    0.000      8.435122      9.17713
  ------------+----------------------------------------------------------------
   whpsleep <- |
       BSIanx |    4.065276    .3890115     10.45    0.000      3.302827     4.827724
        _cons |   -10.25323    3.770224     -2.72    0.007     -17.64274     -2.86373
--------------+----------------------------------------------------------------
 Variance      |
     e.BSIanx |    12.53781    .9306433                    10.84025     14.50119
   e.whpsleep |    732.1075    54.34211                     632.984     846.7534
--------------+----------------------------------------------------------------
LR test of model vs. saturated: chi2(1)   =     12.76, Prob > chi2 = 0.0004
(1 observations with missing values excluded;
 specify option 'method(mlmv)' to use all observations)

Endogenous variables


Observed:  HP2probsoc whpsleep


Exogenous variables


Observed:  crhtw1


Fitting target model:

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

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

----------------+----------------------------------------------------------------
                |                 OIM
                |      Coef.    Std. Err.       z    P>|z|      [95% Con f. Interval]
----------------+----------------------------------------------------------------
 Structural      |
   HP2probsoc <- |
         crhtw1 |    .1108202    .0216995      5.11    0.000      .0682899     .1533505
          _cons |     .189671    .0206738      9.17    0.000       .149151     .2301909
  --------------+----------------------------------------------------------------
   whpsleep <-   |
     HP2probsoc |    32.64489    3.642603      8.96    0.000      25.50552     39.78426
          _cons |    19.60917    1.646923     11.91    0.000      16.38126     22.83708
----------------+----------------------------------------------------------------
 Variance        |
   e.HP2probsoc |    .1517024    .0112759                    .1311363     .1754938
     e.whpsleep |    781.1586    58.06306                    675.2583     903.6672
----------------+----------------------------------------------------------------
LR test of model vs. saturated: chi2(1)   =      0.52, Prob > chi2 = 0.4705

Endogenous variables


Observed:  HP2probsoc whpsleep


Exogenous variables


Observed:  crhtw2


Fitting target model:

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

Structural equation model                       Number of obs      =       363
 Estimation method = ml
 Log likelihood    = -2355.9971

----------------+----------------------------------------------------------------
                |                 OIM
                |      Coef.    Std. Err.       z    P>|z|      [95% Con f. Interval]
----------------+----------------------------------------------------------------
 Structural      |
   HP2probsoc <- |
         crhtw2 |    .1387378    .0231912      5.98    0.000      .0932838     .1841918
          _cons |    .1800158    .0205636      8.75    0.000      .1397119     .2203198
  --------------+----------------------------------------------------------------
   whpsleep <-   |
     HP2probsoc |    32.66925    3.636612      8.98    0.000      25.54162     39.79687
          _cons |    19.58481    1.641949     11.93    0.000      16.36665     22.80297
----------------+----------------------------------------------------------------
 Variance        |
   e.HP2probsoc |     .147734    .0109658                    .1277316     .1708687
     e.whpsleep |    779.1427    57.83339                    673.6509     901.1542
----------------+----------------------------------------------------------------
LR test of model vs. saturated: chi2(1)   =      7.55, Prob > chi2 = 0.0060

Endogenous variables


Observed:  HP2probsoc whpsleep


Exogenous variables


Observed:  crhtw3


Fitting target model:

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

Structural equation model                       Number of obs      =       363
 Estimation method = ml
 Log likelihood    = -2364.6787

----------------+----------------------------------------------------------------
                |                 OIM
                |      Coef.    Std. Err.       z    P>|z|      [95% Con f. Interval]
----------------+----------------------------------------------------------------
 Structural      |
   HP2probsoc <- |
         crhtw3 |    .1338891    .0227001      5.90    0.000      .0893977     .1783804
          _cons |     .180817    .0205733      8.79    0.000       .140494       .22114
  --------------+----------------------------------------------------------------
   whpsleep <-   |
     HP2probsoc |    32.66925    3.636612      8.98    0.000      25.54162     39.79687
          _cons |    19.58481    1.641949     11.93    0.000      16.36665     22.80297
----------------+----------------------------------------------------------------
 Variance        |
   e.HP2probsoc |    .1481053    .0109934                    .1280527     .1712982
     e.whpsleep |    779.1427    57.83339                    673.6509     901.1542
----------------+----------------------------------------------------------------
LR test of model vs. saturated: chi2(1)   =     10.95, Prob > chi2 = 0.0009
(1 observations with missing values excluded;
 specify option 'method(mlmv)' to use all observations)

Endogenous variables


Observed:  BSIposymp whpsleep


Exogenous variables


Observed:  crhtw1


Fitting target model:

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

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

---------------+----------------------------------------------------------------
               |                 OIM
               |      Coef.    Std. Err.       z    P>|z|      [95% Con f. Interval]
---------------+----------------------------------------------------------------
 Structural     |
   BSIposymp <- |
        crhtw1 |    1.618649    1.561776      1.04    0.300     -1.442375     4.679673
         _cons |    86.38402    1.487954     58.06    0.000      83.46769     89.30036
  -------------+----------------------------------------------------------------
   whpsleep <-  |
     BSIposymp |    .6363545    .0471879     13.49    0.000      .5438679     .7288411
         _cons |   -28.82552    4.295816     -6.71    0.000     -37.24516    -20.40587
---------------+----------------------------------------------------------------
 Variance       |
   e.BSIposymp |     785.831    58.41036                    679.2972     909.0724
    e.whpsleep |    635.3099    47.22221                     549.182     734.9452
---------------+----------------------------------------------------------------
LR test of model vs. saturated: chi2(1)   =      6.93, Prob > chi2 = 0.0085

Endogenous variables


Observed:  BSIposymp whpsleep


Exogenous variables


Observed:  crhtw2


Fitting target model:

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

Structural equation model                       Number of obs      =       363
 Estimation method = ml
 Log likelihood    =  -3858.806

---------------+----------------------------------------------------------------
               |                 OIM
               |      Coef.    Std. Err.       z    P>|z|      [95% Con f. Interval]
---------------+----------------------------------------------------------------
 Structural     |
   BSIposymp <- |
        crhtw2 |    9.803972    1.614914      6.07    0.000      6.638798     12.96915
         _cons |    84.82216    1.431941     59.24    0.000       82.0156     87.62871
  -------------+----------------------------------------------------------------
   whpsleep <-  |
     BSIposymp |    .6354575    .0470361     13.51    0.000      .5432683     .7276466
         _cons |   -28.72679    4.278099     -6.71    0.000     -37.11171    -20.34187
---------------+----------------------------------------------------------------
 Variance       |
   e.BSIposymp |    716.3602    53.17324                    619.3688     828.5401
    e.whpsleep |    633.7207    47.03916                    547.9183     732.9595
---------------+----------------------------------------------------------------
LR test of model vs. saturated: chi2(1)   =      3.48, Prob > chi2 = 0.0620

Endogenous variables


Observed:  BSIposymp whpsleep


Exogenous variables


Observed:  crhtw3


Fitting target model:

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

Structural equation model                       Number of obs      =       363
 Estimation method = ml
 Log likelihood    = -3861.1852

---------------+----------------------------------------------------------------
               |                 OIM
               |      Coef.    Std. Err.       z    P>|z|      [95% Con f. Interval]
---------------+----------------------------------------------------------------
 Structural     |
   BSIposymp <- |
        crhtw3 |    10.97906    1.553505      7.07    0.000      7.934242     14.02387
         _cons |     84.6176    1.407959     60.10    0.000      81.85805     87.37715
  -------------+----------------------------------------------------------------
   whpsleep <-  |
     BSIposymp |    .6354575    .0470361     13.51    0.000      .5432683     .7276466
         _cons |   -28.72679    4.278099     -6.71    0.000     -37.11171    -20.34187
---------------+----------------------------------------------------------------
 Variance       |
   e.BSIposymp |     693.651    51.48761                    599.7344     802.2748
    e.whpsleep |    633.7207    47.03916                    547.9183     732.9595
---------------+----------------------------------------------------------------
LR test of model vs. saturated: chi2(1)   =      3.41, Prob > chi2 = 0.0649
(1 observations with missing values excluded;
 specify option 'method(mlmv)' to use all observations)

Endogenous variables


Observed:  BSItotal whpsleep


Exogenous variables


Observed:  crhtw1


Fitting target model:

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

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

--------------+----------------------------------------------------------------
              |                 OIM
              |      Coef.    Std. Err.       z    P>|z|      [95% Con f. Interval]
--------------+----------------------------------------------------------------
 Structural    |
   BSItotal <- |
       crhtw1 |    1.656929    1.641141      1.01    0.313     -1.559648     4.873506
        _cons |    91.61374    1.563568     58.59    0.000       88.5492     94.67827
  ------------+----------------------------------------------------------------
   whpsleep <- |
     BSItotal |    .5924032    .0453939     13.05    0.000      .5034328     .6813735
        _cons |   -28.12047    4.378499     -6.42    0.000     -36.70217    -19.53877
--------------+----------------------------------------------------------------
 Variance      |
   e.BSItotal |     867.728    64.49772                    750.0916     1003.813
   e.whpsleep |    649.0948    48.24684                    561.0982      750.892
--------------+----------------------------------------------------------------
LR test of model vs. saturated: chi2(1)   =      6.96, Prob > chi2 = 0.0083

Endogenous variables


Observed:  BSItotal whpsleep


Exogenous variables


Observed:  crhtw2


Fitting target model:

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

Structural equation model                       Number of obs      =       363
 Estimation method = ml
 Log likelihood    = -3880.7778

--------------+----------------------------------------------------------------
              |                 OIM
              |      Coef.    Std. Err.       z    P>|z|      [95% Con f. Interval]
--------------+----------------------------------------------------------------
 Structural    |
   BSItotal <- |
       crhtw2 |     10.2622    1.697409      6.05    0.000      6.935335     13.58906
        _cons |    89.97482    1.505089     59.78    0.000       87.0249     92.92474
  ------------+----------------------------------------------------------------
   whpsleep <- |
     BSItotal |    .5916426    .0452491     13.08    0.000       .502956     .6803291
        _cons |   -28.03162    4.360618     -6.43    0.000     -36.57828    -19.48497
--------------+----------------------------------------------------------------
 Variance      |
   e.BSItotal |    791.4173     58.7445                    684.2636      915.351
   e.whpsleep |    647.4379    48.05735                    559.7782     748.8248
--------------+----------------------------------------------------------------
LR test of model vs. saturated: chi2(1)   =      3.82, Prob > chi2 = 0.0506

Endogenous variables


Observed:  BSItotal whpsleep


Exogenous variables


Observed:  crhtw3


Fitting target model:

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

Structural equation model                       Number of obs      =       363
 Estimation method = ml
 Log likelihood    = -3883.0889

--------------+----------------------------------------------------------------
              |                 OIM
              |      Coef.    Std. Err.       z    P>|z|      [95% Con f. Interval]
--------------+----------------------------------------------------------------
 Structural    |
   BSItotal <- |
       crhtw3 |    11.51934    1.632556      7.06    0.000      8.319587     14.71909
        _cons |    89.75604    1.479604     60.66    0.000      86.85606     92.65601
  ------------+----------------------------------------------------------------
   whpsleep <- |
     BSItotal |    .5916426    .0452491     13.08    0.000       .502956     .6803291
        _cons |   -28.03162    4.360618     -6.43    0.000     -36.57828    -19.48497
--------------+----------------------------------------------------------------
 Variance      |
   e.BSItotal |    766.0413    56.86092                    662.3234     886.0012
   e.whpsleep |    647.4379    48.05735                    559.7782     748.8248
--------------+----------------------------------------------------------------
LR test of model vs. saturated: chi2(1)   =      3.77, Prob > chi2 = 0.0522
(1 observations with missing values excluded;
 specify option 'method(mlmv)' to use all observations)

Endogenous variables


Observed:  BSIdep whpsleep


Exogenous variables


Observed:  crhtw1


Fitting target model:

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

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

--------------+----------------------------------------------------------------
              |                 OIM
              |      Coef.    Std. Err.       z    P>|z|      [95% Con f. Interval]
--------------+----------------------------------------------------------------
 Structural    |
   BSIdep <-   |
       crhtw1 |    .3880135    .2084177      1.86    0.063     -.0204776     .7965047
        _cons |    9.614106    .1985662     48.42    0.000      9.224923     10.00329
  ------------+----------------------------------------------------------------
   whpsleep <- |
       BSIdep |    3.444293    .3922347      8.78    0.000      2.675527     4.213058
        _cons |   -7.009225    4.067821     -1.72    0.085     -14.98201     .9635566
--------------+----------------------------------------------------------------
 Variance      |
     e.BSIdep |    13.99461    1.040211                    12.09739     16.18938
   e.whpsleep |    786.8643    58.48716                    680.1905     910.2678
--------------+----------------------------------------------------------------
LR test of model vs. saturated: chi2(1)   =      4.71, Prob > chi2 = 0.0299

Endogenous variables


Observed:  BSIdep whpsleep


Exogenous variables


Observed:  crhtw2


Fitting target model:

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

Structural equation model                       Number of obs      =       363
 Estimation method = ml
 Log likelihood    =  -3169.858

--------------+----------------------------------------------------------------
              |                 OIM
              |      Coef.    Std. Err.       z    P>|z|      [95% Con f. Interval]
--------------+----------------------------------------------------------------
 Structural    |
   BSIdep <-   |
       crhtw2 |    1.222397    .2175138      5.62    0.000      .7960774     1.648716
        _cons |    9.445589     .192869     48.97    0.000      9.067572     9.823605
  ------------+----------------------------------------------------------------
   whpsleep <- |
       BSIdep |    3.445067    .3911827      8.81    0.000      2.678363     4.211771
        _cons |   -7.019698    4.053191     -1.73    0.083     -14.96381     .9244117
--------------+----------------------------------------------------------------
 Variance      |
     e.BSIdep |    12.99589    .9646457                    11.23632     15.03101
   e.whpsleep |    784.6999    58.24589                    678.4557     907.5816
--------------+----------------------------------------------------------------
LR test of model vs. saturated: chi2(1)   =      8.55, Prob > chi2 = 0.0035

Endogenous variables


Observed:  BSIdep whpsleep


Exogenous variables


Observed:  crhtw3


Fitting target model:

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

Structural equation model                       Number of obs      =       363
 Estimation method = ml
 Log likelihood    = -3173.3492

--------------+----------------------------------------------------------------
              |                 OIM
              |      Coef.    Std. Err.       z    P>|z|      [95% Con f. Interval]
--------------+----------------------------------------------------------------
 Structural    |
   BSIdep <-   |
       crhtw3 |    1.360398    .2098845      6.48    0.000      .9490322     1.771764
        _cons |    9.421549    .1902207     49.53    0.000      9.048723     9.794375
  ------------+----------------------------------------------------------------
   whpsleep <- |
       BSIdep |    3.445067    .3911827      8.81    0.000      2.678363     4.211771
        _cons |   -7.019698    4.053191     -1.73    0.083     -14.96381     .9244117
--------------+----------------------------------------------------------------
 Variance      |
     e.BSIdep |    12.66125    .9398062                    10.94698     14.64397
   e.whpsleep |    784.6999    58.24589                    678.4557     907.5816
--------------+----------------------------------------------------------------
LR test of model vs. saturated: chi2(1)   =      9.76, Prob > chi2 = 0.0018
(1 observations with missing values excluded;
 specify option 'method(mlmv)' to use all observations)

Endogenous variables


Observed:  BFdepx1 whpsleep


Exogenous variables


Observed:  crhtw1


Fitting target model:

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

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

--------------+----------------------------------------------------------------
              |                 OIM
              |      Coef.    Std. Err.       z    P>|z|      [95% Con f. Interval]
--------------+----------------------------------------------------------------
 Structural    |
   BFdepx1 <-  |
       crhtw1 |    .7787754    .9676886      0.80    0.421      -1.11786      2.67541
        _cons |    15.42812    .9219482     16.73    0.000      13.62114     17.23511
  ------------+----------------------------------------------------------------
   whpsleep <- |
      BFdepx1 |    1.100099    .0733544     15.00    0.000      .9563268     1.243871
        _cons |    9.195953    1.710066      5.38    0.000      5.844285     12.54762
--------------+----------------------------------------------------------------
 Variance      |
    e.BFdepx1 |    301.6915    22.42456                    260.7917     349.0056
   e.whpsleep |    588.7085    43.75836                    508.8983     681.0353
--------------+----------------------------------------------------------------
LR test of model vs. saturated: chi2(1)   =      8.20, Prob > chi2 = 0.0042

Endogenous variables


Observed:  BFdepx1 whpsleep


Exogenous variables


Observed:  crhtw2


Fitting target model:

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

Structural equation model                       Number of obs      =       363
 Estimation method = ml
 Log likelihood    = -3683.5325

--------------+----------------------------------------------------------------
              |                 OIM
              |      Coef.    Std. Err.       z    P>|z|      [95% Con f. Interval]
--------------+----------------------------------------------------------------
 Structural    |
   BFdepx1 <-  |
       crhtw2 |    3.186106    1.035232      3.08    0.002      1.157088     5.215124
        _cons |    14.94148    .9179385     16.28    0.000      13.14235      16.7406
  ------------+----------------------------------------------------------------
   whpsleep <- |
      BFdepx1 |    1.099622    .0731748     15.03    0.000      .9562021     1.243042
        _cons |     9.21263    1.703528      5.41    0.000      5.873776     12.55148
--------------+----------------------------------------------------------------
 Variance      |
    e.BFdepx1 |    294.3799    21.85093                    254.5225      340.479
   e.whpsleep |     587.118    43.57998                    507.6253     679.0589
--------------+----------------------------------------------------------------
LR test of model vs. saturated: chi2(1)   =     14.39, Prob > chi2 = 0.0001

Endogenous variables


Observed:  BFdepx1 whpsleep


Exogenous variables


Observed:  crhtw3


Fitting target model:

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

Structural equation model                       Number of obs      =       363
 Estimation method = ml
 Log likelihood    =   -3688.42

--------------+----------------------------------------------------------------
              |                 OIM
              |      Coef.    Std. Err.       z    P>|z|      [95% Con f. Interval]
--------------+----------------------------------------------------------------
 Structural    |
   BFdepx1 <-  |
       crhtw3 |    4.059266    1.002772      4.05    0.000       2.09387     6.024662
        _cons |    14.79046    .9088233     16.27    0.000       13.0092     16.57172
  ------------+----------------------------------------------------------------
   whpsleep <- |
      BFdepx1 |    1.099622    .0731748     15.03    0.000      .9562021     1.243042
        _cons |     9.21263    1.703528      5.41    0.000      5.873776     12.55148
--------------+----------------------------------------------------------------
 Variance      |
    e.BFdepx1 |    289.0146    21.45268                    249.8836     334.2735
   e.whpsleep |     587.118    43.57998                    507.6253     679.0589
--------------+----------------------------------------------------------------
LR test of model vs. saturated: chi2(1)   =     13.43, Prob > chi2 = 0.0002
(1 observations with missing values excluded;
 specify option 'method(mlmv)' to use all observations)

Endogenous variables


Observed:  icdx3nr5 whpsleep


Exogenous variables


Observed:  crhtw1


Fitting target model:

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

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

--------------+----------------------------------------------------------------
              |                 OIM
              |      Coef.    Std. Err.       z    P>|z|      [95% Con f. Interval]
--------------+----------------------------------------------------------------
 Structural    |
   icdx3nr5 <- |
       crhtw1 |   -.0070949    .0107359     -0.66    0.509     -.0281368     .0139471
        _cons |    .0396183    .0102284      3.87    0.000       .019571     .0596656
  ------------+----------------------------------------------------------------
   whpsleep <- |
     icdx3nr5 |    15.56897     8.38152      1.86    0.063     -.8585064     31.99645
        _cons |    25.68032    1.648286     15.58    0.000      22.44974      28.9109
--------------+----------------------------------------------------------------
 Variance      |
   e.icdx3nr5 |    .0371336    .0027601                    .0320994     .0429572
   e.whpsleep |    945.4623    70.27566                    817.2876     1093.739
--------------+----------------------------------------------------------------
LR test of model vs. saturated: chi2(1)   =      7.99, Prob > chi2 = 0.0047

Endogenous variables


Observed:  icdx3nr5 whpsleep


Exogenous variables


Observed:  crhtw2


Fitting target model:

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

Structural equation model                       Number of obs      =       363
 Estimation method = ml
 Log likelihood    = -2139.3735

--------------+----------------------------------------------------------------
              |                 OIM
              |      Coef.    Std. Err.       z    P>|z|      [95% Con f. Interval]
--------------+----------------------------------------------------------------
 Structural    |
   icdx3nr5 <- |
       crhtw2 |    .0108512    .0116046      0.94    0.350     -.0118935     .0335958
        _cons |    .0367028    .0102898      3.57    0.000      .0165352     .0568705
  ------------+----------------------------------------------------------------
   whpsleep <- |
     icdx3nr5 |    15.60654    8.371598      1.86    0.062     -.8014958     32.01457
        _cons |    25.64275    1.644065     15.60    0.000      22.42044     28.86506
--------------+----------------------------------------------------------------
 Variance      |
   e.icdx3nr5 |    .0369909    .0027457                    .0319826     .0427836
   e.whpsleep |    943.3299    70.02051                     815.608     1091.053
--------------+----------------------------------------------------------------
LR test of model vs. saturated: chi2(1)   =     23.09, Prob > chi2 = 0.0000

Endogenous variables


Observed:  icdx3nr5 whpsleep


Exogenous variables


Observed:  crhtw3


Fitting target model:

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

Structural equation model                       Number of obs      =       363
 Estimation method = ml
 Log likelihood    = -2147.6116

--------------+----------------------------------------------------------------
              |                 OIM
              |      Coef.    Std. Err.       z    P>|z|      [95% Con f. Interval]
--------------+----------------------------------------------------------------
 Structural    |
   icdx3nr5 <- |
       crhtw3 |    .0104594     .011345      0.92    0.357     -.0117764     .0326951
        _cons |    .0367676    .0102821      3.58    0.000      .0166151     .0569202
  ------------+----------------------------------------------------------------
   whpsleep <- |
     icdx3nr5 |    15.60654    8.371598      1.86    0.062     -.8014958     32.01457
        _cons |    25.64275    1.644065     15.60    0.000      22.42044     28.86506
--------------+----------------------------------------------------------------
 Variance      |
   e.icdx3nr5 |    .0369934    .0027459                    .0319847     .0427865
   e.whpsleep |    943.3299    70.02051                     815.608     1091.053
--------------+----------------------------------------------------------------
LR test of model vs. saturated: chi2(1)   =     27.90, Prob > chi2 = 0.0000
(1 observations with missing values excluded;
 specify option 'method(mlmv)' to use all observations)

Endogenous variables


Observed:  icdx4nr6 whpsleep


Exogenous variables


Observed:  crhtw1


Fitting target model:

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

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

--------------+----------------------------------------------------------------
              |                 OIM
              |      Coef.    Std. Err.       z    P>|z|      [95% Con f. Interval]
--------------+----------------------------------------------------------------
 Structural    |
   icdx4nr6 <- |
       crhtw1 |   -.0082753    .0076597     -1.08    0.280     -.0232879     .0067374
        _cons |    .0204384    .0072976      2.80    0.005      .0061353     .0347414
  ------------+----------------------------------------------------------------
   whpsleep <- |
     icdx4nr6 |    25.54438    11.71493      2.18    0.029      2.583539     48.50522
        _cons |    25.78848     1.62905     15.83    0.000       22.5956     28.98136
--------------+----------------------------------------------------------------
 Variance      |
   e.icdx4nr6 |    .0189022     .001405                    .0163396     .0218666
   e.whpsleep |    942.1004    70.02577                    814.3815     1089.849
--------------+----------------------------------------------------------------
LR test of model vs. saturated: chi2(1)   =      8.38, Prob > chi2 = 0.0038

Endogenous variables


Observed:  icdx4nr6 whpsleep


Exogenous variables


Observed:  crhtw2


Fitting target model:

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

Structural equation model                       Number of obs      =       363
 Estimation method = ml
 Log likelihood    = -2016.2402

--------------+----------------------------------------------------------------
              |                 OIM
              |      Coef.    Std. Err.       z    P>|z|      [95% Con f. Interval]
--------------+----------------------------------------------------------------
 Structural    |
   icdx4nr6 <- |
       crhtw2 |    .0099703    .0082811      1.20    0.229     -.0062602     .0262009
        _cons |    .0175704    .0073428      2.39    0.017      .0031788      .031962
  ------------+----------------------------------------------------------------
   whpsleep <- |
     icdx4nr6 |    25.58151    11.70145      2.19    0.029      2.647084     48.51593
        _cons |    25.75135    1.624933     15.85    0.000      22.56654     28.93616
--------------+----------------------------------------------------------------
 Variance      |
   e.icdx4nr6 |    .0188367    .0013982                    .0162863     .0217864
   e.whpsleep |    939.9851    69.77224                    812.7161     1087.184
--------------+----------------------------------------------------------------
LR test of model vs. saturated: chi2(1)   =     22.75, Prob > chi2 = 0.0000

Endogenous variables


Observed:  icdx4nr6 whpsleep


Exogenous variables


Observed:  crhtw3


Fitting target model:

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

Structural equation model                       Number of obs      =       363
 Estimation method = ml
 Log likelihood    = -2024.1593

--------------+----------------------------------------------------------------
              |                 OIM
              |      Coef.    Std. Err.       z    P>|z|      [95% Con f. Interval]
--------------+----------------------------------------------------------------
 Structural    |
   icdx4nr6 <- |
       crhtw3 |    .0116273    .0080887      1.44    0.151     -.0042262     .0274808
        _cons |    .0172829    .0073308      2.36    0.018      .0029147     .0316511
  ------------+----------------------------------------------------------------
   whpsleep <- |
     icdx4nr6 |    25.58151    11.70145      2.19    0.029      2.647084     48.51593
        _cons |    25.75135    1.624933     15.85    0.000      22.56654     28.93616
--------------+----------------------------------------------------------------
 Variance      |
   e.icdx4nr6 |    .0188048    .0013958                    .0162588     .0217496
   e.whpsleep |    939.9851    69.77224                    812.7161     1087.184
--------------+----------------------------------------------------------------
LR test of model vs. saturated: chi2(1)   =     27.27, Prob > chi2 = 0.0000
(1 observations with missing values excluded;
 specify option 'method(mlmv)' to use all observations)

Endogenous variables


Observed:  icdx5nr4 whpsleep


Exogenous variables


Observed:  crhtw1


Fitting target model:

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

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

--------------+----------------------------------------------------------------
              |                 OIM
              |      Coef.    Std. Err.       z    P>|z|      [95% Con f. Interval]
--------------+----------------------------------------------------------------
 Structural    |
   icdx5nr4 <- |
       crhtw1 |   -.0038018    .0029173     -1.30    0.193     -.0095196      .001916
        _cons |    .0032684    .0027794      1.18    0.240     -.0021791      .008716
  ------------+----------------------------------------------------------------
   whpsleep <- |
     icdx5nr4 |    74.92454    30.68568      2.44    0.015      14.78172     135.0674
        _cons |    26.07546    1.612803     16.17    0.000      22.91442     29.23649
--------------+----------------------------------------------------------------
 Variance      |
   e.icdx5nr4 |    .0027419    .0002038                    .0023702      .003172
   e.whpsleep |    939.0095    69.79603                    811.7096     1086.274
--------------+----------------------------------------------------------------
LR test of model vs. saturated: chi2(1)   =      8.68, Prob > chi2 = 0.0032

Endogenous variables


Observed:  icdx5nr4 whpsleep


Exogenous variables


Observed:  crhtw2


Fitting target model:

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

Structural equation model                       Number of obs      =       363
 Estimation method = ml
 Log likelihood    = -1665.5852

--------------+----------------------------------------------------------------
              |                 OIM
              |      Coef.    Std. Err.       z    P>|z|      [95% Con f. Interval]
--------------+----------------------------------------------------------------
 Structural    |
   icdx5nr4 <- |
       crhtw2 |    .0035669     .003157      1.13    0.259     -.0026206     .0097544
        _cons |    .0021419    .0027993      0.77    0.444     -.0033446     .0076283
  ------------+----------------------------------------------------------------
   whpsleep <- |
     icdx5nr4 |    74.96185    30.65146      2.45    0.014      14.88609     135.0376
        _cons |    26.03815    1.608784     16.18    0.000      22.88499     29.19131
--------------+----------------------------------------------------------------
 Variance      |
   e.icdx5nr4 |    .0027376    .0002032                    .0023669     .0031663
   e.whpsleep |    936.9238    69.54501                    810.0693     1083.643
--------------+----------------------------------------------------------------
LR test of model vs. saturated: chi2(1)   =     22.75, Prob > chi2 = 0.0000

Endogenous variables


Observed:  icdx5nr4 whpsleep


Exogenous variables


Observed:  crhtw3


Fitting target model:

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

Structural equation model                       Number of obs      =       363
 Estimation method = ml
 Log likelihood    = -1673.8799

--------------+----------------------------------------------------------------
              |                 OIM
              |      Coef.    Std. Err.       z    P>|z|      [95% Con f. Interval]
--------------+----------------------------------------------------------------
 Structural    |
   icdx5nr4 <- |
       crhtw3 |     .003294    .0030868      1.07    0.286     -.0027561      .009344
        _cons |     .002188    .0027976      0.78    0.434     -.0032952     .0076712
  ------------+----------------------------------------------------------------
   whpsleep <- |
     icdx5nr4 |    74.96185    30.65146      2.45    0.014      14.88609     135.0376
        _cons |    26.03815    1.608784     16.18    0.000      22.88499     29.19131
--------------+----------------------------------------------------------------
 Variance      |
   e.icdx5nr4 |    .0027386    .0002033                    .0023678     .0031675
   e.whpsleep |    936.9238    69.54501                    810.0693     1083.643
--------------+----------------------------------------------------------------
LR test of model vs. saturated: chi2(1)   =     27.60, Prob > chi2 = 0.0000
(1 observations with missing values excluded;
 specify option 'method(mlmv)' to use all observations)

Endogenous variables


Observed:  icdx6nr5 whpsleep


Exogenous variables


Observed:  crhtw1


Fitting target model:

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

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

--------------+----------------------------------------------------------------
              |                 OIM
              |      Coef.    Std. Err.       z    P>|z|      [95% Con f. Interval]
--------------+----------------------------------------------------------------
 Structural    |
   icdx6nr5 <- |
       crhtw1 |   -.0010146    .0041293     -0.25    0.806     -.0091079     .0070786
        _cons |    .0056599    .0039341      1.44    0.150     -.0020508     .0133706
  ------------+----------------------------------------------------------------
   whpsleep <- |
     icdx6nr5 |    75.13267     21.5475      3.49    0.000      32.90035      117.365
        _cons |    25.86733    1.601613     16.15    0.000      22.72823     29.00644
--------------+----------------------------------------------------------------
 Variance      |
   e.icdx6nr5 |    .0054934    .0004083                    .0047487      .006355
   e.whpsleep |     923.459    68.64017                    798.2673     1068.285
--------------+----------------------------------------------------------------
LR test of model vs. saturated: chi2(1)   =      8.06, Prob > chi2 = 0.0045

Endogenous variables


Observed:  icdx6nr5 whpsleep


Exogenous variables


Observed:  crhtw2


Fitting target model:

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

Structural equation model                       Number of obs      =       363
 Estimation method = ml
 Log likelihood    = -1788.1372

--------------+----------------------------------------------------------------
              |                 OIM
              |      Coef.    Std. Err.       z    P>|z|      [95% Con f. Interval]
--------------+----------------------------------------------------------------
 Structural    |
   icdx6nr5 <- |
       crhtw2 |    -.003772    .0044619     -0.85    0.398     -.0125171     .0049732
        _cons |    .0061578    .0039563      1.56    0.120     -.0015965     .0139121
  ------------+----------------------------------------------------------------
   whpsleep <- |
     icdx6nr5 |     75.1695    21.52331      3.49    0.000      32.98459     117.3544
        _cons |     25.8305     1.59761     16.17    0.000      22.69924     28.96176
--------------+----------------------------------------------------------------
 Variance      |
   e.icdx6nr5 |    .0054685    .0004059                    .0047281     .0063249
   e.whpsleep |    921.4008    68.39278                     796.648      1065.69
--------------+----------------------------------------------------------------
LR test of model vs. saturated: chi2(1)   =     26.25, Prob > chi2 = 0.0000

Endogenous variables


Observed:  icdx6nr5 whpsleep


Exogenous variables


Observed:  crhtw3


Fitting target model:

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

Structural equation model                       Number of obs      =       363
 Estimation method = ml
 Log likelihood    = -1796.2697

--------------+----------------------------------------------------------------
              |                 OIM
              |      Coef.    Std. Err.       z    P>|z|      [95% Con f. Interval]
--------------+----------------------------------------------------------------
 Structural    |
   icdx6nr5 <- |
       crhtw3 |   -.0041419    .0043608     -0.95    0.342     -.0126888     .0044051
        _cons |    .0062224    .0039522      1.57    0.115     -.0015239     .0139686
  ------------+----------------------------------------------------------------
   whpsleep <- |
     icdx6nr5 |     75.1695    21.52331      3.49    0.000      32.98459     117.3544
        _cons |     25.8305     1.59761     16.17    0.000      22.69924     28.96176
--------------+----------------------------------------------------------------
 Variance      |
   e.icdx6nr5 |    .0054657    .0004057                    .0047257     .0063216
   e.whpsleep |    921.4008    68.39278                     796.648      1065.69
--------------+----------------------------------------------------------------
LR test of model vs. saturated: chi2(1)   =     31.69, Prob > chi2 = 0.0000

.
. des BFdepx1 BSIdep BSItotal BSIposymp HP2probsoc BSIanx WHPpain age

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


BFdepx1           float    %9.0g                         max(0, WHPer + 3.14933e-007)


BSIdep            double   %9.0g                        Basic symptom inventory Depression subscale


BSItotal          double   %9.0g                        Basic symptom inventory total scale score


BSIposymp         double   %9.0g                        Brief Symptom inventory positive symptom total subscale


HP2probsoc        double   %9.0g         hp2fmt         Hlth profile Pt2: Hlth causing probs with social life


BSIanx            double   %9.0g                        Basic symptom inventory Anxiety subscale


WHPpain           double   %9.0g                        Wtd Health Profile Pain Pt 1 subscale


age               double   %8.0g                    *   Respondent's age


.
. * Possible female sleep indirect effects ----------------------------------------
. *-+- Possible female indirect effects through
. * wave 1 HP2probsoc BSIanx age
. * wave 2 BFdepx1 BSIdep BSItotal BSIposymp HP2probsoc BSIanx WHPpain age
. * wave 3 BFdepx1 BSIdep BSItotal BSIposymp HP2probsoc BSIanx WHPpain age
. *------------------------------------------------------------------------------
.
.
.
. *-+- Wave one chornobyl related health threat and female sleep relationship
.