Parameter estimation in high dimensional Gaussian distributions
Computation
2011-05-30 v1 Numerical Analysis
Abstract
In order to compute the log-likelihood for high dimensional spatial Gaussian models, it is necessary to compute the determinant of the large, sparse, symmetric positive definite precision matrix, Q. Traditional methods for evaluating the log-likelihood for very large models may fail due to the massive memory requirements. We present a novel approach for evaluating such likelihoods when the matrix-vector product, Qv, is fast to compute. In this approach we utilise matrix functions, Krylov subspaces, and probing vectors to construct an iterative method for computing the log-likelihood.
Cite
@article{arxiv.1105.5256,
title = {Parameter estimation in high dimensional Gaussian distributions},
author = {Erlend Aune and Daniel P. Simpson},
journal= {arXiv preprint arXiv:1105.5256},
year = {2011}
}
Comments
1 table