Active Slices for Sliced Stein Discrepancy

Abstract

Sliced Stein discrepancy (SSD) and its kernelized variants have demonstrated promising successes in goodness-of-fit tests and model learning in high dimensions. Despite their theoretical elegance, their empirical performance depends crucially on the search of optimal slicing directions to discriminate between two distributions. Unfortunately, previous gradient-based optimisation approaches for this task return sub-optimal results: they are computationally expensive, sensitive to initialization, and they lack theoretical guarantees for convergence. We address these issues in two steps. First, we provide theoretical results stating that the requirement of using optimal slicing directions in the kernelized version of SSD can be relaxed, validating the resulting discrepancy with finite random slicing directions. Second, given that good slicing directions are crucial for practical performance, we propose a fast algorithm for finding such slicing directions based on ideas of active sub-space construction and spectral decomposition. Experiments on goodness-of-fit tests and model learning show that our approach achieves both improved performance and faster convergence. Especially, we demonstrate a 14-80x speed-up in goodness-of-fit tests when comparing with gradient-based alternatives.

arXiv comment22 pages, 7 figures, International Conference on Machine Learning (ICML) 2021

References in corpus (18)

• Adam: A Method for Stochastic Optimization
• Generative Modeling by Estimating Gradients of the Data Distribution
• Stein Variational Gradient Descent: A General Purpose Bayesian Inference Algorithm
• Measuring Sample Quality with Stein's Method
• A Kernel Test of Goodness of Fit
• Measuring Sample Quality with Kernels
• Density Estimation in Infinite Dimensional Exponential Families
• Sliced Score Matching: A Scalable Approach to Density and Score Estimation
• Gradient Estimators for Implicit Models
• A Spectral Approach to Gradient Estimation for Implicit Distributions
• A Linear-Time Kernel Goodness-of-Fit Test
• Projected Stein Variational Newton: A Fast and Scalable Bayesian Inference Method in High Dimensions
• Learning the Stein Discrepancy for Training and Evaluating Energy-Based Models without Sampling
• Random Feature Stein Discrepancies
• Statistical and Topological Properties of Sliced Probability Divergences
• Stein Neural Sampler
• Kernelized Complete Conditional Stein Discrepancy
• Nonparametric Score Estimators