-------------------------------------------------------------------------------
  log type:  smcl
 opened on:  16 Apr 2026, 14:06:25

. do "C:\Users\sabri\AppData\Local\Temp\STD71b0_00000b.tmp"

. 
. use "C:\Users\sabri\Dropbox\ReUP\Emotionen-Experiment\Data\reup_reduced_SRM.d
> ta", clear
(Welle 4: Die Ukraine und die Europäische Union in Zeiten des Krieges (SID: 148
> 999))

. **# Table 1
. tab wave2

      wave2 |      Freq.     Percent        Cum.
------------+-----------------------------------
          1 |      1,980      100.00      100.00
------------+-----------------------------------
      Total |      1,980      100.00

. tab wave3

      wave3 |      Freq.     Percent        Cum.
------------+-----------------------------------
          1 |      1,542      100.00      100.00
------------+-----------------------------------
      Total |      1,542      100.00

. tab wave4

      wave4 |      Freq.     Percent        Cum.
------------+-----------------------------------
          1 |      1,649      100.00      100.00
------------+-----------------------------------
      Total |      1,649      100.00

. tab wave4a

     wave4a |      Freq.     Percent        Cum.
------------+-----------------------------------
          0 |     16,845       93.62       93.62
          1 |      1,147        6.38      100.00
------------+-----------------------------------
      Total |     17,992      100.00

. 
. 
. tab emo_exp_w2 if wave2 == 1

{if(is_empty |
(randnumber) |
, rand(1,5), |
randnumber.N |
       AOK)} |      Freq.     Percent        Cum.
-------------+-----------------------------------
control cond |      1,582       79.90       79.90
    exp cond |        398       20.10      100.00
-------------+-----------------------------------
       Total |      1,980      100.00

. tab emo_exp_w3 if wave3 == 1

{if(is_empty |
(rand1bild), |
  rand(1,2), |
rand1bild.NA |
        OK)} |      Freq.     Percent        Cum.
-------------+-----------------------------------
control cond |        791       51.30       51.30
    exp cond |        751       48.70      100.00
-------------+-----------------------------------
       Total |      1,542      100.00

. tab emo_exp_w4a if wave4a == 1 & wave4 == 1

 emo_exp_w4a |      Freq.     Percent        Cum.
-------------+-----------------------------------
control cond |        548       51.75       51.75
    exp cond |        511       48.25      100.00
-------------+-----------------------------------
       Total |      1,059      100.00

. tab emo_exp_w4b if wave4a == 0 & wave4 == 1

 emo_exp_w4b |      Freq.     Percent        Cum.
-------------+-----------------------------------
control cond |        283       48.79       48.79
    exp cond |        297       51.21      100.00
-------------+-----------------------------------
       Total |        580      100.00

. 
. **# Table SI
. 
. ** Surveystats
. sum zeit_w2 if wave2 == 1, det

                           zeit_w2
-------------------------------------------------------------
      Percentiles      Smallest
 1%     .9608334          .2455
 5%     14.62167       .3348333
10%     16.78158       .3578333       Obs               1,980
25%     20.57925       .4023333       Sum of wgt.       1,980

50%     26.22958                      Mean           33.68195
                        Largest       Std. dev.       35.4356
75%        35.19         493.45
90%     49.81008         496.81       Variance       1255.682
95%     67.93733       533.2267       Skewness        8.16264
99%     188.8333       591.6583       Kurtosis       97.11025

. sum zeit_w3 if wave3 == 1, det

                           zeit_w3
-------------------------------------------------------------
      Percentiles      Smallest
 1%     6.904666       4.940167
 5%        8.207         5.3355
10%     9.300833         5.3355       Obs               1,542
25%      11.1985         5.3355       Sum of wgt.       1,542

50%     14.59408                      Mean           19.03242
                        Largest       Std. dev.      22.22968
75%       19.538         245.74
90%       27.433       271.6433       Variance       494.1586
95%     39.78783       271.6433       Skewness        9.18386
99%     111.7033       434.5583       Kurtosis       122.9736

. sum zeit_w4 if wave4 == 1, det

                           zeit_w4
-------------------------------------------------------------
      Percentiles      Smallest
 1%       5.4135       .2073333
 5%     7.639333       .2283333
10%     8.723499       .2606667       Obs               1,649
25%       11.111       .2878333       Sum of wgt.       1,649

50%     14.86117                      Mean           19.94633
                        Largest       Std. dev.       22.4414
75%     21.48117       217.2283
90%     32.59367        260.125       Variance       503.6166
95%     44.93666       291.9667       Skewness       8.333175
99%     108.4402       432.0016       Kurtosis       108.6139

. 
. tab intmig1_w2 if wave2 == 1,

  Befinden Sie |
sich gerade in |
 der Ukraine?  |      Freq.     Percent        Cum.
---------------+-----------------------------------
         1. Ja |      1,481       75.72       75.72
       2. Nein |        474       24.23       99.95
98. Weiß nicht |          1        0.05      100.00
---------------+-----------------------------------
         Total |      1,956      100.00

. tab intmig1_w3 if wave3 == 1,

   Befinden |
   Sie sich |
  gerade in |
        der |
  Ukraine?  |      Freq.     Percent        Cum.
------------+-----------------------------------
          1 |      1,178       76.64       76.64
          2 |        357       23.23       99.87
         98 |          2        0.13      100.00
------------+-----------------------------------
      Total |      1,537      100.00

. tab IntMig1 if wave4 == 1,

   Befinden |
   Sie sich |
  gerade in |
        der |
  Ukraine?  |      Freq.     Percent        Cum.
------------+-----------------------------------
         Ja |      1,263       77.34       77.34
       Nein |        368       22.54       99.88
 Weiß nicht |          2        0.12      100.00
------------+-----------------------------------
      Total |      1,633      100.00

. 
. 
. ** W2 Statistics for Table SI
. tab male if wave2 == 1

       male |      Freq.     Percent        Cum.
------------+-----------------------------------
          0 |      1,018       51.41       51.41
          1 |        962       48.59      100.00
------------+-----------------------------------
      Total |      1,980      100.00

. tab ukrainian if wave2 == 1

  ukrainian |      Freq.     Percent        Cum.
------------+-----------------------------------
          0 |      1,026       51.82       51.82
          1 |        954       48.18      100.00
------------+-----------------------------------
      Total |      1,980      100.00

. mean age if wave2==1

Mean estimation                          Number of obs = 1,945

--------------------------------------------------------------
             |       Mean   Std. err.     [95% conf. interval]
-------------+------------------------------------------------
         age |   39.87455   .3071233      39.27222    40.47688
--------------------------------------------------------------

. tab university if wave2 == 1

 university |      Freq.     Percent        Cum.
------------+-----------------------------------
          0 |        393       19.85       19.85
          1 |      1,587       80.15      100.00
------------+-----------------------------------
      Total |      1,980      100.00

. 
. ** W3 Statistics for Table SI
. tab male if wave3 == 1

       male |      Freq.     Percent        Cum.
------------+-----------------------------------
          0 |        764       49.55       49.55
          1 |        778       50.45      100.00
------------+-----------------------------------
      Total |      1,542      100.00

. tab ukrainian if wave3 == 1

  ukrainian |      Freq.     Percent        Cum.
------------+-----------------------------------
          0 |        813       52.72       52.72
          1 |        729       47.28      100.00
------------+-----------------------------------
      Total |      1,542      100.00

. mean age if wave3==1

Mean estimation                          Number of obs = 1,519

--------------------------------------------------------------
             |       Mean   Std. err.     [95% conf. interval]
-------------+------------------------------------------------
         age |   39.87097   .3432065      39.19776    40.54418
--------------------------------------------------------------

. tab university if wave3 == 1

 university |      Freq.     Percent        Cum.
------------+-----------------------------------
          0 |        306       19.84       19.84
          1 |      1,236       80.16      100.00
------------+-----------------------------------
      Total |      1,542      100.00

. 
. ** W4 Statistics for Table SI
. tab male if wave4 == 1

       male |      Freq.     Percent        Cum.
------------+-----------------------------------
          0 |        825       50.03       50.03
          1 |        824       49.97      100.00
------------+-----------------------------------
      Total |      1,649      100.00

. tab ukrainian if wave4 == 1

  ukrainian |      Freq.     Percent        Cum.
------------+-----------------------------------
          0 |        859       52.09       52.09
          1 |        790       47.91      100.00
------------+-----------------------------------
      Total |      1,649      100.00

. mean age if wave4==1

Mean estimation                          Number of obs = 1,612

--------------------------------------------------------------
             |       Mean   Std. err.     [95% conf. interval]
-------------+------------------------------------------------
         age |   40.90881   .3336379       40.2544    41.56322
--------------------------------------------------------------

. tab university if wave4 == 1

 university |      Freq.     Percent        Cum.
------------+-----------------------------------
          0 |        327       19.83       19.83
          1 |      1,322       80.17      100.00
------------+-----------------------------------
      Total |      1,649      100.00

. univar age if wave4 == 1
                                        -------------- Quantiles --------------
Variable       n     Mean     S.D.      Min      .25      Mdn      .75      Max
-------------------------------------------------------------------------------
     age    1612    40.91    13.40    18.00    30.00    40.00    50.00    79.00
-------------------------------------------------------------------------------

. 
. 
. 
. 
. **# Analyses Table 2
. 
. * Study 1
. ttest emo2_3_w2 if len_w2 >= 20 & wave2 == 1, by(emo_exp_w2)

Two-sample t test with equal variances
------------------------------------------------------------------------------
   Group |     Obs        Mean    Std. err.   Std. dev.   [95% conf. interval]
---------+--------------------------------------------------------------------
 control |   1,140    4.410526    .0515451    1.740363    4.309392     4.51166
exp cond |     256    4.804688    .1091855    1.746969    4.589667    5.019708
---------+--------------------------------------------------------------------
Combined |   1,396    4.482808     .046774    1.747623    4.391053    4.574563
---------+--------------------------------------------------------------------
    diff |           -.3941612    .1204514               -.6304468   -.1578756
------------------------------------------------------------------------------
    diff = mean(control) - mean(exp cond)                         t =  -3.2724
