.
. *------ H8: Physical ability as a function of crhtw
. /*
> EQ( 4) Modelling WHPpa 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 3 observation(s) with missing values from the sample
>
>                   Coefficient  Std.Error  t-value  t-prob Part.R^2
> BSIsoma               1.13862     0.1624     7.01  0.0000   0.1321
> BSIdep               0.785203     0.2309     3.40  0.0008   0.0346
> BSIhos               -1.04403     0.2655    -3.93  0.0001   0.0457
> childw1              -1.70882     0.8279    -2.06  0.0398   0.0130
> emplw33               7.24315      1.888     3.84  0.0002   0.0436
> deaw1                 3.83888      1.257     3.05  0.0025   0.0280
> illw3                 2.76037     0.7702     3.58  0.0004   0.0382
> shfamw1            -0.0363219    0.01900    -1.91  0.0568   0.0112
> shhousw1            0.0943166    0.02406     3.92  0.0001   0.0454
> shrelaw2           -0.0351518    0.02429    -1.45  0.1488   0.0064
> sufamw2             0.0588756    0.01588     3.71  0.0002   0.0408
> suchrw1              0.147311    0.05851     2.52  0.0123   0.0192
> phlthw3             -0.106420    0.02165    -4.91  0.0000   0.0696
> physdisagw2         -0.141342    0.04688    -3.02  0.0028   0.0274
>
> sigma                 10.9813  RSS                38949.9482
> log-likelihood       -1278.55
> no. of observations       337  no. of parameters          14
> mean(WHPpa)            9.5427  se(WHPpa)             14.5571
> When the log-likelihood constant is NOT included:
> AIC                   4.83304  SC                    4.99173
> HQ                    4.89629  FPE                   125.598
> When the log-likelihood constant is included:
> AIC                   7.67091  SC                    7.82961
> HQ                    7.73417  FPE                   2145.14
>
> Normality test:   Chi^2(2)  =   79.671 [0.0000]**
> Hetero test:      F(27,309) =   1.2388 [0.1963]
> Hetero-X test:    F(105,231)=   1.8860 [0.0000]**
> RESET23 test:     F(2,321)  =   7.4221 [0.0007]**
>
> Robust standard errors
>               Coefficients           SE        HACSE         HCSE        JHCSE
> BSIsoma             1.1386      0.16239      0.18603      0.17424      0.18438
> BSIdep             0.78520      0.23090      0.27351      0.25808      0.27164
> BSIhos             -1.0440      0.26549      0.30466      0.29082      0.30433
> childw1            -1.7088      0.82790      0.83367      0.79466      0.82170
> emplw33             7.2432       1.8879       2.2931       2.2673       2.3631
> deaw1               3.8389       1.2574       1.5782       1.6798       1.7874
> illw3               2.7604      0.77019       1.0096      0.90800      0.97151
> shfamw1          -0.036322     0.019002     0.020292     0.020705     0.021557
> shhousw1          0.094317     0.024061     0.035886     0.033831     0.035250
> shrelaw2         -0.035152     0.024287     0.029372     0.027854     0.029034
> sufamw2           0.058876     0.015879     0.016071     0.015954     0.016687
> suchrw1            0.14731     0.058513      0.10492      0.10099      0.12569
> phlthw3           -0.10642     0.021654     0.018411     0.019472     0.020527
> physdisagw2       -0.14134     0.046878     0.039933     0.040963     0.043172
>
>               Coefficients         t-SE      t-HACSE       t-HCSE      t-JHCSE
> BSIsoma             1.1386       7.0118       6.1207       6.5349       6.1755
> BSIdep             0.78520       3.4006       2.8708       3.0425       2.8906
> BSIhos             -1.0440      -3.9325      -3.4269      -3.5900      -3.4306
> childw1            -1.7088      -2.0640      -2.0498      -2.1504      -2.0796
> emplw33             7.2432       3.8366       3.1586       3.1946       3.0651
> deaw1               3.8389       3.0529       2.4324       2.2853       2.1478
> illw3               2.7604       3.5840       2.7343       3.0401       2.8413
> shfamw1          -0.036322      -1.9115      -1.7900      -1.7543      -1.6849
> shhousw1          0.094317       3.9199       2.6282       2.7878       2.6756
> shrelaw2         -0.035152      -1.4474      -1.1968      -1.2620      -1.2107
> sufamw2           0.058876       3.7079       3.6636       3.6902       3.5282
> suchrw1            0.14731       2.5176       1.4040       1.4586       1.1720
> phlthw3           -0.10642      -4.9145      -5.7803      -5.4652      -5.1844
> physdisagw2       -0.14134      -3.0151      -3.5395      -3.4505      -3.2739
>
> WHPpa =  + 1.14*BSIsoma + 0.785*BSIdep - 1.04*BSIhos - 1.71*childw1
> (SE)      (0.162)        (0.231)        (0.265)       (0.828)
>          + 7.24*emplw33 + 3.84*deaw1 + 2.76*illw3 - 0.0363*shfamw1
>           (1.89)         (1.26)       (0.77)       (0.019)
>          + 0.0943*shhousw1 - 0.0352*shrelaw2 + 0.0589*sufamw2
>           (0.0241)          (0.0243)          (0.0159)
>          + 0.147*suchrw1 - 0.106*phlthw3 - 0.141*physdisagw2
>           (0.0585)        (0.0217)        (0.0469)
>
> /*
> EQ( 2) Modelling WHPpa 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
>
>                   Coefficient  Std.Error  t-value  t-prob Part.R^2
> BSIsoma               1.33609     0.1591     8.40  0.0000   0.1742
> BSIdep                1.16750     0.2295     5.09  0.0000   0.0719
> BSIhos               -1.46692     0.2616    -5.61  0.0000   0.0860
> deaw1                 5.22227      1.287     4.06  0.0001   0.0470
> sufamw2             0.0634838    0.01584     4.01  0.0001   0.0459
> phlthw3             -0.128363    0.02032    -6.32  0.0000   0.1067
>
> sigma                 11.7675  RSS                46250.2562
> log-likelihood       -1317.63
> no. of observations       340  no. of parameters           6
> mean(WHPpa)           9.54859  se(WHPpa)             14.5562
> When the log-likelihood constant is NOT included:
> AIC                   4.94817  SC                    5.01574
> HQ                    4.97509  FPE                   140.917
> When the log-likelihood constant is included:
> AIC                   7.78605  SC                    7.85362
> HQ                    7.81297  FPE                   2406.80
>
> Normality test:   Chi^2(2)  =   100.71 [0.0000]**
> Hetero test:      F(12,327) =   1.3791 [0.1740]
> Hetero-X test:    F(27,312) =   1.0351 [0.4202]
> RESET23 test:     F(2,332)  =   3.7470 [0.0246]*
>
> */
>
>
> cap gen whppa = WHPpa
> // candidate illness variables for Physical ability
> /* candiate illnesses  for males
>
> foreach i in 7 {
> foreach j in 1 5 6 {
> foreach k in 3 {
> foreach m in 2 {
> foreach n in 4 {
> regress whppa icdx`i'nr1-icdx`i'nr28 if gender==1
> regress whppa icdx`j'nr1-icdx`j'nr18 if gender==1
> regress whppa icdx`k'nr1-icdx`k'nr11 if gender==1
> regress whppa icdx`m'nr1-icdx`m'nr8 if gender==1
> regress whppa icdx`n'nr1-icdx`n'nr20 if gender==1
> }
> }
> }
> }
> }
>
> */
>
>
> sjlog using H8PAtest, replace
> *-------Full male Physical ability crhtw model  **************/
. des icdx1nr5 icdx1nr18 icdx2nr4 icdx3nr3 icdx3nr7 icdx4nr8-icdx4nr10 icdx4nr12  ///
>   icdx5nr7 icdx5nr10 icdx6nr7 icdx7nr1

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


