capture program drop xx
program xx, rclass
args N alpha d0
quietly {
        drop _all
        set obs `N'
        g treat = (_n > 0.5 * `N')
        g attend = runiform() > (.6 - `d0' * treat)
        reg attend treat
        return scalar p2 = abs(_b[treat] / _se[treat]) < invnormal(1 - `alpha' / 2)
}
end


qui simulate p2 = r(p2) , reps(1000) : xx 200 .05 0.1
sum
di 1 - r(mean)

qui simulate p2 = r(p2) , reps(1000) : xx 200 .10 0.1
sum
di 1 - r(mean)

qui simulate p2 = r(p2) , reps(1000) : xx 200 .01 0.1
sum
di 1 - r(mean)

qui simulate p2 = r(p2) , reps(1000) : xx 400 .05 0.1
sum
di 1 - r(mean)

万分感谢