.
. set more off

. title "Female PTSD model set 1"


*******************************************************************************
*******************************************************************************
*****                                                                     *****
*****                                                                     *****
*****                      Female PTSD model set 1                        *****
*****                                                                     *****
*****                                                                     *****
*****                                            11 Aug 2012    15:01:19  *****
*******************************************************************************
*******************************************************************************


.
. regress PTSDw1 cumdose1 injselfr fdferw1 goferw1 injothr healthef ///
>    woman crhtw1 if gender==2, vce(cluster id)

Linear regression                                      Number of obs =     362
                                                       F(  8,   361) =    3.66
                                                       Prob > F      =  0.0004
                                                       R-squared     =  0.1000
                                                       Root MSE      =  31.913

(Std. Err. adjusted for 362 clusters in id)
-------------+----------------------------------------------------------------
             |               Robust
       PTSDw1 |      Coef.    Std. Err.       t    P>|t|      [95% Con f. Interval]
-------------+----------------------------------------------------------------
    cumdose1 |    12.45539    4.949922      2.52    0.012      2.721084     22.18969
    injselfr |    6.649535    4.492536      1.48    0.140     -2.185293     15.48436
     fdferw1 |    .2499567     .072792      3.43    0.001      .1068071     .3931062
     goferw1 |   -.1324596    .0771075     -1.72    0.087     -.2840959     .0191767
     injothr |   -3.806116    5.160595     -0.74    0.461     -13.95472      6.34249
    healthef |    .0422295    .0692399      0.61    0.542     -.0939348     .1783938
       woman |   -.0234387    .0786745     -0.30    0.766     -.1781566     .1312793
      crhtw1 |   -.5906547    1.912303     -0.31    0.758     -4.351308     3.169999
       _cons |    5.616458    5.039679      1.11    0.266     -4.294358     15.52727
-------------+----------------------------------------------------------------

. eststo PTSDwave1f

.
. regress PTSDw2 cumdose2 injselfr fdferw1 goferw1 injothr healthef ///
>    woman crhtw1 crhtw2 if gender==2, vce(cluster id)

Linear regression                                      Number of obs =     362
                                                       F(  9,   361) =    3.27
                                                       Prob > F      =  0.0008
                                                       R-squared     =  0.0531
                                                       Root MSE      =  8.1486

(Std. Err. adjusted for 362 clusters in id)
-------------+----------------------------------------------------------------
             |               Robust
       PTSDw2 |      Coef.    Std. Err.       t    P>|t|      [95% Con f. Interval]
-------------+----------------------------------------------------------------
    cumdose2 |    .4203283    .1918823      2.19    0.029      .0429807     .7976759
    injselfr |    -.451494    1.408818     -0.32    0.749     -3.222014     2.319026
     fdferw1 |     .032874    .0355397      0.92    0.356     -.0370169     .1027648
     goferw1 |   -.0267401    .0391623     -0.68    0.495     -.1037549     .0502748
     injothr |    2.264203    1.094003      2.07    0.039      .1127849     4.415622
    healthef |    .0291624    .0098886      2.95    0.003      .0097159     .0486089
       woman |    .0122656    .0180713      0.68    0.498     -.0232726     .0478039
      crhtw1 |    .4326842     .553304      0.78    0.435     -.6554196     1.520788
      crhtw2 |   -.1472218    .6196786     -0.24    0.812     -1.365855     1.071411
       _cons |   -1.751567    1.312371     -1.33    0.183      -4.33242     .8292854
-------------+----------------------------------------------------------------

. eststo PTSDwave2f

.
. regress MiPTSD cumdose3 injselfr fdferw1 goferw1 injothr healthef ///
>   woman crhtw1 crhtw2 crhtw3 if gender==2, vce(cluster id)

Linear regression                                      Number of obs =     362
                                                       F( 10,   361) =    9.36
                                                       Prob > F      =  0.0000
                                                       R-squared     =  0.2281
                                                       Root MSE      =  10.734

(Std. Err. adjusted for 362 clusters in id)
-------------+----------------------------------------------------------------
             |               Robust
       MiPTSD |      Coef.    Std. Err.       t    P>|t|      [95% Con f. Interval]
