{"doi":"10.1007/s13253-021-00485-9","title":"Greater Than the Sum of its Parts: Computationally Flexible Bayesian Hierarchical Modeling","abstract":"<jats:title>Abstract</jats:title><jats:p>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<jats:italic>a posteriori</jats:italic>(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.</jats:p>","journal":"Journal of Agricultural, Biological and Environmental Statistics","year":2022,"id":629717,"datarank":0.29188652235829704,"base_score":1.9459101490553132,"endowment":1.9459101490553132,"self_citation_contribution":0.29188652235829704,"citation_network_contribution":0.0,"self_endowment_contribution":0.29188652235829704,"citer_contribution":0.0,"corpus_percentile":null,"corpus_rank":null,"citation_count":6,"citer_count":0,"citers_with_citation_signal":0,"citers_with_endowment":0,"datacite_reuse_total":0,"is_dataset":false,"is_dataset_confidence":null,"is_data_producer":false,"deposit_databanks":null,"is_oa":false,"file_count":0,"downloads":0,"has_version_chain":false,"published_date":null,"fair_score":null,"fair_percentile":null,"algorithm_id":"datarank_citation_only_1hop_v6","ranking_scope":"data_only","authors":[{"id":1631002,"name":"Brian M. Brost","orcid":null,"position":1,"is_corresponding":false},{"id":1404373,"name":"Mevin B. Hooten","orcid":"0000-0002-1614-723X","position":2,"is_corresponding":false},{"id":1631000,"name":"Devin S. Johnson","orcid":"0000-0002-0068-7098","position":0,"is_corresponding":false}],"reference_count":0,"raw_metadata":{"has_enrichment":true,"resolved":true,"title":"Greater Than the Sum of its Parts: Computationally Flexible Bayesian Hierarchical Modeling","abstract":"<jats:title>Abstract</jats:title><jats:p>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<jats:italic>a posteriori</jats:italic>(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.</jats:p>","is_dataset_classified":null,"base_score":1.9459101490553132,"endowment":1.9459101490553132,"datacite_reuse_total":0,"file_count":0,"downloads":0,"views":0,"has_version_chain":false,"is_dataset":false,"is_oa":false,"pmid":"19767382","pmcid":null,"openalex_id":"https://openalex.org/W3093969086","authors":[],"funders":[{"funder_name":"National Science Foundation","grant_id":"1927177","title":"Collaborative research:  RoL:  Using reaction norms to link genomic and phenotypic variation with regional-scale population responses to environmental change"}],"total_grants":1,"fwci":0.6593,"citation_percentile":0.65929284,"influential_citations":0,"citation_trend":[{"year":2022,"count":2},{"year":2025,"count":1},{"year":2026,"count":3}],"oa_status":"hybrid","license":"cc-by","oa_locations":[{"url":"https://link.springer.com/content/pdf/10.1007/s13253-021-00485-9.pdf","host_type":"journal"},{"url":"https://link.springer.com/content/pdf/10.1007/s13253-021-00485-9.pdf","host_type":"publisher"},{"url":"https://link.springer.com/article/10.1007/s13253-021-00485-9/fulltext.html","host_type":"publisher"},{"url":"https://doi.org/10.1007/s13253-021-00485-9","host_type":"journal"},{"url":"http://link.springer.com/10.1007/s13253-021-00485-9","host_type":"repository"},{"url":"https://dx.doi.org/10.48550/arxiv.2010.12568","host_type":""},{"url":"http://arxiv.org/abs/2010.12568","host_type":""},{"url":"https://zbmath.org/7570509","host_type":""}],"fields_of_study":["Statistical Methods and Bayesian Inference","Economic and Environmental Valuation","Bayesian Methods and Mixture Models","01 natural sciences","0101 mathematics"],"mesh_terms":[],"keywords":["Estimator","Hierarchical database model","Bayesian probability","Bayesian inference","Mathematics","Markov chain Monte Carlo","Random effects model","Algorithm","Mixed model","Computer science","Graphical model","Multivariate statistics","Statistics","Data mining","FOS: Computer and information sciences","Statistics - Computation","meta-analysis","Methodology (stat.ME)","integrated population model","Applications of statistics to environmental and related topics","multistage estimation","approximation","linear mixed model","Bayesian hierarchical model","Statistics - Methodology","Computation (stat.CO)","Bayesian MAP estimation"],"sdg_mappings":[],"linked_datasets":[],"clinical_trials":[],"software_tools":[],"database_accessions":[],"source":"live","citation_network_status":"fetched"},"created_at":"2026-08-05T19:30:12.702148Z","pmid":null,"pmcid":null,"fwci":null,"citation_percentile":null,"influential_citations":0,"oa_status":null,"license":null,"views":0,"total_file_size_bytes":0,"version_count":0,"fair_f":null,"fair_a":null,"fair_i":null,"fair_r":null,"fair_zscore":null,"fair_rationale":null,"fair_model":null,"fair_agent_version":null,"fair_fulltext_source":null,"fair_has_llm":null,"fair_computed_at":null,"clinical_trials":[],"software_tools":[],"db_accessions":[],"linked_datasets":[],"topics":[]}