zbMATH — the first resource for mathematics

Constructing priors that penalize the complexity of Gaussian random fields. (English) Zbl 07095889
Summary: Priors are important for achieving proper posteriors with physically meaningful covariance structures for Gaussian random fields (GRFs) since the likelihood typically only provides limited information about the covariance structure under in-fill asymptotics. We extend the recent penalized complexity prior framework and develop a principled joint prior for the range and the marginal variance of one-dimensional, two-dimensional, and three-dimensional Matérn GRFs with fixed smoothness. The prior is weakly informative and penalizes complexity by shrinking the range toward infinity and the marginal variance toward zero. We propose guidelines for selecting the hyper parameters, and a simulation study shows that the new prior provides a principled alternative to reference priors that can leverage prior knowledge to achieve shorter credible intervals while maintaining good coverage.
We extend the prior to a nonstationary GRF parameterized through local ranges and marginal standard deviations, and introduce a scheme for selecting the hyper parameters based on the coverage of the parameters when fitting simulated stationary data. The approach is applied to a dataset of annual precipitation in southern Norway and the scheme for selecting the hyper parameters leads to conservative estimates of nonstationarity and improved predictive performance over the stationary model.

62M40 Random fields; image analysis
62F15 Bayesian inference
PDF BibTeX Cite
Full Text: DOI
[1] Anderes, E., On the Consistent Separation of Scale and Variance for Gaussian Random Fields, The Annals of Statistics, 38, 870-893, (2010) · Zbl 1204.60041
[2] Berger, J. O.; Bernardo, J. M.; Sun, D., The Formal Definition of Reference Priors, The Annals of Statistics, 37, 905-938, (2009) · Zbl 1162.62013
[3] Berger, J. O.; De Oliveira, V.; Sansó, B., Objective Bayesian Analysis of Spatially Correlated Data, Journal of the American Statistical Association, 96, 1361-1374, (2001) · Zbl 1051.62095
[4] Chen, M.-H.; Shao, Q.-M., Monte Carlo Estimation of Bayesian Credible and HPD Intervals, Journal of Computational and Graphical Statistics, 8, 69-92, (1999)
[5] De Oliveira, V., Objective Bayesian Analysis of Spatial Data With Measurement Error, The Canadian Journal of Statistics/La Revue Canadienne de Statistique, 35, 283-301, (2007) · Zbl 1129.62086
[6] Ingebrigtsen, R.; Lindgren, F.; Steinsland, I., Spatial Models with Explanatory Variables in the Dependence Structure, Spatial Statistics, 8, 20-38, (2014)
[7] Ingebrigtsen, R.; Lindgren, F.; Steinsland, I.; Martino, S., Estimation of a Non-stationary Model for Annual Precipitation in Southern Norway Using Replicates of the Spatial Field, Spatial Statistics, 14, 338-364, (2015)
[8] Kazianka, H., Objective Bayesian Analysis of Geometrically Anisotropic Spatial Data, Journal of Agricultural, Biological, and Environmental Statistics, 18, 514-537, (2013) · Zbl 1303.62080
[9] Kazianka, H.; Pilz, J., Objective Bayesian Analysis of Spatial Data with Uncertain Nugget and Range Parameters, Canadian Journal of Statistics, 40, 304-327, (2012) · Zbl 1348.62082
[10] Klein, N.; Kneib, T., Scale-Dependent Priors for Variance Parameters in Structured Additive Distributional Regression, Bayesian Analysis, 11, 1071-1106, (2016) · Zbl 1357.62115
[11] Lindgren, F.; Rue, H.; Lindström, J., An Explicit Link Between Gaussian Fields and Gaussian Markov Random Fields: The Stochastic Partial Differential Equation Approach,, Journal of the Royal Statistical Society, 73, 423-498, (2011) · Zbl 1274.62360
[12] Paulo, R., Default Priors for Gaussian Processes, The Annals of Statistics, 33, 556-582, (2005) · Zbl 1069.62030
[13] Rue, H.; Martino, S.; Chopin, N., Approximate Bayesian Inference for Latent Gaussian Models by Using Integrated Nested Laplace Approximations,, Journal of the Royal Statistical Society, 71, 319-392, (2009) · Zbl 1248.62156
[14] Simpson, D.; Rue, H.; Riebler, A.; Martins, T. G.; Sørbye, S. H., Penalising Model Component Complexity: A Principled, Practical Approach to Constructing Priors, Statistical Science, 32, 1-28, (2017) · Zbl 1442.62060
[15] Stein, M. L., Interpolation of Spatial Data: Some Theory for Kriging, (1999), New York: Springer, New York · Zbl 0924.62100
[16] van der Vaart, A. W.; van Zanten, J. H., Adaptive Bayesian Estimation Using a Gaussian Random Field With Inverse Gamma Bandwidth, The Annals of Statistics, 37, 2655-2675, (2009) · Zbl 1173.62021
[17] Warnes, J.; Ripley, B., Problems with Likelihood Estimation of Covariance Functions of Spatial Gaussian Processes, Biometrika, 74, 640-642, (1987) · Zbl 0628.62087
[18] Ying, Z., Asymptotic Properties of a Maximum Likelihood Estimator With Data From a Gaussian Process, Journal of Multivariate Analysis, 36, 280-296, (1991) · Zbl 0733.62091
[19] Zhang, H., Inconsistent Estimation and Asymptotically Equal Interpolations in Model-Based Geostatistics, Journal of the American Statistical Association, 99, 250-261, (2004) · Zbl 1089.62538
This reference list is based on information provided by the publisher or from digital mathematics libraries. Its items are heuristically matched to zbMATH identifiers and may contain data conversion errors. It attempts to reflect the references listed in the original paper as accurately as possible without claiming the completeness or perfect precision of the matching.