We apply a new simulation method that solves the multidimensional probability integrals that arise in maximum likelihood estimation of a broad class of limited dependent variable models. The simulation method has four key features: the simulated choice probabilities are unbiased; they are a...