A full scale approximation of covariance functions for large spatial data sets

Huiyan Sang, Jianhua Z. Huang

Research output: Contribution to journalArticlepeer-review

160 Scopus citations


Gaussian process models have been widely used in spatial statistics but face tremendous computational challenges for very large data sets. The model fitting and spatial prediction of such models typically require O(n 3) operations for a data set of size n. Various approximations of the covariance functions have been introduced to reduce the computational cost. However, most existing approximations cannot simultaneously capture both the large- and the small-scale spatial dependence. A new approximation scheme is developed to provide a high quality approximation to the covariance function at both the large and the small spatial scales. The new approximation is the summation of two parts: a reduced rank covariance and a compactly supported covariance obtained by tapering the covariance of the residual of the reduced rank approximation. Whereas the former part mainly captures the large-scale spatial variation, the latter part captures the small-scale, local variation that is unexplained by the former part. By combining the reduced rank representation and sparse matrix techniques, our approach allows for efficient computation for maximum likelihood estimation, spatial prediction and Bayesian inference. We illustrate the new approach with simulated and real data sets. © 2011 Royal Statistical Society.
Original languageEnglish (US)
Pages (from-to)111-132
Number of pages22
JournalJournal of the Royal Statistical Society: Series B (Statistical Methodology)
Issue number1
StatePublished - Oct 10 2011
Externally publishedYes


Dive into the research topics of 'A full scale approximation of covariance functions for large spatial data sets'. Together they form a unique fingerprint.

Cite this