Abstract
This paper deals with an empirical Bayes approach for spatial prediction of a Gaussian random field. In fact, we estimate the hyperparameters of the prior distribution by using the maximum likelihood method. In order to maximize the marginal distribution of the data, the EM algorithm is used. Since this algorithm requires the evaluation of analytically intractable and high dimensionally integrals, a Monte Carlo method based on discretizing parameter space, is proposed to estimate the relevant integrals. Then, the approach is illustrated by its application to a spatial data set. Finally, we compare the predictive performance of this approach with the reference prior method.
Similar content being viewed by others
References
Abramowitz M, Stegun I (1965) Handbook of mathematical functions. Dover, New York
Banerjee S (2005) On Geodetic distance computations in spatial modeling. Biometrics 61: 617–625
Banerjee S, Carlin BP, Gelfand AE (2004) Hierarchical modeling and analysis for spatial data. Chapman and Hall/CRC, Boca Raton
Berger JO, De Oliveira V, Sanso B (2001) Objective Bayesian analysis of spatially correlated data. J Am Stat Assoc 96: 1361–1374
Carlin BP, Louis TA (2001) Bayes and Empirical Bayes Methods for Data Analysis. Chapman and Hall, Boca Raton
Chiles JP, Delfiner P (1999) Geostatistics; modelling spatial uncertainty. Wiley, New York
Christensen OF (2004) Monte Carlo maximum likelihood in model-based geostatistics. J Computat Graphical Stat 13: 702–718
Cressie N (1993) Statistics for spatial data. Wiley, New York
Cowles MK (2003) Efficient model-fitting and model-comparison for high-dimensional Bayesian geostatistical models. J Stat Plann Inference 111: 221–239
Dempster AP, Laird NM, Rubin DB (1977) Maximum likelihood estimation from incomplete data via the EM algorithm. J Roy Stat Soc B 96: 1361–1374
De Oliveira V (2005) Bayesian inference and prediction of Gaussian random fields based on censored data. J Computat Graphical Stat 14: 95–115
De Oliveira V (2000) Bayesian prediction of clipped Gaussian random fields. Computat Stat Data Anal 34: 299–314
De Oliveira V, Kedem B, Short DA (1997) Bayesian prediction of transformed Gaussian random fields. J Am Stat Assoc 92: 1442–1433
Diggle PJ, Tawn JA, Moyeed RA (1998) Model based geostatistics. Appl Stat 47: 299–350
Gelfand AE, Dey DK, Chang H (1992) Model determination using predictive distributions with implementation via sampling-based methods. Bayesian Stat 4: 147–167
Heagerty PJ, Lele SR (1998) A composite likelihood approach to binary data in space. J Am Stat Assoc 93: 1099–1111
Meilijson I (1989) A fast improvement to the EM algorithm on its own terms. J Roy Stat Soc B 51: 127–138
McLachlan G, Krishnan T (1997) The EM algorithm and extensions. Wiley, New York
Stein M (1999) Spatial interpolation, some theory for kriging. Springer, New York
Palacios MB, Steel MJF (2006) Non-Gaussian Bayesian geostatistical modeling. J Am Stat Assoc 101: 604–618
Tanner MA (1996) Tools for statistical inference: methods for the exploration of posterior distributions and likelihood functions, 3rd edn. Springer, New York
Author information
Authors and Affiliations
Corresponding author
Rights and permissions
About this article
Cite this article
Khaledi, M.J., Rivaz, F. Empirical Bayes spatial prediction using a Monte Carlo EM algorithm. Stat Methods Appl 18, 35–47 (2009). https://doi.org/10.1007/s10260-007-0081-5
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s10260-007-0081-5