skip to main content
article

Parallel algorithms for Markov chain Monte Carlo methods in latent spatial Gaussian models

Published: 01 August 2004 Publication History

Abstract

Markov chain Monte Carlo (MCMC) implementations of Bayesian inference for latent spatial Gaussian models are very computationally intensive, and restrictions on storage and computation time are limiting their application to large problems. Here we propose various parallel MCMC algorithms for such models. The algorithms' performance is discussed with respect to a simulation study, which demonstrates the increase in speed with which the algorithms explore the posterior distribution as a function of the number of processors. We also discuss how feasible problem size is increased by use of these algorithms.

References

[1]
Blackford L.S., Ghol J., Cleary A.E., D'Azevedo E., Demmel J., Dhillon I., Dongarra J., Hammarling S., Henry G., Petitet A., Stanley K., Walker D., and Whaley R.C. 1997. ScaLAPACK Users' Guide, Society for Industrial & Applied Mathematics, Philadelphia.
[2]
Bradford R. and Thomas A. 1996. Markov chain Monte Carlo methods for family trees using parallel processor. Statistics and Computing 6: 67-75.
[3]
Cressie, N. 1993. Statistics for Spatial Data, 2nd edition. Wiley, New York.
[4]
Cybenko G. and Allen T.G. 1991. Multidimensional binary partitions: Distributed data structures for spatial partitioning. International Journal of Control 54: 1335-1352.
[5]
Diggle P.J., Tawn J.A., and Moyeed R.A. 1998. Model-based geostatistics (with discussion). Applied Statistics 47: 299- 350.
[6]
Duff I.S. and van der Vorst H.A. 1999. Developments and trends in the parallel solution of linear systems. Parallel Computing 25: 1931- 1970.
[7]
Gallivan K.A., Plemmons R.J., and Sameh A.H. 1990. Parallel algorithms for dense linear algebra computations. SIAM Review 32: 54-135.
[8]
Geyer C.J. and Thompson E.A. 1995. Annealing Markov chain Monte Carlo with applications to ancestral inference. J. Amer. Statist. Assoc. 90: 909-920.
[9]
Heath M.T., Ng E., and Peyton B.W. 1991. Parallel algorithms for sparse linear systems. SIAM Review 33: 420-460.
[10]
Matérn B. 1960. Spatial Variation, 2nd edition. Springer-Verlag, New York.
[11]
Neal R.M. 2002. Circularly-coupled Markov chain sampling, Technical Report 9910 (revised), Department of Statistics, University of Toronto.
[12]
Robles Sánchez O.D. and Wilson S.P. 2002. A set partitioning algorithm for use with a parallel approach to MCMC for latent Gaussian spatial models. Technical Report 02/01, Department of Statistics, Trinity College Dublin.
[13]
Rosenthal J. 2000. Parallel computing and Monte Carlo algorithms. Far East J. Theor. Stat. 4: 207-236.
[14]
Rue H. 2001. Fast sampling of Gaussian Markov random fields. Journal of the Royal Statistical Society B. 63: 325- 338.
[15]
Rue H. and Tjelmeland H. 2002. Fitting Gaussian Markov random fields to Gaussian fields. Scandanavian Journal of Statistics 29: 30-48.
[16]
Schervish M.J. 1988. Applications of parallel computation to statistical inference. Journal of the American Statistical Associsation 83: 976-983.

Cited By

View all

Recommendations

Comments

Information & Contributors

Information

Published In

cover image Statistics and Computing
Statistics and Computing  Volume 14, Issue 3
August 2004
106 pages

Publisher

Kluwer Academic Publishers

United States

Publication History

Published: 01 August 2004

Author Tags

  1. Bayesian inference
  2. Markov chain Monte Carlo
  3. latent models
  4. linear algebra
  5. parallel algorithms
  6. spatial modelling

Qualifiers

  • Article

Contributors

Other Metrics

Bibliometrics & Citations

Bibliometrics

Article Metrics

  • Downloads (Last 12 months)0
  • Downloads (Last 6 weeks)0
Reflects downloads up to 16 Oct 2024

Other Metrics

Citations

Cited By

View all

View Options

View options

Media

Figures

Other

Tables

Share

Share

Share this Publication link

Share on social media