Rand() doesn't work for NormalInverseGaussian

#1

rand() doesn’t work for NormalInverseGaussian.

x = rand(dnig, 100)

Error: MethodError: no method matching iterate(::NormalInverseGaussian{Float64})

Would somebody help to generate random values?
#2

Please quote your code, and include an MWE.

#3

You should probably file an issue on GitHub, because it looks like that functionality is missing for that distribution. Or you can try to implement it and file a PR.