“…This intermediate correlation is then imposed on uncorrelated, standard-normal variates before they are transformed to non-normality via the third-order polynomial method, and that yields the full multivariate extension of this algorithm. The Fleishman (1978) method and its multivariate extension are by far the most widely used data-generating algorithms in the quantitative behavioural sciences and they have seen widespread use outside these fields as well, as can be seen in Demirtas (2016), Demirtas and Hedeker (2016), Demirtas, Ahmadian, Atis, Can, and Ercan (2016) and Demirtas and Vardar-Acar (2017). The third-order polynomial method is implemented in latent variable software packages such as EQS (Bentler, 1995), the lavaan R package (Rosseel, 2012), the simsem R package (Pornprasertmanit, Miller, & Schoemann, 2012) and the BinOrdNonNor R package (Demirtas, Wang, & Allozi, 2017).…”