Mailing-List: contact issues-help@commons.apache.org; run by ezmlm
Precedence: bulk
Reply-To: issues@commons.apache.org
Date: Mon, 26 Nov 2012 10:03:01 +0000 (UTC)
From: "Thomas Neidhart (JIRA)" <jira@apache.org>
To: issues@commons.apache.org
Message-ID: <866190729.22393.1353924181993.JavaMail.jiratomcat@arcas>
In-Reply-To: <1543210504.21693.1353886738381.JavaMail.jiratomcat@arcas>
Subject: [jira] [Commented] (MATH-909) FDistribution NoBracketingException
 in BrentSolver
MIME-Version: 1.0
Content-Type: text/plain; charset=utf-8
Content-Transfer-Encoding: quoted-printable


    [ https://issues.apache.org/jira/browse/MATH-909?page=3Dcom.atlassian.j=
ira.plugin.system.issuetabpanels:comment-tabpanel&focusedCommentId=3D135037=
03#comment-13503703 ]=20

Thomas Neidhart commented on MATH-909:
--------------------------------------

Ah ok my fault, but adding more regression tests is definitely useful.
Will add the test case later on.
               =20
> FDistribution NoBracketingException in BrentSolver
> --------------------------------------------------
>
>                 Key: MATH-909
>                 URL: https://issues.apache.org/jira/browse/MATH-909
>             Project: Commons Math
>          Issue Type: Bug
>    Affects Versions: 3.0
>            Reporter: Patrick Meyer
>   Original Estimate: 24h
>  Remaining Estimate: 24h
>
> I get an exception when running the code below. the exception is=20
> {code}
> function values at endpoints do not have different signs, endpoints: [0, =
1.002], values: [-0.025, -=E2=88=9E]
> {code}
> The problematic code:
> {code}
> double df1 =3D 10675;
> double df2 =3D 501725;
> FDistribution fDist =3D new FDistribution(df1, df2);
> System.out.println(fDist.inverseCumulativeProbability(0.025));//NoBracket=
ingException
> {code}
> However, R returns the value 0.9733505. The R code is:
> {code}
> qf(p=3D.025, df1=3D10675, df2=3D501725)
> {code}
> I don't know enough about the FDistribution class to know the solution to=
 the exception, but I thought I would report it.

--
This message is automatically generated by JIRA.
If you think it was sent incorrectly, please contact your JIRA administrato=
rs
For more information on JIRA, see: http://www.atlassian.com/software/jira