[Top][All Lists]
[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
Re: [Bug-gsl] Error in gsl_cdf_fdist_Pinv
From: |
Brian Gough |
Subject: |
Re: [Bug-gsl] Error in gsl_cdf_fdist_Pinv |
Date: |
Fri, 10 Nov 2006 13:42:05 +0000 |
User-agent: |
Wanderlust/2.14.0 (Africa) Emacs/21.3 Mule/5.0 (SAKAKI) |
At Fri, 10 Nov 2006 12:08:31 +1100,
Daniel Falster wrote:
>
> The function double gsl_cdf_fdist_Pinv (double P, double nu1, double
> nu2) is in GSL ver1.8 but hasn't been included in earlier versions. It
> returns the inverse of the cumulative distribution function for the
> F-distribution. The arguments nu1 and nu2 are the degrees of freedom
> (1,..., Inf). When nu1 =1, and nu2 greater than ~ 250 the functions
> fails. I have verified this for both windows and unix versions.
> For example,
> for (int y = 1; y<500; y+=20)
> cout<<y<<"
> \t"<<gsl_cdf_fdist_Pinv(0.95,1.0,y)<<endl;
> fails for y > 240 on windows XP.
Hello,
Thanks for your email. Could you send the whole example program (with
#include's etc) and its complete output showing how it fails, so I can
reproduce the problem correctly.
--
best regards,
Brian Gough
(GSL Maintainer)
Network Theory Ltd,
Publishing the GSL Manual - http://www.network-theory.co.uk/gsl/manual/