Abstract
We first review quasi Monte Carlo (QMC) integration for approximating integrals, which we believe is a useful tool often overlooked by statistics researchers. We then present a manually-adaptive extension of QMC for approximating marginal densities when the joint density is known up to a normalization constant. Randomization and a batch-wise approach involving (0,s)-sequences are the cornerstones of our method. By incorporating a variety of graphical diagnostics the method allows the user to adaptively allocate points for joint density function evaluations. Through intelligent allocation of resources to different regions of the marginal space, the method can quickly produce reliable marginal density approximations in moderate dimensions. We demonstrate by examples that adaptive QMC can be a viable alternative to the Metropolis algorithm.
Similar content being viewed by others
Explore related subjects
Discover the latest articles, news and stories from top researchers in related subjects.References
Box, G. E. P. and Jenkins, G. M. (1976) Time Series Analysis Forecasting and Control, San Francisco: Holden Day.
Chib, S. and Greenberg, E. (1994) Bayes inference for regression models with ARMA(p,q) errors. Journal of Econometrics, 64, 183–206.
Chib, S. and Greenberg, E. (1995) Understanding the Metropolis-Hastings algorithm. American Statistician, 49, 327–35.
Cowles, M. K. and Carlin, B. P. (1996) Markov chain monte carlo convergence diagnostics: a comparitive review. Journal of the American Statistical Association, 91, 883–904.
Cranley, R. and Patterson, T. N. L. (1976) Randomization of number theoretic methods for multiple integration. SIAM Journal of Numerical Analysis, 13, 904–14.
Davis, P. J. and Rabinowitz, P. (1984) Methods of Numerical Integration, 2nd edn, Orlando: Academic Press Inc.
Fang, K. T. and Wang, Y. (1993) Number Theoretic Methods in Statistics, London: Chapman & Hall.
Faure, H. (1982) Discré pance de suites associté s à un syteme de numé ration (en dimensions). Acta Arithmetica, 41, 337–351.
Gelman, A. and Rubin, D. B. (1992) Inference from iterative simulation using multiple sequences (with discussion). Statistical Science, 7, 457–72.
Gelman, A., Roberts, G. and Gilks, W. (1995) Efficient Metropolis jumping rules. In J.M. Bernardo, J.O. Berger, A. P. Dawid and A. F. M. Smith (eds) Bayesian Statistics5. New York: Oxford Press.
Geyer, C. (1992) Practical Markov chain Monte Carlo (with discussion). Statistical Science, 7, 473–511.
Hua, L. K. and Wang, Y. (1981) Application of Number Theory to Numerical Analysis, Berlin: Springer.
Leonard, T., Hsu, J. S. J. and Ritter, C. (1994) The Laplacian tapproximation in Bayesian inference. Statistica Sinica, 4, 127–42.
Niederreiter, H. (1978) Quasi-Monte Carlo methods and pseudorandom numbers. Bulletin of the American Mathematical Society, 84, 957–1042.
Niederreiter, H. (1992) Random Number Generation and Quasi-Monte Carlo Methods, Philadelphia: Siam.
Owen, A. (1995) Randomly permuted (t,m,s)-nets and (t,s)-sequences. In H. Niederreiter and J.S. Shiue (eds) Monte Carlo and Quasi-Monte Carlo Methods in Scientific Computing NY: Springer.
Owen, A. (1996) Monte Carlo variance of scrambled equidistribution quadriture. SIAM Journal of Numerical Analysis, to appear.
Ripley, B. (1987) Stochiastic Simulation, New York: Wiley.
Rubinstein, R. Y. (1981) Simulation and the Monte Carlo Method, New York: Wiley.
Shaw, J. E. H. (1988) A quasirandom approach to integration in Bayesian statistics. Annals of Statistics, 16, 895–914.
Tierney, L. (1994) Markov chains for exploring posterior distributions. Annals of Statistics, 22, 1701–62.
Tierney, L., Kass, R. E. and Kadane, J. B. (1989) Approximate marginal densities of nonlinear functions. Biometrika, 76, 425–33.
Author information
Authors and Affiliations
Rights and permissions
About this article
Cite this article
Ostland, M., Yu, B. Exploring quasi Monte Carlo for marginal density approximation. Statistics and Computing 7, 217–228 (1997). https://doi.org/10.1023/A:1018542303861
Issue Date:
DOI: https://doi.org/10.1023/A:1018542303861