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)
万分感谢


雷达卡




京公网安备 11010802022788号







