I can measure a distribution's kurtosis and skewness using the moments package, and simulate a distribution using rnorm(), but I don't know how to apply the two together to simulate a non-normal distribution. Is there a package that can help? Most of the things I looked up haven't been very helpful.
For example, I would like to generate a distribution that has a leptokurtic distribution (say a kurtosis of 5), and a skewness of -.3 . Is there a way to apply that specific kurtosis and skewness to a distribution that is assumed normal?