In generating non-normal data, H. Luo gives the coding to generate non-normal data using Fleishman Power method. www.diva-portal.org/smash/get/diva2:407995/FULLTEXT01.pdf. But this coding just can be used in specific skewness and kurtosis. I would like to simulate the non-normal data for more large values of skewness and kurtosis (for example skewness = 3.5 and kurtosis = -4.0).

my question is how can i do this or are there any package in R that can generate the data? Thanks...

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