Antithetic methods for Gibbs samplers

Chris Holmes, Ajay Jasra

Research output: Contribution to journalArticlepeer-review

3 Scopus citations

Abstract

In this article we propose a modification to the output fromMetropolis-within-Gibbs samplers that can lead to substantial reductions in the variance over standard estimates. The idea is simple: at each time step of the algorithm, introduce an extra sample into the estimate that is negatively correlated with the current sample, the rationale being that this provides a two-sample numerical approximation to a Rao-Blackwellized estimate. As the conditional sampling distribution at each step has already been constructed, the generation of the antithetic sample often requires negligible computational effort. Our method is implementable whenever one subvector of the state can be sampled from its full conditional and the corresponding distribution function may be inverted, or the full conditional has a symmetric density. We demonstrate our approach in the context of logistic regression and hierarchical Poisson models. The data and computer code used in this article are available online. © 2009 American Statistical Association, Institute of Mathematical Statistics, and Interface Foundation of North America.
Original languageEnglish (US)
JournalJournal of Computational and Graphical Statistics
Volume18
Issue number2
DOIs
StatePublished - Dec 1 2009
Externally publishedYes

Bibliographical note

Generated from Scopus record by KAUST IRTS on 2019-11-20

Fingerprint

Dive into the research topics of 'Antithetic methods for Gibbs samplers'. Together they form a unique fingerprint.

Cite this