Interpolating the Trace of the Inverse of Matrix 𝐀 + t 𝐁

by   Siavash Ameli, et al.

We develop heuristic interpolation methods for the function t ↦trace( (𝐀 + t 𝐁)^-1), where the matrices 𝐀 and 𝐁 are symmetric and positive definite and t is a real variable. This function is featured in many applications in statistics, machine learning, and computational physics. The presented interpolation functions are based on the modification of a sharp upper bound that we derive for this function, which is a new trace inequality for matrices. We demonstrate the accuracy and performance of the proposed method with numerical examples, namely, the marginal maximum likelihood estimation for linear Gaussian process regression and the estimation of the regularization parameter of ridge regression with the generalized cross-validation method.


page 1

page 2

page 3

page 4

βˆ™ 05/20/2020

On randomized trace estimates for indefinite matrices with an application to determinants

Randomized trace estimation is a popular and well studied technique that...
βˆ™ 07/14/2023

Comparing Scale Parameter Estimators for Gaussian Process Regression: Cross Validation and Maximum Likelihood

Gaussian process (GP) regression is a Bayesian nonparametric method for ...
βˆ™ 08/10/2018

A new Newton-type inequality and the concavity of a class of k-trace functions

In this paper, we prove a new Newton-type inequality that generalizes Ne...
βˆ™ 03/17/2023

Scalable Physics-based Maximum Likelihood Estimation using Hierarchical Matrices

Physics-based covariance models provide a systematic way to construct co...
βˆ™ 04/23/2015

Regularization-free estimation in trace regression with symmetric positive semidefinite matrices

Over the past few years, trace regression models have received considera...
βˆ™ 03/21/2021

Detecting Label Noise via Leave-One-Out Cross-Validation

We present a simple algorithm for identifying and correcting real-valued...
βˆ™ 02/11/2013

The trace norm constrained matrix-variate Gaussian process for multitask bipartite ranking

We propose a novel hierarchical model for multitask bipartite ranking. T...

Please sign up or login with your details

Forgot password? Click here to reset