IDEAS home Printed from https://ideas.repec.org/a/bla/scjsta/v36y2009i1p1-22.html
   My bibliography  Save this article

Approximate Bayesian Inference in Spatial Generalized Linear Mixed Models

Author

Listed:
  • JO EIDSVIK
  • SARA MARTINO
  • HÅVARD RUE

Abstract

In this paper we propose fast approximate methods for computing posterior marginals in spatial generalized linear mixed models. We consider the common geostatistical case with a high dimensional latent spatial variable and observations at known registration sites. The methods of inference are deterministic, using no simulation-based inference. The first proposed approximation is fast to compute and is 'practically sufficient', meaning that results do not show any bias or dispersion effects that might affect decision making. Our second approximation, an improvement of the first version, is 'practically exact', meaning that one would have to run MCMC simulations for very much longer than is typically done to detect any indication of error in the approximate results. For small-count data the approximations are slightly worse, but still very accurate. Our methods are limited to likelihood functions that give unimodal full conditionals for the latent variable. The methods help to expand the future scope of non-Gaussian geostatistical models as illustrated by applications of model choice, outlier detection and sampling design. The approximations take seconds or minutes of CPU time, in sharp contrast to overnight MCMC runs for solving such problems. Copyright (c) 2008 Board of the Foundation of the Scandinavian Journal of Statistics.

