I’m getting NaN from a multinomial distribution, where I expected 1. Is this a bug or am I doing something wrong? using Distributions dist = Multinomial(2, [0.999,0.001]) pdf(dist,[2,0]) ## 0.998, ok dist = Multinomial(2, [1.0,0.0]) pdf(dist,[2,0]) #NaN, not ok The equivalent limit of the binom…

NaN from multinomial distribution at p=1

improbable22 October 21, 2017, 10:13pm 5

OK, I made issue 675.

I also had a brief look at the code but am way out of my depth!