Just in case anyone is following this, I was looking for a function to generate Monte Carlo Poisson rates given an observation k.
Given the "exact" CI expression I used mu = 1/2*CHIINV(k,2*(n+k)), where u is Uniform(0,1). This gives about 2k df for small values of u, and 2k+2 df for large u. It...