News & Events

Subscribe to email list

Please select the email list(s) to which you wish to subscribe.

You are here

Multi-resolution spatial methods for large data sets (joint EOS,SCAIM,PIMS,STAT)

Tuesday, October 14, 2014 - 11:00
Dr. Doug Nychka, Director of the Institute for Mathematics Applied to Geosciences of the National Center for Atmospheric Research
Michael Smith Labs, Room 102

Spatial data is ubiquitous arising in numerous areas in the geophysical and environmental sciences. A basic problem for statisticians is to estimate complete surfaces from irregular observations or measurements and to quantify the uncertainty in the result. However, standard  statistical methods break when applied to large data sets and so alternative approaches are needed that balance shortcuts in the statistical models for increases in computational efficiency. A useful method expands the surface in a set of compact basis functions and places a Markov random field model on the basis coefficients. The impact is that evaluating the model likelihood and computing spatial predictions is feasible even for tens of thousands of spatial observations on a single computational core (e.g. a laptop). Moreover, by varying the support of the basis functions and the correlations among basis coefficients it is possible to entertain multi-resolution and non-stationary spatial models that mirror the rich covariance structure often found in large geophysical data sets.

See: A multi-resolution Gaussian process model for the analysis of large spatial data sets.

D Nychka, S Bandyopadhyay, D Hammerling, F Lindgren, S Sain (2014) Journal of Computational and Graphical Statistics (In press).