Efficient Covariance Approximations for Large Sparse Precision Matrices
Artikel i vetenskaplig tidskrift, 2018

Published with license by Taylor & Francis Group, LLC. The use of sparse precision (inverse covariance) matrices has become popular because they allow for efficient algorithms for joint inference in high-dimensional models. Many applications require the computation of certain elements of the covariance matrix, such as the marginal variances, which may be nontrivial to obtain when the dimension is large. This article introduces a fast Rao–Blackwellized Monte Carlo sampling-based method for efficiently approximating selected elements of the covariance matrix. The variance and confidence bounds of the approximations can be precisely estimated without additional computational costs. Furthermore, a method that iterates over subdomains is introduced, and is shown to additionally reduce the approximation errors to practically negligible levels in an application on functional magnetic resonance imaging data. Both methods have low memory requirements, which is typically the bottleneck for competing direct methods.

Sparse precision matrix

Stochastic approximation

Selected inversion

Gaussian Markov random fields

Spatial analysis

Författare

Per Sidén

Linköpings universitet

Finn Lindgren

University of Edinburgh

David Bolin

Chalmers, Matematiska vetenskaper, Tillämpad matematik och statistik

Mattias Villani

Linköpings universitet

Journal of Computational and Graphical Statistics

1061-8600 (ISSN) 1537-2715 (eISSN)

Vol. 27 4 898-909

Ämneskategorier

Beräkningsmatematik

Sannolikhetsteori och statistik

Signalbehandling

DOI

10.1080/10618600.2018.1473782

Mer information

Senast uppdaterat

2019-01-08