On the specification of prior distributions for variance components in disease mapping models

In this paper, we consider the problem of specifying priors for the variance components in the Bayesian analysis of the Besag-York-Molli´e model, a model that is popular among epidemiologists for disease mapping. The model encompasses two sets of random effects: one spatially structured to model spatial autocorrelation and the other spatially unstructured to describe residual heterogeneity. In this model, prior specification for variance components is an important problem because these priors maintain their influence on the posterior distributions of relative risks when mapping rare diseases. We propose using generalised inverse Gaussian priors, a broad class of distributions that includes many distributions commonly used as priors in this context, such as inverse gammas. We discuss the prior parameter choice with the aim of balancing the prior weight of the two sets of random effects on total variation and controlling the amount of shrinkage. The suggested prior specification strategy is compared to popular alternatives using a simulation exercise and applications to real data sets.