A multiresolution Gaussian process model for the analysis of large spatial datasets

AMS Citation:
Nychka, D., S. Bandyopadhyay, D. Hammerling, F. Lindgren, and S. Sain, 2015: A multiresolution Gaussian process model for the analysis of large spatial datasets. Journal of Computational and Graphical Statistics, 24, 579-599, doi:10.1080/10618600.2014.914946.
Date:2015-06-16
Resource Type:article
Title:A multiresolution Gaussian process model for the analysis of large spatial datasets
Abstract: We develop a multiresolution model to predict two-dimensional spatial fields based on irregularly spaced observations. The radial basis functions at each level of resolution are constructed using a Wendland compactly supported correlation function with the nodes arranged on a rectangular grid. The grid at each finer level increases by a factor of two and the basis functions are scaled to have a constant overlap. The coefficients associated with the basis functions at each level of resolution are distributed according to a Gaussian Markov random field (GMRF) and take advantage of the fact that the basis is organized as a lattice. Several numerical examples and analytical results establish that this scheme gives a good approximation to standard covariance functions such as the Matérn and also has flexibility to fit more complicated shapes. The other important feature of this model is that it can be applied to statistical inference for large spatial datasets because key matrices in the computations are sparse. The computational efficiency applies to both the evaluation of the likelihood and spatial predictions.
Peer Review:Refereed
Copyright Information:Copyright 2015 Taylor & Francis.
OpenSky citable URL: ark:/85065/d7h13366
Publisher's Version: 10.1080/10618600.2014.914946
Author(s):
  • Doug Nychka - NCAR/UCAR
  • Soutir Bandyopadhyay
  • Dorit Hammerling - NCAR/UCAR
  • Finn Lindgren - NCAR/UCAR
  • Stephan Sain - NCAR/UCAR
  • Random Profile

    SOFT ENG/PROG III

    Recent & Upcoming Visitors