icdx1nr5          double   %8.0g                        icdx1nr==401 hypertension


icdx1nr18         double   %8.0g                        icdx1nr==732 osteochondropathies


icdx2nr4          double   %8.0g                        icdx2nr==401 hypertension


icdx3nr3          double   %8.0g                        icdx3nr==diabetes militus


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


icdx4nr8          double   %8.0g                        icdx4nr==angina pectoris


icdx4nr9          double   %8.0g                        icdx4nr==434.91 crbrl art ocl nos w infarc


icdx4nr10         double   %8.0g                        icdx4nr==varicose veins in legs


icdx4nr12         double   %8.0g                        icdx4nr==gastritis/duodenitis


icdx5nr7          double   %8.0g                        icdx5nr==angina pectoris


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


icdx6nr7          double   %8.0g                        icdx6nr==varicose veins in legs


icdx7nr1          double   %8.0g                        icdx7nr==malignant melanoma skin


.
. regress WHPpa age BSIsoma BSIdep deaw1 physdisagw2 ///
>    illw3 BSIanx medcow1 WHPpain crhtw1-crhtw3 ///
>    icdx1nr5 icdx1nr18 icdx2nr4 icdx3nr3 icdx3nr7 ///
>    icdx4nr8-icdx4nr10 icdx4nr12 icdx5nr7 icdx5nr10 icdx6nr7 icdx7nr1 ///
>    if gender==1, vce(cluster id)

