[Top][All Lists]
[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
Re: [Bug-gsl] possible bug in F distribution computation function gsl_ra
From: |
Brian Gough |
Subject: |
Re: [Bug-gsl] possible bug in F distribution computation function gsl_ran_fdist_pdf |
Date: |
Wed, 24 Feb 2010 19:54:54 +0000 |
User-agent: |
Wanderlust/2.14.0 (Africa) Emacs/22.1 Mule/5.0 (SAKAKI) |
At Tue, 5 Jan 2010 14:55:19 -0500,
Gunjan Verma wrote:
> When I run gsl_ran_fdist_pdf(.012, 8, 249) I get infinity.
> But when I run
>
> gsl_ran_fdist_pdf(.012, 8, 248) I get .000074, which is correct (checked
> against MATLAB and R)
>
> one d.f. change makes this difference - but I don't see why it should from
> the formula for F pdf.
>
> i am running a Dell Latitude E6500 machine with ubuntu 9 linux , 64 bit and
> gcc 4.4.1 . I built the gsl library from source for my system.
I have just committed a fix for this problem.
http://git.savannah.gnu.org/cgit/gsl.git/commit/?id=87c4d916076889709158ff189db1ec0d80211e78
Thanks for the bug report.
--
Brian Gough
GNU Scientific Library -
http://www.gnu.org/software/gsl/
[Prev in Thread] |
Current Thread |
[Next in Thread] |
- Re: [Bug-gsl] possible bug in F distribution computation function gsl_ran_fdist_pdf,
Brian Gough <=