[Top][All Lists]
[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
Re: [Help-gsl] Re: question about dirichlet lnpdf function (potential bu
From: |
Brian Gough |
Subject: |
Re: [Help-gsl] Re: question about dirichlet lnpdf function (potential bug) |
Date: |
Sun, 19 Jul 2009 11:04:35 +0100 |
User-agent: |
Wanderlust/2.14.0 (Africa) Emacs/22.2 Mule/5.0 (SAKAKI) |
At Thu, 16 Jul 2009 16:21:54 -0500,
Rodney Sparapani wrote:
> Now, I'm not so sure. I see in the GSL docs that you allow the
> alphas to be positive including zero. However, GCSR says that
> the alphas>0. What does alpha1=0 mean?
You're right, the function should only be valid for alpha>0. I'll
change the documentation in the manual to be "alpha > 0" rather than
"alpha >= 0". For alpha = 0 a nan is the correct result.
--
Brian