. est clear

. // Zero-Order Female Pain analysis
. regress WHPpain avgcumdosew1 if gender==2, vce(cluster id)

Linear regression                                      Number of obs =     363
                                                       F(  1,   362) =   11.99
                                                       Prob > F      =  0.0006
                                                       R-squared     =  0.0383
                                                       Root MSE      =  21.926

(Std. Err. adjusted for 363 clusters in id)
-------------+----------------------------------------------------------------
             |               Robust
      WHPpain |      Coef.    Std. Err.       t    P>|t|      [95% Con f. Interval]
-------------+----------------------------------------------------------------
avgcumdosew1 |    7.949532    2.295408      3.46    0.001      3.435524     12.46354
       _cons |    15.34436     1.37731     11.14    0.000      12.63583     18.05289
-------------+----------------------------------------------------------------

. eststo wave1

.
. regress WHPpain avgcumdosew2 if gender==2, vce(cluster id)

Linear regression                                      Number of obs =     363
                                                       F(  1,   362) =    5.81
                                                       Prob > F      =  0.0165
                                                       R-squared     =  0.0230
                                                       Root MSE      =  22.099

(Std. Err. adjusted for 363 clusters in id)
-------------+----------------------------------------------------------------
             |               Robust
      WHPpain |      Coef.    Std. Err.       t    P>|t|      [95% Con f. Interval]
-------------+----------------------------------------------------------------
avgcumdosew2 |    2.451314    1.017161      2.41    0.016      .4510282       4.4516
       _cons |    15.81625    1.444076     10.95    0.000      12.97642     18.65608
-------------+----------------------------------------------------------------

. eststo wave2

.
. regress WHPpain avgcumdosew3 if gender==2, vce(cluster id)

Linear regression                                      Number of obs =     363
                                                       F(  1,   362) =    5.19
                                                       Prob > F      =  0.0233
                                                       R-squared     =  0.0199
                                                       Root MSE      =  22.134

(Std. Err. adjusted for 363 clusters in id)
-------------+----------------------------------------------------------------
             |               Robust
      WHPpain |      Coef.    Std. Err.       t    P>|t|      [95% Con f. Interval]
-------------+----------------------------------------------------------------
avgcumdosew3 |    1.793888    .7873128      2.28    0.023       .245607     3.342169
       _cons |    15.85567     1.46564     10.82    0.000      12.97343     18.73791
-------------+----------------------------------------------------------------

. eststo wave3

.
.
. estout  *,      cells(b(star fmt(%9.3f)) se(par) t p)                ///
>                 stats(r2_a bic N, fmt(%9.3f %9.0g) labels(adjR^2))      ///
>                 legend collabels(none) varlabels(_cons Constant)  ///
>                 title("Zero-order regression coefficients" ///
>                  "of reconstructed dose on female pain")

Zero-order regression coefficients of reconstructed dose on female pain
------------------------------------------------------------
                    wave1           wave2           wave3
------------------------------------------------------------
avgcumdosew1        7.950***
                  (2.295)
                    3.463
                    0.001
avgcumdosew2                        2.451*
                                  (1.017)
                                    2.410
                                    0.016
avgcumdosew3                                        1.794*
                                                  (0.787)
                                                    2.278
                                                    0.023
Constant           15.344***       15.816***       15.856***
                  (1.377)         (1.444)         (1.466)
                   11.141          10.953          10.818
                    0.000           0.000           0.000
------------------------------------------------------------
adjR^2              0.036           0.020           0.017
bic              3281.569        3287.288        3288.427
N                     363             363             363
------------------------------------------------------------
* p<0.05, ** p<0.01, *** p<0.001