/* r_nrml generates Normal(0,1) distributed random variable */ program define r_nrml version 2.1 gen %_1 = invnorm(uniform()) end