#6 Potential bug in Gamma.quantile

[Roby Joehanes]

Originally posted by Roland Ewald:
https://sourceforge.net/p/jdistlib/discussion/general/thread/1fccccac/

Sorry to post this here, but the bug tracker does not seem to allow others to create new tickets.

There might be a problem in Gamma.quantile. It seems to work fine for most values, but not always.

For example:

Gamma.quantile(0.42, 617.0, 2.0, true, false)

yields 1206.088240779075, which differs from qgamma in R (3.0.1):

> qgamma(p=0.42,shape=617.0,scale=2.0,lower.tail=TRUE,log.p=FALSE)
> 1223.334

Best regards,
Roland

[Roby Joehanes]

This turns out a bug in Normal.cumulative function. Fixed in version 0.0.8.