I’m fitting Negative Binomial in Bayesian frame work. I want to check the sensitivity of priors in the dispersion parameters. I used two different priors and compared it using the Bayes Factor, and I chose bridgesampling. Unfortunately, I got the following error for many times. I read here about the related error, and they recommend increasing the number of iteration. I tried the number of iteration with 30,000, but I got the same.

FYI – My number of observations is around 6,000.

```
Error in jj[2,]: Subscript out of bounds
In addition: Warning messages:
1: Infinite value in iterative scheme, returning NA.
Try rerunning with more samples.
2: logml could not be estimated within maxiter, rerunning with adjusted starting value.
Estimate might be more variable than usual.
```

Therefore:

1- Would you please recommend something to solve it? Or

2 - Do we have other packages or techniques use to compare the model with different priors?

Regards