Seol mar théacs é seo: Scalable Bayesian methods for the analysis of neuroimaging data