K-Means Requires Exponentially Many Iterations Even in the Plane

Total Page:16

File Type:pdf, Size:1020Kb

K-Means Requires Exponentially Many Iterations Even in the Plane k-means Requires Exponentially Many Iterations Even in the Plane Andrea Vattani University of California, San Diego [email protected] ABSTRACT centers as the means (or centers of mass) of their assigned The k-means algorithm is a well-known method for parti- points. This process of assigning data points and readjusting tioning n points that lie in the d-dimensional space into k centers is repeated until it stabilizes. clusters. Its main features are simplicity and speed in prac- Despite its age, k-means is still very popular today and tice. Theoretically, however, the best known upper bound is considered \by far the most popular clustering algorithm on its running time (i.e. O(nkd)) is, in general, exponential used in scientific and industrial applications", as Berkhin re- in the number of points (when kd = Ω(n= log n)). Recently, marks in his survey on data mining [4]. Its widespread usage Arthur and Vassilvitskii [2] showed a super-polynomial worst- extends over a variety of different areas, such as artificial in- case analysis, improving the best known lower bound from telligence, computational biology, computer graphics, just to p p name a few (see [1, 8]). It is particularly popular because of Ω(n) to 2Ω( n) with a construction in d = Ω( n) dimen- its simplicity and observed speed: as Duda et al. say in their sions. In [2] they also conjectured the existence of super- text on pattern classification [6], \In practice the number of polynomial lower bounds for any d ≥ 2. iterations is much less than the number of samples". Our contribution is twofold: we prove this conjecture and we improve the lower bound, by presenting a simple con- Even if, in practice, speed is recognized as one of k-means' struction in the plane that leads to the exponential lower main qualities (see [11] for empirical studies), on the other bound 2Ω(n). hand there are a few theoretical bounds on its worst-case running time and they do not corroborate this feature. Categories and Subject Descriptors An upper bound of O(kn) can be trivially established since it can be shown that no clustering occurs twice during the F.2.2 [Analysis of Algorithms and Problem Complex- course of the algorithm. In [10], Inaba et al. improved ity]: Nonnumerical Algorithms and Problems|Geometrical this bound to O(nkd) by counting the number of Voronoi problems and computations d partitions of n points in R into k classes. Other bounds are known for some special cases. Namely, Dasgupta [5] General Terms analyzed the case d = 1, proving an upper bound of O(n) Algorithms, Theory when k < 5, and a worst-case lower bound of Ω(n). Later, Har-Peled and Sadri [9], again for the one-dimensional case, showed an upper bound of O(n∆2) where ∆ is the spread Keywords of the point set (i.e. the ratio between the largest and the K-means, Lower bounds smallest pairwise distance), and conjectured that k-means might run in time polynomial in n and ∆ for any d. kd 1. INTRODUCTION The upper bound O(n ) for the general case has not been improved since more than a decade, and this suggests that The k-means method is one of the most widely used algo- it might be not far from the truth. Arthur and Vassilvit- rithms for geometric clustering. It was originally proposed skii [2] showed that k-means can run for super-polynomially by Forgy in 1965 [7] and McQueen in 1967 [13], and is often many iterations, improvingp the best known lower bound known as Lloyd's algorithm [12]. It is a local search algo- Ω( n) from Ω(n) [5]p to 2 . Their contruction lies in a space rithm and partitions n data points into k clusters in this with d = Θ( n) dimensions, and they leave an open ques- way: seeded with k initial cluster centers, it assigns every tion about the performance of k-means for a smaller number data point to its closest center, and then recomputes the new of dimensions d, conjecturing the existence of superpolyno- mial lower bounds when d > 1. Also they show that their construction can be modified to have low spread,p disproving the aforementioned conjecture in [9] for d = Ω( n). Permission to make digital or hard copies of all or part of this work for personal or classroom use is granted without fee provided that copies are A more recent line of work that aims to close the gap not made or distributed for profit or commercial advantage and that copies between practical and theoretical performance makes use bear this notice and the full citation on the first page. To copy otherwise, to of the smoothed analysis introduced by Spielman and Teng republish, to post on servers or to redistribute to lists, requires prior specific [15]. Arthur and Vassilvitskii [3] proved a smoothed upper p permission and/or a fee. O(k) O( k) SCG’09, June 8–10, 2009, Aarhus, Denmark. bound of poly(n ), recently improved to poly(n ) by Copyright 2009 ACM 978-1-60558-501-7/09/06 ...$5.00. Manthey and R¨oglin [14]. 324 1.1 Our result We stress that when k-means runs on our constructions, it In this work we are interested in the performance of k- does not fall into any of these situations, so the lower bound means in a low dimensional space. We said it is conjectured does not exploit these degeneracies. [2] that there exist instances in d dimensions for any d ≥ 2, Our construction uses points that have constant integer for which k-means runs for a super-polynomial number of weights. This means that the data set that k-means will iterations. take in input is actually a multiset, and the center of mass Our main result is a construction in the plane (d = 2) of a cluster Ci (that is Step 2 of k-means) is computed as P P for which k-means requires exponentially many iterations to wxx= wx, where wx is the weight of x. This is x2Ci x2Ci stabilize. Specifically, we present a set of n data points ly- not a restriction since integer weights in the range [1;C] can 2 ing in R , and a set of k = Θ(n) adversarially chosen cluster be simulated by blowing up the size of the data set by at 2 Ω(n) centers in R , for which the algorithm runs for 2 iter- most C: it is enough to replace each point x of weight w with ations. This proves the aforementioned conjecture and, at a set of w distinct points (of unitary weight) whose center the samep time, it also improves the best known lower bound of mass is x, and so close to each other that the behavior of from 2Ω( n) to 2Ω(n). Notice that the exponent is optimal k-means (as well as its number of iterations) is not affected. disregarding logarithmic factor, since the bound for the gen- eral case O(nkd) can be rewritten as 2O(n log n) when d = 2 3. LOWER BOUND and k = Θ(n). For any k = o(n), our lower bound easily Ω(k) In this section we present a construction in the plane for translates to 2 , which, analogously, is almost optimal Ω(n) O(k log n) which k-means requires 2 iterations. We start with some since the upper bound is 2 . high level intuition of the construction, then we give some A common practice for seeding k-means is to choose the definitions explaining the idea behind the construction, and initial centers as a subset of the data points. We show that finally we proceed to the formal proof. even in this case (i.e. cluster centers adversarially chosen At the end of the section, we show a couple of extensions: among the data points), the running time of k-means is still the first one is a modification of our construction so that the exponential. initial set of centers is a subset of the data points, and the Also, using a result in [2], our construction can be modified second one describes how to obtain low spread. to an instance in d = 3 dimensions having low spread for Ω(n) A simple implementation in Python of the lower bound is which k-means requires 2 iterations, which disproves the available at the web address [16]. conjecture of Har-Peled and Sadri [9] for any d ≥ 3. Finally, we observe that our result implies that the smoothed 3.1 High level intuition analysis helps even for a small number of dimensions, since p The idea behind our construction is simple and can be O( k) the best smoothed upper bound is n , while our lower related to the saying \Who watches the watchmen?" (or the Ω(k) 2 bound is 2 which is larger for k = !(log n). In other original Latin phrase \Quis custodiet ipsos custodes?"). words, perturbing each data point and then running k-means Consider a sequence of t watchmen W0;W1;:::;Wt−1.A would improve the performance of the algorithm. \day" of a watchman Wi (i > 0) can be described as follows (see Fig. 1): Wi watches Wi−1, waking it up once it falls 2. THE K-MEANS ALGORITHM asleep, and does so twice; afterwards, Wi falls asleep itself. The k-means algorithm partitions a set X of n points in The watchman W0 instead will simply fall asleep directly d R into k clusters.
Recommended publications
  • The Smoothed Possibility of Social Choice
    The Smoothed Possibility of Social Choice Lirong Xia, RPI, [email protected] Abstract We develop a framework that leverages the smoothed complexity analysis by Spielman and Teng [60] to circumvent paradoxes and impossibility theorems in social choice, motivated by modern applications of social choice powered by AI and ML. For Condrocet’s paradox, we prove that the smoothed likelihood of the paradox either vanishes at an exponential rate as the number of agents in- creases, or does not vanish at all. For the ANR impossibility on the non-existence of voting rules that simultaneously satisfy anonymity, neutrality, and resolvability, we characterize the rate for the impossibility to vanish, to be either polynomially fast or exponentially fast. We also propose a novel easy-to-compute tie-breaking mechanism that optimally preserves anonymity and neutrality for even number of alternatives in natural settings. Our results illustrate the smoothed possibility of social choice—even though the paradox and the impossibility theorem hold in the worst case, they may not be a big concern in practice. 1 Introduction Dealing with paradoxes and impossibility theorems is a major challenge in social choice theory, because “the force and widespread presence of impossibility results generated a consolidated sense of pessimism, and this became a dominant theme in welfare economics and social choice theory in general”, as the eminent economist Amartya Sen commented in his Nobel prize lecture [57]. Many paradoxes and impossibility theorems in social choice are based on worst-case analysis. Take perhaps the earliest one, namely Condorcet’s (voting) paradox [15], for example. Condorcet’s para- dox states that, when there are at least three alternatives, it is impossible for pairwise majority aggregation to be transitive.
    [Show full text]
  • Smoothed Complexity Theory And, Thus, to Embed Smoothed Analysis Into Compu- Tational Complexity
    Smoothed Complexity Theory∗ Markus Bl¨aser1 Bodo Manthey2 1Saarland University, Department of Computer Science, [email protected] 2University of Twente, Department of Applied Mathematics, [email protected] Smoothed analysis is a new way of analyzing algorithms introduced by Spielman and Teng (J. ACM, 2004). Classical methods like worst-case or average-case anal- ysis have accompanying complexity classes, like P and Avg-P, respectively. While worst-case or average-case analysis give us a means to talk about the running time of a particular algorithm, complexity classes allows us to talk about the inherent difficulty of problems. Smoothed analysis is a hybrid of worst-case and average-case analysis and com- pensates some of their drawbacks. Despite its success for the analysis of single algorithms and problems, there is no embedding of smoothed analysis into com- putational complexity theory, which is necessary to classify problems according to their intrinsic difficulty. We propose a framework for smoothed complexity theory, define the relevant classes, and prove some first hardness results (of bounded halting and tiling) and tractability results (binary optimization problems, graph coloring, satisfiability). Furthermore, we discuss extensions and shortcomings of our model and relate it to semi-random models. 1 Introduction The goal of computational complexity theory is to classify computational problems according to their intrinsic difficulty. While the analysis of algorithms is concerned with analyzing, say, the running time of a particular algorithm, complexity theory rather analyses the amount of resources that all algorithms need at least to solve a given problem. Classical complexity classes, like P, reflect worst-case analysis of algorithms.
    [Show full text]
  • Smoothed Analysis of the 2-Opt Heuristic for the TSP: Polynomial Bounds for Gaussian Noise∗
    Smoothed Analysis of the 2-Opt Heuristic for the TSP: Polynomial Bounds for Gaussian Noise∗ Bodo Manthey Rianne Veenstra University of Twente, Enschede, Netherlands, Department of Applied Mathematics [email protected], [email protected] The 2-opt heuristic is a very simple local search heuristic for the traveling sales- man problem. While it usually converges quickly in practice, its running-time can be exponential in the worst case. Englert, R¨oglin,and V¨ocking (Algorithmica, 2014) provided a smoothed analysis in the so-called one-step model in order to explain the performance of 2-opt on d- dimensional Euclidean instances. However, translating their results to the classical model of smoothed analysis, where points are perturbed by Gaussian distributions with standard deviation σ, yields only bounds that are polynomial in n and 1/σd. We prove bounds that are polynomial in n and 1/σ for the smoothed running- time with Gaussian perturbations. In addition, our analysis for Euclidean distances is much simpler than the existing smoothed analysis. 1 2-Opt and Smoothed Analysis The traveling salesman problem (TSP) is one of the classical combinatorial optimization prob- lems. Euclidean TSP is the following variant: given points X ⊆ [0; 1]d, find the shortest Hamiltonian cycle that visits all points in X (also called a tour). Even this restricted variant is NP-hard for d ≥ 2 [21]. We consider Euclidean TSP with Manhattan and Euclidean dis- tances as well as squared Euclidean distances to measure the distances between points. For the former two, there exist polynomial-time approximation schemes (PTAS) [2, 20].
    [Show full text]
  • Smoothed Analysis of the TSP Algorithms
    Smoothed Analysis of the TSP Algorithms Jens Ziegler Geboren am 11. Februar 1989 in Siegburg September 22, 2012 Bachelorarbeit Mathematik Betreuer: Prof. Dr. Nitin Saxena Mathematisches Institut Mathematisch-Naturwissenschaftliche Fakultat¨ der Rheinischen Friedrich-Wilhelms-Universitat¨ Bonn Contents 1 Introduction 2 2 Preliminaries 6 3 2-Opt 8 3.1 Expected Number of 2-Changes on the L1 Metric . 8 3.1.1 Pairs of Linked 2-Changes . 10 3.1.2 Analysing the Pairs of Linked 2-Changes . 12 3.1.3 Proofs of Theorem 1 a) and Theorem 2 a) . 16 3.2 Expected Number of 2-Changes on the L2 Metric . 17 3.2.1 Analysis of a Single 2-Change . 18 3.2.2 Simplified Random Experiments . 29 3.2.3 Analysis of Pairs of Linked 2-Changes . 31 3.2.4 ProofsofTheorem1b)and2b) . 34 3.3 Expected Number of 2-Changes on General Graphs . 36 3.3.1 Definition of Witness Sequences . 36 3.3.2 Improvement to the Tour Made by Witness Sequences 37 3.3.3 Identifying Witness Sequences . 39 3.3.4 ProofofTheorem1c)and2c) . 43 3.4 ProofofTheorem3........................ 44 3.5 SmoothedAnalysisof2-Opt. 46 4 Karp’s Partitioning Scheme 48 4.1 Preliminaries ........................... 48 4.2 NecessaryLemmasandTheorems. 49 4.3 Smoothed Analysis of Karp’s Partitioning Scheme . 52 5 Conclusions and Open Questions 54 1 1 Introduction Smoothed Analysis was first introduced by Daniel A. Spielman and Shang- Hua Teng in 2001. Up until that time, there were only two standard ways of measuring an algorithms complexity: worst-case and average-case-analysis.
    [Show full text]
  • Smoothed Analysis of Local Search Algorithms
    Smoothed Analysis of Local Search Algorithms Bodo Manthey University of Twente, Department of Applied Mathematics Enschede, The Netherlands, [email protected] Abstract. Smoothed analysis is a method for analyzing the perfor- mance of algorithms for which classical worst-case analysis fails to ex- plain the performance observed in practice. Smoothed analysis has been applied to explain the performance of a variety of algorithms in the last years. One particular class of algorithms where smoothed analysis has been used successfully are local search algorithms. We give a survey of smoothed analysis, in particular applied to local search algorithms. 1 Smoothed Analysis 1.1 Motivation The goal of the analysis of algorithms is to provide measures for the performance of algorithms. In this way, it helps to compare algorithms and to understand their behavior. The most commonly used method for the performance of algorithms is worst-case analysis. If an algorithm has a good worst-case performance, then this is a very strong statement and, up to constants and lower order terms, the algorithm should also perform well in practice. However, there are many algorithms that work surprisingly well in practice although they have a very poor worst-case performance. The reason for this is that the worst-case performance can be dominated by a few pathological instances that hardly or never occur in practice. A frequently used alternative to worst-case analysis is average-case analysis. In average-case analysis, the expected performance is measured with respect to some fixed probability distribution. Many algorithms with poor worst-case but good practical performance show a good average-case performance.
    [Show full text]
  • Smoothed Analysis of Algorithms, We Measure the Expected Performance of Algorithms Under Slight Random Perturbations of Worst-Case Inputs
    ICM 2002 Vol. I 597–606 · · Smoothed Analysis of Algorithms Daniel A. Spielman∗ Shang-Hua Teng† Abstract Spielman and Teng [STOC ’01] introduced the smoothed analysis of al- gorithms to provide a framework in which one could explain the success in practice of algorithms and heuristics that could not be understood through the traditional worst-case and average-case analyses. In this talk, we survey some of the smoothed analyses that have been performed. 2000 Mathematics Subject Classification: 65Y20, 68Q17, 68Q25, 90C05. Keywords and Phrases: Smoothed analysis, Simplex method, Condition number, Interior point method, Perceptron algorithm. 1. Introduction The most common theoretical approach to understanding the behavior of al- gorithms is worst-case analysis. In worst-case analysis, one proves a bound on the worst possible performance an algorithm can have. A triumph of the Algo- rithms community has been the proof that many algorithms have good worst-case performance—a strong guarantee that is desirable in many applications. However, there are many algorithms that work exceedingly well in practice, but which are known to perform poorly in the worst-case or lack good worst-case analyses. In an attempt to rectify this discrepancy between theoretical analysis and observed performance, researchers introduced the average-case analysis of algorithms. In average-case analysis, one bounds the expected performance of an algorithm on arXiv:math/0212413v1 [math.OC] 1 Dec 2002 random inputs. While a proof of good average-case performance provides evidence that an algorithm may perform well in practice, it can rarely be understood to ex- plain the good behavior of an algorithm in practice.
    [Show full text]
  • Exploring Parameter Spaces in Coping with Computational Intractability
    Exploring Parameter Spaces in Coping with Computational Intractability Vorgelegt von Diplom-Informatiker Sepp Hartung geboren in Jena, Thüringen Von der Fakultät IV – Elektrotechnik und Informatik der Technischen Universität Berlin zur Erlangung des akademischen Grades Doktor der Naturwissenschaften – Dr. rer. nat. – genehmigte Dissertation Promotionsausschuss: Vorsitzender: Prof. Dr. Stephan Kreutzer Gutachter: Prof. Dr. Rolf Niedermeier Gutachter: Prof. Dr. Pinar Heggernes Gutachter: Prof. Dr. Dieter Kratsch Tag der wissenschaftlichen Aussprache: 06. Dezember 2013 Berlin 2014 D83 Zusammenfassung In dieser Arbeit werden verschiedene Ansätze zur Identifikation von effizient lösbaren Spezialfällen für schwere Berechnungsprobleme untersucht. Dabei liegt der Schwerpunkt der Betrachtungen auf der Ermittlung des Einflusses sogenannter Problemparameter auf die Berechnungsschwere. Für Betrach- tungen dieser Art bietet die parametrisierte Komplexitätstheorie unsere theoretische Basis. In der parametrisierten Komplexitätstheorie wird versucht sogenannte Pro- blemparameter in schweren Berechnungsproblemen (meistens, gehören diese zu den NP-schweren Problemen) zu identifizieren und diese zur Entwicklung von parametrisierten Algorithmen zu nutzen. Parametrisierte Algorithmen sind dadurch gekennzeichnet, dass der exponentielle Laufzeitanteil durch eine Funktion abhängig allein vom Problemparameter beschränkt werden kann. Die Entwicklung von parametrisierten Algorithmen ist stark motiviert durch sogenannte heuristische Algorithmen. Entgegen der theoretisch
    [Show full text]
  • Approximation Algorithms Under the Worst-Case Analysis and the Smoothed Analysis
    Approximation Algorithms under the Worst-case Analysis and the Smoothed Analysis by Weitian Tong A thesis submitted in partial fulfillment of the requirements for the degree of Doctor of Philosophy Department of Computing Science University of Alberta c Weitian Tong, 2015 Abstract How to evaluate the performance of an algorithm is a very important subject in com- puter science, for understanding its applicability, for understanding the problem which it is applied to, and for the development of new ideas that help to improve the ex- isting algorithms. There are two main factors, i.e. the performance measure and the analysis model, that affect the evaluation of an algorithm. In this thesis, we analyse algorithms with approximation ratio as a measure under the worst-case analysis and the smoothed analysis. In particular, several interesting discrete optimization problems from bioinformatics, networking and graph theory are investigated, and currently the best approximation algorithms are designed for them under the classic analysis model | the worst-case analysis and the more \realistic" model | the smoothed analysis. ii Acknowledgements It would not have been possible to write this thesis without the support and encourage- ment of many people, to whom I wish to express my sincere appreciation. Most of all, I thank my supervisor, Dr. Guohui Lin, for his invaluable guidance and support. It is difficult to summarize how many ways I am grateful to Dr. Lin, who is more like a close friend than a supervisor. Of course, I thank him mostly for his enthusiasm and for always leading my research into the right direction while leaving me a lot of freedom at the same time.
    [Show full text]
  • Beyond Worst-Case Analysis Lecture #1: Three Motivating Examples∗
    CS264: Beyond Worst-Case Analysis Lecture #1: Three Motivating Examples∗ Tim Roughgardeny September 22, 2014 1 Preamble We begin in Sections 2{4 with three examples that clearly justify the need for a course like this, and for alternatives to the traditional worst-case analysis of algorithms. Our discussion is brief and informal, to provide motivation and context. Later lectures treat these examples, and many others, in depth. 2 Caching: LRU vs. FIFO Our first example demonstrates how the worst-case analysis of algorithms can fail to differ- entiate between empirically \better" and \worse" solutions. The setup should be familiar from your system courses, if not your algorithms courses. There is a small fast memory (the \cache") and a big slow memory. For example, the former could be an on-chip cache and the latter main memory; or the former could be main memory and the latter a disk. There is a program that periodically issues read and write requests to data stored on \pages" in the big slow memory. In the happy event that the requested page is also in the cache, then it can be accessed directly. If not | on a \page fault" a.k.a. \cache miss" | the requested page needs to be brought into the cache, and this requires evicting an incumbent page. The key algorithmic question is then: which page should be evicted? For example, consider a cache that has room for four pages (Figure 1). Suppose the first four page requests are for a, b, c, and d; these are brought into the (initially empty) cache in succession.
    [Show full text]
  • Professor Shang-Hua Teng
    Professor Shang-Hua Teng Department of Computer Science, University of Southern California 1039 25 Street 3710 McClintock, RTH 505, Los Angeles, CA 90089 Santa Monica, CA 90403 (213) 740-3667, [email protected] (617) 440-4281 EMPLOYMENT: University Professor USC (2017 – ) Seeley G. Mudd Professor (Computer Science and Mathematics) USC (2012 – ) Chair and Seeley G. Mudd Professor (Computer Science) USC (2009 – 2012) Professor (Computer Science) Boston University (2002 – 2009) Research Affiliate Professor (Mathematics) MIT (1999 – Present) Visiting Research Scientist Microsoft Research Asia and New England (2004 – Present) Senior Research Scientist Akamai (1997 – 2009) Visiting Professor (Computer Science) Tsinghua University (2004 – 2008) Professor (Computer Science)[Associate Professor before 1999] UIUC (2000 – 2002) Research Scientist IBM Almaden Research Center (1997 – 1999) Assistant Professor (Computer Science) Univ. of Minnesota (1994 – 1997) Research Scientist Intel Corporation (Summer 1994, 1995) Instructor (Mathematics) MIT (1992 – 1994) Computer Scientist NASA Ames Research Center (Summer 1993) Research Scientist Xerox Palo Alto Research Center (1991 – 1992) EDUCATION: Ph.D. (Computer Science) Carnegie Mellon University (1991) Thesis: “A Unified Geometric Approach to Graph Partitioning.” Advisor: Gary Lee Miller M.S. (Computer Science) University of Southern California (1988) Thesis: “Security, Verifiablity, and Universality in Distributed Computing”. Advised by Leonard Adleman and Ming-Deh Huang B.S. (Computer Science) & B.A. (Electrical
    [Show full text]
  • 1 Introduction
    15-850: Advanced Algorithms CMU, Spring 2017 Lecture #29: Smoothed Complexity Apr 2015, Feb 2016 Lecturer: Anupam Gupta Scribe: Melanie Schmidt, Goran Zuzic 1 Introduction Smoothed analysis originates from a very influential paper of Spielman and Teng [ST01, ST04]. It provides an alternative to worst-case analysis. Assume that we want to analyze an algorithm's cost, e.g., the running time of an algorithm. Let cost(A(I)) be the cost that the algorithm has for input instance I. We will usually group the possible input instances according to their \size" (depending on the problem, this might be the number of jobs, vertices, variables and so on). Let In be the set of all inputs of `size' n. For a worst case analysis, we are now interested in max cost(A(I)); I2In the maximum cost for any input of size n. Consider Figure 29.1 and imagine that all instances of size n are lined up on the x-axis. The blue line could be the running time of an algorithm that is fast for most inputs but very slow for a small number of instances. The worst case is the height of the tallest peak. For the green curve, which could for example be the running time of a dynamic programming algorithm, the worst case is a tight bound since all instances induce the same cost. Figure 29.1: An example of a cost function with outliers. Worst case analysis provides a bound that is true for all instances, but it might be very loose for some or most instances like it happens for the blue curve.
    [Show full text]
  • Patterns, Predictions, and Actions
    PATTERNS, PREDICTIONS, AND ACTIONS A story about machine learning Moritz Hardt and Benjamin Recht arXiv:2102.05242v1 [cs.LG] 10 Feb 2021 Working draft Created: Tue 09 Feb 2021 07:25:15 PM PST Latest version available at https://mlstory.org. For Isaac, Leonora, and Quentin Contents Introduction 9 Ambitions of the 20th century 10 Pattern classification 12 Prediction and action 16 About this book 18 Chapter notes 20 Acknowledgments 21 Decision making 23 Bayesian binary hypothesis testing 25 Canonical cases of likelihood ratio tests 28 Types of errors and successes 29 The Neyman-Pearson Lemma 30 Properties of ROC curves 32 Looking ahead: what if we don’t know the models? 33 Decisions that discriminate 34 Chapter notes 38 Supervised learning 41 Sample versus population 42 Supervised learning 43 A first learning algorithm: The perceptron 45 Chapter notes 53 6 Representations and features 55 Measurement 56 Quantization 57 Template matching 59 Summarization and histograms 60 Nonlinear predictors 61 Chapter notes 72 Optimization 75 Optimization basics 76 Gradient descent 77 Applications to empirical risk minimization 80 Insights on convergence of gradient methods from quadratic functions 82 Stochastic gradient descent 84 Analysis of the stochastic gradient method 90 Implicit convexity 92 Regularization 96 Squared loss methods and other optimization tools 100 Chapter notes 101 Generalization 103 Generalization gap 103 Overparameterization: empirical phenomena 104 Theories of generalization 108 Algorithmic stability 110 Model complexity and uniform
    [Show full text]