Integrated non-factorised variational Bayes and Markov chain Monte Carlo methods for estimation of spatial negative binomial models, as presented in the following paper:
Prateek Bansal, Rico Krueger, Daniel J. Graham: Fast Bayesian Estimation of Spatial Count Data Models.