H0: diff = 0                                     Degrees of freedom =     1394

    Ha: diff < 0                 Ha: diff != 0                 Ha: diff > 0
 Pr(T < t) = 0.0005         Pr(|T| > |t|) = 0.0011          Pr(T > t) = 0.9995

. * Study 2
. ttest emo2_3_w3 if len_w3 >= 20 & wave3 == 1, by(emo_exp_w3)

Two-sample t test with equal variances
------------------------------------------------------------------------------
   Group |     Obs        Mean    Std. err.   Std. dev.   [95% conf. interval]
---------+--------------------------------------------------------------------
 control |     571    4.420315    .0740592    1.769688    4.274853    4.565777
exp cond |     551    4.272232    .0766646    1.799578    4.121641    4.422824
---------+--------------------------------------------------------------------
Combined |   1,122    4.347594    .0532946     1.78517    4.243025    4.452162
---------+--------------------------------------------------------------------
    diff |            .1480829    .1065619               -.0610005    .3571663
------------------------------------------------------------------------------
    diff = mean(control) - mean(exp cond)                         t =   1.3896
H0: diff = 0                                     Degrees of freedom =     1120

    Ha: diff < 0                 Ha: diff != 0                 Ha: diff > 0
 Pr(T < t) = 0.9175         Pr(|T| > |t|) = 0.1649          Pr(T > t) = 0.0825

. * Study 3
. ttest Emo2_SQ003 if wave4==1 & id < 1170 & len_w4 >=20, by(emo_exp_w4a)

Two-sample t test with equal variances
------------------------------------------------------------------------------
   Group |     Obs        Mean    Std. err.   Std. dev.   [95% conf. interval]
---------+--------------------------------------------------------------------
 control |     221    4.384615    .1217171    1.809455    4.144735    4.624496
exp cond |     160     4.84375    .1436679    1.817272    4.560006    5.127494
---------+--------------------------------------------------------------------
Combined |     381    4.577428    .0934728    1.824516    4.393639    4.761216
---------+--------------------------------------------------------------------
    diff |           -.4591346    .1881661               -.8291149   -.0891543
------------------------------------------------------------------------------
    diff = mean(control) - mean(exp cond)                         t =  -2.4400
H0: diff = 0                                     Degrees of freedom =      379

    Ha: diff < 0                 Ha: diff != 0                 Ha: diff > 0
 Pr(T < t) = 0.0076         Pr(|T| > |t|) = 0.0151          Pr(T > t) = 0.9924

. * Study 4
. ttest Emo2_SQ003 if wave4==1 & id >=1170, by(emo_exp_w4b)

Two-sample t test with equal variances
------------------------------------------------------------------------------
   Group |     Obs        Mean    Std. err.   Std. dev.   [95% conf. interval]
---------+--------------------------------------------------------------------
 control |     232    4.456897     .114089    1.737751    4.232109    4.681684
exp cond |     261    4.279693    .1137848     1.83825    4.055636    4.503751
---------+--------------------------------------------------------------------
Combined |     493    4.363083    .0807093    1.792038    4.204506    4.521661
---------+--------------------------------------------------------------------
    diff |            .1772031    .1616655               -.1404385    .4948447
------------------------------------------------------------------------------
    diff = mean(control) - mean(exp cond)                         t =   1.0961
H0: diff = 0                                     Degrees of freedom =      491

    Ha: diff < 0                 Ha: diff != 0                 Ha: diff > 0
 Pr(T < t) = 0.8632         Pr(|T| > |t|) = 0.2736          Pr(T > t) = 0.1368

. 
. 
. 
. 
. **# Analyses Table 2 revision
. 
. * Study 1
. ttest emo2_3_w2  if  wave2 == 1, by(emo_exp_w2)

Two-sample t test with equal variances
------------------------------------------------------------------------------
   Group |     Obs        Mean    Std. err.   Std. dev.   [95% conf. interval]
---------+--------------------------------------------------------------------
 control |   1,301     4.43505     .047784    1.723539    4.341308    4.528792
exp cond |     330    4.730303    .0954427    1.733802    4.542548    4.918058
---------+--------------------------------------------------------------------
Combined |   1,631    4.494788    .0428163    1.729163    4.410808    4.578769
---------+--------------------------------------------------------------------
    diff |           -.2952531    .1063593               -.5038686   -.0866376
------------------------------------------------------------------------------
    diff = mean(control) - mean(exp cond)                         t =  -2.7760
H0: diff = 0                                     Degrees of freedom =     1629

    Ha: diff < 0                 Ha: diff != 0                 Ha: diff > 0
 Pr(T < t) = 0.0028         Pr(|T| > |t|) = 0.0056          Pr(T > t) = 0.9972

. ttest emo2_3_w3  if wave3 == 1, by(emo_exp_w3)

Two-sample t test with equal variances
------------------------------------------------------------------------------
   Group |     Obs        Mean    Std. err.   Std. dev.   [95% conf. interval]
---------+--------------------------------------------------------------------
 control |     664    4.448795    .0681773    1.756805    4.314926    4.582665
exp cond |     627     4.30622    .0723261    1.811044    4.164189    4.448251
---------+--------------------------------------------------------------------
Combined |   1,291    4.379551    .0496538    1.784086     4.28214    4.476962
---------+--------------------------------------------------------------------
    diff |            .1425751    .0993076               -.0522471    .3373973
------------------------------------------------------------------------------
    diff = mean(control) - mean(exp cond)                         t =   1.4357
H0: diff = 0                                     Degrees of freedom =     1289

    Ha: diff < 0                 Ha: diff != 0                 Ha: diff > 0
 Pr(T < t) = 0.9243         Pr(|T| > |t|) = 0.1513          Pr(T > t) = 0.0757

. ttest Emo2_SQ003 if wave4==1 & id < 1170, by(emo_exp_w4a)

Two-sample t test with equal variances
------------------------------------------------------------------------------
   Group |     Obs        Mean    Std. err.   Std. dev.   [95% conf. interval]
---------+--------------------------------------------------------------------
 control |     252    4.325397    .1146427    1.819896    4.099613    4.551181
exp cond |     202    4.747525    .1275085    1.812236    4.496099    4.998951
---------+--------------------------------------------------------------------
Combined |     454    4.513216    .0857266    1.826601    4.344745    4.681687
---------+--------------------------------------------------------------------
    diff |           -.4221279    .1715481               -.7592587   -.0849972
------------------------------------------------------------------------------
    diff = mean(control) - mean(exp cond)                         t =  -2.4607
H0: diff = 0                                     Degrees of freedom =      452

    Ha: diff < 0                 Ha: diff != 0                 Ha: diff > 0
 Pr(T < t) = 0.0071         Pr(|T| > |t|) = 0.0142          Pr(T > t) = 0.9929

. ttest Emo2_SQ003 if wave4==1 & id >=1170, by(emo_exp_w4b)

Two-sample t test with equal variances
------------------------------------------------------------------------------
   Group |     Obs        Mean    Std. err.   Std. dev.   [95% conf. interval]
---------+--------------------------------------------------------------------
 control |     232    4.456897     .114089    1.737751    4.232109    4.681684
exp cond |     261    4.279693    .1137848     1.83825    4.055636    4.503751
---------+--------------------------------------------------------------------
Combined |     493    4.363083    .0807093    1.792038    4.204506    4.521661
---------+--------------------------------------------------------------------
    diff |            .1772031    .1616655               -.1404385    .4948447
------------------------------------------------------------------------------
    diff = mean(control) - mean(exp cond)                         t =   1.0961
H0: diff = 0                                     Degrees of freedom =      491

    Ha: diff < 0                 Ha: diff != 0                 Ha: diff > 0
 Pr(T < t) = 0.8632         Pr(|T| > |t|) = 0.2736          Pr(T > t) = 0.1368

. 
. 
. 
. ** needs sdtest, equal variances yes!
. sdtest emo2_3_w2  if wave2 == 1, by(emo_exp_w2)

Variance ratio test
------------------------------------------------------------------------------
   Group |     Obs        Mean    Std. err.   Std. dev.   [95% conf. interval]
---------+--------------------------------------------------------------------
 control |   1,301     4.43505     .047784    1.723539    4.341308    4.528792
exp cond |     330    4.730303    .0954427    1.733802    4.542548    4.918058
---------+--------------------------------------------------------------------
Combined |   1,631    4.494788    .0428163    1.729163    4.410808    4.578769
------------------------------------------------------------------------------
    ratio = sd(control) / sd(exp cond)                            f =   0.9882
H0: ratio = 1                                   Degrees of freedom = 1300, 329

    Ha: ratio < 1               Ha: ratio != 1                 Ha: ratio > 1
  Pr(F < f) = 0.4390         2*Pr(F < f) = 0.8780           Pr(F > f) = 0.5610

. sdtest emo2_3_w3  if wave3 == 1, by(emo_exp_w3)

Variance ratio test
------------------------------------------------------------------------------
   Group |     Obs        Mean    Std. err.   Std. dev.   [95% conf. interval]
---------+--------------------------------------------------------------------
 control |     664    4.448795    .0681773    1.756805    4.314926    4.582665
exp cond |     627     4.30622    .0723261    1.811044    4.164189    4.448251
---------+--------------------------------------------------------------------
Combined |   1,291    4.379551    .0496538    1.784086     4.28214    4.476962
------------------------------------------------------------------------------
    ratio = sd(control) / sd(exp cond)                            f =   0.9410
H0: ratio = 1                                    Degrees of freedom = 663, 626

    Ha: ratio < 1               Ha: ratio != 1                 Ha: ratio > 1
  Pr(F < f) = 0.2200         2*Pr(F < f) = 0.4400           Pr(F > f) = 0.7800

. sdtest Emo2_SQ003 if wave4a==1, by(emo_exp_w4a)

Variance ratio test
------------------------------------------------------------------------------
   Group |     Obs        Mean    Std. err.   Std. dev.   [95% conf. interval]
---------+--------------------------------------------------------------------
 control |     265    4.316981    .1124963    1.831308    4.095477    4.538485
