Please use this identifier to cite or link to this item:
Ardia, David
Hoogerheide, Lennart
van Dijk, Herman K.
Year of Publication: 
Series/Report no.: 
Tinbergen Institute Discussion Paper No. 09-017/4
This discussion paper led to a publication in 'Computational Statistics & Data Analysis' 56(11), pp. 3398-1414.Important choices for efficient and accurate evaluation of marginal likelihoods by means of Monte Carlo simulation methods are studied for the case of highly non-elliptical posterior distributions. We focus on the situation where one makes use of importance sampling or the independence chain Metropolis-Hastings algorithm for posterior analysis. A comparative analysis is presented of possible advantages and limitations of different simulation techniques; of possible choices of candidate distributions and choices of target or warped target distributions; and finally of numerical standard errors. The importance of a robust and flexible estimation strategy is demonstrated where the complete posterior distribution is explored. In this respect, the adaptive mixture of Student-t distributions of Hoogerheide et al.(2007) works particularly well. Given an appropriately yet quickly tuned candidate, straightforward importance sampling provides the most efficient estimator of the marginal likelihood in the cases investigated in this paper, which include a non-linear regression model of Ritter and Tanner (1992) and a conditional normal distribution of Gelman and Meng (1991). A poor choice of candidate density may lead to a huge loss of efficiency where the numerical standard error may be highly unreliable.
marginal likelihood
Bayes factor
importance sampling
Markov chain Monte Carlo
bridge sampling
adaptive mixture of Student-t distributions
Document Type: 
Working Paper

Files in This Item:
574.02 kB

Items in EconStor are protected by copyright, with all rights reserved, unless otherwise indicated.