Bayesian spatial prediction

This paper presents a complete Bayesian methodology for analyzing spatial data, one which employs proper priors and features diagnostic methods in the Bayesian spatial setting. The spatial covariance structure is modeled using a rich class of covariance functions for Gaussian random fields. A general class of priors for trend, scale, and structural covariance parameters is considered. In particular, we obtain analytic results that allow easy computation of the predictive distribution for an arbitrary prior on the parameters of the covariance function using importance sampling. The computations, as well as model diagnostics and sensitivity analysis, are illustrated with a set of precipitation data.