Linear regression                                      Number of obs =     338
                                                       F( 21,   337) =       .
                                                       Prob > F      =       .
                                                       R-squared     =  0.6024
                                                       Root MSE      =  9.5478

(Std. Err. adjusted for 338 clusters in id)
-------------+----------------------------------------------------------------
             |               Robust
        WHPpa |      Coef.    Std. Err.       t    P>|t|      [95% Con f. Interval]
-------------+----------------------------------------------------------------
         age |    .0545287    .0524031      1.04    0.299     -.0485496      .157607
     BSIsoma |    .5628542    .1857164      3.03    0.003      .1975448     .9281636
      BSIdep |    .3040762    .3011014      1.01    0.313     -.2881988     .8963511
       deaw1 |    2.266718    1.255357      1.81    0.072      -.202604      4.73604
 physdisagw2 |   -.1155609    .0393136     -2.94    0.004     -.1928918      -.03823
       illw3 |    1.565851    .8296602      1.89    0.060     -.0661142     3.197816
      BSIanx |   -.7255457     .307856     -2.36    0.019     -1.331107    -.1199843
     medcow1 |     .540166    .2349823      2.30    0.022      .0779493     1.002383
     WHPpain |    .4783801    .0540109      8.86    0.000      .3721391     .5846211
      crhtw1 |    1.008249    1.301731      0.77    0.439     -1.552293     3.568791
      crhtw2 |    1.057135    2.049181      0.52    0.606     -2.973662     5.087932
      crhtw3 |   -1.433471    1.804705     -0.79    0.428     -4.983377     2.116434
    icdx1nr5 |   -3.954297    2.596418     -1.52    0.129     -9.061525      1.15293
   icdx1nr18 |    3.780637    3.399621      1.11    0.267     -2.906512     10.46779
    icdx2nr4 |   -.7822253    3.536503     -0.22    0.825     -7.738627     6.174177
    icdx3nr3 |    11.29403    4.892569      2.31    0.022      1.670208     20.91785
    icdx3nr7 |    10.79598    12.74224      0.85    0.397     -14.26837     35.86034
    icdx4nr8 |    2.570385    4.893629      0.53    0.600     -7.055522     12.19629
    icdx4nr9 |   -17.53347    6.224572     -2.82    0.005     -29.77738     -5.28956
   icdx4nr10 |   -2.445988    3.073787     -0.80    0.427     -8.492214     3.600239
   icdx4nr12 |    16.12243    3.274253      4.92    0.000       9.68188     22.56298
    icdx5nr7 |    7.716505    3.671209      2.10    0.036      .4951339     14.93788
   icdx5nr10 |      14.423     2.06132      7.00    0.000      10.36833     18.47768
    icdx6nr7 |    22.84422    1.781839     12.82    0.000      19.33929     26.34914
    icdx7nr1 |    15.19774    2.458806      6.18    0.000       10.3612     20.03428
       _cons |   -2.796892    3.241522     -0.86    0.389     -9.173058     3.579273
