At Sun, 23 Oct 2011 15:54:11 +0100, Yan Zhou wrote: > > The problem is that gsl fails for not too small degree of freedoms. > I think this is reported and fixed before > (http://lists.gnu.org/archive/html/bug-gsl/2008-04/msg00036.html). The > GSL version I was using is 1.15, compiled with Clang 3.0 (Apple) and > the `make check` test passed successfully when installing GSL. It > appears that the same problem remains. > > And this time I don't see a simple fix for this problem as in the > original post. Maybe this is a a flow in the algorithm itself rather > the implementation and I would suggest taking a look into R's pgamma > implementation, which seems to be more robust.
Thanks for the bug report, I've added a test case for it. -- Brian Gough _______________________________________________ Bug-gsl mailing list [email protected] https://lists.gnu.org/mailman/listinfo/bug-gsl
