Calculation of multivariate normal probabilities by simulation, with applications to maximum simulated likelihood estimation

We discuss methods for calculating multivariate normal probabilities by simulation and two new Stata programs for this purpose: mdraws for deriving draws from the standard uniform density using either Halton or pseudorandom sequences, and an egen function, mvnp(), for calculating the probabilities themselves. Several illustrations show how the programs may be used for maximum simulated likelihood estimation.