-------------+----------------------------------------------------------------

.    di e(r2_a)
.5705122

.
. *-------Trimmed male Physical ability crhtw model  **************/
. des icdx1nr5 icdx1nr18 icdx3nr3  icdx4nr9 icdx4nr12  ///
>   icdx5nr7 icdx5nr10 icdx6nr7 icdx7nr1

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


icdx1nr5          double   %8.0g                        icdx1nr==401 hypertension


icdx1nr18         double   %8.0g                        icdx1nr==732 osteochondropathies


icdx3nr3          double   %8.0g                        icdx3nr==diabetes militus


icdx4nr9          double   %8.0g                        icdx4nr==434.91 crbrl art ocl nos w infarc


icdx4nr12         double   %8.0g                        icdx4nr==gastritis/duodenitis


icdx5nr7          double   %8.0g                        icdx5nr==angina pectoris


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


icdx6nr7          double   %8.0g                        icdx6nr==varicose veins in legs


icdx7nr1          double   %8.0g                        icdx7nr==malignant melanoma skin


.
. regress WHPpa age BSIsoma BSIdep deaw1 physdisagw2 ///
>    illw3 BSIanx medcow1 WHPpain crhtw1-crhtw3 ///
>     icdx3nr3 ///
>     icdx4nr9 icdx4nr12 icdx5nr7 icdx5nr10 icdx6nr7 icdx7nr1 ///
>    if gender==1, vce(cluster id)

Linear regression                                      Number of obs =     338
                                                       F( 15,   337) =       .
                                                       Prob > F      =       .
                                                       R-squared     =  0.5878
                                                       Root MSE      =  9.6286

(Std. Err. adjusted for 338 clusters in id)
-------------+----------------------------------------------------------------
             |               Robust
        WHPpa |      Coef.    Std. Err.       t    P>|t|      [95% Con f. Interval]
-------------+----------------------------------------------------------------
         age |    .0578442    .0597399      0.97    0.334     -.0596658     .1753543
     BSIsoma |    .5294401    .1869968      2.83    0.005      .1616121      .897268
      BSIdep |    .3052377    .3041135      1.00    0.316     -.2929623     .9034376
       deaw1 |    2.045296    1.306777      1.57    0.118     -.5251717     4.615763
 physdisagw2 |   -.1115784    .0384099     -2.90    0.004     -.1871319     -.036025
       illw3 |    1.827375    .8479948      2.15    0.032      .1593456     3.495405
      BSIanx |   -.7640778    .3060775     -2.50    0.013     -1.366141    -.1620147
     medcow1 |     .521033    .2245599      2.32    0.021      .0793172     .9627487
     WHPpain |    .4871103    .0573704      8.49    0.000      .3742612     .5999594
      crhtw1 |    1.401308    1.365234      1.03    0.305     -1.284147     4.086763
      crhtw2 |    .3348103    1.972859      0.17    0.865      -3.54586     4.215481
      crhtw3 |   -1.045175    1.796972     -0.58    0.561     -4.579869      2.48952
    icdx3nr3 |    10.71167    4.781384      2.24    0.026      1.306554     20.11679
    icdx4nr9 |   -17.75787    4.863508     -3.65    0.000     -27.32453    -8.191217
   icdx4nr12 |    16.40801    3.340722      4.91    0.000      9.836713      22.9793
    icdx5nr7 |    6.313072    3.440606      1.83    0.067     -.4546969     13.08084
   icdx5nr10 |    14.78105    1.392747     10.61    0.000      12.04148     17.52063
    icdx6nr7 |     22.5554    1.848801     12.20    0.000      18.91875     26.19204
    icdx7nr1 |    14.53545     2.90797      5.00    0.000      8.815386     20.25551
       _cons |   -2.409054    3.293573     -0.73    0.465     -8.887604     4.069497
