`neg_binomial_2_lcdf()` "grad_2F1: k (internal counter) exceeded 100000 iterations"

Thanks for the quick reply. Here is an example that gives the error:

y=1 
mu=0.000519101
phi=1.70972

I get the error about half the times I fit my data (less often with the simulated data above). As soon as I get the error, that chain crashes.

Interestingly, I don’t get it with BRMS, which can also fit this model. I guess the priors used there avoid it.