In a working paper, Luo (2011) analyzes *generating non-normal data*, using *Fleishman power method*. But this approach just can be used for specific skewness and kurtosis. I would like to simulate non-normal data for larger values of skewness and kurtosis (for example, skewness = 3.5 and kurtosis = -4.0). My **questions** are: 1) how can I do this; 2) are there any `R` package that can generate such data; 3) how can I generate non-normal data for residuals of the regression model?

**References**

Luo, H. (2011). *Generation of non-normal data – A study of Fleishman’s power method.* [Working paper]. Uppsala University. Retrieved from www.diva-portal.org/smash/get/diva2:407995/FULLTEXT01.pdf