exp cond |     215    4.734884    .1231592    1.805868    4.492123    4.977644
---------+--------------------------------------------------------------------
Combined |     480    4.504167    .0835243    1.829925    4.340047    4.668286
------------------------------------------------------------------------------
    ratio = sd(control) / sd(exp cond)                            f =   1.0284
H0: ratio = 1                                    Degrees of freedom = 264, 214

    Ha: ratio < 1               Ha: ratio != 1                 Ha: ratio > 1
  Pr(F < f) = 0.5832         2*Pr(F > f) = 0.8335           Pr(F > f) = 0.4168

. sdtest Emo2_SQ003 if wave4a==0, by(emo_exp_w4b)

Variance ratio test
------------------------------------------------------------------------------
   Group |     Obs        Mean    Std. err.   Std. dev.   [95% conf. interval]
---------+--------------------------------------------------------------------
 control |     235     4.46383    .1133435    1.737523    4.240526    4.687134
exp cond |     273     4.25641    .1114325    1.841167     4.03703     4.47579
---------+--------------------------------------------------------------------
Combined |     508    4.352362     .079649    1.795198    4.195879    4.508845
------------------------------------------------------------------------------
    ratio = sd(control) / sd(exp cond)                            f =   0.8906
H0: ratio = 1                                    Degrees of freedom = 234, 272

    Ha: ratio < 1               Ha: ratio != 1                 Ha: ratio > 1
  Pr(F < f) = 0.1805         2*Pr(F < f) = 0.3611           Pr(F > f) = 0.8195

. 
. 
. ** treatment and item-non response 
. tab1 break*

-> tabulation of break_w2  

   break_w2 |      Freq.     Percent        Cum.
------------+-----------------------------------
          0 |     17,964       99.78       99.78
          1 |         39        0.22      100.00
------------+-----------------------------------
      Total |     18,003      100.00

-> tabulation of break_w3  

   break_w3 |      Freq.     Percent        Cum.
------------+-----------------------------------
          0 |     17,981       99.88       99.88
          1 |         22        0.12      100.00
------------+-----------------------------------
      Total |     18,003      100.00

-> tabulation of break_w4a  

  break_w4a |      Freq.     Percent        Cum.
------------+-----------------------------------
          0 |     17,961       99.77       99.77
          1 |         42        0.23      100.00
------------+-----------------------------------
      Total |     18,003      100.00

-> tabulation of break_w4b  

  break_w4b |      Freq.     Percent        Cum.
------------+-----------------------------------
          0 |     17,993       99.94       99.94
          1 |         10        0.06      100.00
------------+-----------------------------------
      Total |     18,003      100.00

. tab mis_em_w2 if wave2 == 1

  mis_em_w2 |      Freq.     Percent        Cum.
------------+-----------------------------------
          0 |      1,631       82.37       82.37
          1 |        349       17.63      100.00
------------+-----------------------------------
      Total |      1,980      100.00

. tab mis_em_w3 if wave3 == 1

  mis_em_w3 |      Freq.     Percent        Cum.
------------+-----------------------------------
          0 |      1,291       83.72       83.72
          1 |        251       16.28      100.00
------------+-----------------------------------
      Total |      1,542      100.00

. tab mis_em_w4a if randmani==2 & wave4a == 1

 mis_em_w4a |      Freq.     Percent        Cum.
------------+-----------------------------------
          0 |        464       84.21       84.21
          1 |         87       15.79      100.00
------------+-----------------------------------
      Total |        551      100.00

. tab mis_em_w4b if wave4a == 0

 mis_em_w4b |      Freq.     Percent        Cum.
------------+-----------------------------------
          0 |        508       78.03       78.03
          1 |        143       21.97      100.00
------------+-----------------------------------
      Total |        651      100.00

. 
. 
. ttest mis_em_w2  if wave2 == 1, by(emo_exp_w2)

Two-sample t test with equal variances
------------------------------------------------------------------------------
   Group |     Obs        Mean    Std. err.   Std. dev.   [95% conf. interval]
---------+--------------------------------------------------------------------
 control |   1,582    .1776233    .0096121    .3823161    .1587694    .1964771
exp cond |     398    .1708543      .01889    .3768553    .1337172    .2079913
---------+--------------------------------------------------------------------
Combined |   1,980    .1762626    .0085655    .3811397    .1594643    .1930609
---------+--------------------------------------------------------------------
    diff |             .006769    .0213782               -.0351571    .0486951
------------------------------------------------------------------------------
    diff = mean(control) - mean(exp cond)                         t =   0.3166
H0: diff = 0                                     Degrees of freedom =     1978

    Ha: diff < 0                 Ha: diff != 0                 Ha: diff > 0
 Pr(T < t) = 0.6242         Pr(|T| > |t|) = 0.7516          Pr(T > t) = 0.3758

. ttest mis_em_w3  if wave3 == 1, by(emo_exp_w3)

Two-sample t test with equal variances
------------------------------------------------------------------------------
   Group |     Obs        Mean    Std. err.   Std. dev.   [95% conf. interval]
---------+--------------------------------------------------------------------
 control |     791    .1605563    .0130616    .3673534    .1349167    .1861958
exp cond |     751    .1651132    .0135573    .3715301    .1384984     .191728
---------+--------------------------------------------------------------------
Combined |   1,542    .1627756     .009404    .3692806    .1443296    .1812217
---------+--------------------------------------------------------------------
    diff |           -.0045569    .0188202               -.0414728    .0323589
------------------------------------------------------------------------------
    diff = mean(control) - mean(exp cond)                         t =  -0.2421
H0: diff = 0                                     Degrees of freedom =     1540

    Ha: diff < 0                 Ha: diff != 0                 Ha: diff > 0
 Pr(T < t) = 0.4044         Pr(|T| > |t|) = 0.8087          Pr(T > t) = 0.5956

. ttest mis_em_w4a if wave4a==1, by(emo_exp_w4a)

Two-sample t test with equal variances
------------------------------------------------------------------------------
   Group |     Obs        Mean    Std. err.   Std. dev.   [95% conf. interval]
---------+--------------------------------------------------------------------
 control |     304    .1480263    .0204015    .3557116    .1078798    .1881728
exp cond |     247    .1700405    .0239517    .3764311    .1228638    .2172171
---------+--------------------------------------------------------------------
Combined |     551    .1578947    .0155484    .3649736    .1273532    .1884362
---------+--------------------------------------------------------------------
    diff |           -.0220142    .0312789               -.0834551    .0394268
------------------------------------------------------------------------------
    diff = mean(control) - mean(exp cond)                         t =  -0.7038
H0: diff = 0                                     Degrees of freedom =      549

    Ha: diff < 0                 Ha: diff != 0                 Ha: diff > 0
 Pr(T < t) = 0.2409         Pr(|T| > |t|) = 0.4819          Pr(T > t) = 0.7591

. ttest mis_em_w4b if wave4a==0, by(emo_exp_w4b)

Two-sample t test with equal variances
------------------------------------------------------------------------------
   Group |     Obs        Mean    Std. err.   Std. dev.   [95% conf. interval]
---------+--------------------------------------------------------------------
 control |     315    .2539683    .0245643    .4359722    .2056369    .3022996
exp cond |     336       .1875     .021325    .3908945    .1455521    .2294479
---------+--------------------------------------------------------------------
Combined |     651    .2196621    .0162391     .414336    .1877746    .2515495
---------+--------------------------------------------------------------------
    diff |            .0664683    .0324153                .0028167    .1301198
------------------------------------------------------------------------------
    diff = mean(control) - mean(exp cond)                         t =   2.0505
H0: diff = 0                                     Degrees of freedom =      649

    Ha: diff < 0                 Ha: diff != 0                 Ha: diff > 0
 Pr(T < t) = 0.9796         Pr(|T| > |t|) = 0.0407          Pr(T > t) = 0.0204

. 
. 
. logit mis_em_w2 (c.emo_exp_w2)## (c.age i.male i.university ost)   if wave2 =
> = 1

Iteration 0:  Log likelihood = -910.61394  
Iteration 1:  Log likelihood = -905.17512  
Iteration 2:  Log likelihood = -905.10867  
Iteration 3:  Log likelihood = -905.10866  

Logistic regression                                     Number of obs =  1,945
                                                        LR chi2(9)    =  11.01
                                                        Prob > chi2   = 0.2750
Log likelihood = -905.10866                             Pseudo R2     = 0.0060

-------------------------------------------------------------------------------
    mis_em_w2 | Coefficient  Std. err.      z    P>|z|     [95% conf. interval]
--------------+----------------------------------------------------------------
   emo_exp_w2 |  -.1222981   .5768084    -0.21   0.832    -1.252822    1.008226
          age |   .0037507   .0140694     0.27   0.790    -.0238247    .0313262
       1.male |  -.0966866   .3830389    -0.25   0.801     -.847429    .6540558
 1.university |  -.1574798   .4698334    -0.34   0.737    -1.078336    .7633767
        1.ost |  -.8337653   .6523271    -1.28   0.201    -2.112303    .4447723
              |
 c.emo_exp_w2#|
        c.age |    .000286   .0112735     0.03   0.980    -.0218097    .0223817
              |
         male#|
 c.emo_exp_w2 |
           1  |   .2087565     .30446     0.69   0.493    -.3879741    .8054871
              |
   university#|
 c.emo_exp_w2 |
           1  |  -.1487144   .3746763    -0.40   0.691    -.8830665    .5856377
              |
          ost#|
 c.emo_exp_w2 |
           1  |   .7353015   .4888669     1.50   0.133      -.22286    1.693463
              |
        _cons |  -1.364162   .7107745    -1.92   0.055    -2.757255    .0289303
-------------------------------------------------------------------------------

. testparm c.emo_exp_w2#c.age c.emo_exp_w2#i.male c.emo_exp_w2#i.university c.e
> mo_exp_w2#i.ost

 ( 1)  [mis_em_w2]c.emo_exp_w2#c.age = 0
 ( 2)  [mis_em_w2]1.male#c.emo_exp_w2 = 0
 ( 3)  [mis_em_w2]1.university#c.emo_exp_w2 = 0
 ( 4)  [mis_em_w2]1.ost#c.emo_exp_w2 = 0

           chi2(  4) =    3.00
         Prob > chi2 =    0.5585

