Source author record

Wai Ming Tai

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

4works
6topics
4close collaborators

Published work

4 publications

preprint2022arXiv

Optimal Coreset for Gaussian Kernel Density Estimation

Given a point set P⊂RdP\subset \mathbb{R}^d, the kernel density estimate of PP is defined as G‾P(x)=1∣P∣∑p∈Pe−∥x−p∥2 \overline{\mathcal{G}}_P(x) = \frac{1}{\left|P\right|}\sum_{p\in P}e^{-\left\lVert x-p \right\rVert^2} for any x∈Rdx\in\mathbb{R}^d. We study how to construct a small subset QQ of PP such that the kernel density estimate of PP is approximated by the kernel density estimate of QQ. This subset QQ is called a coreset. The main technique in this work is constructing a ±1\pm 1 coloring on the point set PP by discrepancy theory and we leverage Banaszczyk's Theorem. When d>1d>1 is a constant, our construction gives a coreset of size O(1ε)O\left(\frac{1}{\varepsilon}\right) as opposed to the best-known result of O(1εlog⁡1ε)O\left(\frac{1}{\varepsilon}\sqrt{\log\frac{1}{\varepsilon}}\right). It is the first result to give a breakthrough on the barrier of log⁡\sqrt{\log} factor even when d=2d=2.

preprint2022arXiv

Optimal estimation of Gaussian DAG models

We study the optimal sample complexity of learning a Gaussian directed acyclic graph (DAG) from observational data. Our main results establish the minimax optimal sample complexity for learning the structure of a linear Gaussian DAG model in two settings of interest: 1) Under equal variances without knowledge of the true ordering, and 2) For general linear models given knowledge of the ordering. In both cases the sample complexity is n≍qlog⁡(d/q)n\asymp q\log(d/q), where qq is the maximum number of parents and dd is the number of nodes. We further make comparisons with the classical problem of learning (undirected) Gaussian graphical models, showing that under the equal variance assumption, these two problems share the same optimal sample complexity. In other words, at least for Gaussian models with equal error variances, learning a directed graphical model is statistically no more difficult than learning an undirected graphical model. Our results also extend to more general identification assumptions as well as subgaussian errors.

preprint2020arXiv

The GaussianSketch for Almost Relative Error Kernel Distance

We introduce two versions of a new sketch for approximately embedding the Gaussian kernel into Euclidean inner product space. These work by truncating infinite expansions of the Gaussian kernel, and carefully invoking the RecursiveTensorSketch [Ahle et al. SODA 2020]. After providing concentration and approximation properties of these sketches, we use them to approximate the kernel distance between points sets. These sketches yield almost (1+ε)(1+\varepsilon)-relative error, but with a small additive αα term. In the first variants the dependence on 1/α1/α is poly-logarithmic, but has higher degree of polynomial dependence on the original dimension dd. In the second variant, the dependence on 1/α1/α is still poly-logarithmic, but the dependence on dd is linear.

preprint2014arXiv

Tracking the Frequency Moments at All Times

The traditional requirement for a randomized streaming algorithm is just {\em one-shot}, i.e., algorithm should be correct (within the stated $\eps$-error bound) at the end of the stream. In this paper, we study the {\em tracking} problem, where the output should be correct at all times. The standard approach for solving the tracking problem is to run O(log⁡m)O(\log m) independent instances of the one-shot algorithm and apply the union bound to all mm time instances. In this paper, we study if this standard approach can be improved, for the classical frequency moment problem. We show that for the FpF_p problem for any 1<p≤21 < p \le 2, we actually only need O(log⁡log⁡m+log⁡n)O(\log \log m + \log n) copies to achieve the tracking guarantee in the cash register model, where nn is the universe size. Meanwhile, we present a lower bound of Ω(log⁡mlog⁡log⁡m)Ω(\log m \log\log m) bits for all linear sketches achieving this guarantee. This shows that our upper bound is tight when n=(log⁡m)O(1)n=(\log m)^{O(1)}. We also present an Ω(log⁡2m)Ω(\log^2 m) lower bound in the turnstile model, showing that the standard approach by using the union bound is essentially optimal.

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.