There’s another approx from an Intel guy:
In Stan you can always take the Bayes approach, sample from an Inverse Gaussian and then use this for the Poisson distribution as input. Then test your data and if appropriate optimize, otherwise move on. Just my 2 cents.