event thumbnail image
Workshops

Statistical Leverage and Improved Matrix Algorithms

author: Michael Mahoney, Department of Computer Science, Stanford University

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.

You might be experiencing some problems with Your Video player.
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.

Link this page

Would you like to put a link to this lecture on your homepage?
Go ahead! Copy the HTML snippet !

Write your own review or comment:

make sure you have javascript enabled or clear this field: