Source author record

Paul Valiant

Explore this author’s publications, collaborators and research topics. This catalog record has not been claimed by its author.

ResearcherUnclaimed source record

Catalog footprint

Research overview

9works
9topics
4close collaborators

Published work

9 publications

preprint2022arXiv

Finite-Sample Maximum Likelihood Estimation of Location

We consider 1-dimensional location estimation, where we estimate a parameter λλ from nn samples λ+ηiλ+ η_i, with each ηiη_i drawn i.i.d. from a known distribution ff. For fixed ff the maximum-likelihood estimate (MLE) is well-known to be optimal in the limit as n→∞n \to \infty: it is asymptotically normal with variance matching the Cramér-Rao lower bound of 1nI\frac{1}{n\mathcal{I}}, where I\mathcal{I} is the Fisher information of ff. However, this bound does not hold for finite nn, or when ff varies with nn. We show for arbitrary ff and nn that one can recover a similar theory based on the Fisher information of a smoothed version of ff, where the smoothing radius decays with nn.

preprint2021arXiv

Uncertainty about Uncertainty: Optimal Adaptive Algorithms for Estimating Mixtures of Unknown Coins

Given a mixture between two populations of coins, "positive" coins that each have -- unknown and potentially different -- bias ≥12+Δ\geq\frac{1}{2}+Δ and "negative" coins with bias ≤12−Δ\leq\frac{1}{2}-Δ, we consider the task of estimating the fraction ρρ of positive coins to within additive error εε. We achieve an upper and lower bound of Θ(ρε2Δ2log⁡1δ)Θ(\fracρ{ε^2Δ^2}\log\frac{1}δ) samples for a 1−δ1-δ probability of success, where crucially, our lower bound applies to all fully-adaptive algorithms. Thus, our sample complexity bounds have tight dependence for every relevant problem parameter. A crucial component of our lower bound proof is a decomposition lemma (see Lemmas 17 and 18) showing how to assemble partially-adaptive bounds into a fully-adaptive bound, which may be of independent interest: though we invoke it for the special case of Bernoulli random variables (coins), it applies to general distributions. We present simulation results to demonstrate the practical efficacy of our approach for realistic problem parameters for crowdsourcing applications, focusing on the "rare events" regime where ρρ is small. The fine-grained adaptive flavor of both our algorithm and lower bound contrasts with much previous work in distributional testing and learning.

preprint2020arXiv

Implicit regularization for deep neural networks driven by an Ornstein-Uhlenbeck like process

We consider networks, trained via stochastic gradient descent to minimize ℓ2\ell_2 loss, with the training labels perturbed by independent noise at each iteration. We characterize the behavior of the training dynamics near any parameter vector that achieves zero training error, in terms of an implicit regularization term corresponding to the sum over the data points, of the squared ℓ2\ell_2 norm of the gradient of the model with respect to the parameter vector, evaluated at each data point. This holds for networks of any connectivity, width, depth, and choice of activation function. We interpret this implicit regularization term for three simple settings: matrix sensing, two layer ReLU networks trained on one-dimensional data, and two layer networks with sigmoid activations trained on a single datapoint. For these settings, we show why this new and general implicit regularization effect drives the networks towards "simple" models.

preprint2016arXiv

Information Theoretically Secure Databases

We introduce the notion of a database system that is information theoretically "Secure In Between Accesses"--a database system with the properties that 1) users can efficiently access their data, and 2) while a user is not accessing their data, the user's information is information theoretically secure to malicious agents, provided that certain requirements on the maintenance of the database are realized. We stress that the security guarantee is information theoretic and everlasting: it relies neither on unproved hardness assumptions, nor on the assumption that the adversary is computationally or storage bounded. We propose a realization of such a database system and prove that a user's stored information, in between times when it is being legitimately accessed, is information theoretically secure both to adversaries who interact with the database in the prescribed manner, as well as to adversaries who have installed a virus that has access to the entire database and communicates with the adversary. The central idea behind our design is the construction of a "re-randomizing database" that periodically changes the internal representation of the information that is being stored. To ensure security, these remappings of the representation of the data must be made sufficiently often in comparison to the amount of information that is being communicated from the database between remappings and the amount of local memory in the database that a virus may preserve during the remappings. The core of the proof of the security guarantee is the following communication/data tradeoff for the problem of learning sparse parities from uniformly random nn-bit examples: any algorithm that can learn a parity of size kk with probability at least pp and extracts at most rr bits of information from each example, must see at least p⋅(nr)k/2ckp\cdot \left(\frac{n}{r}\right)^{k/2} c_k examples.

preprint2016arXiv

Optimizing Star-Convex Functions

We introduce a polynomial time algorithm for optimizing the class of star-convex functions, under no restrictions except boundedness on a region about the origin, and Lebesgue measurability. The algorithm's performance is polynomial in the requested number of digits of accuracy, contrasting with the previous best known algorithm of Nesterov and Polyak that has exponential dependence, and that further requires Lipschitz second differentiability of the function, but has milder dependence on the dimension of the domain. Star-convex functions constitute a rich class of functions generalizing convex functions to new parameter regimes, and which confound standard variants of gradient descent; more generally, we construct a family of star-convex functions where gradient-based algorithms provably give no information about the location of the global optimum. We introduce a new randomized algorithm for finding cutting planes based only on function evaluations, where, counterintuitively, the algorithm must look outside the feasible region to discover the structure of the star-convex function that lets it compute the next cut of the feasible region. We emphasize that the class of star-convex functions we consider is as unrestricted as possible: the class of Lebesgue measurable star-convex functions has theoretical appeal, introducing to the domain of polynomial-time algorithms a huge class with many interesting pathologies. We view our results as a step forward in understanding the scope of optimization techniques beyond the garden of convex optimization and local gradient-based methods.

