ResearchTrend.AI
  • Papers
  • Communities
  • Events
  • Blog
  • Pricing
Papers
Communities
Social Events
Terms and Conditions
Pricing
Parameter LabParameter LabTwitterGitHubLinkedInBlueskyYoutube

© 2025 ResearchTrend.AI, All rights reserved.

  1. Home
  2. Papers
  3. 0804.0671
84
92

A theoretical comparison of the data augmentation, marginal augmentation and PX-DA algorithms

4 April 2008
J. Hobert
Dobrin Marchev
ArXiv (abs)PDFHTML
Abstract

The data augmentation (DA) algorithm is a widely used Markov chain Monte Carlo (MCMC) algorithm that is based on a Markov transition density of the form p(x∣x′)=∫YfX∣Y(x∣y)fY∣X(y∣x′)dyp(x|x')=\int_{\mathsf{Y}}f_{X|Y}(x|y)f_{Y|X}(y|x') dyp(x∣x′)=∫Y​fX∣Y​(x∣y)fY∣X​(y∣x′)dy, where fX∣Yf_{X|Y}fX∣Y​ and fY∣Xf_{Y|X}fY∣X​ are conditional densities. The PX-DA and marginal augmentation algorithms of Liu and Wu [J. Amer. Statist. Assoc. 94 (1999) 1264--1274] and Meng and van Dyk [Biometrika 86 (1999) 301--320] are alternatives to DA that often converge much faster and are only slightly more computationally demanding. The transition densities of these alternative algorithms can be written in the form pR(x∣x′)=∫Y∫YfX∣Y(x∣y′)R(y,dy′)fY∣X(y∣x′)dyp_R(x|x')=\int_{\mathsf{Y}}\int _{\mathsf{Y}}f_{X|Y}(x|y')R(y,dy')f_{Y|X}(y|x') dypR​(x∣x′)=∫Y​∫Y​fX∣Y​(x∣y′)R(y,dy′)fY∣X​(y∣x′)dy, where RRR is a Markov transition function on Y\mathsf{Y}Y. We prove that when RRR satisfies certain conditions, the MCMC algorithm driven by pRp_RpR​ is at least as good as that driven by ppp in terms of performance in the central limit theorem and in the operator norm sense. These results are brought to bear on a theoretical comparison of the DA, PX-DA and marginal augmentation algorithms. Our focus is on situations where the group structure exploited by Liu and Wu is available. We show that the PX-DA algorithm based on Haar measure is at least as good as any PX-DA algorithm constructed using a proper prior on the group.

View on arXiv
Comments on this paper