> 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 _______________________________________________ Numpy-discussion mailing list Numpy-discussion@scipy.org http://projects.scipy.org/mailman/listinfo/numpy-discussion