[Top][All Lists]
[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
Re: [Bug-gsl] bug in gsl_ran_multinomial_lnpdf function
From: |
Brian Gough |
Subject: |
Re: [Bug-gsl] bug in gsl_ran_multinomial_lnpdf function |
Date: |
Thu, 20 Sep 2007 16:04:31 +0100 |
User-agent: |
Wanderlust/2.14.0 (Africa) Emacs/22.1 Mule/5.0 (SAKAKI) |
At Tue, 18 Sep 2007 16:48:54 -0400,
mingxi wrote:
> I encountered a bug in gsl_ran_multinomial_lnpdf function.
> The case this function does not consider is the degenerated cases.
> Suppose I have input size_t K=3, p[]={0.5,0.5,0}; n[]={2,3,0}. This
> function will not return correct results.
>
> In this case, we can just ignore the last class probability, which is
> p[]={0.5,0.5}, n[]={2,3}.
Thanks for the bug report. I see the problem, and will correct that
in the next release.
--
Brian Gough