NIMBLE for Bayesian Disease Mapping

Abstract

This tutorial describes the basic implementation of Bayesian hierarchical models for spatial health data using the R package nimble. To quote the nimble R description:A system for writing hierarchical statistical models largely compatible with 'BUGS' and 'JAGS', writing nimbleFunctions to operate models and do basic R-style math, and compiling both models and nimbleFunctions via custom-generated C++. 'NIMBLE' includes default methods for MCMC, particle filtering, Monte Carlo Expectation Maximization, and some other tools. The nimbleFunction system makes it easy to do things like implement new MCMC samplers from R, customize the assignment of samplers to different parts of a model from R, and compile the new samplers automatically via C++ alongside the samplers 'NIMBLE' provides.Examples of the use of the package for a small range of Bayesian Disease Mapping (BDM) models is explored and focus on different approaches to model fitting and analysis are discussed. Examples of publicly available small area health data is used throughout.