. 
. logit mis_em_w3 c.emo_exp_w3## (c.age i.male i.university ost)     if wave3 =
> = 1

Iteration 0:  Log likelihood = -677.65516  
Iteration 1:  Log likelihood = -671.14931  
Iteration 2:  Log likelihood = -671.03979  
Iteration 3:  Log likelihood = -671.03952  
Iteration 4:  Log likelihood = -671.03952  

Logistic regression                                     Number of obs =  1,519
                                                        LR chi2(9)    =  13.23
                                                        Prob > chi2   = 0.1524
Log likelihood = -671.03952                             Pseudo R2     = 0.0098

-------------------------------------------------------------------------------
    mis_em_w3 | Coefficient  Std. err.      z    P>|z|     [95% conf. interval]
--------------+----------------------------------------------------------------
   emo_exp_w3 |  -.2281208   .5436107    -0.42   0.675    -1.293578    .8373367
          age |  -.0153917    .016694    -0.92   0.357    -.0481114    .0173279
       1.male |   .2085921   .4448267     0.47   0.639    -.6632522    1.080436
 1.university |   1.061147   .6104294     1.74   0.082    -.1352729    2.257566
        1.ost |  -1.806654    1.01735    -1.78   0.076    -3.800623     .187316
              |
 c.emo_exp_w3#|
        c.age |   .0157054   .0104644     1.50   0.133    -.0048044    .0362153
              |
         male#|
 c.emo_exp_w3 |
           1  |  -.0006925   .2833275    -0.00   0.998    -.5560042    .5546193
              |
   university#|
 c.emo_exp_w3 |
           1  |  -.5882297   .3815902    -1.54   0.123    -1.336133    .1596733
              |
          ost#|
 c.emo_exp_w3 |
           1  |   1.018913   .5823823     1.75   0.080    -.1225355    2.160361
              |
        _cons |  -1.868017   .8576128    -2.18   0.029    -3.548907    -.187127
-------------------------------------------------------------------------------

. testparm c.emo_exp_w3#c.age c.emo_exp_w3#i.male c.emo_exp_w3#i.university c.e
> mo_exp_w3#i.ost

 ( 1)  [mis_em_w3]c.emo_exp_w3#c.age = 0
 ( 2)  [mis_em_w3]1.male#c.emo_exp_w3 = 0
 ( 3)  [mis_em_w3]1.university#c.emo_exp_w3 = 0
 ( 4)  [mis_em_w3]1.ost#c.emo_exp_w3 = 0

           chi2(  4) =    6.66
         Prob > chi2 =    0.1547

. 
. logit mis_em_w4a c.emo_exp_w4a## (c.age i.male i.university ost)   if wave4a 
> == 1

Iteration 0:  Log likelihood =  -238.4142  
Iteration 1:  Log likelihood = -233.79344  
Iteration 2:  Log likelihood = -233.68103  
Iteration 3:  Log likelihood = -233.68074  
Iteration 4:  Log likelihood = -233.68074  

Logistic regression                                     Number of obs =    540
                                                        LR chi2(9)    =   9.47
                                                        Prob > chi2   = 0.3953
Log likelihood = -233.68074                             Pseudo R2     = 0.0199

-------------------------------------------------------------------------------
   mis_em_w4a | Coefficient  Std. err.      z    P>|z|     [95% conf. interval]
--------------+----------------------------------------------------------------
  emo_exp_w4a |  -.1373119   .8890171    -0.15   0.877    -1.879753     1.60513
          age |   .0232107   .0270815     0.86   0.391    -.0298681    .0762896
       1.male |  -.4339933   .7491939    -0.58   0.562    -1.902386      1.0344
 1.university |  -1.317353    .869918    -1.51   0.130    -3.022361    .3876545
        1.ost |   1.147798   1.383063     0.83   0.407    -1.562956    3.858552
              |
           c. |
  emo_exp_w4a#|
        c.age |  -.0070178   .0179069    -0.39   0.695    -.0421148    .0280791
              |
         male#|
c.emo_exp_w4a |
           1  |    .427326   .4814505     0.89   0.375    -.5162996    1.370952
              |
   university#|
c.emo_exp_w4a |
           1  |   .5942808   .5727004     1.04   0.299    -.5281912    1.716753
              |
          ost#|
c.emo_exp_w4a |
           1  |  -1.086765   .9568981    -1.14   0.256    -2.962251     .788721
              |
        _cons |  -1.707098   1.380287    -1.24   0.216     -4.41241    .9982141
-------------------------------------------------------------------------------

. testparm c.emo_exp_w4a#c.age c.emo_exp_w4a#i.male c.emo_exp_w4a#i.university 
> c.emo_exp_w4a#i.ost

 ( 1)  [mis_em_w4a]c.emo_exp_w4a#c.age = 0
 ( 2)  [mis_em_w4a]1.male#c.emo_exp_w4a = 0
 ( 3)  [mis_em_w4a]1.university#c.emo_exp_w4a = 0
 ( 4)  [mis_em_w4a]1.ost#c.emo_exp_w4a = 0

           chi2(  4) =    2.91
         Prob > chi2 =    0.5724

. 
. logit mis_em_w4b c.emo_exp_w4b## (c.age i.male i.university ost)   if wave4a 
> == 0

Iteration 0:  Log likelihood = -339.48396  
Iteration 1:  Log likelihood = -334.78464  
Iteration 2:  Log likelihood = -334.72504  
Iteration 3:  Log likelihood = -334.72501  

Logistic regression                                     Number of obs =    643
                                                        LR chi2(9)    =   9.52
                                                        Prob > chi2   = 0.3909
Log likelihood = -334.72501                             Pseudo R2     = 0.0140

-------------------------------------------------------------------------------
   mis_em_w4b | Coefficient  Std. err.      z    P>|z|     [95% conf. interval]
--------------+----------------------------------------------------------------
  emo_exp_w4b |  -.3592619   .7535814    -0.48   0.634    -1.836254    1.117731
          age |   .0014218   .0218203     0.07   0.948    -.0413452    .0441888
       1.male |  -.6930778   .6019342    -1.15   0.250    -1.872847    .4866916
 1.university |   .2119488   .8041339     0.26   0.792    -1.364125    1.788022
        1.ost |  -1.931887   1.387262    -1.39   0.164     -4.65087    .7870958
              |
           c. |
  emo_exp_w4b#|
        c.age |   -.005983   .0144427    -0.41   0.679    -.0342902    .0223243
              |
         male#|
c.emo_exp_w4b |
           1  |   .3696473    .388123     0.95   0.341    -.3910598    1.130354
              |
   university#|
c.emo_exp_w4b |
           1  |   -.055645    .536651    -0.10   0.917    -1.107462    .9961717
              |
          ost#|
c.emo_exp_w4b |
           1  |   1.057754   .8348865     1.27   0.205    -.5785939    2.694101
              |
        _cons |  -.4298106     1.1451    -0.38   0.707    -2.674164    1.814543
-------------------------------------------------------------------------------

. testparm c.emo_exp_w4b#c.age c.emo_exp_w4b#i.male c.emo_exp_w4b#i.university 
> c.emo_exp_w4b#i.ost

 ( 1)  [mis_em_w4b]c.emo_exp_w4b#c.age = 0
 ( 2)  [mis_em_w4b]1.male#c.emo_exp_w4b = 0
 ( 3)  [mis_em_w4b]1.university#c.emo_exp_w4b = 0
 ( 4)  [mis_em_w4b]1.ost#c.emo_exp_w4b = 0

           chi2(  4) =    2.98
         Prob > chi2 =    0.5616

. 
. ** treatment and complete dropouts
. tab break_w2 if wave2_all == 1

   break_w2 |      Freq.     Percent        Cum.
------------+-----------------------------------
          0 |      2,277       98.32       98.32
          1 |         39        1.68      100.00
------------+-----------------------------------
      Total |      2,316      100.00

. tab break_w3 if wave3_all == 1

   break_w3 |      Freq.     Percent        Cum.
------------+-----------------------------------
          0 |      1,668       98.70       98.70
          1 |         22        1.30      100.00
------------+-----------------------------------
      Total |      1,690      100.00

. tab break_w4a if wave4_all == 1

  break_w4a |      Freq.     Percent        Cum.
------------+-----------------------------------
          0 |      1,767       97.68       97.68
          1 |         42        2.32      100.00
------------+-----------------------------------
      Total |      1,809      100.00

. tab break_w4b if wave4_all == 1

  break_w4b |      Freq.     Percent        Cum.
------------+-----------------------------------
          0 |      1,799       99.45       99.45
          1 |         10        0.55      100.00
------------+-----------------------------------
      Total |      1,809      100.00

.  
. ttest break_w2 if wave2_all == 1, by(emo_exp_w2)

Two-sample t test with equal variances
------------------------------------------------------------------------------
   Group |     Obs        Mean    Std. err.   Std. dev.   [95% conf. interval]
---------+--------------------------------------------------------------------
 control |   1,834    .0125409    .0025992    .1113121    .0074432    .0176386
exp cond |     482     .033195    .0081683    .1793316     .017145     .049245
---------+--------------------------------------------------------------------
Combined |   2,316    .0168394    .0026742    .1286972    .0115952    .0220835
---------+--------------------------------------------------------------------
    diff |           -.0206541    .0065748               -.0335473   -.0077609
------------------------------------------------------------------------------
    diff = mean(control) - mean(exp cond)                         t =  -3.1414
H0: diff = 0                                     Degrees of freedom =     2314

    Ha: diff < 0                 Ha: diff != 0                 Ha: diff > 0
 Pr(T < t) = 0.0009         Pr(|T| > |t|) = 0.0017          Pr(T > t) = 0.9991

. ttest break_w3 if wave3_all == 1, by(emo_exp_w3)

Two-sample t test with equal variances
------------------------------------------------------------------------------
   Group |     Obs        Mean    Std. err.   Std. dev.   [95% conf. interval]
---------+--------------------------------------------------------------------
 control |     866    .0103926    .0034481    .1014716    .0036249    .0171603
