This study develops a multivariate, nonnormal density function that can accurately and separately account for skewness, kurtosis, heteroskedasticity, and the correlation among the random variables of interest. The statistical attributes of the underlying random variables and correlation processes are examined. The potential applications of this modeling tool are discussed and exemplified by analyzing and simulating Corn Belt corn, soybean, and wheat yields. While corn and soybean yields are found to be skewed and kurtotic and exhibit different variances through time, wheat yields appear normal but also heteroskedastic. A strong correlation is detected between corn and soybean yields. Copyright 1997, Oxford University Press.