[Numpy-discussion] multinomial question

Alan G Isaac aisaac@american....
Thu Dec 6 07:19:40 CST 2007


> Alan G Isaac wrote:
>> I would think that 
>>     multinomial(1,prob,size=ntrials).sum(axis=0)
>> would be equivalent to 
>>     multinomial(ntrials,prob)
>> but the first gives a surprising result.  (See below.)
>> Explanation? 


On Wed, 05 Dec 2007, Robert Kern apparently wrote:
> Pretty much anyone who derives their binomial distribution algorithm from 
> http://www.unc.edu/~gfish/fcmc.html is also wrong. 
> SVN now has a bound such that CDF(bound) is within 1e-16 (or so) of 1.0. 


Thanks!
Alan





More information about the Numpy-discussion mailing list