exp cond |     824    .0157767    .0043437    .1246863    .0072508    .0243026
---------+--------------------------------------------------------------------
Combined |   1,690    .0130178    .0027581    .1133838    .0076081    .0184274
---------+--------------------------------------------------------------------
    diff |           -.0053841     .005518               -.0162068    .0054387
------------------------------------------------------------------------------
    diff = mean(control) - mean(exp cond)                         t =  -0.9757
H0: diff = 0                                     Degrees of freedom =     1688

    Ha: diff < 0                 Ha: diff != 0                 Ha: diff > 0
 Pr(T < t) = 0.1647         Pr(|T| > |t|) = 0.3293          Pr(T > t) = 0.8353

. ttest break_w4a if wave4_all == 1, by(emo_exp_w4a)

Two-sample t test with equal variances
------------------------------------------------------------------------------
   Group |     Obs        Mean    Std. err.   Std. dev.   [95% conf. interval]
---------+--------------------------------------------------------------------
 control |     586    .0170648    .0053547    .1296237    .0065481    .0275816
exp cond |     561    .0409982    .0083791     .198463    .0245399    .0574566
---------+--------------------------------------------------------------------
Combined |   1,147    .0287707    .0049379    .1672344    .0190823    .0384591
---------+--------------------------------------------------------------------
    diff |           -.0239334    .0098572               -.0432735   -.0045933
------------------------------------------------------------------------------
    diff = mean(control) - mean(exp cond)                         t =  -2.4280
H0: diff = 0                                     Degrees of freedom =     1145

    Ha: diff < 0                 Ha: diff != 0                 Ha: diff > 0
 Pr(T < t) = 0.0077         Pr(|T| > |t|) = 0.0153          Pr(T > t) = 0.9923

. ttest break_w4b if wave4_all == 1, by(emo_exp_w4b)

Two-sample t test with equal variances
------------------------------------------------------------------------------
   Group |     Obs        Mean    Std. err.   Std. dev.   [95% conf. interval]
---------+--------------------------------------------------------------------
 control |     315           0           0           0           0           0
exp cond |     336    .0297619    .0092843    .1701832    .0114991    .0480247
---------+--------------------------------------------------------------------
Combined |     651     .015361    .0048238    .1230784    .0058888    .0248331
---------+--------------------------------------------------------------------
    diff |           -.0297619    .0095892               -.0485915   -.0109323
------------------------------------------------------------------------------
    diff = mean(control) - mean(exp cond)                         t =  -3.1037
H0: diff = 0                                     Degrees of freedom =      649

    Ha: diff < 0                 Ha: diff != 0                 Ha: diff > 0
 Pr(T < t) = 0.0010         Pr(|T| > |t|) = 0.0020          Pr(T > t) = 0.9990

. 
. 
. 
. logit break_w2 (c.emo_exp_w2)## (c.age i.male i.university ost)   if wave2_al
> l == 1

Iteration 0:  Log likelihood = -193.18252  
Iteration 1:  Log likelihood = -191.09138  
Iteration 2:  Log likelihood = -186.03636  
Iteration 3:  Log likelihood = -186.01397  
Iteration 4:  Log likelihood = -186.01396  

Logistic regression                                     Number of obs =  2,275
                                                        LR chi2(9)    =  14.34
                                                        Prob > chi2   = 0.1108
Log likelihood = -186.01396                             Pseudo R2     = 0.0371

-------------------------------------------------------------------------------
     break_w2 | Coefficient  Std. err.      z    P>|z|     [95% conf. interval]
--------------+----------------------------------------------------------------
   emo_exp_w2 |   1.102784   1.224044     0.90   0.368    -1.296299    3.501867
          age |  -.0322365   .0416224    -0.77   0.439     -.113815     .049342
       1.male |  -.7005187   1.011053    -0.69   0.488    -2.682146    1.281109
 1.university |   .7548169   1.292267     0.58   0.559     -1.77798    3.287613
        1.ost |   .3790779   1.832746     0.21   0.836    -3.213037    3.971193
              |
 c.emo_exp_w2#|
        c.age |   .0044131    .027268     0.16   0.871    -.0490311    .0578573
              |
         male#|
 c.emo_exp_w2 |
           1  |   .5210024   .6757913     0.77   0.441    -.8035241    1.845529
              |
   university#|
 c.emo_exp_w2 |
           1  |  -.5477546   .8283074    -0.66   0.508    -2.171207    1.075698
              |
          ost#|
 c.emo_exp_w2 |
           1  |  -.4194615    1.29263    -0.32   0.746    -2.952969    2.114046
              |
        _cons |  -4.538838   1.831435    -2.48   0.013    -8.128385   -.9492917
-------------------------------------------------------------------------------

. testparm c.emo_exp_w2#c.age c.emo_exp_w2#i.male c.emo_exp_w2#i.university c.e
> mo_exp_w2#i.ost

 ( 1)  [break_w2]c.emo_exp_w2#c.age = 0
 ( 2)  [break_w2]1.male#c.emo_exp_w2 = 0
 ( 3)  [break_w2]1.university#c.emo_exp_w2 = 0
 ( 4)  [break_w2]1.ost#c.emo_exp_w2 = 0

           chi2(  4) =    1.25
         Prob > chi2 =    0.8700

. 
. logit break_w3 c.emo_exp_w3## (c.age i.male i.university ost)     if wave3_al
> l == 1

Iteration 0:  Log likelihood = -112.65039  
Iteration 1:  Log likelihood = -107.64224  
Iteration 2:  Log likelihood = -105.75846  
Iteration 3:  Log likelihood = -105.68455  
Iteration 4:  Log likelihood = -105.66691  
Iteration 5:  Log likelihood = -105.66259  
Iteration 6:  Log likelihood = -105.66175  
Iteration 7:  Log likelihood = -105.66161  
Iteration 8:  Log likelihood = -105.66158  
Iteration 9:  Log likelihood = -105.66157  

Logistic regression                                     Number of obs =  1,661
                                                        LR chi2(9)    =  13.98
                                                        Prob > chi2   = 0.1231
Log likelihood = -105.66157                             Pseudo R2     = 0.0620

-------------------------------------------------------------------------------
     break_w3 | Coefficient  Std. err.      z    P>|z|     [95% conf. interval]
--------------+----------------------------------------------------------------
   emo_exp_w3 |   1.083764   1.646646     0.66   0.510    -2.143603    4.311131
          age |  -.0702245   .0746884    -0.94   0.347    -.2166112    .0761622
       1.male |   2.029445   1.609962     1.26   0.207    -1.126022    5.184913
 1.university |   .5947714   1.860486     0.32   0.749    -3.051714    4.241257
        1.ost |  -26.28627    1507.07    -0.02   0.986    -2980.089    2927.517
              |
 c.emo_exp_w3#|
        c.age |   .0113165   .0439074     0.26   0.797    -.0747404    .0973734
              |
         male#|
 c.emo_exp_w3 |
           1  |  -1.512121   .9628741    -1.57   0.116    -3.399319     .375078
              |
   university#|
 c.emo_exp_w3 |
           1  |  -.3753659   1.107415    -0.34   0.735    -2.545859    1.795127
              |
          ost#|
 c.emo_exp_w3 |
           1  |   13.36325   753.5353     0.02   0.986    -1463.539    1490.265
              |
        _cons |  -4.070928   2.790796    -1.46   0.145    -9.540789    1.398932
-------------------------------------------------------------------------------
Note: 15 failures and 0 successes completely determined.

. testparm c.emo_exp_w3#c.age c.emo_exp_w3#i.male c.emo_exp_w3#i.university c.e
> mo_exp_w3#i.ost

 ( 1)  [break_w3]c.emo_exp_w3#c.age = 0
 ( 2)  [break_w3]1.male#c.emo_exp_w3 = 0
 ( 3)  [break_w3]1.university#c.emo_exp_w3 = 0
 ( 4)  [break_w3]1.ost#c.emo_exp_w3 = 0

           chi2(  4) =    2.58
         Prob > chi2 =    0.6307

. 
. logit break_w4a c.emo_exp_w4a## (c.age i.male i.university ost)   if wave4_al
> l == 1

Iteration 0:  Log likelihood = -141.85106  
Iteration 1:  Log likelihood =  -135.3805  
Iteration 2:  Log likelihood =  -133.2473  
Iteration 3:  Log likelihood = -133.22472  
Iteration 4:  Log likelihood = -133.22467  
Iteration 5:  Log likelihood = -133.22467  

Logistic regression                                     Number of obs =  1,123
                                                        LR chi2(9)    =  17.25
                                                        Prob > chi2   = 0.0449
Log likelihood = -133.22467                             Pseudo R2     = 0.0608

-------------------------------------------------------------------------------
    break_w4a | Coefficient  Std. err.      z    P>|z|     [95% conf. interval]
--------------+----------------------------------------------------------------
  emo_exp_w4a |   1.985295   1.615058     1.23   0.219     -1.18016    5.150751
          age |   .0228292   .0575079     0.40   0.691    -.0898842    .1355426
       1.male |  -2.843401   1.669593    -1.70   0.089    -6.115744    .4289415
 1.university |   .8985556   2.232183     0.40   0.687    -3.476443    5.273554
        1.ost |    .377813   2.239508     0.17   0.866    -4.011541    4.767167
              |
           c. |
  emo_exp_w4a#|
        c.age |  -.0382847   .0345508    -1.11   0.268    -.1060031    .0294337
              |
         male#|
c.emo_exp_w4a |
           1  |   1.482538   .9249752     1.60   0.109    -.3303798    3.295456
              |
   university#|
c.emo_exp_w4a |
           1  |  -.4358168   1.215059    -0.36   0.720    -2.817289    1.945655
              |
          ost#|
c.emo_exp_w4a |
           1  |  -.1187448    1.25422    -0.09   0.925    -2.576971    2.339481
              |
        _cons |  -5.375878   2.875969    -1.87   0.062    -11.01267    .2609173
-------------------------------------------------------------------------------

