Hi, I've started using Julia for some machine learning tasks and love it so far. Currently I need to sample from a generalized inverse Gaussian distribution , where p=1. However, it seems the Distribution is not implemented in the Distributions.jl package. Does anyone know if there exists an efficient implementation? If not, maybe someone knows an easy way how to implement a sampler myself? Of course I could try to use rejection sampling with a gamma distribution as proposal distribution, but maybe there is a better way.
This post has NOT been accepted by the mailing list yet.
I have been working with Bayesian Quantile regression model and needed to code GIG sampler in Julia. I'm not sure if you have already coded yourself but if you need it, I can send it to you so please let me know :)