# Randomised regression

January 14, 2017 — December 1, 2020

Tackling your regression, by using random embeddings of the predictors (and/or predictions?). Usually this means using low-dimensional projections, to reduce the dimensionality of a high dimensional regression. In this case it is not far from compressed sensing, except in how we handle noise. In this linear model case, this is of course random linear algebra, and may be a randomised matrix factorisation. You can do it the other way and prject something into a higher dimensional space, which is a popular trick for kernel approximation.

I am especially interested in seeing how this might be useful for dependent data, especially time series.

Brian McWilliams, Gabriel Krummenacher and Mario Lučić, Randomized Linear Regression: A brief overview and recent results. Gabriel implemented some of the algorithms mentioned, e.g.

- Subsampled Randomized Fourier Transform.
- SRHT: Uses the Subsampled Randomized Hadamard Transform (SRHT), equivalent to leverage based sampling.
- aIWS, aRWS: Samples based on approximated statistical influence.
- Uluru: SRHT with bias correction.

Martin Wainright, Statistics meets Optimization: Randomization and approximation for high-dimensional problems.

In the modern era of high-dimensional data, the interface between mathematical statistics and optimization has become an increasingly vibrant area of research. In this course, we provide some vignettes into this interface, including the following topics:

Dimensionality reduction via random projection. The naive idea of projecting high-dimensional data to a randomly chosen low-dimensional space is remarkably effective. We discuss the classical Johnson-Lindenstrauss lemma, as well as various modern variants that provide computationally-efficient embeddings with strong guarantees.

When is it possible to quickly obtain approximate solutions of large-scale convex programs? In practice, methods based on randomized projection can work very well, and arguments based on convex analysis and concentration of measure provide a rigorous underpinning to these observations.

Optimization problems with some form of nonconvexity arise frequently in statistical settings – for instance, in problems with latent variables, combinatorial constraints, or rank constraints. Nonconvex programs are known to be intractable in a complexity-theoretic sense, but the random ensembles arising in statistics are not adversarially constructed. Under what conditions is it possible to make rigorous guarantees about the behavior of simple iterative algorithms for such problems? We develop some general theory for addressing these questions, exploiting tools from both optimization theory and empirical process theory.

## 1 References

*arXiv:1702.05327 [Cs, Math, Stat]*.

*arXiv:1703.00864 [Stat]*.

*IEEE Transactions on Electronic Computers*.

*Advances in Neural Information Processing Systems 28*. NIPS’15.

*Advances in Neural Information Processing Systems*.

*arXiv:1705.08664 [Cs, Stat]*.

*arXiv:2007.07383 [Physics, Stat]*.

*arXiv:1706.07180 [Cs, Math, Stat]*.

*arXiv:1605.07913 [Math, Stat]*.

*arXiv:1406.3469 [Stat]*.

*Advances in Neural Information Processing Systems 29*.

*Randomized Algorithms for Matrices and Data*.

*Advances in Neural Information Processing Systems*.

*Advances in Neural Information Processing Systems*.

*2008 46th Annual Allerton Conference on Communication, Control, and Computing*.

*Wiley Interdisciplinary Reviews: Data Mining and Knowledge Discovery*.

*Advances in Neural Information Processing Systems 29*.

*arXiv:1702.05181 [Cs, Stat]*.

*arXiv:1701.05325 [Math, Stat]*.

*arXiv:1702.01166 [Stat]*.

*arXiv:1701.00481 [Stat]*.