ABSTRACT
The family of polynomial-normal distributions includes a number of widely used distributions, such as the Gram–Charlier and Edgeworth distributions. In this note, we present three simple algorithms, (i) CDF Inversion, (ii) Acceptance–Rejection, (iii) and Ratio–of–Uniforms, for simulating variates from a polynomial-normal distribution. Details on the efficiency of the Acceptance–Rejection and the Ratio–of–Uniforms algorithms and a comparison across the various implementations are provided.
Notes
1 The MatLab code used is available upon request. An effort was made to optimize the code, in particular in terms of vectorization.