-------------+----------------------------------------------------------------

.    di e(r2_a)
.56320566

. *-+-  No male threat physical ability relationship found
.
.
.
.
.
.
.
. *------  female Physical ability
. /* testing for candiate illnesses for female model
>
> foreach i in 7 {
> foreach j in 1 5 6 {
> foreach k in 3 {
> foreach m in 2 {
> foreach n in 4 {
> regress whppa icdx`i'nr1-icdx`i'nr28 if gender==2
> regress whppa icdx`j'nr1-icdx`j'nr18 if gender==2
> regress whppa icdx`k'nr1-icdx`k'nr11 if gender==2
> regress whppa icdx`m'nr1-icdx`m'nr8 if gender==2
> regress whppa icdx`n'nr1-icdx`n'nr20 if gender==2
> }
> }
> }
> }
> }
> */
.
. // Full female Physical ability model with candidate illnesses
. des icdx1nr12 icdx4nr5 icdx4nr12 icdx5nr4 icdx5nr5 icdx6nr5 icdx6nr6 icdx6nr12 ///
>   icdx6nr15 icdx7nr10 icdx7nr19 icdx7nr26

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


icdx1nr12         double   %8.0g                        icdx1nr==531 gastric ulcer


icdx4nr5          double   %8.0g                        icdx4nr==rheum fev w/o hrt involv


icdx4nr12         double   %8.0g                        icdx4nr==gastritis/duodenitis


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


icdx5nr5          double   %8.0g                        icdx5nr==hypertension


icdx6nr5          double   %8.0g                        icdx6nr==hypertension


icdx6nr6          double   %8.0g                        icdx6nr==acute myocardial infarct


icdx6nr12         double   %8.0g                        icdx6nr==575.1 cholecystitis


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


icdx7nr10         double   %8.0g                        icdx7nr==acute myocardial infarct


icdx7nr19         double   %8.0g                        icdx7nr==oth gallbladder disordr


icdx7nr26         double   %8.0g                        icdx7nr==intervertebral disc dis*


. regress WHPpa age marrw14 marrw15 CSprbslv injselfr ///
>   BSIanx illw3  medcow1 sepaw2 inc2w1 inc2w3 inc4w3 WHPpain ///
>   crhtw1 crhtw2 crhtw3  icdx1nr12  icdx1nr12 icdx4nr5 icdx4nr12 icdx5nr4 ///
>     icdx5nr5 icdx6nr5 icdx6nr6 icdx6nr12 ///
>   icdx6nr15 icdx7nr10 icdx7nr19 icdx7nr26 ///
>   if gender==2, vce(cluster id)
note: icdx1nr12 omitted because of collinearity
note: icdx6nr15 omitted because of collinearity

Linear regression                                      Number of obs =     360
                                                       F( 20,   359) =       .
                                                       Prob > F      =       .
                                                       R-squared     =  0.6941
                                                       Root MSE      =   12.35

(Std. Err. adjusted for 360 clusters in id)
-------------+----------------------------------------------------------------
             |               Robust
        WHPpa |      Coef.    Std. Err.       t    P>|t|      [95% Con f. Interval]