. testparm c.emo_exp_w4a#c.age c.emo_exp_w4a#i.male c.emo_exp_w4a#i.university 
> c.emo_exp_w4a#i.ost

 ( 1)  [break_w4a]c.emo_exp_w4a#c.age = 0
 ( 2)  [break_w4a]1.male#c.emo_exp_w4a = 0
 ( 3)  [break_w4a]1.university#c.emo_exp_w4a = 0
 ( 4)  [break_w4a]1.ost#c.emo_exp_w4a = 0

           chi2(  4) =    4.43
         Prob > chi2 =    0.3509

. 
. *not enough cases
. *logit break_w4b c.emo_exp_w4b## (c.age i.male i.university ost)   if wave4_a
> ll == 0
. *testparm c.emo_exp_w4b#c.age c.emo_exp_w4b#i.male c.emo_exp_w4b#i.university
>  c.emo_exp_w4b#i.ost
. 
. foreach var in emo_exp_w2 emo_exp_w3 emo_exp_w4a emo_exp_w4b{
  2.         ttest age, by(`var') 
  3.         ttest male, by(`var') 
  4.         ttest university, by(`var') 
  5.         ttest ost, by(`var') 
  6.         di `var'
  7. }

Two-sample t test with equal variances
------------------------------------------------------------------------------
   Group |     Obs        Mean    Std. err.   Std. dev.   [95% conf. interval]
---------+--------------------------------------------------------------------
 control |   1,801    39.48251     .320693    13.60963    38.85354    40.11148
exp cond |     474    40.13502    .6160889    13.41321    38.92441    41.34563
---------+--------------------------------------------------------------------
Combined |   2,275    39.61846    .2844751     13.5686     39.0606    40.17632
---------+--------------------------------------------------------------------
    diff |           -.6525114    .7004743               -2.026147    .7211244
------------------------------------------------------------------------------
    diff = mean(control) - mean(exp cond)                         t =  -0.9315
H0: diff = 0                                     Degrees of freedom =     2273

    Ha: diff < 0                 Ha: diff != 0                 Ha: diff > 0
 Pr(T < t) = 0.1758         Pr(|T| > |t|) = 0.3517          Pr(T > t) = 0.8242

Two-sample t test with equal variances
------------------------------------------------------------------------------
   Group |     Obs        Mean    Std. err.   Std. dev.   [95% conf. interval]
---------+--------------------------------------------------------------------
 control |   1,834    .4961832    .0116782    .5001218    .4732792    .5190872
exp cond |     482    .4585062    .0227194     .498793    .4138647    .5031478
---------+--------------------------------------------------------------------
Combined |   2,316     .488342    .0103891     .499972    .4679691    .5087148
---------+--------------------------------------------------------------------
    diff |             .037677    .0255848               -.0124946    .0878485
------------------------------------------------------------------------------
    diff = mean(control) - mean(exp cond)                         t =   1.4726
H0: diff = 0                                     Degrees of freedom =     2314

    Ha: diff < 0                 Ha: diff != 0                 Ha: diff > 0
 Pr(T < t) = 0.9295         Pr(|T| > |t|) = 0.1410          Pr(T > t) = 0.0705

Two-sample t test with equal variances
------------------------------------------------------------------------------
   Group |     Obs        Mean    Std. err.   Std. dev.   [95% conf. interval]
---------+--------------------------------------------------------------------
 control |   1,834    .7928026    .0094666    .4054087    .7742362     .811369
exp cond |     482    .8091286    .0179187    .3933962      .77392    .8443372
---------+--------------------------------------------------------------------
Combined |   2,316    .7962003    .0083722    .4029087    .7797826    .8126181
---------+--------------------------------------------------------------------
    diff |            -.016326    .0206247               -.0567709    .0241188
------------------------------------------------------------------------------
    diff = mean(control) - mean(exp cond)                         t =  -0.7916
H0: diff = 0                                     Degrees of freedom =     2314

    Ha: diff < 0                 Ha: diff != 0                 Ha: diff > 0
 Pr(T < t) = 0.2143         Pr(|T| > |t|) = 0.4287          Pr(T > t) = 0.7857

Two-sample t test with equal variances
------------------------------------------------------------------------------
   Group |     Obs        Mean    Std. err.   Std. dev.   [95% conf. interval]
---------+--------------------------------------------------------------------
 control |   1,834    .0839695    .0064779    .2774177    .0712646    .0966743
exp cond |     482    .0850622    .0127201    .2792641    .0600683    .1100561
---------+--------------------------------------------------------------------
Combined |   2,316    .0841969    .0057713    .2777428    .0728794    .0955143
---------+--------------------------------------------------------------------
    diff |           -.0010928    .0142194               -.0289769    .0267914
------------------------------------------------------------------------------
    diff = mean(control) - mean(exp cond)                         t =  -0.0769
H0: diff = 0                                     Degrees of freedom =     2314

    Ha: diff < 0                 Ha: diff != 0                 Ha: diff > 0
 Pr(T < t) = 0.4694         Pr(|T| > |t|) = 0.9387          Pr(T > t) = 0.5306
.

Two-sample t test with equal variances
------------------------------------------------------------------------------
   Group |     Obs        Mean    Std. err.   Std. dev.   [95% conf. interval]
---------+--------------------------------------------------------------------
 control |     852    39.08568     .455271    13.28892     38.1921    39.97927
exp cond |     809    39.99011    .4731658    13.45822    39.06133    40.91889
---------+--------------------------------------------------------------------
Combined |   1,661    39.52619    .3281841    13.37526    38.88249    40.16989
---------+--------------------------------------------------------------------
    diff |           -.9044305    .6564106               -2.191911    .3830499
------------------------------------------------------------------------------
    diff = mean(control) - mean(exp cond)                         t =  -1.3778
H0: diff = 0                                     Degrees of freedom =     1659

    Ha: diff < 0                 Ha: diff != 0                 Ha: diff > 0
 Pr(T < t) = 0.0842         Pr(|T| > |t|) = 0.1684          Pr(T > t) = 0.9158

Two-sample t test with equal variances
------------------------------------------------------------------------------
   Group |     Obs        Mean    Std. err.   Std. dev.   [95% conf. interval]
---------+--------------------------------------------------------------------
 control |     866    .4907621    .0169976    .5002035    .4574007    .5241235
exp cond |     824    .5194175    .0174158    .4999263    .4852329     .553602
---------+--------------------------------------------------------------------
Combined |   1,690    .5047337    .0121657    .5001256    .4808724    .5285951
---------+--------------------------------------------------------------------
    diff |           -.0286554    .0243361               -.0763874    .0190767
------------------------------------------------------------------------------
    diff = mean(control) - mean(exp cond)                         t =  -1.1775
H0: diff = 0                                     Degrees of freedom =     1688

    Ha: diff < 0                 Ha: diff != 0                 Ha: diff > 0
 Pr(T < t) = 0.1196         Pr(|T| > |t|) = 0.2392          Pr(T > t) = 0.8804

Two-sample t test with equal variances
------------------------------------------------------------------------------
   Group |     Obs        Mean    Std. err.   Std. dev.   [95% conf. interval]
---------+--------------------------------------------------------------------
 control |     866    .7875289    .0139083    .4092928    .7602308    .8148269
exp cond |     824    .8082524    .0137227    .3939147    .7813169     .835188
---------+--------------------------------------------------------------------
Combined |   1,690    .7976331    .0097759    .4018832     .778459    .8168073
---------+--------------------------------------------------------------------
    diff |           -.0207236    .0195571               -.0590823    .0176352
------------------------------------------------------------------------------
    diff = mean(control) - mean(exp cond)                         t =  -1.0596
H0: diff = 0                                     Degrees of freedom =     1688

    Ha: diff < 0                 Ha: diff != 0                 Ha: diff > 0
 Pr(T < t) = 0.1447         Pr(|T| > |t|) = 0.2895          Pr(T > t) = 0.8553

Two-sample t test with equal variances
------------------------------------------------------------------------------
   Group |     Obs        Mean    Std. err.   Std. dev.   [95% conf. interval]
---------+--------------------------------------------------------------------
 control |     866    .0681293    .0085672    .2521133    .0513145    .0849442
exp cond |     824    .0910194    .0100264    .2878114    .0713391    .1106997
---------+--------------------------------------------------------------------
Combined |   1,690    .0792899    .0065744    .2702707    .0663951    .0921848
---------+--------------------------------------------------------------------
    diff |           -.0228901    .0131449               -.0486722     .002892
------------------------------------------------------------------------------
    diff = mean(control) - mean(exp cond)                         t =  -1.7414
H0: diff = 0                                     Degrees of freedom =     1688

    Ha: diff < 0                 Ha: diff != 0                 Ha: diff > 0
 Pr(T < t) = 0.0409         Pr(|T| > |t|) = 0.0818          Pr(T > t) = 0.9591
.

Two-sample t test with equal variances
------------------------------------------------------------------------------
   Group |     Obs        Mean    Std. err.   Std. dev.   [95% conf. interval]
---------+--------------------------------------------------------------------
 control |     577    40.31023    .5664796    13.60731    39.19761    41.42284
exp cond |     546     39.4359    .5508222    12.87087     38.3539    40.51789
---------+--------------------------------------------------------------------
Combined |   1,123    39.88513    .3955599    13.25568    39.10901    40.66125
---------+--------------------------------------------------------------------
    diff |            .8743279    .7913436               -.6783535    2.427009
------------------------------------------------------------------------------
    diff = mean(control) - mean(exp cond)                         t =   1.1049
H0: diff = 0                                     Degrees of freedom =     1121

    Ha: diff < 0                 Ha: diff != 0                 Ha: diff > 0
 Pr(T < t) = 0.8653         Pr(|T| > |t|) = 0.2695          Pr(T > t) = 0.1347

Two-sample t test with equal variances
------------------------------------------------------------------------------
   Group |     Obs        Mean    Std. err.   Std. dev.   [95% conf. interval]
---------+--------------------------------------------------------------------
 control |     586     .496587     .020672    .5004155    .4559867    .5371874
exp cond |     561    .5062389    .0211272    .5004073    .4647406    .5477371
---------+--------------------------------------------------------------------
Combined |   1,147    .5013078    .0147699    .5002164    .4723288    .5302867
---------+--------------------------------------------------------------------
    diff |           -.0096518    .0295583               -.0676462    .0483426
------------------------------------------------------------------------------
    diff = mean(control) - mean(exp cond)                         t =  -0.3265
H0: diff = 0                                     Degrees of freedom =     1145

    Ha: diff < 0                 Ha: diff != 0                 Ha: diff > 0
 Pr(T < t) = 0.3720         Pr(|T| > |t|) = 0.7441          Pr(T > t) = 0.6280

Two-sample t test with equal variances
------------------------------------------------------------------------------
   Group |     Obs        Mean    Std. err.   Std. dev.   [95% conf. interval]
---------+--------------------------------------------------------------------
 control |     586    .8054608    .0163662    .3961838    .7733171    .8376044
exp cond |     561    .7914439    .0171683    .4066389    .7577217     .825166
---------+--------------------------------------------------------------------
Combined |   1,147    .7986051    .0118467    .4012174    .7753614    .8218487
---------+--------------------------------------------------------------------
    diff |            .0140169    .0237058               -.0324948    .0605286
------------------------------------------------------------------------------
    diff = mean(control) - mean(exp cond)                         t =   0.5913
H0: diff = 0                                     Degrees of freedom =     1145

    Ha: diff < 0                 Ha: diff != 0                 Ha: diff > 0
 Pr(T < t) = 0.7228         Pr(|T| > |t|) = 0.5544          Pr(T > t) = 0.2772

Two-sample t test with equal variances
------------------------------------------------------------------------------
   Group |     Obs        Mean    Std. err.   Std. dev.   [95% conf. interval]
---------+--------------------------------------------------------------------
 control |     586    .0733788     .010781    .2609801    .0522047     .094553
exp cond |     561    .1016043    .0127672    .3023968    .0765268    .1266817
---------+--------------------------------------------------------------------
Combined |   1,147     .087184    .0083333    .2822275    .0708337    .1035342
---------+--------------------------------------------------------------------
    diff |           -.0282254     .016657               -.0609071    .0044562
------------------------------------------------------------------------------
    diff = mean(control) - mean(exp cond)                         t =  -1.6945
H0: diff = 0                                     Degrees of freedom =     1145

    Ha: diff < 0                 Ha: diff != 0                 Ha: diff > 0
 Pr(T < t) = 0.0452         Pr(|T| > |t|) = 0.0904          Pr(T > t) = 0.9548
2

Two-sample t test with equal variances
------------------------------------------------------------------------------
   Group |     Obs        Mean    Std. err.   Std. dev.   [95% conf. interval]
---------+--------------------------------------------------------------------
 control |     309     41.9288    .8080394    14.20404    40.33883    43.51878
exp cond |     334    41.58383    .7343039    13.41989    40.13937    43.02829
---------+--------------------------------------------------------------------
Combined |     643    41.74961    .5439255    13.79256    40.68152     42.8177
---------+--------------------------------------------------------------------
    diff |            .3449703    1.089438               -1.794328    2.484268
------------------------------------------------------------------------------
    diff = mean(control) - mean(exp cond)                         t =   0.3166
H0: diff = 0                                     Degrees of freedom =      641

    Ha: diff < 0                 Ha: diff != 0                 Ha: diff > 0
 Pr(T < t) = 0.6242         Pr(|T| > |t|) = 0.7516          Pr(T > t) = 0.3758

Two-sample t test with equal variances
------------------------------------------------------------------------------
   Group |     Obs        Mean    Std. err.   Std. dev.   [95% conf. interval]
---------+--------------------------------------------------------------------
 control |     315    .5206349    .0281926    .5003689    .4651646    .5761052
exp cond |     336    .4761905    .0272869    .5001777    .4225152    .5298658
---------+--------------------------------------------------------------------
Combined |     651    .4976959    .0196114    .5003792    .4591865    .5362052
---------+--------------------------------------------------------------------
    diff |            .0444444    .0392347               -.0325978    .1214867
------------------------------------------------------------------------------
    diff = mean(control) - mean(exp cond)                         t =   1.1328
H0: diff = 0                                     Degrees of freedom =      649

    Ha: diff < 0                 Ha: diff != 0                 Ha: diff > 0
 Pr(T < t) = 0.8711         Pr(|T| > |t|) = 0.2577          Pr(T > t) = 0.1289

Two-sample t test with equal variances
------------------------------------------------------------------------------
   Group |     Obs        Mean    Std. err.   Std. dev.   [95% conf. interval]
---------+--------------------------------------------------------------------
 control |     315    .7873016    .0230934    .4098672    .7418642     .832739
exp cond |     336    .8511905    .0194449    .3564314    .8129409      .88944
---------+--------------------------------------------------------------------
Combined |     651    .8202765      .01506    .3842522    .7907043    .8498487
---------+--------------------------------------------------------------------
    diff |           -.0638889    .0300545               -.1229047   -.0048731
------------------------------------------------------------------------------
    diff = mean(control) - mean(exp cond)                         t =  -2.1258
H0: diff = 0                                     Degrees of freedom =      649

    Ha: diff < 0                 Ha: diff != 0                 Ha: diff > 0
 Pr(T < t) = 0.0169         Pr(|T| > |t|) = 0.0339          Pr(T > t) = 0.9831

Two-sample t test with equal variances
------------------------------------------------------------------------------
   Group |     Obs        Mean    Std. err.   Std. dev.   [95% conf. interval]
---------+--------------------------------------------------------------------
 control |     315    .0698413    .0143837    .2552849    .0415407    .0981418
exp cond |     336    .0654762     .013515    .2477331    .0388913    .0920611
---------+--------------------------------------------------------------------
Combined |     651    .0675883    .0098465    .2512312    .0482535    .0869232
---------+--------------------------------------------------------------------
    diff |            .0043651    .0197177               -.0343532    .0430833
------------------------------------------------------------------------------
    diff = mean(control) - mean(exp cond)                         t =   0.2214
H0: diff = 0                                     Degrees of freedom =      649

    Ha: diff < 0                 Ha: diff != 0                 Ha: diff > 0
 Pr(T < t) = 0.5876         Pr(|T| > |t|) = 0.8249          Pr(T > t) = 0.4124
.

. 
. 
. 
. tabstat age male university ost, by(emo_exp_w2) stat(mean) format(%9.2f)

Summary statistics: Mean
Group variable: emo_exp_w2 ({if(is_empty(randnumber), rand(1,5), randnumber.NAO
> K)})

  emo_exp_w2 |       age      male  univer~y       ost
-------------+----------------------------------------
control cond |     39.48      0.50      0.79      0.08
    exp cond |     40.14      0.46      0.81      0.09
-------------+----------------------------------------
       Total |     39.62      0.49      0.80      0.08
------------------------------------------------------

. tabstat age male university ost, by(emo_exp_w3) stat(mean) format(%9.2f)

Summary statistics: Mean
Group variable: emo_exp_w3 ({if(is_empty(rand1bild), rand(1,2), rand1bild.NAOK)
> })

  emo_exp_w3 |       age      male  univer~y       ost
-------------+----------------------------------------
control cond |     39.09      0.49      0.79      0.07
    exp cond |     39.99      0.52      0.81      0.09
-------------+----------------------------------------
       Total |     39.53      0.50      0.80      0.08
------------------------------------------------------

. tabstat age male university ost, by(emo_exp_w4a) stat(mean) format(%9.2f)

Summary statistics: Mean
Group variable: emo_exp_w4a 

 emo_exp_w4a |       age      male  univer~y       ost
-------------+----------------------------------------
control cond |     40.31      0.50      0.81      0.07
    exp cond |     39.44      0.51      0.79      0.10
-------------+----------------------------------------
       Total |     39.89      0.50      0.80      0.09
------------------------------------------------------

. tabstat age male university ost, by(emo_exp_w4b) stat(mean) format(%9.2f)

Summary statistics: Mean
Group variable: emo_exp_w4b 

 emo_exp_w4b |       age      male  univer~y       ost
-------------+----------------------------------------
control cond |     41.93      0.52      0.79      0.07
    exp cond |     41.58      0.48      0.85      0.07
-------------+----------------------------------------
       Total |     41.75      0.50      0.82      0.07
------------------------------------------------------

. 
. 
. 
. ** Consort diagram
. recode len_w2(0/19=0)(20/50000=1), gen(l_w2)
(1,786 differences between len_w2 and l_w2)

. recode len_w3(0/19=0)(20/50000=1), gen(l_w3)
(1,564 differences between len_w3 and l_w3)

. recode len_w4(0/19=0)(20/50000=1), gen(l_w4)
(1,146 differences between len_w4 and l_w4)

. 
. tab l_w2 mis_em_w2 if wave2==1 & emo_exp_w2==1

 RECODE of |       mis_em_w2
    len_w2 |         0          1 |     Total
-----------+----------------------+----------
         0 |       161         70 |       231 
         1 |     1,140        211 |     1,351 
-----------+----------------------+----------
     Total |     1,301        281 |     1,582 

. tab l_w2 mis_em_w2 if wave2==1 & emo_exp_w2==2

 RECODE of |       mis_em_w2
    len_w2 |         0          1 |     Total
-----------+----------------------+----------
         0 |        74         23 |        97 
         1 |       256         45 |       301 
-----------+----------------------+----------
     Total |       330         68 |       398 

. 
. 
. tab l_w3 mis_em_w3 if wave3==1 & emo_exp_w3==1

 RECODE of |       mis_em_w3
    len_w3 |         0          1 |     Total
-----------+----------------------+----------
         0 |        93         30 |       123 
         1 |       571         97 |       668 
-----------+----------------------+----------
     Total |       664        127 |       791 

. tab l_w3 mis_em_w3 if wave3==1 & emo_exp_w3==2

 RECODE of |       mis_em_w3
    len_w3 |         0          1 |     Total
-----------+----------------------+----------
         0 |        76         17 |        93 
         1 |       551        107 |       658 
-----------+----------------------+----------
     Total |       627        124 |       751 

. 
. tab lastpage  if id<1170

     Letzte |
      Seite |      Freq.     Percent        Cum.
------------+-----------------------------------
          2 |          2        0.18        0.18
          3 |         10        0.90        1.07
          4 |          3        0.27        1.34
          5 |          2        0.18        1.52
          7 |          1        0.09        1.61
          9 |          1        0.09        1.70
         10 |          1        0.09        1.79
         12 |          2        0.18        1.97
         19 |          3        0.27        2.24
         20 |         28        2.51        4.74
         22 |          2        0.18        4.92
         23 |          3        0.27        5.19
         24 |          5        0.45        5.64
         25 |          1        0.09        5.73
         27 |          3        0.27        6.00
         29 |          2        0.18        6.18
         31 |          2        0.18        6.36
         34 |          1        0.09        6.45
         35 |          1        0.09        6.54
         36 |          2        0.18        6.71
         37 |          3        0.27        6.98
         39 |      1,039       93.02      100.00
------------+-----------------------------------
      Total |      1,117      100.00

. tab emo_exp_w4a  if id<1170

 emo_exp_w4a |      Freq.     Percent        Cum.
-------------+-----------------------------------
control cond |        576       51.57       51.57
    exp cond |        541       48.43      100.00
-------------+-----------------------------------
       Total |      1,117      100.00

. tab emo_exp_w4a  if id<1170 & wave4 == 1

 emo_exp_w4a |      Freq.     Percent        Cum.
-------------+-----------------------------------
control cond |        540       51.97       51.97
    exp cond |        499       48.03      100.00
-------------+-----------------------------------
       Total |      1,039      100.00

. tab randmani  if id<1170 & wave4 == 1

{if(is_empt |
y(randmani) |
          , |
 rand(1,1), |
randmani.NA |
       OK)} |      Freq.     Percent        Cum.
------------+-----------------------------------
          1 |        523       50.34       50.34
          2 |        516       49.66      100.00
------------+-----------------------------------
      Total |      1,039      100.00

. tab randmani  if id<1170 & wave4 == 1 & emo_exp_w4a

{if(is_empt |
y(randmani) |
          , |
 rand(1,1), |
randmani.NA |
       OK)} |      Freq.     Percent        Cum.
------------+-----------------------------------
          1 |        523       50.34       50.34
          2 |        516       49.66      100.00
------------+-----------------------------------
      Total |      1,039      100.00

. tab randmani  if id<1170 & wave4 == 1 & emo_exp_w4a==1

{if(is_empt |
y(randmani) |
          , |
 rand(1,1), |
randmani.NA |
       OK)} |      Freq.     Percent        Cum.
------------+-----------------------------------
          1 |        256       47.41       47.41
          2 |        284       52.59      100.00
------------+-----------------------------------
      Total |        540      100.00

. tab randmani  if id<1170 & wave4 == 1 & emo_exp_w4a==2

{if(is_empt |
y(randmani) |
          , |
 rand(1,1), |
randmani.NA |
       OK)} |      Freq.     Percent        Cum.
------------+-----------------------------------
          1 |        267       53.51       53.51
          2 |        232       46.49      100.00
------------+-----------------------------------
      Total |        499      100.00

. 
. tab emo_exp_w4a  if id>=1170 & wave4 == 1

 emo_exp_w4a |      Freq.     Percent        Cum.
-------------+-----------------------------------
control cond |          8       40.00       40.00
    exp cond |         12       60.00      100.00
-------------+-----------------------------------
       Total |         20      100.00

. tab lastpage emo_exp_w4b if id >=1170 & wave4==1

    Letzte |      emo_exp_w4b
     Seite | control c   exp cond |     Total
-----------+----------------------+----------
        39 |         4          1 |         5 
        40 |       279        296 |       575 
-----------+----------------------+----------
     Total |       283        297 |       580 

. 
. tab l_w4 mis_em_w4a if wave4==1 & id < 1170 & emo_exp_w4a==1

 RECODE of |      mis_em_w4a
    len_w4 |         0          1 |     Total
-----------+----------------------+----------
         0 |        31         10 |        41 
         1 |       221         22 |       243 
-----------+----------------------+----------
     Total |       252         32 |       284 

. tab l_w4 mis_em_w4a if wave4==1 & id < 1170 &  emo_exp_w4a==2

 RECODE of |      mis_em_w4a
    len_w4 |         0          1 |     Total
-----------+----------------------+----------
         0 |        42          8 |        50 
         1 |       160         22 |       182 
-----------+----------------------+----------
     Total |       202         30 |       232 

. 
. 
. 
. tab mis_em_w4b if id>=1170 & wave4==1&emo_exp_w4b==1

 mis_em_w4b |      Freq.     Percent        Cum.
------------+-----------------------------------
          0 |        232       81.98       81.98
          1 |         51       18.02      100.00
------------+-----------------------------------
      Total |        283      100.00

. tab mis_em_w4b if id>=1170 & wave4==1&emo_exp_w4b==2

 mis_em_w4b |      Freq.     Percent        Cum.
------------+-----------------------------------
          0 |        261       87.88       87.88
          1 |         36       12.12      100.00
------------+-----------------------------------
      Total |        297      100.00

. 
. 
. ** Manski Bounds
. 


. 
. 
. * Note: Define variables (dv=y, iv=d01)
. 
. capture drop y d01

. 
. *Wave 1:
. gen y = emo2_3_w2
(16,333 missing values generated)

. gen d01 = emo_exp_w2 - 1
(15,687 missing values generated)

. manski_ate_mts  y d01, yL(1) yU(7)
Manski extreme bounds:
  Lower =    -3.194
  Upper =     2.806
MTS-only bounds:
  Lower =     0.306
  Upper =     2.806

. 
. *Wave 2:
. replace y = emo2_3_w3
(1,998 real changes made, 935 to missing)

. replace d01 = emo_exp_w3 - 1
(2,204 real changes made, 1,108 to missing)

. manski_ate_mts  y d01, yL(1) yU(7)
Manski extreme bounds:
  Lower =    -3.073
  Upper =     2.927
MTS-only bounds:
  Lower =    -0.127
  Upper =     2.927

. 
. *Wave 3:
. replace y = Emo2_SQ003
(1,670 real changes made, 832 to missing)

. replace d01 = emo_exp_w4a - 1
(1,673 real changes made, 896 to missing)

. manski_ate_mts  y d01, yL(1) yU(7)
Manski extreme bounds:
  Lower =    -2.803
  Upper =     3.197
MTS-only bounds:
  Lower =     0.418
  Upper =     3.197

. 
. *Wave 4:
. replace y = Emo2_SQ003
(0 real changes made)

. replace d01 = emo_exp_w4b - 1
(1,798 real changes made, 1,147 to missing)

. manski_ate_mts  y d01, yL(1) yU(7)
Manski extreme bounds:
  Lower =    -3.092
  Upper =     2.908
MTS-only bounds:
  Lower =    -0.207
  Upper =     2.908

. 
. 
. 
. manski_ci, outcome(emo2_3_w2) treat(d01) dis
Inferring maximum/minimum values from outcome data...

    Variable |        Obs        Mean    Std. dev.       Min        Max
-------------+---------------------------------------------------------
   emo2_3_w2 |      1,670    4.502395    1.734555          1          7
(16,333 missing values generated)
(172 real changes made)
(206 real changes made)
WORST CASE
logit break_w4a c.emo_exp_w4a## (c.age i.male i.university ost)   if wave4_all 
> == 1

      Source |       SS           df       MS      Number of obs   =       651
-------------+----------------------------------   F(1, 649)       =    568.88
       Model |  1983.26683         1  1983.26683   Prob > F        =    0.0000
    Residual |  2262.57956       649  3.48625511   R-squared       =    0.4671
-------------+----------------------------------   Adj R-squared   =    0.4663
       Total |  4245.84639       650  6.53207137   Root MSE        =    1.8672

------------------------------------------------------------------------------
    __000002 | Coefficient  Std. err.      t    P>|t|     [95% conf. interval]
-------------+----------------------------------------------------------------
         d01 |  -3.492659   .1464351   -23.85   0.000    -3.780202   -3.205115
       _cons |    5.88254   .1052021    55.92   0.000     5.675962    6.089117
------------------------------------------------------------------------------
WITHOUT BOUNDS
regress __000002 d01 ,  level(95)

      Source |       SS           df       MS      Number of obs   =       273
-------------+----------------------------------   F(1, 271)       =      0.07
       Model |  .197435897         1  .197435897   Prob > F        =    0.7911
    Residual |  760.930769       271  2.80786262   R-squared       =    0.0003
-------------+----------------------------------   Adj R-squared   =   -0.0034
       Total |  761.128205       272  2.79826546   Root MSE        =    1.6757

------------------------------------------------------------------------------
   emo2_3_w2 | Coefficient  Std. err.      t    P>|t|     [95% conf. interval]
-------------+----------------------------------------------------------------
         d01 |   .0538462   .2030624     0.27   0.791    -.3459343    .4536266
       _cons |   4.538462   .1401264    32.39   0.000     4.262587    4.814336
------------------------------------------------------------------------------
(16,333 missing values generated)
(206 real changes made)
(172 real changes made)
BEST CASE
regress emo2_3_w2 d01 ,  level(95)

      Source |       SS           df       MS      Number of obs   =       651
-------------+----------------------------------   F(1, 649)       =    574.71
       Model |  1948.71736         1  1948.71736   Prob > F        =    0.0000
    Residual |   2200.6129       649  3.39077488   R-squared       =    0.4696
-------------+----------------------------------   Adj R-squared   =    0.4688
       Total |  4149.33026       650  6.38358502   Root MSE        =    1.8414

------------------------------------------------------------------------------
    __000004 | Coefficient  Std. err.      t    P>|t|     [95% conf. interval]
-------------+----------------------------------------------------------------
         d01 |   3.462103   .1444159    23.97   0.000     3.178524    3.745682
       _cons |   2.606349   .1037515    25.12   0.000      2.40262    2.810078
------------------------------------------------------------------------------

e(manski_ci)[3,1]
                     d01
        coef   .05384615
95lowerbound  -3.7802025
95upperbound    3.745682

. 
. 
. 
. 
. 
. 
. 
end of do-file

. log close
      name:  <unnamed>
       log:  C:\Users\sabri\Dropbox\ReUP\Emotionen-Experiment\srm.smcl
  log type:  smcl
 closed on:  16 Apr 2026, 14:06:50
-------------------------------------------------------------------------------
