Going-Greek/2810Model.do at main · HigherEdData/Going-Greek · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
*** do 2800Model

 ** Model 1

    mi estimate, ni(20) dots post: regress $Y $X1
    estimates store Model_1

    forvalues i = 1(1)20 {
        qui mi xeq `i': regress $Y $X1
        local R`i' = e(r2)
    }
    local R_m1 = (`R1'+`R2'+`R3'+`R4'+`R5'+`R6'+`R7'+`R8'+`R9'+`R10'+`R11'+`R12'+`R13'+`R14'+`R15'+`R16'+`R17'+`R18'+`R19'+`R20')/20

 ** Model 2

    mi estimate, ni(20) dots post: regress $Y $X1 i.$X2 $Select
    estimates store Model_2

    forvalues i = 1(1)20 {
        qui mi xeq `i': regress $Y $X1 i.$X2 $Select
        local R`i' = e(r2)
    }
    local R_m2 = (`R1'+`R2'+`R3'+`R4'+`R5'+`R6'+`R7'+`R8'+`R9'+`R10'+`R11'+`R12'+`R13'+`R14'+`R15'+`R16'+`R17'+`R18'+`R19'+`R20')/20

 ** Model 3

    mi estimate, ni(20) dots post saving(Model_3, replace): regress $Y c.$X1#i1.$X2 c.$X1#i2.$X2 c.$X1#i3.$X2 i.$X2 $Select
    estimates store Model_3
    qui mi estimate (diff: _b[c.$X1#i1.$X2]-_b[c.$X1#i2.$X2]) using Model_3, nocoef
    mi testtransform diff
    qui mi estimate (diff: _b[c.$X1#i1.$X2]-_b[c.$X1#i3.$X2]) using Model_3, nocoef
    mi testtransform diff
    qui mi estimate (diff: _b[c.$X1#i2.$X2]-_b[c.$X1#i3.$X2]) using Model_3, nocoef
    mi testtransform diff

    forvalues i = 1(1)20 {
        qui mi xeq `i': regress $Y c.$X1#i1.$X2 c.$X1#i2.$X2 c.$X1#i3.$X2 i.$X2 $Select
        local R`i' = e(r2)
    }
    local R_m3 = (`R1'+`R2'+`R3'+`R4'+`R5'+`R6'+`R7'+`R8'+`R9'+`R10'+`R11'+`R12'+`R13'+`R14'+`R15'+`R16'+`R17'+`R18'+`R19'+`R20')/20

 ** Model 4

    mi estimate, ni(20) dots post saving(Model_4, replace): regress $Y c.$X1#i1.$X2 c.$X1#i2.$X2 c.$X1#i3.$X2 i.$X2 $Select $iControls
    estimates store Model_4
    qui mi estimate (diff: _b[c.$X1#i1.$X2]-_b[c.$X1#i2.$X2]) using Model_4, nocoef
    mi testtransform diff
    qui mi estimate (diff: _b[c.$X1#i1.$X2]-_b[c.$X1#i3.$X2]) using Model_4, nocoef
    mi testtransform diff
    qui mi estimate (diff: _b[c.$X1#i2.$X2]-_b[c.$X1#i3.$X2]) using Model_4, nocoef
    mi testtransform diff

    forvalues i = 1(1)20 {
        qui mi xeq `i': regress $Y c.$X1#i1.$X2 c.$X1#i2.$X2 c.$X1#i3.$X2 i.$X2 $Select $iControls
        local R`i' = e(r2)
    }
    local R_m4 = (`R1'+`R2'+`R3'+`R4'+`R5'+`R6'+`R7'+`R8'+`R9'+`R10'+`R11'+`R12'+`R13'+`R14'+`R15'+`R16'+`R17'+`R18'+`R19'+`R20')/20


    // estimates table Model_1 Model_2 Model_3, b(%9.3f) star(.01 .05 .10) stats(N F_mi p_mi) title($title)
    estimates table Model_1 Model_2 Model_3 Model_4, b(%9.3f) star(.001 .01 .05) stats(N F_mi p_mi) title($title)
    di _col(6) %7s "R2" _col(14) "|" _col(20) %4.3f `R_m1' _col(35) %4.3f `R_m2' _col(50) %4.3f `R_m3' _col(65) %4.3f `R_m4'