Estimating basis functions in massive fields under the spatial mixed effects model
- National Security Directorate, Pacific Northwest National Laboratory Richland Washington USA, Department of Statistics North Carolina State University Raleigh North Carolina USA
- Department of Statistics Iowa State University Ames Iowa
Abstract Spatial prediction is commonly achieved under the assumption of a Gaussian random field by obtaining maximum likelihood estimates of parameters, and then using the kriging equations to arrive at predicted values. For massive datasets, fixed rank kriging using the expectation–maximization algorithm for estimation has been proposed as an alternative to the usual but computationally prohibitive kriging method. The method reduces computation cost of estimation by redefining the spatial process as a linear combination of basis functions and spatial random effects. A disadvantage of this method is that it imposes constraints on the relationship between the observed locations and the knots. We develop an alternative method that utilizes the spatial mixed effects model, but allows for additional flexibility by estimating the range of the spatial dependence between the observations and the knots via an alternating expectation conditional maximization algorithm. Experiments show that our methodology improves estimation without sacrificing prediction accuracy while also minimizing the additional computational burden of extra parameter estimation. The methodology is applied to a temperature dataset archived by the United States National Climate Data Center, with improved results over previous methodology.
- Sponsoring Organization:
- USDOE
- OSTI ID:
- 1811549
- Journal Information:
- Statistical Analysis and Data Mining, Journal Name: Statistical Analysis and Data Mining Journal Issue: 5 Vol. 14; ISSN 1932-1864
- Publisher:
- Wiley Blackwell (John Wiley & Sons)Copyright Statement
- Country of Publication:
- United States
- Language:
- English
Similar Records
Empirical Bayes regionalization techniques for estimation in ground water models
End-to-end learning of multiple sequence alignments with differentiable Smith–Waterman