generate N random numbers from a skew normal distribution using numpy

Question

I need a function in python to return N random numbers from a skew normal distribution. The skew needs to be taken as a parameter.

e.g. my current use is

x = numpy.random.randn(1000)

and the ideal function would be e.g.

x = randn_skew(1000, skew=0.7)

Solution needs to conform with: python version 2.7, numpy v.1.9

A similar answer is here: skew normal distribution in scipy However this generates a PDF not the random numbers.

user_6288518 · Answer

from scipy.stats import skewnorm
a=10
data= skewnorm.rvs(a, size=1000)

Here, a is a parameter which you can refer to: https://docs.scipy.org/doc/scipy/reference/generated/scipy.stats.skewnorm.html