Random orthogonal matrix simulation with exact means, covariances, and multivariate skewness

We develop a simulation algorithm that generates multivariate samples with exact means, covariances, and multivariate skewness. If required for financial applications, absence of arbitrage can be ensured. Potential applications include the simulation of risk factors for the risk management of financial institutions. We use the Kollo measure of multivariate skewness, which is more informative for these applications than the Mardia skewness previously used in this context.