Antithetic Coupling for Perfect Sampling

SUMMARY This paper reports some initial investigations of the use of antithetic variates in perfect sampling. A simple random walk example is presented to illustrate the key ingredients of antithetic coupling for perfect sampling as well as its potential benefit. A key step in implementing antithetic coupling is to generate random variates that are negatively associated, a stronger condition than negative correlation as it requires that the variates remain non-positively correlated after any (component-wise) monotone transformations have been applied. For , this step is typically trivial (e.g., by taking and , where ) and it constitutes much of the common use of antithetic variates in Monte Carlo simulation. Our emphasis is on because we have observed some general gains in going beyond the commonly used pair of antithetic variates. We discuss several ways of generating negatively associated random variates for arbitrary , and our comparison generally favors Iterative Latin Hypercube Sampling.

[1]  A. Owen A Central Limit Theorem for Latin Hypercube Sampling , 1992 .

[2]  Richard J. Beckman,et al.  A Comparison of Three Methods for Selecting Values of Input Variables in the Analysis of Output From a Computer Code , 2000, Technometrics.

[3]  Jørund Gåsemyr,et al.  Antithetic coupling of two Gibbs sampler chains , 2000 .

[4]  David Bruce Wilson,et al.  Exact sampling with coupled Markov chains and applications to statistical mechanics , 1996, Random Struct. Algorithms.

[5]  J. Hammersley,et al.  A new Monte Carlo technique: antithetic variates , 1956, Mathematical Proceedings of the Cambridge Philosophical Society.

[6]  C. Holmes,et al.  Perfect simulation for Bayesian curve and surface ttingC , 1999 .

[7]  ANTITHETIC VARIATES FOR MONTE CARLO ESTIMATION OF PROBABILITIES , 1984 .

[8]  M. Stein Large sample properties of simulations using latin hypercube sampling , 1987 .

[9]  Nils I. Arvidsen,et al.  Variance reduction through negative correlation, a simulation study , 1982 .

[10]  David B. Wilson How to couple from the past using a read-once source of randomness , 2000 .

[11]  David Bruce Wilson Layered Multishift Coupling for use in Perfect Sampling Algorithms (with a primer on CFTP) , 1999 .

[12]  D. Murdoch Exact Sampling for Bayesian Inference: Unbounded State Spaces , 2000 .

[13]  Geoff K. Nicholls,et al.  Perfect simulation for sample-based inference , 1999 .

[14]  P. Green,et al.  Exact Sampling from a Continuous State Space , 1998 .

[15]  David Bruce Wilson,et al.  How to couple from the past using a read-once source of randomness , 1999, Random Struct. Algorithms.

[16]  K. Joag-dev,et al.  Negative Association of Random Variables with Applications , 1983 .