ENH: add predict_prob to poisson

Adds a method to predict the probability of a certain count to poisson results. Compare to Stata's `predict varname, n(#)` and predprob from R's pscl. Needs tests but works.