@BillT wrote:
to calculate the cumulative probabilities, I added a line before the 'output' statement of the skellam datastep: 'cp+p;'
This should be fine. If in doubt, just increase the value of macro variable c. As soon as the largest computed cp values are equal to 1 (within the accuracy limits), you see that the omitted lower tail of the distribution (i.e., the sum of the p values for x<floor(m-&c*s)) was negligible.
... View more