Hi,

I am trying to improvise an example from help of
Code:
simulate
command as shown below

Code:
     program define lnsim, rclass
        version 17
        syntax [, obs(integer 1) mu(real 0) sigma(real 1) ]
        drop _all
        set obs `obs'
        tempvar z
        gen `z' = exp(rnormal(`mu',`sigma'))
        summarize `z'
        return scalar mean = r(mean)
        return scalar Var  = r(Var)
    end
The goal is to have the program return percentiles and not the mean and variance.

Any ideas?.

Thanks in advance