Statistical Leverage and Improved Matrix Algorithms
Description
Given an m x n matrix A and a rank parameter k, define the leverage of the i-th row of A to be the i-th diagonal element of the projection matrix onto the span of the top k left singular vectors of A. In this case, "high leverage" rows have a disproportionately large amount of the "mass" in the top singular vectors. Historically, this statistical concept (and generalizations of it) has found extensive applications in diagnostic regression analysis. Recently, this concept has also been central in the development of improved randomized algorithms for several fundamental matrix problems that have broad applications in machine learning and data analysis. Two examples of the use of statistical leverage for improved worst-case analysis of matrix algorithms will be described. The first problem is the least squares approximation problem, in which there are n constraints and d variables. Classical algorithms, dating back to Gauss and Legendre, use O(nd2) time. We describe a randomized algorithm that uses only O(n d log d) time to compute a relative-error, i.e., 1+/-epsilon, approximation. The second problem is the problem of selecting a "good" set of exactly k columns from an m x n matrix, and the algorithm of Gu and Eisenstat provides the best previously existing result. We describe a two-stage algorithm that improves on their result. Recent applications of statistical leverage ideas in modern large-scale machine learning and data analysis will also be briefly described. This concept has proven to be particularly fruitful in large data applications where modeling decisions regarding what computations to perform are made for computational reasons, as opposed to having any realistic hope that the statistical assumptions implicit in those computations are satisfied by the data.
| Slides | |
| 0:00 | Statistical Leverage and Improved Matrix Algorithm |
| 2:12 | Least Squares (LS) Approximation |
| 4:16 | Many applications of this! |
| 4:53 | Exact solution to LS Approximation |
| 8:41 | Modeling with Least Squares |
| 9:27 | Statistical Issues and Regression Diagnostics |
| 13:58 | Hat Matrix and Regression Diagnostics |
| 14:22 | Statistical Leverage and Large Internet Data |
| 17:15 | Overview |
| 19:00 | Original (expensive) sampling algorithm |
| 25:34 | A "fast" LS sampling algorithm |
| 27:19 | A structural lemma |
| 28:13 | Randomized Hadamard preprocessing |
| 29:28 | Fast LS via sparse projection |
| 29:32 | Overview |
| 30:42 | Column Subset Selection Problem (CSSP) |
| 33:05 | A lower bound for CSS problem |
| 33:48 | Prior work in NLA |
| 37:24 | Working on p(k,n): 1965 - today |
| 37:47 | Theoretical computer science contributions |
| 39:20 | Prior work in TCS |
| 40:08 | The strongest Frobenius norm bound |
| 40:54 | Subspace sampling probabilities |
| 41:23 | Prior work bridging NLA/TCS |
| 42:10 | A hybrid two-stage algorithm |
| 46:52 | Comparison: spectral norm |
| 48:10 | Comparison: Frobenius norm |
| 48:52 | Overview |
| 49:25 | Empirical Evaluation: Data Sets |
| 52:22 | Empirical Evaluation: Algorithms |
| 52:50 | S&P 500 Financial Data |
| 54:37 | TechTC Term-document data |
| 54:42 | TechTC Term-document data |
| 55:07 | TechTC Term-document data |
| 57:42 | DNA HapMap SNP data |
| 58:32 | DNA HapMap SNA data |
| 59:50 | Conclusion |
Lecture rating
| People found this lecture: | ||
| Worth seeing | ||
| because it is: | ||
| Valuable and informative | ||
| Well presented | ||
| Easily understandable | ||
| Acceptably recorded | ||
| You need to login to cast your vote. | ||
Report a problem or upload files
If you have found a problem with this lecture or would like to send us extra material, articles, exercises, etc., please use our ticket system to describe your request and upload the data.Enter your e-mail into the 'Cc' field, and we will keep you updated with your request's status.
Related content
Link this page
Would you like to put a link to this lecture on your homepage?Go ahead! Copy the HTML snippet !





