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.
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?