Evaluating marginal likelihood with thermodynamic integration method and comparison with several other numerical methods
Abstract
Abstract Evaluating marginal likelihood is the most critical and computationally expensive task, when conducting Bayesian model averaging to quantify parametric and model uncertainties. The evaluation is commonly done by using Laplace approximations to evaluate semianalytical expressions of the marginal likelihood or by using Monte Carlo (MC) methods to evaluate arithmetic or harmonic mean of a joint likelihood function. This study introduces a new MC method, i.e., thermodynamic integration, which has not been attempted in environmental modeling. Instead of using samples only from prior parameter space (as in arithmetic mean evaluation) or posterior parameter space (as in harmonic mean evaluation), the thermodynamic integration method uses samples generated gradually from the prior to posterior parameter space. This is done through a path sampling that conducts Markov chain Monte Carlo simulation with different power coefficient values applied to the joint likelihood function. The thermodynamic integration method is evaluated using three analytical functions by comparing the method with two variants of the Laplace approximation method and three MC methods, including the nested sampling method that is recently introduced into environmental modeling. The thermodynamic integration method outperforms the other methods in terms of their accuracy, convergence, and consistency. The thermodynamic integration method is also appliedmore »
- Authors:
-
- Hefei Univ. of Technology, Hefei (China); Florida State Univ., Tallahassee, FL (United States)
- Florida State Univ., Tallahassee, FL (United States)
- Nanjing Univ., Nanjing (China)
- Oak Ridge National Lab. (ORNL), Oak Ridge, TN (United States)
- Hefei Univ. of Technology, Hefei (China)
- Publication Date:
- Research Org.:
- Oak Ridge National Laboratory (ORNL), Oak Ridge, TN (United States)
- Sponsoring Org.:
- USDOE
- OSTI Identifier:
- 1326496
- Alternate Identifier(s):
- OSTI ID: 1402380
- Grant/Contract Number:
- AC05-00OR22725; DEāSC0008272
- Resource Type:
- Accepted Manuscript
- Journal Name:
- Water Resources Research
- Additional Journal Information:
- Journal Volume: 52; Journal Issue: 2; Journal ID: ISSN 0043-1397
- Publisher:
- American Geophysical Union (AGU)
- Country of Publication:
- United States
- Language:
- English
- Subject:
- 97 MATHEMATICS AND COMPUTING; Bayesian model averaging; arithmetic mean; harmonic mean; nested sampling; model uncertainty; Markov chain Monte Carlo
Citation Formats
Liu, Peigui, Elshall, Ahmed S., Ye, Ming, Beerli, Peter, Zeng, Xiankui, Lu, Dan, and Tao, Yuezan. Evaluating marginal likelihood with thermodynamic integration method and comparison with several other numerical methods. United States: N. p., 2016.
Web. doi:10.1002/2014WR016718.
Liu, Peigui, Elshall, Ahmed S., Ye, Ming, Beerli, Peter, Zeng, Xiankui, Lu, Dan, & Tao, Yuezan. Evaluating marginal likelihood with thermodynamic integration method and comparison with several other numerical methods. United States. https://doi.org/10.1002/2014WR016718
Liu, Peigui, Elshall, Ahmed S., Ye, Ming, Beerli, Peter, Zeng, Xiankui, Lu, Dan, and Tao, Yuezan. Fri .
"Evaluating marginal likelihood with thermodynamic integration method and comparison with several other numerical methods". United States. https://doi.org/10.1002/2014WR016718. https://www.osti.gov/servlets/purl/1326496.
@article{osti_1326496,
title = {Evaluating marginal likelihood with thermodynamic integration method and comparison with several other numerical methods},
author = {Liu, Peigui and Elshall, Ahmed S. and Ye, Ming and Beerli, Peter and Zeng, Xiankui and Lu, Dan and Tao, Yuezan},
abstractNote = {Abstract Evaluating marginal likelihood is the most critical and computationally expensive task, when conducting Bayesian model averaging to quantify parametric and model uncertainties. The evaluation is commonly done by using Laplace approximations to evaluate semianalytical expressions of the marginal likelihood or by using Monte Carlo (MC) methods to evaluate arithmetic or harmonic mean of a joint likelihood function. This study introduces a new MC method, i.e., thermodynamic integration, which has not been attempted in environmental modeling. Instead of using samples only from prior parameter space (as in arithmetic mean evaluation) or posterior parameter space (as in harmonic mean evaluation), the thermodynamic integration method uses samples generated gradually from the prior to posterior parameter space. This is done through a path sampling that conducts Markov chain Monte Carlo simulation with different power coefficient values applied to the joint likelihood function. The thermodynamic integration method is evaluated using three analytical functions by comparing the method with two variants of the Laplace approximation method and three MC methods, including the nested sampling method that is recently introduced into environmental modeling. The thermodynamic integration method outperforms the other methods in terms of their accuracy, convergence, and consistency. The thermodynamic integration method is also applied to a synthetic case of groundwater modeling with four alternative models. The application shows that model probabilities obtained using the thermodynamic integration method improves predictive performance of Bayesian model averaging. The thermodynamic integration method is mathematically rigorous, and its MC implementation is computationally general for a wide range of environmental problems.},
doi = {10.1002/2014WR016718},
journal = {Water Resources Research},
number = 2,
volume = 52,
place = {United States},
year = {Fri Feb 05 00:00:00 EST 2016},
month = {Fri Feb 05 00:00:00 EST 2016}
}
Web of Science
Works referenced in this record:
Bayes Factors
journal, June 1995
- Kass, Robert E.; Raftery, Adrian E.
- Journal of the American Statistical Association, Vol. 90, Issue 430
Maximum likelihood Bayesian averaging of spatial variability models in unsaturated fractured tuff: MAXIMUM LIKELIHOOD BAYESIAN MODEL AVERAGING
journal, May 2004
- Ye, Ming; Neuman, Shlomo P.; Meyer, Philip D.
- Water Resources Research, Vol. 40, Issue 5
Comment on āInverse groundwater modeling for hydraulic conductivity estimation using Bayesian model averaging and variance windowā by Frank T.-C. Tsai and Xiaobao Li: COMMENTARY
journal, February 2010
- Ye, Ming; Lu, Dan; Neuman, Shlomo P.
- Water Resources Research, Vol. 46, Issue 2
Nested Sampling
conference, January 2004
- Skilling, John
- BAYESIAN INFERENCE AND MAXIMUM ENTROPY METHODS IN SCIENCE AND ENGINEERING: 24th International Workshop on Bayesian Inference and Maximum Entropy Methods in Science and Engineering, AIP Conference Proceedings
Marginal likelihood estimation via power posteriors
journal, July 2008
- Friel, N.; Pettitt, A. N.
- Journal of the Royal Statistical Society: Series B (Statistical Methodology), Vol. 70, Issue 3
Simulating normalizing constants: from importance sampling to bridge sampling to path sampling
journal, May 1998
- Gelman, Andrew; Meng, Xiao-Li
- Statistical Science, Vol. 13, Issue 2
Nested sampling algorithm for subsurface flow model selection, uncertainty quantification, and nonlinear calibration: NESTED SAMPLING FOR SUBSURFACE FLOW MODELS
journal, December 2013
- Elsheikh, A. H.; Wheeler, M. F.; Hoteit, I.
- Water Resources Research, Vol. 49, Issue 12
Treatment of input uncertainty in hydrologic modeling: Doing hydrology backward with Markov chain Monte Carlo simulation: FORCING DATA ERROR USING MCMC SAMPLING
journal, December 2008
- Vrugt, Jasper A.; ter Braak, Cajo J. F.; Clark, Martyn P.
- Water Resources Research, Vol. 44, Issue 12
Conceptual model uncertainty in groundwater modeling: Combining generalized likelihood uncertainty estimation and Bayesian model averaging: EVALUATING CONCEPTUAL MODEL
journal, December 2008
- Rojas, Rodrigo; Feyen, Luc; Dassargues, Alain
- Water Resources Research, Vol. 44, Issue 12
Hydrological model selection: A Bayesian alternative: HYDROLOGICAL MODEL SELECTION
journal, October 2005
- Marshall, Lucy; Nott, David; Sharma, Ashish
- Water Resources Research, Vol. 41, Issue 10
Review of surrogate modeling in water resources: REVIEW
journal, July 2012
- Razavi, Saman; Tolson, Bryan A.; Burn, Donald H.
- Water Resources Research, Vol. 48, Issue 7
Assessment of parametric uncertainty for groundwater reactive transport modeling
journal, May 2014
- Shi, Xiaoqing; Ye, Ming; Curtis, Gary P.
- Water Resources Research, Vol. 50, Issue 5
A multimodel data assimilation framework via the ensemble Kalman filter
journal, May 2014
- Xue, Liang; Zhang, Dongxiao
- Water Resources Research, Vol. 50, Issue 5
Estimation of Aquifer Parameters Under Transient and Steady State Conditions: 1. Maximum Likelihood Method Incorporating Prior Information
journal, February 1986
- Carrera, Jesus; Neuman, Shlomo P.
- Water Resources Research, Vol. 22, Issue 2
Assessing five evolving microbial enzyme models against field measurements from a semiarid savannah-What are the mechanisms of soil respiration pulses?
journal, September 2014
- Zhang, Xia; Niu, Guo-Yue; Elshall, Ahmed S.
- Geophysical Research Letters, Vol. 41, Issue 18
Bayesian analysis of data-worth considering model and parameter uncertainties
journal, February 2012
- Neuman, Shlomo P.; Xue, Liang; Ye, Ming
- Advances in Water Resources, Vol. 36
A Model-Averaging Method for Assessing Groundwater Conceptual Model Uncertainty
journal, August 2010
- Ye, Ming; Pohlmann, Karl F.; Chapman, Jenny B.
- Ground Water, Vol. 48, Issue 5
Analysis of regression confidence intervals and Bayesian credible intervals for uncertainty quantification: REGRESSION CONFIDENCE AND BAYESIAN CREDIBLE INTERVALS
journal, September 2012
- Lu, Dan; Ye, Ming; Hill, Mary C.
- Water Resources Research, Vol. 48, Issue 9
Computing the Bayes Factor from a Markov Chain Monte Carlo Simulation of the Posterior Distribution
journal, September 2012
- Weinberg, Martin D.
- Bayesian Analysis, Vol. 7, Issue 3
Hierarchical Bayesian model averaging for hydrostratigraphic modeling: Uncertainty segregation and comparative evaluation: HIERARCHICAL BAYESIAN MODEL AVERAGING
journal, September 2013
- Tsai, Frank T. -C.; Elshall, Ahmed S.
- Water Resources Research, Vol. 49, Issue 9
An adaptive sparse-grid high-order stochastic collocation method for Bayesian inference in groundwater reactive transport modeling: Sparse-Grid Method for Bayesian Inference
journal, October 2013
- Zhang, Guannan; Lu, Dan; Ye, Ming
- Water Resources Research, Vol. 49, Issue 10
Efficient posterior exploration of a high-dimensional groundwater model from two-stage Markov chain Monte Carlo simulation and polynomial chaos expansion: Speeding up MCMC Simulation of a Groundwater Model
journal, May 2013
- Laloy, Eric; Rogiers, Bart; Vrugt, Jasper A.
- Water Resources Research, Vol. 49, Issue 5
Estimating the evidence - a review
journal, January 2012
- Friel, Nial; Wyse, Jason
- Statistica Neerlandica, Vol. 66, Issue 3
Model selection on solid ground: Rigorous comparison of nine ways to evaluate Bayesian model evidence
journal, December 2014
- Schƶniger, Anneli; Wƶhling, Thomas; Samaniego, Luis
- Water Resources Research, Vol. 50, Issue 12
Postaudit evaluation of conceptual model uncertainty for a glacial aquifer groundwater flow and contaminant transport model
journal, January 2010
- Lemke, Lawrence D.; Cypher, Joseph A.
- Hydrogeology Journal, Vol. 18, Issue 4
High-dimensional posterior exploration of hydrologic models using multiple-try DREAM (ZS) and high-performance computing : EFFICIENT MCMC FOR HIGH-DIMENSIONAL PROBLEMS
journal, January 2012
- Laloy, Eric; Vrugt, Jasper A.
- Water Resources Research, Vol. 48, Issue 1
Unified Framework to Evaluate Panmixia and Migration Direction Among Multiple Sampling Locations
journal, February 2010
- Beerli, Peter; Palczewski, Michal
- Genetics, Vol. 185, Issue 1
Computing Bayes Factors Using Thermodynamic Integration
journal, April 2006
- Lartillot, Nicolas; Philippe, HervƩ
- Systematic Biology, Vol. 55, Issue 2
Hydrologic data assimilation using particle Markov chain Monte Carlo simulation: Theory, concepts and applications
journal, January 2013
- Vrugt, Jasper A.; ter Braak, Cajo J. F.; Diks, Cees G. H.
- Advances in Water Resources, Vol. 51
Maximum likelihood Bayesian averaging of uncertain model predictions
journal, November 2003
- Neuman, S. P.
- Stochastic Environmental Research and Risk Assessment (SERRA), Vol. 17, Issue 5
On model selection criteria in multimodel analysis: ON MODEL SELECTION CRITERIA IN MULTIMODEL ANALYSIS
journal, March 2008
- Ye, Ming; Meyer, Philip D.; Neuman, Shlomo P.
- Water Resources Research, Vol. 44, Issue 3
A Bayesian approach for inverse modeling, data assimilation, and conditional simulation of spatial random fields
journal, October 2010
- Rubin, Yoram; Chen, Xingyuan; Murakami, Haruko
- Water Resources Research, Vol. 46, Issue 10
Constructive epistemic modeling of groundwater flow with geological structure and boundary condition uncertainty under the Bayesian paradigm
journal, September 2014
- Elshall, Ahmed S.; Tsai, Frank T. -C.
- Journal of Hydrology, Vol. 517
Marginal Likelihood From the MetropolisāHastings Output
journal, March 2001
- Chib, Siddhartha; Jeliazkov, Ivan
- Journal of the American Statistical Association, Vol. 96, Issue 453
Improving Marginal Likelihood Estimation for Bayesian Phylogenetic Model Selection
journal, December 2010
- Xie, Wangang; Lewis, Paul O.; Fan, Yu
- Systematic Biology, Vol. 60, Issue 2
Markov Chain Monte Carlo Methods for Computing Bayes Factors: A Comparative Review
journal, September 2001
- Han, Cong; Carlin, Bradley P.
- Journal of the American Statistical Association, Vol. 96, Issue 455
Multimodel Bayesian analysis of data-worth applied to unsaturated fractured tuffs
journal, January 2012
- Lu, Dan; Ye, Ming; Neuman, Shlomo P.
- Advances in Water Resources, Vol. 35
Model complexity control for hydrologic prediction: MODEL COMPLEXITY CONTROL
journal, August 2008
- Schoups, G.; van de Giesen, N. C.; Savenije, H. H. G.
- Water Resources Research, Vol. 44, Issue 12
Using MCMC chain outputs to efficiently estimate Bayes factors
journal, October 2011
- Morey, Richard D.; Rouder, Jeffrey N.; Pratte, Michael S.
- Journal of Mathematical Psychology, Vol. 55, Issue 5
Dependence of Bayesian Model Selection Criteria and Fisher Information Matrix on Sample Size
journal, October 2011
- Lu, Dan; Ye, Ming; Neuman, Shlomo P.
- Mathematical Geosciences, Vol. 43, Issue 8
Effects of error covariance structure on estimation of model averaging weights and predictive performance: EFFECTS OF ERROR COVARIANCE STRUCTURE ON MODEL AVERAGING
journal, September 2013
- Lu, Dan; Ye, Ming; Meyer, Philip D.
- Water Resources Research, Vol. 49, Issue 9
David Draper and E. I. George, and a rejoinder by the authors
journal, November 1999
- Volinsky, Chris T.; Raftery, Adrian E.; Madigan, David
- Statistical Science, Vol. 14, Issue 4
Markov Chain Sampling Methods for Dirichlet Process Mixture Models
journal, June 2000
- Neal, Radford M.
- Journal of Computational and Graphical Statistics, Vol. 9, Issue 2
Works referencing / citing this record:
Hydrogeological Model Selection Among Complex Spatial Priors
journal, August 2019
- Brunetti, C.; Bianchi, M.; Pirot, G.
- Water Resources Research, Vol. 55, Issue 8
Relative model score: a scoring rule for evaluating ensemble simulations with application to microbial soil respiration modeling
journal, August 2018
- Elshall, Ahmed S.; Ye, Ming; Pei, Yongzhen
- Stochastic Environmental Research and Risk Assessment, Vol. 32, Issue 10