On Monte Carlo computation of posterior expectations with uncertainty

Wei Wei*, Wenxin Jiang

*Corresponding author for this work

Research output: Contribution to journalArticlepeer-review

Abstract

In this paper, we study computation of the range of posterior expectations that arise from robust Bayesian statistics. We compute supremum and infimum of the posterior expectations, when allowing uncertainty for the choice of the likelihood function, or uncertainty for the choice of the prior distribution. In the standard approach of sensitivity analysis, posterior statistics is computed a multiple number of times for each choice of the uncertainty scenarios, which might involve heavy computation due to running Monte Carlo sampling many times. Our paper proposes a more efficient computational method that only requires one Monte Carlo sample for all possible choices of the uncertainty scenarios. The proposed computational method involves three steps (with the mnemonic PSI): (Prior step.) Introduce an auxiliary hyperprior distribution on a parameter λ that indexes the uncertainty. (Sampling step.) For any parameter of interest h, we derive a sample of (h, λ) from the joint posterior distribution given the observed data, using a Monte Carlo method. 1 (Inference step.) Based on this posterior sample of (h, λ), we estimate the range of posterior expectations {infλ E(h|λ, data), supλ E(h|λ, data)} (and similarly the range for any posterior quantile).

Original languageEnglish (US)
Pages (from-to)2038-2049
Number of pages12
JournalJournal of Statistical Computation and Simulation
Volume87
Issue number10
DOIs
StatePublished - Jul 3 2017

Keywords

  • Bayes
  • Expectation
  • Monte Carlo
  • posterior
  • prior
  • uncertainty

ASJC Scopus subject areas

  • Statistics and Probability
  • Modeling and Simulation
  • Statistics, Probability and Uncertainty
  • Applied Mathematics

Fingerprint Dive into the research topics of 'On Monte Carlo computation of posterior expectations with uncertainty'. Together they form a unique fingerprint.

Cite this