maybe you should assert sign(d) == sign(p) (or whatever the correct julia equivalent is?)
(sorry if this is dumb - i was just looking at the code out of curiousity and don't know the maths) andrew On Friday, 24 July 2015 09:29:27 UTC-3, [email protected] wrote: > > Here is a generalized gamma distribution. > > https://github.com/jlapeyre/GenGammaDist.jl > > https://en.wikipedia.org/wiki/Generalized_gamma_distribution > > Only rand() and mean() are implemented, and again no tests. This depends > on Distributions. The authors of Distributions may want to include > something like this. The parameters may be chosen to get the Amoroso > distribution (inverse generalized gamma distribution) mentioned in the > Wikipedia article. > > John > > > >
