Importance Tempering

dc.creatorGramacy, Robert B.
dc.creatorSamworth, Richard J.
dc.creatorKing, Ruth
dc.date2007-07-28
dc.date2008-11-03
dc.date.accessioned2026-07-07T10:14:21Z
dc.date.available2026-07-07T10:14:21Z
dc.descriptionSimulated tempering (ST) is an established Markov chain Monte Carlo (MCMC) method for sampling from a multimodal density $π(θ)$. Typically, ST involves introducing an auxiliary variable $k$ taking values in a finite subset of $[0,1]$ and indexing a set of tempered distributions, say $π_k(θ) \propto π(θ)^k$. In this case, small values of $k$ encourage better mixing, but samples from $π$ are only obtained when the joint chain for $(θ,k)$ reaches $k=1$. However, the entire chain can be used to estimate expectations under $π$ of functions of interest, provided that importance sampling (IS) weights are calculated. Unfortunately this method, which we call importance tempering (IT), can disappoint. This is partly because the most immediately obvious implementation is naïve and can lead to high variance estimators. We derive a new optimal method for combining multiple IS estimators and prove that the resulting estimator has a highly desirable property related to the notion of effective sample size. We briefly report on the success of the optimal combination in two modelling scenarios requiring reversible-jump MCMC, where the naïve approach fails.
dc.description16 pages, 2 tables, significantly shortened from version 4 in response to referee comments, to appear in Statistics and Computing
dc.identifierhttps://arxiv.org/abs/0707.4242
dc.identifierhttp://arxiv.org/abs/0707.4242
dc.identifier.urihttp://salesiana.dossiersoluciones.com/handle/123456789/172852
dc.subjectComputation
dc.subjectApplications
dc.titleImportance Tempering
dc.typetext

Files

Collections