clear
set obs 20000
gen random = invnorm(uniform())
sum random , detail
return list // basic stats are stored in r( )

gen dummy = 1 if ( random > r(mean) + r(sd) ) | ( random < r(mean) - r(sd) )
replace dummy = 0 if dummy == .

sum dummy
Anonymous | 01.23.08 - 6:19 pm | #

I LOVE YOU ANONYMOUS.
Belle Lettre | Homepage | 01.23.08 - 6:40 pm | #