A Spectral Series Approach to High-Dimensional Nonparametric Regression

Abstract

A key question in modern statistics is how to make fast and reliable inferences for complex, high-dimensional data. While there has been much interest in sparse techniques, current methods do not generalize well to data with nonlinear structure. In this work, we present an orthogonal series estimator for predictors that are complex aggregate objects, such as natural images, galaxy spectra, trajectories, and movies. Our series approach ties together ideas from kernel machine learning, and Fourier methods. We expand the unknown regression on the data in terms of the eigenfunctions of a kernel-based operator, and we take advantage of orthogonality of the basis with respect to the underlying data distribution, P, to speed up computations and tuning of parameters. If the kernel is appropriately chosen, then the eigenfunctions adapt to the intrinsic geometry and dimension of the data. We provide theoretical guarantees for a radial kernel with varying bandwidth, and we relate smoothness of the regression function with respect to P to sparsity in the eigenbasis. Finally, using simulated and real-world data, we systematically compare the performance of the spectral series approach with classical kernel smoothing, k-nearest neighbors regression, kernel ridge regression, and state-of-the-art manifold and local regression methods.

References in corpus (4)

• Local polynomial regression on unknown manifolds
• Rodeo: Sparse, greedy nonparametric regression
• Regression on manifolds: Estimation of the exterior derivative
• A Spectral Series Approach to High-Dimensional Nonparametric Regression

Cited by in corpus (11)

• Cryo-EM reconstruction of continuous heterogeneity by Laplacian spectral volumes
• Conditional Density Estimation Tools in Python and R with Applications to Photometric Redshifts and Likelihood-Free Cosmological Inference
• A Spectral Series Approach to High-Dimensional Nonparametric Regression
• ABC-CDE: Towards Approximate Bayesian Computation with Complex High-Dimensional Data and Limited Simulations
• Manifold learning with arbitrary norms
• Learning Sparse Nonparametric DAGs
• Spectral decomposition of atomic structures in heterogeneous cryo-EM
• Multiscale regression on unknown manifolds
• Product Manifold Learning
• A Tale of Two Bases: Local-Nonlocal Regularization on Image Patches with Convolution Framelets
• Minimax Optimal Regression over Sobolev Spaces via Laplacian Eigenmaps on Neighborhood Graphs