You can run a Monte Carlo estimate of a Poisson
parameter.
It is the Bayes estimate of the mean q of a Poisson distribution assuming squared error loss
and a heirarchical
prior of gamma distributions. See the Poisson gamma heirarchy in Lehmann and Casella, Theory of
Point Estimation,
for details.