-------------+----------------------------------------------------------------
         age |     .342642    .0586434      5.84    0.000      .2273142     .4579697
     marrw14 |     20.0629    17.61741      1.14    0.256     -14.58339      54.7092
     marrw15 |    6.572265    10.32968      0.64    0.525     -13.74203     26.88656
    CSprbslv |   -.2490141    .1556114     -1.60    0.110     -.5550386     .0570105
    injselfr |    3.927676    1.402029      2.80    0.005      1.170455     6.684897
      BSIanx |   -.5985432    .2248685     -2.66    0.008     -1.040768    -.1563182
       illw3 |     1.64225    .7453057      2.20    0.028      .1765363     3.107964
     medcow1 |    .6925395    .2106715      3.29    0.001      .2782343     1.106845
      sepaw2 |   -9.833499    2.856635     -3.44    0.001     -15.45134    -4.215658
      inc2w1 |   -1.478954     1.50891     -0.98    0.328     -4.446368     1.488459
      inc2w3 |   -2.827452    1.413466     -2.00    0.046     -5.607166    -.0477392
      inc4w3 |   -8.222078    3.466326     -2.37    0.018     -15.03893    -1.405222
     WHPpain |    .5732947    .0405816     14.13    0.000      .4934872     .6531023
      crhtw1 |    1.179102    1.153757      1.02    0.307     -1.089869     3.448074
      crhtw2 |     2.46974    2.156825      1.15    0.253     -1.771859     6.711339
      crhtw3 |   -3.161378    1.911385     -1.65    0.099     -6.920295     .5975404
   icdx1nr12 |    6.717806    6.058974      1.11    0.268     -5.197735     18.63335
   icdx1nr12 |           0   (omitted)
    icdx4nr5 |   -.6352011    5.472252     -0.12    0.908      -11.3969      10.1265
   icdx4nr12 |    .7603838    3.107902      0.24    0.807     -5.351597     6.872364
    icdx5nr4 |    29.49139    3.320407      8.88    0.000       22.9615     36.02128
    icdx5nr5 |    5.528182    11.70114      0.47    0.637      -17.4832     28.53956
    icdx6nr5 |       9.906    2.129735      4.65    0.000      5.717677     14.09432
    icdx6nr6 |    -18.4459    4.310556     -4.28    0.000     -26.92301    -9.968787
   icdx6nr12 |    14.59511    5.834292      2.50    0.013      3.121424     26.06879
   icdx6nr15 |           0   (omitted)
   icdx7nr10 |    14.50939    10.50235      1.38    0.168     -6.144459     35.16325
   icdx7nr19 |    43.74917    12.49369      3.50    0.001      19.17915      68.3192
   icdx7nr26 |    14.34909    12.02749      1.19    0.234     -9.304108     38.00228
       _cons |    -1.96652    4.347224     -0.45    0.651     -10.51574     6.582704
-------------+----------------------------------------------------------------

.
. // Trimmed female Physical ability model with candidate illnesses
. des  WHPpa age marrw14 marrw15 CSprbslv injselfr ///
>                 BSIanx illw3  medcow1 sepaw2 inc2w1 inc2w3 inc4w3 WHPpain ///
>                 crhtw1 crhtw2 crhtw3    icdx5nr4 ///
>         icdx6nr5 icdx6nr6 icdx6nr12 ///
>                 icdx7nr19

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


WHPpa             double   %9.0g                        Wtd Health Profile Physical Ability Pt 1 Subscale


age               double   %8.0g                    *   Respondent's age


marrw14           byte     %8.0g                        marrw1==4. separated


marrw15           byte     %8.0g                        marrw1==5. divorced


CSprbslv          double   %9.0g                        Coping Problem Solving Subscale


injselfr          double   %9.0g         dum            Were u injured because of Chornobyl acc in 1986?


BSIanx            double   %9.0g                        Basic symptom inventory Anxiety subscale


illw3             double   %8.0g                        Total number of illnesses experienced in time period 1996-NOW


medcow1           double   %8.0g                        number of medical visits  for a medical condition per year 1976-1986


sepaw2            double   %8.0g                        Total number of separations experienced in time period 1987-1996


inc2w1            double   %15.0g        LABJ           Income is just sufficient for basic neccessities in 1986


