Mailing-List: contact issues-help@commons.apache.org; run by ezmlm
Precedence: bulk
Reply-To: issues@commons.apache.org
Date: Wed, 13 Apr 2016 14:47:26 +0000 (UTC)
From: "Gilles (JIRA)" <jira@apache.org>
To: issues@commons.apache.org
Message-ID: <JIRA.12958453.1460550435000.219318.1460558846219@Atlassian.JIRA>
In-Reply-To: <JIRA.12958453.1460550435000@Atlassian.JIRA>
References: <JIRA.12958453.1460550435000@Atlassian.JIRA>
 <JIRA.12958453.1460550435702@arcas>
Subject: [jira] [Commented] (MATH-1356) HypergeometricDistribution
 probability give NaN result
MIME-Version: 1.0
Content-Type: text/plain; charset=utf-8
Content-Transfer-Encoding: 7bit


    [ https://issues.apache.org/jira/browse/MATH-1356?page=com.atlassian.jira.plugin.system.issuetabpanels:comment-tabpanel&focusedCommentId=15239357#comment-15239357 ] 

Gilles commented on MATH-1356:
------------------------------

Bug(s) indeed:
* failing to check precondition(s) that lead to computing "0 * infinity", and/or
* failing to give the correct result in partcular cases (that throw off the implementation)

The NaN is created in {{SaddleExpansion.logBinomialProbability}} and {{SaddleExpansion.getDeviancePart}}.
Both perform a computation like the above.

Case
{noformat}
double probIT = hgd.probability(0);
{noformat}
also produces NaN.


> HypergeometricDistribution probability give NaN result
> ------------------------------------------------------
>
>                 Key: MATH-1356
>                 URL: https://issues.apache.org/jira/browse/MATH-1356
>             Project: Commons Math
>          Issue Type: Bug
>    Affects Versions: 3.6.1
>         Environment: Windows
>            Reporter: Lacroix Thomas
>            Priority: Minor
>
> Hi,
> Unless I am mistaken the HypergeometricDistribution probability method returns NaN for the following cases :
> HypergeometricDistribution hgd = new HypergeometricDistribution(11,11,1);
> double probIT = hgd.probability(1);
> HypergeometricDistribution hgd = new HypergeometricDistribution(11,11,11);
> double probIT = hgd.probability(11);
> I think it should return 1.0
> Thanks,
> Thomas


--
This message was sent by Atlassian JIRA
(v6.3.4#6332)