Further properties of random orthogonal matrix simulation

Random orthogonal matrix (ROM) simulation is a very fast procedure for generating multivariate random samples that always have exactly the same mean, covariance and Mardia multivariate skewness and kurtosis. This paper investigates how the properties of parametric, data‐specific and deterministic ROM simulations are influenced by the choice of orthogonal matrix. Specifically, we consider how cyclic and general permutation matrices alter their time‐series properties, and how three classes of rotation matrices – upper Hessenberg, Cayley, and exponential – influence both the unconditional moments of the marginal distributions and the behaviour of skewness when samples are concatenated. We also perform an experiment which demonstrates that parametric ROM simulation can be hundreds of times faster than equivalent Monte Carlo simulation.