Stata users have access to two easy-to-use implementations of Bayesian inference: Stata's native bayesmh command and StataStan, which calls the general Bayesian engine, Stan. We compare these implementations on two important models for education research: the Rasch model and the hierarchical Rasch model. StataStan fits a more general range of models than can be fit by bayesmh and uses a superior sampling algorithm, that is, Hamiltonian Monte Carlo using the no-U-turn sampler. Furthermore, StataStan can run in parallel on multiple CPU cores, regardless of the flavor of Stata. Given these advantages and given that Stan is open source and can be run directly from Stata do-files, we recommend that Stata users interested in Bayesian methods consider using StataStan.
HoffmanM. D., and GelmanA.2014. The no-U-turn sampler: Adaptively setting path lengths in Hamiltonian Monte Carlo. Journal of Machine Learning Research15: 1593–1623.
5.
KucukelbirA., RanganathR., GelmanA., and BleiD. M.2015. Automatic variational inference in Stan. ArXiv Working Paper No. arXiv:1506.03431. https://arxiv.org/abs/1506.03431.
6.
NealR. M.2011. MCMC using Hamiltonian dynamics. In Handbook of Markov Chain Monte Carlo, ed. BrooksS., GelmanA., JonesG. L., and MengX.-L., 113–162. Boca Raton, FL: Chapman & Hall/CRC.
7.
NocedalJ., and WrightS. J.2006. Numerical Optimization. 2nd ed. Berlin: Springer.