Suggested Citation

  • Jo Eidsvik & Sara Martino & Håvard Rue, 2009. "Approximate Bayesian Inference in Spatial Generalized Linear Mixed Models," Scandinavian Journal of Statistics, Danish Society for Theoretical Statistics;Finnish Statistical Society;Norwegian Statistical Association;Swedish Statistical Association, vol. 36(1), pages 1-22.
  • Handle: RePEc:bla:scjsta:v:36:y:2009:i:1:p:1-22
    as

    Download full text from publisher

    File URL: http://www.blackwell-synergy.com/doi/abs/10.1111/j.1467-9469.2008.00621.x
    File Function: link to full text
    Download Restriction: Access to full text is restricted to subscribers.

    As the access to this document is restricted, you may want to search for a different version of it.

    References listed on IDEAS

    as
    1. Steinsland, Ingelin, 2007. "Parallel exact sampling and evaluation of Gaussian Markov random fields," Computational Statistics & Data Analysis, Elsevier, vol. 51(6), pages 2969-2981, March.
    2. Yulia Gel & Adrian E. Raftery & Tilmann Gneiting, 2004. "Calibrated Probabilistic Mesoscale Weather Field Forecasting: The Geostatistical Output Perturbation Method," Journal of the American Statistical Association, American Statistical Association, vol. 99, pages 575-583, January.
    3. Peter Diggle & Søren Lophaven, 2006. "Bayesian Geostatistical Design," Scandinavian Journal of Statistics, Danish Society for Theoretical Statistics;Finnish Statistical Society;Norwegian Statistical Association;Swedish Statistical Association, vol. 33(1), pages 53-64.
    4. Skaug, Hans J. & Fournier, David A., 2006. "Automatic approximation of the marginal likelihood in non-Gaussian hierarchical models," Computational Statistics & Data Analysis, Elsevier, vol. 51(2), pages 699-709, November.
    5. Soubeyrand, Samuel & Chadoeuf, Joel, 2007. "Residual-based specification of a hidden random field included in a hierarchical spatial model," Computational Statistics & Data Analysis, Elsevier, vol. 51(12), pages 6404-6422, August.
    6. J. G. Booth & J. P. Hobert, 1999. "Maximizing generalized linear mixed model likelihoods with an automated Monte Carlo EM algorithm," Journal of the Royal Statistical Society Series B, Royal Statistical Society, vol. 61(1), pages 265-285.
    7. Ainsworth, L.M. & Dean, C.B., 2006. "Approximate inference for disease mapping," Computational Statistics & Data Analysis, Elsevier, vol. 50(10), pages 2552-2570, June.
    8. Hao Zhang, 2002. "On Estimation and Prediction for Spatial Generalized Linear Mixed Models," Biometrics, The International Biometric Society, vol. 58(1), pages 129-136, March.
    9. Paciorek, Christopher J., 2007. "Computational techniques for spatial logistic regression with large data sets," Computational Statistics & Data Analysis, Elsevier, vol. 51(8), pages 3631-3653, May.
    Full references (including those not matched with items on IDEAS)

    Citations

    Citations are extracted by the CitEc Project, subscribe to its RSS feed for this item.
    as


    Cited by:

    1. Pierrette Chagneau & Frédéric Mortier & Nicolas Picard & Jean-Noël Bacro, 2011. "A Hierarchical Bayesian Model for Spatial Prediction of Multivariate Non-Gaussian Random Fields," Biometrics, The International Biometric Society, vol. 67(1), pages 97-105, March.
    2. Håvard Rue & Sara Martino & Nicolas Chopin, 2009. "Approximate Bayesian inference for latent Gaussian models by using integrated nested Laplace approximations," Journal of the Royal Statistical Society Series B, Royal Statistical Society, vol. 71(2), pages 319-392.
    3. Martins, Thiago G. & Simpson, Daniel & Lindgren, Finn & Rue, Håvard, 2013. "Bayesian computing with INLA: New features," Computational Statistics & Data Analysis, Elsevier, vol. 67(C), pages 68-83.
    4. Eidsvik, Jo & Finley, Andrew O. & Banerjee, Sudipto & Rue, Håvard, 2012. "Approximate Bayesian inference for large spatial datasets using predictive process models," Computational Statistics & Data Analysis, Elsevier, vol. 56(6), pages 1362-1380.
    5. De Oliveira, Victor, 2013. "Hierarchical Poisson models for spatial count data," Journal of Multivariate Analysis, Elsevier, vol. 122(C), pages 393-408.
    6. Hosseini, Fatemeh & Eidsvik, Jo & Mohammadzadeh, Mohsen, 2011. "Approximate Bayesian inference in spatial GLMM with skew normal latent variables," Computational Statistics & Data Analysis, Elsevier, vol. 55(4), pages 1791-1806, April.
    7. Baghishani, Hossein & Mohammadzadeh, Mohsen, 2012. "Asymptotic normality of posterior distributions for generalized linear mixed models," Journal of Multivariate Analysis, Elsevier, vol. 111(C), pages 66-77.
    8. Higgs, Megan Dailey & Hoeting, Jennifer A., 2010. "A clipped latent variable model for spatially correlated ordered categorical data," Computational Statistics & Data Analysis, Elsevier, vol. 54(8), pages 1999-2011, August.
    9. Baghishani, Hossein & Mohammadzadeh, Mohsen, 2011. "A data cloning algorithm for computing maximum likelihood estimates in spatial generalized linear mixed models," Computational Statistics & Data Analysis, Elsevier, vol. 55(4), pages 1748-1759, April.

    More about this item

    Statistics

    Access and download statistics

    Corrections

    All material on this site has been provided by the respective publishers and authors. You can help correct errors and omissions. When requesting a correction, please mention this item's handle: RePEc:bla:scjsta:v:36:y:2009:i:1:p:1-22. See general information about how to correct material in RePEc.

    For technical questions regarding this item, or to correct its authors, title, abstract, bibliographic or download information, contact: (Wiley-Blackwell Digital Licensing) or (Christopher F. Baum). General contact details of provider: http://www.blackwellpublishing.com/journal.asp?ref=0303-6898 .

    If you have authored this item and are not yet registered with RePEc, we encourage you to do it here. This allows to link your profile to this item. It also allows you to accept potential citations to this item that we are uncertain about.

    If CitEc recognized a reference but did not link an item in RePEc to it, you can help with this form .

    If you know of missing items citing this one, you can help us creating those links by adding the relevant references in the same way as above, for each refering item. If you are a registered author of this item, you may also want to check the "citations" tab in your RePEc Author Service profile, as there may be some citations waiting for confirmation.

    Please note that corrections may take a couple of weeks to filter through the various RePEc services.

    IDEAS is a RePEc service hosted by the Research Division of the Federal Reserve Bank of St. Louis . RePEc uses bibliographic data supplied by the respective publishers.