Access

You are not currently logged in.

Access your personal account or get JSTOR access through your library or other institution:

login

Log in to your personal account or through your institution.

If You Use a Screen Reader

This content is available through Read Online (Free) program, which relies on page scans. Since scans are not currently available to screen readers, please contact JSTOR User Support for access. We'll provide a PDF copy for your screen reader.

Inference from Iterative Simulation Using Multiple Sequences

Andrew Gelman and Donald B. Rubin
Statistical Science
Vol. 7, No. 4 (Nov., 1992), pp. 457-472
Stable URL: http://www.jstor.org/stable/2246093
Page Count: 16
  • Read Online (Free)
  • Subscribe ($19.50)
  • Cite this Item
Since scans are not currently available to screen readers, please contact JSTOR User Support for access. We'll provide a PDF copy for your screen reader.
Inference from Iterative Simulation Using Multiple Sequences
Preview not available

Abstract

The Gibbs sampler, the algorithm of Metropolis and similar iterative simulation methods are potentially very helpful for summarizing multivariate distributions. Used naively, however, iterative simulation can give misleading answers. Our methods are simple and generally applicable to the output of any iterative simulation; they are designed for researchers primarily interested in the science underlying the data and models they are analyzing, rather than for researchers interested in the probability theory underlying the iterative simulations themselves. Our recommended strategy is to use several independent sequences, with starting points sampled from an overdispersed distribution. At each step of the iterative simulation, we obtain, for each univariate estimand of interest, a distributional estimate and an estimate of how much sharper the distributional estimate might become if the simulations were continued indefinitely. Because our focus is on applied inference for Bayesian posterior distributions in real problems, which often tend toward normality after transformations and marginalization, we derive our results as normal-theory approximations to exact Bayesian inference, conditional on the observed simulations. The methods are illustrated on a random-effects mixture model applied to experimental measurements of reaction times of normal and schizophrenic patients.

Page Thumbnails

  • Thumbnail: Page 
457
    457
  • Thumbnail: Page 
458
    458
  • Thumbnail: Page 
459
    459
  • Thumbnail: Page 
460
    460
  • Thumbnail: Page 
461
    461
  • Thumbnail: Page 
462
    462
  • Thumbnail: Page 
463
    463
  • Thumbnail: Page 
464
    464
  • Thumbnail: Page 
465
    465
  • Thumbnail: Page 
466
    466
  • Thumbnail: Page 
467
    467
  • Thumbnail: Page 
468
    468
  • Thumbnail: Page 
469
    469
  • Thumbnail: Page 
470
    470
  • Thumbnail: Page 
471
    471
  • Thumbnail: Page 
472
    472