Georgia Institute of Technology School of Mathematics | Georgia Institute of Technology | Atlanta, GA

Random effects models are commonly used to measure genetic variance-covariance matrices of quantitative phenotypic traits. The population eigenvalues of these matrices describe the evolutionary response to selection. However, they may be difficult to estimate from limited samples when the number of traits is large. In this talk, I will present several results describing the eigenvalues of classical MANOVA estimators of these matrices, including dispersion of the bulk eigenvalue distribution, bias and aliasing of large "spike" eigenvalues, and distributional limits of eigenvalues at the spectral edges. I will then discuss a new procedure that uses these results to obtain better estimates of the large population eigenvalues when there are many traits, and a Tracy-Widom test for detecting true principal components in these models. The theoretical results extend proof techniques in random matrix theory and free probability, which I will also briefly describe.This is joint work with Iain Johnstone, Yi Sun, Mark Blows, and Emma Hine.

In this talk, I will discuss some examples of sparse signal detection problems in the context of binary outcomes. These will be motivated by examples from next generation sequencing association studies, understanding heterogeneities in large scale networks, and exploring opinion distributions over networks. Moreover, these examples will serve as templates to explore interesting phase transitions present in such studies. In particular, these phase transitions will be aimed at revealing a difference between studies with possibly dependent binary outcomes and Gaussian outcomes. The theoretical developments will be further complemented with numerical results.

High-dimensional data arise in many fields of contemporary science and introduce new challenges in statistical learning due to the well-known curse of dimensionality. Many data sets in image analysis and signal processing are in a high-dimensional space but exhibit a low-dimensional structure. We are interested in building efficient representations of these data for the purpose of compression and inference, and giving performance guarantees that are only cursed by the intrinsic dimension of data. Specifically, in the setting where a data set in $R^D$ consists of samples from a probability measure concentrated on or near an unknown $d$-dimensional manifold with $d$ much smaller than $D$, we consider two sets of problems: low-dimensional geometric approximation to the manifold and regression of a function on the manifold. In the first case we construct multiscale low-dimensional empirical approximations to the manifold and give finite-sample performance guarantees. In the second case we exploit these empirical geometric approximations of the manifold to construct multiscale approximations to the function. We prove finite-sample guarantees showing that we attain the same learning rates as if the function was defined on a Euclidean domain of dimension $d$. In both cases our approximations can adapt to the regularity of the manifold or the function even when this varies at different scales or locations. All algorithms have complexity $C n\log (n)$ where $n$ is the number of samples, and the constant $C$ is linear in $D$ and exponential in $d$.

Asymptotic equivalence between two statistical models means that they have the same asymptotic (large sample) properties with respect to all decision problems with bounded loss. In nonparametric (infinite-dimensional) statistical models, asymptotic equivalence has been found to be useful since it can allow one to derive certain results by studying simpler models. One of the key results in this area is Nussbaum’s theorem, which states that nonparametric density estimation is asymptotically equivalent to a Gaussian shift model, provided that the densities are smooth enough and uniformly bounded away from zero.We will review the notion of asymptotic equivalence and existing results, before presenting recent work on the extent to which one can relax the assumption of being bounded away from zero. We further derive the optimal (Le Cam) distance between these models, which quantifies how close they are for finite-samples. As an application, we also consider Poisson intensity estimation with low count data. This is joint work with Johannes Schmidt-Hieber.