Greater than the sum of its parts: Computationally flexible Bayesian hierarchical modeling

Journal of Agricultural, Biological and Environmental Statistics
By: , and 

Links

Abstract

We propose a multistage method for making inference at all levels of a Bayesian hierarchical model (BHM) using natural data partitions to increase efficiency by allowing computations to take place in parallel form using software that is most appropriate for each data partition. The full hierarchical model is then approximated by the product of independent normal distributions for the data component of the model. In the second stage, the Bayesian maximum a posteriori (MAP) estimator is found by maximizing the approximated posterior density with respect to the parameters. If the parameters of the model can be represented as normally distributed random effects, then the second-stage optimization is equivalent to fitting a multivariate normal linear mixed model. We consider a third stage that updates the estimates of distinct parameters for each data partition based on the results of the second stage. The method is demonstrated with two ecological data sets and models, a generalized linear mixed effects model (GLMM) and an integrated population model (IPM). The multistage results were compared to estimates from models fit in single stages to the entire data set. In both cases, multistage results were very similar to a full MCMC analysis. Supplementary materials accompanying this paper appear online.

Publication type Article
Publication Subtype Journal Article
Title Greater than the sum of its parts: Computationally flexible Bayesian hierarchical modeling
Series title Journal of Agricultural, Biological and Environmental Statistics
DOI 10.1007/s13253-021-00485-9
Edition 382
Volume 27
Year Published 2022
Language English
Publisher Springer
Contributing office(s) Coop Res Unit Seattle
Description 19 p.
First page 400
Google Analytic Metrics Metrics page
Additional publication details