Generating random samples in Stata is very straightforward if the distribution drawn from is uniform or normal. With any other distribution, an inverse method can be used; but even in this case, the user is limited to the built-in functions. For any other distribution functions, their inverse must be derived analytically or numerical methods must be used if analytical derivation of the inverse function is tedious or impossible. In this article, I introduce a command that generates a random sample from any user-specified distribution function using numeric methods that make this command very generic.
AvramidisA. N., and WilsonJ. R.1994. A flexible method for estimating inverse distribution functions in simulation experiments. ORSA Journal on Computing6: 342–355.
2.
GentleJ. E.2003. Random Number Generation and Monte Carlo Methods. 2nd ed. New York: Springer.
3.
LindD. A., MarchalW. G., and WathenS. A.2008. Basic Statistics for Business and Economics. 6th ed. New York: McGraw–Hill.
4.
UlrichG., and WatsonL. T.1987. A method for computer generation of variates from arbitrary continuous distributions. SIAM Journal on Scientific Computing8: 185–197.
