Improving Pearson's chi-squared test: hypothesis testing of distributions -- optimally

Abstract

Pearson's chi-squared test, from 1900, is the standard statistical tool for "hypothesis testing on distributions": namely, given samples from an unknown distribution Q that may or may not equal a hypothesis distribution P, we want to return "yes" if P=Q and "no" if P is far from Q. While the chi-squared test is easy to use, it has been known for a while that it is not "data efficient", it does not make the best use of its data. Precisely, for accuracy ε and confidence δ, and given n samples from the unknown distribution Q, a tester should return "yes" with probability >1-δ when P=Q, and "no" with probability >1-δ when |P-Q|>ε. The challenge is to find a tester with the best tradeoff between ε, δ, and n. We introduce a new tester, efficiently computable and easy to use, which we hope will replace the chi-squared tester in practical use. Our tester is found via a new non-convex optimization framework that essentially seeks to "find the tester whose Chernoff bounds on its performance are as good as possible". This tester is 1+o(1) optimal, in that the number of samples n needed by the tester is within 1+o(1) factor of the samples needed by any tester, even non-linear testers (for the setting: accuracy ε, confidence δ, and hypothesis P). We complement this algorithmic framework with matching lower bounds saying, essentially, that "our tester is instance-optimal, even to 1+o(1) factors, to the degree that Chernoff bounds are tight". Our overall non-convex optimization framework extends well beyond the current problem and is of independent interest.

0

Turn this paper into a full lesson

ArcXiv compiles a staged curriculum from this paper: 8-12 lessons across beginner → advanced, synthesised section guides, visuals, flashcards, a quiz, exercises, and on-demand deep dives per section. Grounded in the abstract, never invented.

Discussion (0)

Sign in to join the discussion.

Loading comments…