-------------+----------------------------------------------------------------
    cumdose3 |    .2722351    .5222604      0.52    0.603     -.7548197      1.29929
    injselfr |     4.03679    1.563304      2.58    0.010      .9624637     7.111116
     fdferw1 |     .065294    .0275992      2.37    0.019      .0110186     .1195695
     goferw1 |   -.0231637    .0305715     -0.76    0.449     -.0832842     .0369569
     injothr |    1.352724    1.903614      0.71    0.478     -2.390842      5.09629
    healthef |    .0770909    .0226322      3.41    0.001      .0325833     .1215985
       woman |   -.0254135    .0269797     -0.94    0.347     -.0784706     .0276436
      crhtw1 |   -4.012422    .9683655     -4.14    0.000     -5.916768    -2.108076
      crhtw2 |    3.552389    2.319281      1.53    0.126      -1.00861     8.113388
      crhtw3 |      .69108    2.069107      0.33    0.739     -3.377937     4.760097
       _cons |    39.64146    2.103862     18.84    0.000       35.5041     43.77883
-------------+----------------------------------------------------------------

. eststo MiPTSDwave3f

.
.
.
. estout PTSDwave1f PTSDwave2f MiPTSDwave3f, cells(b(star fmt(3)) ///
>   se(par fmt(3)) t(par fmt(2)) p(par fmt(3))) ///
>   stats(r2 r2_a bic N)  starlevels(# 0.10 * 0.05 ** 0.01 *** 0.001 **** 0.0001) ///
>   title("Female PTSD regression analysis set 1" ///
>   "with cluster robust standard errors" ) legend

Female PTSD regression analysis set 1 with cluster robust standard errors
---------------------------------------------------------------
               PTSDwave1f       PTSDwave2f     MiPTSDwave3f
                 b/se/t/p         b/se/t/p         b/se/t/p
---------------------------------------------------------------
cumdose1           12.455*
                  (4.950)
                   (2.52)
                  (0.012)
injselfr            6.650           -0.451            4.037*
                  (4.493)          (1.409)          (1.563)
                   (1.48)          (-0.32)           (2.58)
                  (0.140)          (0.749)          (0.010)
fdferw1             0.250***         0.033            0.065*
                  (0.073)          (0.036)          (0.028)
                   (3.43)           (0.92)           (2.37)
                  (0.001)          (0.356)          (0.019)
goferw1            -0.132#          -0.027           -0.023
                  (0.077)          (0.039)          (0.031)
                  (-1.72)          (-0.68)          (-0.76)
                  (0.087)          (0.495)          (0.449)
injothr            -3.806            2.264*           1.353
                  (5.161)          (1.094)          (1.904)
                  (-0.74)           (2.07)           (0.71)
                  (0.461)          (0.039)          (0.478)
healthef            0.042            0.029**          0.077***
                  (0.069)          (0.010)          (0.023)
                   (0.61)           (2.95)           (3.41)
                  (0.542)          (0.003)          (0.001)
woman              -0.023            0.012           -0.025
                  (0.079)          (0.018)          (0.027)
                  (-0.30)           (0.68)          (-0.94)
                  (0.766)          (0.498)          (0.347)
crhtw1             -0.591            0.433           -4.012****
                  (1.912)          (0.553)          (0.968)
                  (-0.31)           (0.78)          (-4.14)
                  (0.758)          (0.435)          (0.000)
cumdose2                             0.420*
                                   (0.192)
                                    (2.19)
                                   (0.029)
crhtw2                              -0.147            3.552
                                   (0.620)          (2.319)
                                   (-0.24)           (1.53)
                                   (0.812)          (0.126)
cumdose3                                              0.272
                                                    (0.522)
                                                     (0.52)
                                                    (0.603)
crhtw3                                                0.691
                                                    (2.069)
                                                     (0.33)
                                                    (0.739)
_cons               5.616           -1.752           39.641****
                  (5.040)          (1.312)          (2.104)
                   (1.11)          (-1.33)          (18.84)
                  (0.266)          (0.183)          (0.000)
---------------------------------------------------------------
r2                  0.100            0.053            0.228
r2_a                0.080            0.029            0.206
bic              3578.455         2594.925         2799.284
N                 362.000          362.000          362.000
---------------------------------------------------------------
# p<0.10, * p<0.05, ** p<0.01, *** p<0.001, **** p<0.0001

.
. title "Female PTSD model set 2"


*******************************************************************************
*******************************************************************************
*****                                                                     *****
*****                                                                     *****
*****                      Female PTSD model set 2                        *****
*****                                                                     *****
*****                                                                     *****
*****                                            11 Aug 2012    15:01:19  *****
*******************************************************************************
*******************************************************************************


.
. regress MiPTSD cumdose1 injselfr fdferw1 goferw1 injothr healthef ///
>    woman crhtw1 if gender==2, vce(cluster id)

Linear regression                                      Number of obs =     362
                                                       F(  8,   361) =   11.13
                                                       Prob > F      =  0.0000
                                                       R-squared     =  0.2101
                                                       Root MSE      =  10.828

(Std. Err. adjusted for 362 clusters in id)
-------------+----------------------------------------------------------------
             |               Robust
       MiPTSD |      Coef.    Std. Err.       t    P>|t|      [95% Con f. Interval]
-------------+----------------------------------------------------------------
    cumdose1 |    2.944881     1.35381      2.18    0.030      .2825361     5.607226
    injselfr |    5.507291    1.588811      3.47    0.001      2.382804     8.631779
     fdferw1 |    .0704175     .027756      2.54    0.012      .0158337     .1250012
     goferw1 |   -.0258756    .0304101     -0.85    0.395     -.0856788     .0339276
     injothr |    2.616343     1.91557      1.37    0.173     -1.150735      6.38342
    healthef |    .0950178    .0239927      3.96    0.000      .0478347     .1422009
       woman |   -.0066442    .0273885     -0.24    0.808     -.0605053     .0472168
      crhtw1 |   -2.031846    .7086662     -2.87    0.004     -3.425478    -.6382132
       _cons |    34.80656    1.843515     18.88    0.000      31.18118     38.43194
-------------+----------------------------------------------------------------

. eststo MiPTSDwave1f

.
. regress MiPTSD cumdose2 injselfr fdferw1 goferw1 injothr healthef ///
>    woman crhtw1 crhtw2 if gender==2, vce(cluster id)

Linear regression                                      Number of obs =     362
                                                       F(  9,   361) =   10.22
                                                       Prob > F      =  0.0000
                                                       R-squared     =  0.2292
                                                       Root MSE      =  10.711

(Std. Err. adjusted for 362 clusters in id)
-------------+----------------------------------------------------------------
             |               Robust
       MiPTSD |      Coef.    Std. Err.       t    P>|t|      [95% Con f. Interval]
-------------+----------------------------------------------------------------
    cumdose2 |    .4660976    .6264163      0.74    0.457     -.7657859     1.697981
    injselfr |    4.091384    1.568473      2.61    0.009      1.006892     7.175877
     fdferw1 |    .0656592    .0274735      2.39    0.017       .011631     .1196874
     goferw1 |   -.0238586     .030318     -0.79    0.432     -.0834808     .0357635
     injothr |    1.448933    1.864944      0.78    0.438     -2.218586     5.116451
    healthef |    .0778882    .0227057      3.43    0.001      .0332362     .1225402
       woman |   -.0242043    .0268647     -0.90    0.368     -.0770352     .0286266
      crhtw1 |   -4.069685    .9334801     -4.36    0.000     -5.905427    -2.233943
      crhtw2 |    4.174644    1.063426      3.93    0.000      2.083355     6.265932
       _cons |    39.32613     2.00585     19.61    0.000      35.38152     43.27075
-------------+----------------------------------------------------------------

. eststo MiPTSDwave2f

.
. regress MiPTSD cumdose3 injselfr fdferw1 goferw1 injothr healthef ///
>   woman crhtw1 crhtw2 crhtw3 if gender==2, vce(cluster id)

Linear regression                                      Number of obs =     362
                                                       F( 10,   361) =    9.36
                                                       Prob > F      =  0.0000
                                                       R-squared     =  0.2281
                                                       Root MSE      =  10.734

(Std. Err. adjusted for 362 clusters in id)
-------------+----------------------------------------------------------------
             |               Robust
       MiPTSD |      Coef.    Std. Err.       t    P>|t|      [95% Con f. Interval]
-------------+----------------------------------------------------------------
    cumdose3 |    .2722351    .5222604      0.52    0.603     -.7548197      1.29929
    injselfr |     4.03679    1.563304      2.58    0.010      .9624637     7.111116
     fdferw1 |     .065294    .0275992      2.37    0.019      .0110186     .1195695
     goferw1 |   -.0231637    .0305715     -0.76    0.449     -.0832842     .0369569
     injothr |    1.352724    1.903614      0.71    0.478     -2.390842      5.09629
    healthef |    .0770909    .0226322      3.41    0.001      .0325833     .1215985
       woman |   -.0254135    .0269797     -0.94    0.347     -.0784706     .0276436
      crhtw1 |   -4.012422    .9683655     -4.14    0.000     -5.916768    -2.108076
      crhtw2 |    3.552389    2.319281      1.53    0.126      -1.00861     8.113388
      crhtw3 |      .69108    2.069107      0.33    0.739     -3.377937     4.760097
       _cons |    39.64146    2.103862     18.84    0.000       35.5041     43.77883
-------------+----------------------------------------------------------------

. eststo MiPTSDwave3f

.
.
. estout MiPTSDwave1f MiPTSDwave2f MiPTSDwave3f, ///
>   cells(b(star fmt(3)) se(par fmt(3)) t(par fmt(2)) p(par fmt(3))) ///
>   stats(r2 r2_a bic N)  starlevels(# 0.10 * 0.05 ** 0.01 *** 0.001 **** 0.0001) ///
>   title("Female PTSD regression analysis set 1" ///
>   "with cluster robust standard errors" ) legend

Female PTSD regression analysis set 1 with cluster robust standard errors
---------------------------------------------------------------
             MiPTSDwave1f     MiPTSDwave2f     MiPTSDwave3f
                 b/se/t/p         b/se/t/p         b/se/t/p
---------------------------------------------------------------
cumdose1            2.945*
                  (1.354)
                   (2.18)
                  (0.030)
injselfr            5.507***         4.091**          4.037*
                  (1.589)          (1.568)          (1.563)
                   (3.47)           (2.61)           (2.58)
                  (0.001)          (0.009)          (0.010)
fdferw1             0.070*           0.066*           0.065*
                  (0.028)          (0.027)          (0.028)
                   (2.54)           (2.39)           (2.37)
                  (0.012)          (0.017)          (0.019)
goferw1            -0.026           -0.024           -0.023
                  (0.030)          (0.030)          (0.031)
                  (-0.85)          (-0.79)          (-0.76)
                  (0.395)          (0.432)          (0.449)
injothr             2.616            1.449            1.353
                  (1.916)          (1.865)          (1.904)
                   (1.37)           (0.78)           (0.71)
                  (0.173)          (0.438)          (0.478)
healthef            0.095****        0.078***         0.077***
                  (0.024)          (0.023)          (0.023)
                   (3.96)           (3.43)           (3.41)
                  (0.000)          (0.001)          (0.001)
woman              -0.007           -0.024           -0.025
                  (0.027)          (0.027)          (0.027)
                  (-0.24)          (-0.90)          (-0.94)
                  (0.808)          (0.368)          (0.347)
crhtw1             -2.032**         -4.070****       -4.012****
                  (0.709)          (0.933)          (0.968)
                  (-2.87)          (-4.36)          (-4.14)
                  (0.004)          (0.000)          (0.000)
cumdose2                             0.466
                                   (0.626)
                                    (0.74)
                                   (0.457)
crhtw2                               4.175***         3.552
                                   (1.063)          (2.319)
                                    (3.93)           (1.53)
                                   (0.000)          (0.126)
cumdose3                                              0.272
                                                    (0.522)
                                                     (0.52)
                                                    (0.603)
crhtw3                                                0.691
                                                    (2.069)
                                                     (0.33)
                                                    (0.739)
_cons              34.807****       39.326****       39.641****
                  (1.844)          (2.006)          (2.104)
                  (18.88)          (19.61)          (18.84)
                  (0.000)          (0.000)          (0.000)
---------------------------------------------------------------
r2                  0.210            0.229            0.228
r2_a                0.192            0.209            0.206
bic              2795.864         2792.906         2799.284
N                 362.000          362.000          362.000
---------------------------------------------------------------
# p<0.10, * p<0.05, ** p<0.01, *** p<0.001, **** p<0.0001