Fourier analysis perspective for sufficient dimension reduction problem

Abstract

A theory of sufficient dimension reduction (SDR) is developed from an optimizational perspective. In our formulation of the problem, instead of dealing with raw data, we assume that our ground truth includes a mapping f: Rn→ Rm and a probability distribution function p over Rn, both given analytically. We formulate SDR as a problem of finding a function g: Rk→ Rm and a matrix P∈ Rk× n such that E x p( x) | f( x) - g(P x)|2 is minimal. It turns out that the latter problem allows a reformulation in the dual space, i.e. instead of searching for g(P x) we suggest searching for its Fourier transform. First, we characterize all tempered distributions that can serve as the Fourier transform of such functions. The reformulation in the dual space can be interpreted as a problem of finding a k-dimensional linear subspace S and a tempered distribution t supported in S such that t is "close" in a certain sense to the Fourier transform of f. Instead of optimizing over generalized functions with a k-dimensional support, we suggest minimizing over ordinary functions but with an additional term R that penalizes a strong distortion of the support from any k-dimensional linear subspace. For a specific case of R, we develop an algorithm that can be formulated for functions given in the initial form as well as for their Fourier transforms. Eventually, we report results of numerical experiments with a discretized version of the latter algorithm.

0

Turn this paper into a lesson

ArcXiv compiles a structured reading guide from this paper's metadata: plain-English importance, contributions, prerequisite concepts, which sections to read first, flashcards, and a quiz. Grounded in the abstract, never invented.

Discussion (0)

Sign in to join the discussion.

Loading comments…