961
Views
23
CrossRef citations to date
0
Altmetric
Original Articles

Deterministic Sampling of Expensive Posteriors Using Minimum Energy Designs

, , , &
Pages 297-308 | Received 17 Dec 2017, Accepted 16 Nov 2018, Published online: 22 Mar 2019
 

Abstract

Markov chain Monte Carlo ( MCMC) methods require a large number of samples to approximate a posterior distribution, which can be costly when the likelihood or prior is expensive to evaluate. The number of samples can be reduced if we can avoid repeated samples and those that are close to each other. This is the idea behind deterministic sampling methods such as quasi-Monte Carlo (QMC). However, the existing QMC methods aim at sampling from a uniform hypercube, which can miss the high probability regions of the posterior distribution and thus the approximation can be poor. Minimum energy design (MinED) is a recently proposed deterministic sampling method, which makes use of the posterior evaluations to obtain a weighted space-filling design in the region of interest. However, the existing implementation of MinED is inefficient because it requires several global optimizations and thus numerous evaluations of the posterior. In this article, we develop an efficient algorithm that can generate MinED samples with few posterior evaluations. We also make several improvements to the MinED criterion to make it perform better in high dimensions. The advantages of MinED over MCMC and QMC are illustrated using an example of calibrating a friction drilling process.

Additional information

Funding

This research is supported by a U.S. National Science Foundation grant DMS-1712642.

Reprints and Corporate Permissions

Please note: Selecting permissions does not provide access to the full text of the article, please see our help page How do I view content?

To request a reprint or corporate permissions for this article, please click on the relevant link below:

Academic Permissions

Please note: Selecting permissions does not provide access to the full text of the article, please see our help page How do I view content?

Obtain permissions instantly via Rightslink by clicking on the button below:

If you are unable to obtain permissions via Rightslink, please complete and submit this Permissions form. For more information, please visit our Permissions help page.