inc2w3            double   %15.0g        LABJ           Income is just sufficient for basic neccessities NOW


inc4w3            double   %15.0g        LABJ           Income allows to comfortably afford luxury items NOW


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


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


crhtw2            float    %9.0g                        Chornobyl related health threat: wave 2 alpha = .8219889682935094


crhtw3            float    %9.0g                        Chornobyl related health threat: wave 3 alpha = .8347477221944793


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


icdx6nr5          double   %8.0g                        icdx6nr==hypertension


icdx6nr6          double   %8.0g                        icdx6nr==acute myocardial infarct


icdx6nr12         double   %8.0g                        icdx6nr==575.1 cholecystitis


icdx7nr19         double   %8.0g                        icdx7nr==oth gallbladder disordr


.
. regress WHPpa age marrw14 marrw15 CSprbslv injselfr ///
>         BSIanx illw3  medcow1 sepaw2 inc2w1 inc2w3 inc4w3 WHPpain ///
>         crhtw1 crhtw2 crhtw3    icdx5nr4 ///
>     icdx6nr5 icdx6nr6 icdx6nr12 ///
>         icdx7nr19 if gender==2, vce(cluster id)

Linear regression                                      Number of obs =     360
                                                       F( 17,   359) =       .
                                                       Prob > F      =       .
                                                       R-squared     =  0.6892
                                                       Root MSE      =  12.339

(Std. Err. adjusted for 360 clusters in id)
-------------+----------------------------------------------------------------
             |               Robust
        WHPpa |      Coef.    Std. Err.       t    P>|t|      [95% Con f. Interval]
-------------+----------------------------------------------------------------
         age |    .3482028     .059034      5.90    0.000      .2321069     .4642986
     marrw14 |    19.62981    17.55158      1.12    0.264     -14.88702     54.14664
     marrw15 |    7.810871    9.178617      0.85    0.395     -10.23974     25.86148
    CSprbslv |   -.2572574    .1528363     -1.68    0.093     -.5578243     .0433096
    injselfr |    3.826369    1.403084      2.73    0.007      1.067073     6.585665
      BSIanx |   -.6087282    .2229695     -2.73    0.007     -1.047219    -.1702377
       illw3 |    1.756534    .7445011      2.36    0.019      .2924022     3.220665
     medcow1 |    .6841996    .2016887      3.39    0.001      .2875597     1.080839
      sepaw2 |   -9.869247    2.905035     -3.40    0.001     -15.58227    -4.156222
      inc2w1 |   -1.056651    1.497261     -0.71    0.481     -4.001156     1.887854
      inc2w3 |   -3.014879    1.407928     -2.14    0.033     -5.783703     -.246056
      inc4w3 |   -8.317145    3.472502     -2.40    0.017     -15.14615    -1.488143
     WHPpain |    .5860507    .0407996     14.36    0.000      .5058145     .6662869
      crhtw1 |    1.085233    1.211396      0.90    0.371      -1.29709     3.467557
      crhtw2 |    2.663147    2.260917      1.18    0.240     -1.783158     7.109451
      crhtw3 |   -3.222642    1.933385     -1.67    0.096     -7.024826     .5795418
    icdx5nr4 |    28.21617    3.251006      8.68    0.000      21.82276     34.60957
    icdx6nr5 |    16.23052    6.727014      2.41    0.016      3.001214     29.45982
    icdx6nr6 |   -19.87407    4.324977     -4.60    0.000     -28.37955     -11.3686
   icdx6nr12 |    13.41546     5.62233      2.39    0.018      2.358615      24.4723
   icdx7nr19 |    49.72662     4.86081     10.23    0.000      40.16738     59.28586
       _cons |   -2.032377    4.308756     -0.47    0.637     -10.50595     6.441197
-------------+----------------------------------------------------------------

.  di e(r2_a)
.66987567

.
. *-+-  No female threat physical ability relationships significant
.
.