Search All of the Math Forum:
Views expressed in these public forums are not endorsed by
NCTM or The Math Forum.



Re: Simulating iid Poisson variates
Posted:
Jun 19, 1996 3:56 PM


In article <m0uWReL000Iv0C@one.reed.edu>, jones@one.reed.edu wrote: [I wrote]: > >I need to simulate iid Poisson variates for a Monte Carlo > >calculation. I can do it by brute force; are there any > >more elegant or convenient ways? > I'm not sure what you mean by "brute force" here, perhaps you are already aware of the socalled inversetransform method, ie. generate u ~ U(0,1), let p_i = Pr(X = i), then let the poisson rv X be 0 if u < p_0 1 if p_0 <= u < p_0 + p_1 2 if p_0 + p_1 <= u < p_0 + p_1 + p_2 etc. if E(X) is small, this won't require too many comparisons. if E(X) is large then it would be quicker to start the search near the mean. depending on the context, it might be useful to use a table of the values of the CDF, or to make use of the fact that
p_k = (lambda/k)*p_{k1} to compute the probabilites on the fly.
Albyn, good to hear from you!
Yes, obvious schemes like inverting the cumulative distribution, and acceptancerejection methods, were what I had in mind as "brute force." I was, rather, looking to see if there were more elegant or faster schemes available.
Thanks to all who have replied either on the net or by Email. Further ideas will be gratefully received.
Bill



