Randomized Learning of the Second-Moment Matrix of a Smooth Function

Abstract

Consider an open set D⊂eqRn, equipped with a probability measure μ. An important characteristic of a smooth function f:D→R is its second-moment matrix μ:=∫ ∇ f(x) ∇ f(x)* μ(dx) ∈Rn× n, where ∇ f(x)∈Rn is the gradient of f(·) at x∈D and * stands for transpose. For instance, the span of the leading r eigenvectors of μ forms an active subspace of f(·), which contains the directions along which f(·) changes the most and is of particular interest in ridge approximation. In this work, we propose a simple algorithm for estimating μ from random point evaluations of f(·) without imposing any structural assumptions on μ. Theoretical guarantees for this algorithm are established with the aid of the same technical tools that have proved valuable in the context of covariance matrix estimation from partial measurements.

0

Discussion (0)

Sign in to join the discussion.

Loading comments…