A Singular Woodbury and Pseudo-Determinant Matrix Identities and
Application to Gaussian Process Regression
- URL: http://arxiv.org/abs/2207.08038v3
- Date: Mon, 24 Apr 2023 23:12:51 GMT
- Title: A Singular Woodbury and Pseudo-Determinant Matrix Identities and
Application to Gaussian Process Regression
- Authors: Siavash Ameli, Shawn C. Shadden
- Abstract summary: We study a matrix that arises from a singular form of the Woodbury matrix identity.
We present generalized inverse and pseudo-determinant identities for this matrix.
We extend the definition of the precision matrix to the Bott-Duffin inverse of the covariance matrix.
- Score: 1.5002438468152661
- License: http://arxiv.org/licenses/nonexclusive-distrib/1.0/
- Abstract: We study a matrix that arises from a singular form of the Woodbury matrix
identity. We present generalized inverse and pseudo-determinant identities for
this matrix, which have direct applications for Gaussian process regression,
specifically its likelihood representation and precision matrix. We extend the
definition of the precision matrix to the Bott-Duffin inverse of the covariance
matrix, preserving properties related to conditional independence, conditional
precision, and marginal precision. We also provide an efficient algorithm and
numerical analysis for the presented determinant identities and demonstrate
their advantages under specific conditions relevant to computing
log-determinant terms in likelihood functions of Gaussian process regression.
Related papers
- Resolvent-based quantum phase estimation: Towards estimation of parametrized eigenvalues [0.0]
We propose a novel approach for estimating the eigenvalues of non-normal matrices based on the matrix resolvent formalism.
We construct the first efficient algorithm for estimating the phases of the unit-norm eigenvalues of a given non-unitary matrix.
We then construct an efficient algorithm for estimating the real eigenvalues of a given non-Hermitian matrix.
arXiv Detail & Related papers (2024-10-07T08:51:05Z) - Multiresolution kernel matrix algebra [0.0]
We show the compression of kernel matrices by means of samplets produces optimally sparse matrices in a certain S-format.
The inverse of a kernel matrix (if it exists) is compressible in the S-format as well.
The matrix algebra is justified mathematically by pseudo differential calculus.
arXiv Detail & Related papers (2022-11-21T17:50:22Z) - Manifold Gaussian Variational Bayes on the Precision Matrix [70.44024861252554]
We propose an optimization algorithm for Variational Inference (VI) in complex models.
We develop an efficient algorithm for Gaussian Variational Inference whose updates satisfy the positive definite constraint on the variational covariance matrix.
Due to its black-box nature, MGVBP stands as a ready-to-use solution for VI in complex models.
arXiv Detail & Related papers (2022-10-26T10:12:31Z) - On confidence intervals for precision matrices and the
eigendecomposition of covariance matrices [20.20416580970697]
This paper tackles the challenge of computing confidence bounds on the individual entries of eigenvectors of a covariance matrix of fixed dimension.
We derive a method to bound the entries of the inverse covariance matrix, the so-called precision matrix.
As an application of these results, we demonstrate a new statistical test, which allows us to test for non-zero values of the precision matrix.
arXiv Detail & Related papers (2022-08-25T10:12:53Z) - Quantum algorithms for matrix operations and linear systems of equations [65.62256987706128]
We propose quantum algorithms for matrix operations using the "Sender-Receiver" model.
These quantum protocols can be used as subroutines in other quantum schemes.
arXiv Detail & Related papers (2022-02-10T08:12:20Z) - Symplectic decomposition from submatrix determinants [0.0]
An important theorem in Gaussian quantum information tells us that we can diagonalise the covariance matrix of any Gaussian state via a symplectic transformation.
Inspired by a recently presented technique for finding the eigenvectors of a Hermitian matrix from certain submatrix eigenvalues, we derive a similar method for finding the diagonalising symplectic from certain submatrix determinants.
arXiv Detail & Related papers (2021-08-11T18:00:03Z) - Robust 1-bit Compressive Sensing with Partial Gaussian Circulant
Matrices and Generative Priors [54.936314353063494]
We provide recovery guarantees for a correlation-based optimization algorithm for robust 1-bit compressive sensing.
We make use of a practical iterative algorithm, and perform numerical experiments on image datasets to corroborate our results.
arXiv Detail & Related papers (2021-08-08T05:28:06Z) - Variance Reduction for Matrix Computations with Applications to Gaussian
Processes [0.0]
We focus on variance reduction for matrix computations via matrix factorization.
We show how computing the square root factorization of the matrix can achieve in some important cases arbitrarily better performance.
arXiv Detail & Related papers (2021-06-28T10:41:22Z) - Non-PSD Matrix Sketching with Applications to Regression and
Optimization [56.730993511802865]
We present dimensionality reduction methods for non-PSD and square-roots" matrices.
We show how these techniques can be used for multiple downstream tasks.
arXiv Detail & Related papers (2021-06-16T04:07:48Z) - Understanding Implicit Regularization in Over-Parameterized Single Index
Model [55.41685740015095]
We design regularization-free algorithms for the high-dimensional single index model.
We provide theoretical guarantees for the induced implicit regularization phenomenon.
arXiv Detail & Related papers (2020-07-16T13:27:47Z) - Relative Error Bound Analysis for Nuclear Norm Regularized Matrix Completion [101.83262280224729]
We develop a relative error bound for nuclear norm regularized matrix completion.
We derive a relative upper bound for recovering the best low-rank approximation of the unknown matrix.
arXiv Detail & Related papers (2015-04-26T13:12:16Z)
This list is automatically generated from the titles and abstracts of the papers in this site.
This site does not guarantee the quality of this site (including all information) and is not responsible for any consequences.