[
https://issues.apache.org/jira/browse/MATH-753?page=com.atlassian.jira.plugin.system.issuetabpanels:comment-tabpanel&focusedCommentId=13216391#comment-13216391
]
Gilles commented on MATH-753:
-----------------------------
bq. wait until 3.1 to include this fix,
+1
bq. include this fix right now, but keep the ticket open, in order to
thoroughly investigate accurracy issues in 3.1.
-1
Let's freeze the "trunk" as proposed by Luc. This is an enhancement that can
wait a few days, the more so that you suspect that there could be side-effects.
> Better implementation for the gamma distribution density function
> -----------------------------------------------------------------
>
> Key: MATH-753
> URL: https://issues.apache.org/jira/browse/MATH-753
> Project: Commons Math
> Issue Type: Improvement
> Affects Versions: 2.2
> Reporter: Francesco Strino
> Priority: Minor
> Labels: improvement, stability
> Fix For: 2.2
>
>
> The way the density of the gamma distribution function is estimated can be
> improved.
> It's much more stable to calculate first the log of the density and then
> exponentiate, otherwise the function returns NaN for high values of the
> parameters alpha and beta.
> It would be sufficient to change the public double density(double x) function
> at line 204 in the file
> org.apache.commons.math.distribution.GammaDistributionImpl as follows:
> return Math.exp(Math.log( x )*(alpha-1) - Math.log(beta)*alpha - x/beta -
> Gamma.logGamma(alpha));
> In order to improve performance, log(beta) and Gamma.logGamma(alpha) could
> also be precomputed and stored during initialization.
--
This message is automatically generated by JIRA.
If you think it was sent incorrectly, please contact your JIRA administrators:
https://issues.apache.org/jira/secure/ContactAdministrators!default.jspa
For more information on JIRA, see: http://www.atlassian.com/software/jira
