TY - JOUR
T1 - Likelihood approximation with hierarchical matrices for large spatial datasets
AU - Litvinenko, Alexander
AU - Sun, Ying
AU - Genton, Marc G.
AU - Keyes, David E.
N1 - KAUST Repository Item: Exported on 2020-10-01
Acknowledgements: The research reported in this publication was supported by funding from King Abdullah University of Science and Technology (KAUST). Additionally, we would like to express our special thanks of gratitude to Ronald Kriemann (for the HLIBPro software library) as well as to Lars Grasedyck and Steffen Börm for the HLIB software library. Alexander Litvinenko was supported by the Bayesian Computational Statistics & Modeling group (KAUST), the SRI-UQ group (KAUST), the Extreme Computing Research Center (KAUST), and the Alexander von Humboldt Foundation.
PY - 2019/2/12
Y1 - 2019/2/12
N2 - The unknown parameters (variance, smoothness, and covariance length) of a spatial covariance function can be estimated by maximizing the joint Gaussian log-likelihood function. To overcome cubic complexity in the linear algebra, the discretized covariance function is approximated in the hierarchical (H-) matrix format. The H-matrix format has a log-linear computational cost and O(knlogn) storage, where the rank k is a small integer, and n is the number of locations. The H-matrix technique can approximate general covariance matrices (also inhomogeneous) discretized on a fairly general mesh that is not necessarily axes-parallel, and neither the covariance matrix itself nor its inverse has to be sparse. It is investigated how the H-matrix approximation error influences the estimated parameters. Numerical examples with Monte Carlo simulations, where the true values of the unknown parameters are given, and an application to soil moisture data with unknown parameters are presented. The C, C++ codes and data are freely available.
AB - The unknown parameters (variance, smoothness, and covariance length) of a spatial covariance function can be estimated by maximizing the joint Gaussian log-likelihood function. To overcome cubic complexity in the linear algebra, the discretized covariance function is approximated in the hierarchical (H-) matrix format. The H-matrix format has a log-linear computational cost and O(knlogn) storage, where the rank k is a small integer, and n is the number of locations. The H-matrix technique can approximate general covariance matrices (also inhomogeneous) discretized on a fairly general mesh that is not necessarily axes-parallel, and neither the covariance matrix itself nor its inverse has to be sparse. It is investigated how the H-matrix approximation error influences the estimated parameters. Numerical examples with Monte Carlo simulations, where the true values of the unknown parameters are given, and an application to soil moisture data with unknown parameters are presented. The C, C++ codes and data are freely available.
UR - http://hdl.handle.net/10754/625430
UR - https://www.sciencedirect.com/science/article/pii/S0167947319300374
UR - http://www.scopus.com/inward/record.url?scp=85062680500&partnerID=8YFLogxK
U2 - 10.1016/j.csda.2019.02.002
DO - 10.1016/j.csda.2019.02.002
M3 - Article
SN - 0167-9473
VL - 137
SP - 115
EP - 132
JO - Computational Statistics & Data Analysis
JF - Computational Statistics & Data Analysis
ER -