preprint2015arXiv

Eroding dipoles and vorticity growth for Euler flows in R3 \scriptstyle{\mathbb{R}}^3 I. Axisymmetric flow without swirl

A review of analyses based upon anti-parallel vortex structures suggests that structurally stable vortex structures with eroding circulation may offer a path to the study of rapid vorticity growth in solutions of Euler's equations in R3 \scriptstyle{\mathbb{R}}^3. We examine here the possible formation of such a structure in axisymmetric flow without swirl, leading to maximal growth of vorticity as t4/3t^{4/3}. Our study suggests that the optimizing flow giving the t4/3t^{4/3} growth mimics an exact solution of Euler's equations representing an eroding toroidal vortex dipole which locally conserves kinetic energy. The dipole cross-section is a perturbation of the classical Sadovskii dipole having piecewise constant vorticity, which breaks the symmetry of closed streamlines. The structure of this perturbed Sadovskii dipole is analyzed asymptotically at large times, and its predicted properties are verified numerically.

preprint2015arXiv

Instance Optimal Learning

We consider the following basic learning task: given independent draws from an unknown distribution over a discrete support, output an approximation of the distribution that is as accurate as possible in ℓ1\ell_1 distance (i.e. total variation or statistical distance). Perhaps surprisingly, it is often possible to "de-noise" the empirical distribution of the samples to return an approximation of the true distribution that is significantly more accurate than the empirical distribution, without relying on any prior assumptions on the distribution. We present an instance optimal learning algorithm which optimally performs this de-noising for every distribution for which such a de-noising is possible. More formally, given nn independent draws from a distribution pp, our algorithm returns a labelled vector whose expected distance from pp is equal to the minimum possible expected error that could be obtained by any algorithm that knows the true unlabeled vector of probabilities of distribution pp and simply needs to assign labels, up to an additive subconstant term that is independent of pp and goes to zero as nn gets large. One conceptual implication of this result is that for large samples, Bayesian assumptions on the "shape" or bounds on the tail probabilities of a distribution over discrete support are not helpful for the task of learning the distribution. As a consequence of our techniques, we also show that given a set of nn samples from an arbitrary distribution, one can accurately estimate the expected number of distinct elements that will be observed in a sample of any size up to nlog⁡nn \log n. This sort of extrapolation is practically relevant, particularly to domains such as genomics where it is important to understand how much more might be discovered given larger sample sizes, and we are optimistic that our approach is practically viable.

preprint2013arXiv

Optimal Algorithms for Testing Closeness of Discrete Distributions

We study the question of closeness testing for two discrete distributions. More precisely, given samples from two distributions pp and qq over an nn-element set, we wish to distinguish whether p=qp=q versus pp is at least $\eps$-far from qq, in either ℓ1\ell_1 or ℓ2\ell_2 distance. Batu et al. gave the first sub-linear time algorithms for these problems, which matched the lower bounds of Valiant up to a logarithmic factor in nn, and a polynomial factor of $\eps.$ In this work, we present simple (and new) testers for both the ℓ1\ell_1 and ℓ2\ell_2 settings, with sample complexity that is information-theoretically optimal, to constant factors, both in the dependence on nn, and the dependence on $\eps$; for the ℓ1\ell_1 testing problem we establish that the sample complexity is $Θ(\max\{n^{2/3}/\eps^{4/3}, n^{1/2}/\eps^2 \}).$

preprint2011arXiv

Testing kk-Modal Distributions: Optimal Algorithms via Reductions

We give highly efficient algorithms, and almost matching lower bounds, for a range of basic statistical problems that involve testing and estimating the L_1 distance between two k-modal distributions pp and qq over the discrete domain {1,…,n}\{1,\dots,n\}. More precisely, we consider the following four problems: given sample access to an unknown k-modal distribution pp, Testing identity to a known or unknown distribution: 1. Determine whether p=qp = q (for an explicitly given k-modal distribution qq) versus pp is $\eps$-far from qq; 2. Determine whether p=qp=q (where qq is available via sample access) versus pp is $\eps$-far from qq; Estimating L1L_1 distance ("tolerant testing'') against a known or unknown distribution: 3. Approximate dTV(p,q)d_{TV}(p,q) to within additive $\eps$ where qq is an explicitly given k-modal distribution qq; 4. Approximate dTV(p,q)d_{TV}(p,q) to within additive $\eps$ where qq is available via sample access. For each of these four problems we give sub-logarithmic sample algorithms, that we show are tight up to additive $\poly(k)$ and multiplicative $\polylog\log n+\polylog k$ factors. Thus our bounds significantly improve the previous results of \cite{BKR:04}, which were for testing identity of distributions (items (1) and (2) above) in the special cases k=0 (monotone distributions) and k=1 (unimodal distributions) and required O((log⁡n)3)O((\log n)^3) samples. As our main conceptual contribution, we introduce a new reduction-based approach for distribution-testing problems that lets us obtain all the above results in a unified way. Roughly speaking, this approach enables us to transform various distribution testing problems for k-modal distributions over {1,…,n}\{1,\dots,n\} to the corresponding distribution testing problems for unrestricted distributions over a much smaller domain {1,…,ℓ}\{1,\dots,\ell\} where ℓ=O(klog⁡n).\ell = O(k \log n).

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.