/* r_ltnrml generates "long-tail" Normal random variable */ program define r_ltnrml version 2.1 quietly { gen %_1 = invnorm(uniform()) replace u = 1.25*u if abs(u)>1.5 & uniform()>.5 } end