Research

My research sits at the intersection of statistical theory, machine learning, and computation. I develop methods that remain interpretable and statistically principled while scaling to modern high-dimensional and structured datasets.

High-dimensional statistics and implicit regularization

I study estimation and prediction when the number of parameters is comparable to or larger than the sample size. This includes implicit regularization from optimization algorithms, spectral methods, principal component regression, and procedures whose finite-time behavior can differ sharply from their limiting solutions.

Research Highlight · arXiv preprint

Beyond Negative-Ridge Endpoints: Mixed-Sign Spectral Regularization via Negative-Shifted Gradient Descent

Finite-time negative-shifted gradient descent crosses a spectral barrier that constrains every stable negative-ridge endpoint, enabling head anti-shrinkage and lower-spectrum control in one path.

Paper arXiv Code

Why high-dimensional prediction behaves differently

When the number of parameters exceeds the sample size, fitting the training data does not by itself determine test performance. Prediction depends on the sample covariance spectrum: a small leading head can carry the recoverable signal, while many weak directions determine how the fitted inverse behaves.

Notation. Let \(X\in\mathbb R^{n\times p}\) and \(y=X\beta^\star+\varepsilon\) denote the observed design and response, and let \(\Sigma=\mathbb E(xx^\top)\) be the population feature covariance. In the population eigenbasis, split the observed design as \(X=[X_H\;X_T]\), where \(H\) is the signal-bearing head and \(T\) is the weak tail. We use \(\lambda\) for population covariance eigenvalues—\(\lambda_h\) for the head and \(\lambda_T\) for the flat tail—and \(\mu\) for eigenvalues of the observed Gram matrix.

\[K=\frac{XX^\top}{n}=K_H+K_T,\qquad K_H=\frac{X_HX_H^\top}{n},\quad K_T=\frac{X_TX_T^\top}{n}\approx aI_n,\quad a=\frac{\operatorname{tr}(\Sigma_T)}{n}.\]

A high-effective-rank tail can therefore be weak coordinate by coordinate yet create a substantial sample-space floor. The head is learned through \(K_H+aI_n\), so ridgeless regression experiences implicit positive shrinkage on signal-bearing directions. Negative ridge has the corrective sign, but its stable endpoint is limited by a pole and a tail-heavy filter shape.

Start with the common-spike model. If a flat tail has level \(\lambda_T\) and aspect ratio \(\gamma_T=d_T/n\), its many weak directions create the effective floor

Implicit floor\[a=\gamma_T\lambda_T.\]
Geometry-suggested shift\[\nu_\star=a+\lambda_h.\]
MP lower-edge heuristic\[\widehat\mu_{\min}^{+}\approx a-2\sqrt{a\lambda_T}.\]

The shift \(\nu_\star\) therefore lies beyond the admissible stable-endpoint range \(0\leq\nu<\widehat\mu_{\min}^{+}\). This is the Marchenko–Pastur pole barrier.

  1. 1
    Weak tail directions create an implicit ridge floor.

    Collectively, the high-effective-rank tail shrinks the signal-bearing head as if a positive ridge of size \(a\) were present.

  2. 2
    Stable negative ridge meets a spectral barrier.

    The endpoint filter \(A_\nu(\mu)=\mu/(\mu-\nu)\) has a pole and amplifies smaller stable eigenvalues most—the wrong shape for head correction with tail control.

  3. 3
    Finite time removes the endpoint restriction.

    Stopping before convergence makes the singularity removable and lets the displacement \(f_{\nu,t}(\mu)-1\) cross sign once, on a data-adaptive leading spectral prefix.

Figure A. Endpoint versus finite path. The orange curve is the limiting envelope of all admissible stable negative-ridge endpoints. The blue finite path remains smooth at its selected signed level; increasing \(\nu\) can lift its response on large, signal-rich eigenvalues above the endpoint envelope, while stopping time controls lower-spectrum exposure.

Schematic normalization: orange is the admissible endpoint envelope as \(\nu\uparrow\widehat\mu_{\min}^{+}\); blue uses the selected \((\nu,t)\).

Admissible stable-endpoint envelope and supercritical finite-time spectral filter An interactive spectral-filter plot ordered from small eigenvalues on the left to large eigenvalues on the right. The orange curve is the envelope of all admissible stable negative-ridge endpoints and rises toward the stable-endpoint barrier. The blue finite-time path is smooth at its selected signed level and, for supercritical choices, can exceed the orange envelope on large signal-rich eigenvalues. tail control head anti-shrinkage ridgeless f = 1 endpoint barrier μ̂⁺min removable μ = ν crossover small eigenvalues large eigenvalues empirical eigenvalue μ → filter

\[f_{\nu,t}(\mu)=\frac{\mu}{\mu-\nu}\left\{1-e^{-t(\mu-\nu)}\right\},\qquad f_{\nu,t}(\nu)=\nu t.\]

Figure B. The Marchenko–Pastur barrier. The population correction lies beyond the stable endpoint interval, while a finite-time path can operate there without taking an unstable inverse.
Marchenko–Pastur endpoint barrier A horizontal spectral scale marks the stable endpoint interval below the smallest positive empirical eigenvalue, the Marchenko–Pastur bulk, the tail floor a, and nu star equal to a plus lambda h beyond the stable interval. A blue arrow indicates that finite-time dynamics can use nu star. admissible stable endpoints empirical MP bulk μ̂⁺min ≈ a − 2√(aλT) a ν★ = a + λh finite-time path no endpoint convergence required

Networks and relational data

Many modern datasets describe relationships rather than independent observations. My work develops latent-space, covariate-assisted, dynamic, and spatial models for networks, with an emphasis on uncertainty quantification, robustness, and interpretable structure.

Scalable Bayesian inference

I develop variational and approximate Bayesian methods for models where conventional posterior computation is too costly. The goal is to pair computational scalability with rigorous statistical guarantees.

Multivariate and nonparametric Bayes

I am also interested in shrinkage priors, graphical models, dimension reduction, and Bayesian models that adapt to complex dependence without imposing unnecessarily rigid parametric structure.

Current support

My research on Graph Machine Learning is supported by Microsoft Research (co-PI; 2026–present).