A method for generating multivariate nonnormal distributions with specified intercorrelations and marginal means, variances, skews, and kurtoses is proposed. As an example, the method is applied to the generation of simulated scores on three psychological tests administered to a single group of individuals.