Thanks for looking into it! Basically I’m just looking to draw from the distribution, having specified the relevant parameters. For example:
mydist = Normal(0, 1)
rand(mydist)
Ideally I would love to be able to just plug in SkewNormal for Normal. But I guess all I really need is the cdf and the already-existing ability to draw from Uniform(0,1), correct?