English

Likelihood Approximation With Hierarchical Matrices For Large Spatial Datasets

Computation 2018-09-13 v2

Abstract

We use available measurements to estimate the unknown parameters (variance, smoothness parameter, and covariance length) of a covariance function by maximizing the joint Gaussian log-likelihood function. To overcome cubic complexity in the linear algebra, we approximate the discretized covariance function in the hierarchical (H-) matrix format. The H-matrix format has a log-linear computational cost and storage O(knlogn)O(kn \log n), where the rank kk is a small integer and nn is the number of locations. The H-matrix technique allows us to work with general covariance matrices in an efficient way, since H-matrices can approximate inhomogeneous covariance functions, with a fairly general mesh that is not necessarily axes-parallel, and neither the covariance matrix itself nor its inverse have to be sparse. We demonstrate our method with Monte Carlo simulations and an application to soil moisture data. The C, C++ codes and data are freely available.

Keywords

Cite

@article{arxiv.1709.04419,
  title  = {Likelihood Approximation With Hierarchical Matrices For Large Spatial Datasets},
  author = {Alexander Litvinenko and Ying Sun and Marc G. Genton and David Keyes},
  journal= {arXiv preprint arXiv:1709.04419},
  year   = {2018}
}

Comments

23 pages, 24 figures, 3 tables, version after the second major revision

R2 v1 2026-06-22T21:42:08.536Z