Global and Local Two-Sample Tests via Regression

Abstract

Two-sample testing is a fundamental problem in statistics. Despite its long history, there has been renewed interest in this problem with the advent of high-dimensional and complex data. Specifically, in the machine learning literature, there have been recent methodological developments such as classification accuracy tests. The goal of this work is to present a regression approach to comparing multivariate distributions of complex data. Depending on the chosen regression model, our framework can efficiently handle different types of variables and various structures in the data, with competitive power under many practical scenarios. Whereas previous work has been largely limited to global tests which conceal much of the local information, our approach naturally leads to a local two-sample testing framework in which we identify local differences between multivariate distributions with statistical confidence. We demonstrate the efficacy of our approach both theoretically and empirically, under some well-known parametric and nonparametric regression methods. Our proposed methods are applied to simulated data as well as a challenging astronomy data set to assess their practical usefulness.

References in corpus (12)

• The Relationship Between Stellar Light Distributions of Galaxies and their Formation Histories
• A New Non-Parametric Approach to Galaxy Morphological Classification
• Analysis of a Random Forests Model
• Galaxy Morphology and Star Formation in the Illustris Simulation at z=0
• Adaptive Concentration of Regression Trees, with Application to Random Forests
• Local polynomial regression on unknown manifolds
• New Image Statistics for Detecting Disturbed Galaxy Morphologies at High Redshift
• k-NN Regression Adapts to Local Intrinsic Dimension
• Adaptivity and Computation-Statistics Tradeoffs for Kernel and Distance based High Dimensional Two Sample Testing
• Comparing two samples by penalized logistic regression
• The Classification Permutation Test: A Nonparametric Test for Equality of Multivariate Distributions
• Rates of convergence for nearest neighbor estimators with the smoother regression function