Saddlepoint approximation for sums of negative binomials in Stan

Hi, I just wanted to say that this post is great!
In my case (Using integration to fit the difference of gamma distributed random variable) , I actually found the saddle point analytically (well wolfram alpha did), and then it obviously goes super fast :), and it was a much more stable than the real approximate solution by integration.

3 Likes