Two centuries of tools, and three lamps that lit at once in 2006–2012
Show a machine many examples and let it find the rules. Many of the tools behind this came out of statistics, bee classification, parole decisions in prisons, research on sending signals, and psychological tests. In the seven years from 2006, deep neural networks (deep learning) suddenly grew powerful because GPUs and huge amounts of data joined that toolbox.
012006 · The start of deep learning
Three conditions came together
What was new in 2006: an idea, or the conditions?
Neural networks with many stacked layers could not learn well with backpropagation (Chapter 8) alone. In 2006 Geoffrey Hinton and his colleagues showed a way to : first, with no correct answers, the network learns the features of the data (pre-training), and then the whole thing is finished off with backpropagation. The error rate on the handwritten-digit set MNIST was 1.25%.
A Science paper from the same year said it had been . It named speed of computation, amount of data, and a good starting point for the weights, and went on:
In the figure on the right, switch the lamps for the three conditions on and off to compare reconstructions of images of vertical and horizontal bars. When the starting weights are set by pre-training, the network learns one layer at a time with no correct answers, and then the whole thing is finished off with backpropagation (the effects of the lamps are an example made for this site).
Δwij =η(⟨vihj⟩data − ⟨vihj⟩reconstruction)
v, h
neurons of the visible layer (pixels) and the hidden layer
w
weight
⟨ ⟩
average (the weights are corrected by the gap between data and reconstruction; the method of the restricted Boltzmann machine)
What you think you knowThe popular story:In 2006, Hinton invented a brand-new idea called deep learning The original source says:The paper itself said it had been obvious since the 1980s, and what was new was mainly . In the same year, Yoshua Bengio and colleagues published a similar method.
Why one layer at a time
If you train a deep network by backpropagation from random starting weights, the signal that travels backward tends to weaken at every layer, or the network tends to settle on a poor solution. Starting from weights that had already grasped the features of the data one layer at a time made the adjustment of the whole network go well.
It later turned out that, with smarter ways of choosing the starting weights and with the activation function ReLU, networks can be trained without pre-training (see the spread in this chapter, "Keeping the signal alive in deep stacks" (深く重ねても、信号が消えないように)).
Sources & confidence (1)Three conditions came together
28 July 2006Verified in the original
Hinton and colleagues present a way to train a deep neural network “one layer at a time”. It is seen as an important step toward
Neural networks with many stacked layers could not be trained well by backpropagation alone. Hinton and colleagues showed a method in which, first, one layer at a time and without being given correct answers, the network learns the of the data (pre-training), and then the whole network is adjusted by backpropagation. It reached an error rate of 1.25% on the handwritten digits of MNIST (the Neural Computation paper), and it was also better than earlier methods at representing images and documents with a small number of values (the Science paper).
The popular story:“In 2006 Hinton invented the new idea of deep learning,” people tend to say The original source says:The paper itself (Science) writes that it has been for the job of representing data with a small number of values. This needs computers that are fast enough, data that is large enough, and initial weights close enough to a good solution, and it says What was new was mainly what the paper calls . In the same year, 2006, Bengio and colleagues also published a similar method (reference 34 in the Nobel committee's explanation).
What the figure shows:What was new in 2006 was not the idea of deep networks, but 3 conditions coming together: fast computation, plentiful data, and good starting weights.
Reconstruction error (it shrinks as learning proceeds)
This is the same mechanism as "Try it", written as Python code. Press ▶ Run and the result appears below. Change the numbers and run it again.
u09_rbm_layer.py
021763, 1961 · Bayes
From the result back to the cause
Can you classify the field of a document from clue words alone?
Bayes' theorem is the idea of recalculating how likely a cause or hypothesis is after you obtain evidence. The paper was read at the Royal Society in England on December 23, 1763. Its author, Thomas Bayes, had already died; it was his friend Price who found it among his papers and sent it in.
In 1961, M. E. Maron of the RAND Corporation used this idea in an experiment that sorted paper abstracts into 32 fields. He assumed that, once the field is fixed, words occur independently of one another, and simply multiplied the probabilities of the individual words (today called simple Bayes, or naive Bayes). On 85 new abstracts that contained two or more clue words, the top-ranked field was correct for 44 (51.8%), and for 68 (80%) if the top three fields were allowed.
In the figure on the right, pressing the clue words changes the probabilities of the fields by multiplication (the words and probabilities are an example made for this site).
P(field | words) ∝ P(field) × Π P(word | field)
P(A | B)
probability that A happens given condition B
Π
multiply together across the words
∝
proportional to (at the end, divide so that the total is 1)
What you think you knowThe popular story:Naive (simple) Bayes is a method that naively believes words are independent of each other The original source says:Maron himself wrote that the independence assumption is clearly not correct, and that he adopted it to make the calculation simple. The paper contains neither the word "naive" nor the name "Bayes".
Training score versus test score
Maron counted words in 260 of the 405 abstracts from 1959 (90 clue words), and kept the other 145 completely unseen during preparation, to use for testing. In the paper's figures, the top-rank accuracy was 84.6% on the documents used for training and 51.8% on new documents. If you measure only on the data you learned from, you overestimate the performance. This idea of “testing on separate data” (分けて確かめる) comes up again and again until the end of this chapter.
Sources & confidence (2)From the result back to the cause
23 December 1763Verified in the original
Bayes’s essay: in 1763 his friend Price finds the manuscript among his papers and sends it to the Royal Society
Bayes’ theorem is the idea of recalculating how likely a cause or hypothesis is after seeing a result (evidence). It is the foundation of naive Bayes classification. On December 23, 1763, a paper was read at the Royal Society in the form of a letter from Price to the member John Canton. It was , which Price had found among the papers of his late friend Bayes, and it was printed in the society’s journal Philosophical Transactions, volume 53 (from page 370). According to Price’s letter, Bayes’s aim was to find a method for judging the probability that an event will happen, when all we know is how many times it has happened and how many times it has failed, and nothing else.
The popular story:People tend to say “Bayes published Bayes’ theorem in 1763” The original source says:Bayes himself was already dead by then (the title says “the late”). It was his friend Price who found the essay among his papers, sent it to the Royal Society and brought it to the public.
Naive Bayes classifies documents — in 1961, Maron runs an experiment guessing the field of a paper from clue words by probability. He himself wrote that the assumption that words are independent is not correct
Naive Bayes classification is a method that, from the clue words contained in a document, uses Bayes' theorem to calculate the probability that the document belongs to each field, and puts it in the field with the highest probability. It assumes that once the field is fixed, the words appear independently of one another, so the calculation can be done by simply multiplying the probabilities of the individual words. In his 1961 paper , M. E. Maron (The RAND Corporation) ran an experiment using 405 abstracts from 1959 issues of the computer journal IRE Transactions on Electronic Computers. The 260 abstracts from the March and June issues served as training data, used to count the relationship between words and fields. The 145 abstracts from the September issue were left completely unseen during preparation and were used later as test data. There were 32 fields and 90 clue words. Among the 85 test abstracts that contained two or more clue words, the top-ranked field was correct for 44 (51.8%), and the correct field was within the top three for 68 (80%).
The popular story:People tend to think “Naive Bayes is a method that naively believed that words are unrelated to each other” The original source says:in his 1961 paper Maron himself wrote that the assumption that words are independent of each other is clearly not correct, but that he makes it to make the calculation easier (at some cost in accuracy). It was a shortcut chosen with full knowledge. Neither the word naive nor the name Bayes appears in the paper; he calls this calculation reasoning backward from evidence to hypothesis.
What we still don’t know:This 1961 experiment separated training data from test data by leaving the 145 test abstracts completely unseen during preparation. The paper's numbers also show that the accuracy on new documents (51.8%) is much lower than the accuracy on the documents used for training (84.6% with the top-ranked field correct).
What the figure shows:Multiplying word probabilities for each subject area can identify a document's subject from clue words. It works reasonably well even with the false assumption that the words are unrelated to one another.
This is the same mechanism as "Try it", written as Python code. Press ▶ Run and the result appears below. Change the numbers and run it again.
u09_naive_bayes.py
031805 · Least squares
Square the errors, add them, minimize
When measured points scatter, which line “fits” best?
In 1805, Adrien-Marie Legendre in Paris published the method of least squares in an appendix to a book on calculating the orbits of comets. When there are more equations than unknowns, you cannot make every error zero. So he wrote that the most general, accurate and easy approach is to square each error, add them up, and choose the values that make that total smallest. He also explained that this balances the errors against each other and keeps extreme errors from dragging the result around.
In the figure on the right, the errors are drawn as real squares. Move the points or change the slope of the line (the intercept adjusts automatically) and try to make the total area smaller. Press "Least squares (to the valley bottom)" and the line found from the formula in one step is drawn. You fit the numbers in the formula to the data. Many of the methods in this chapter start from this form.
S=Σ(yi − (a xi +b))² made as small as possible a=Σ(xi − x̄)(yi − ȳ) ÷ Σ(xi − x̄)²、b= ȳ − a x̄
(xi, yi)
the observed data points
a, b
slope and intercept of the line
x̄, ȳ
means of x and y
What you think you knowThe popular story:Least squares was devised by Gauss The original source says:Legendre was the first to print and publish it (1805). In his 1809 book Gauss wrote that he had used it since 1795, and Stephen Stigler (1981) concluded that Gauss was probably first but that it is not conclusive. Robert Adrain in America also published it in 1808–1809.
The most famous priority dispute
Stigler (1981) calls this the most famous priority dispute over a discovery in the history of statistics. In later recollections Gauss sometimes says 1794, so even the year wavers.
Sources & confidence (1)Square the errors, add them, minimize
1805–1809Verified in the originalSources disagree
Least squares: Legendre announces it in the appendix of a book in 1805, and Gauss claims in 1809 that he had used it since 1795 (the most famous priority dispute in the history of statistics)
The method of least squares is a way of choosing the constants of a formula (for a straight line, the slope and intercept) so that, when you square each error between the measured values and the values calculated from the formula and add them all up, the total is as small as possible. In 1805 Legendre (A. M. Legendre) published a book in Paris on determining the orbits of comets, . In its appendix (on the method of least squares) he wrote that when there are more equations than unknowns the errors cannot all be made zero, and that making the sum of the squares of the errors as small as possible is the most general, most accurate and most practical choice. He also noted that this balances the errors against one another and guards against the effect of extreme errors. In Section 186 of his book on the motion of the heavenly bodies, , published in Hamburg in 1809, Gauss (C. F. Gauss) wrote that the principle of minimizing the sum of the squares of the differences was one he and others had used since 1795, and that Legendre had recently published it too. Stigler, who re-examined the history in 1981, calls this the most famous priority dispute in the history of statistics.
The popular story:People tend to say “Gauss came up with least squares” The original source says:The one who published it in print first was Legendre (in the appendix of his 1805 book). Gauss wrote in his 1809 book that he had used it since 1795. But Stigler, who re-examined the evidence in 1981, concluded that Gauss probably did discover it first but did not succeed in communicating it to the people of his time, and he adds that even this is not conclusive. In addition, the American Adrain also published the same method at the end of 1808 or the beginning of 1809, so three people are named, not two.
What we still don’t know:Whether Gauss really used it in 1795 is not settled. In Gauss’s later testimony the year is sometimes given as 1794, so the date is uncertain. Stigler writes that it might be possible to check using data from the French meridian survey.
3 sources and how they were checked
Where sources disagree:Year of publication of Legendre’s book: the cover says AN XIII − 1805. Gauss cites it in his 1809 book as Paris 1806 (archive.org also has an 1806 edition with a supplement). The year Gauss began using the method: 1795 in his 1809 book, and, according to Stigler (1981), sometimes 1794 in his later testimony.
What the figure shows:The straight line that fits scattered points is chosen to minimize the sum of the squared differences.
Sum of squared errors as the slope a changes (the bottom of the valley is the least squares)
This is the same mechanism as "Try it", written as Python code. Press ▶ Run and the result appears below. Change the numbers and run it again.
u09_least_squares.py
041838, 1958 · The S-shaped curve
From population formula to probability
How do you build a prediction that always stays between 0 and 1?
In 1838, Pierre François Verhulst in Belgium wrote down how a population grows as a formula. It takes the doubling push mp and subtracts a brake np² that gets stronger as the population grows. The solution curve is S-shaped, and the population approaches an upper limit m/n. In 1845 he named this curve logistique (he did not give a reason).
The same S-shaped curve later came to be used to predict probabilities. For data whose answer is “yes (1)” or “no (0)”, the method that predicts the probability of a 1 from explanatory variables is logistic regression. A straight line can predict below 0 or above 1, but an S-shaped curve always stays between 0 and 1. In 1958, David Cox of the University of London published a paper on regression for 0-and-1 data.
In the figure on the right, in "Population (1838)" move m and n, and in "Yes/No (1958)" press "Learn one step", and watch the S-curve fit itself to the points.
dp/dt=mp − np² P(y=1)=1 ÷(1+e−(ax+b))
p, m, n
the population, the growth push, the strength of the brake
a, b
numbers set by learning
e
about 2.718 (the base of natural logarithms)
What you think you knowThe popular story:Logistic regression was devised by Cox in 1958 The original source says:In a 1994 interview Cox himself said he did not know who was first and that Cornfield might have done it. In 1944, Berkson had used the logistic function and coined the word .
The shape of the brake cannot be fixed
Verhulst compared his formula with population data from France, Belgium, Essex and Russia. He wrote that even if the brake is changed to np³ or n log p, the error against the observations is almost the same, yet the value of the upper limit changes a lot. With little data you cannot settle which formula is correct. This problem leads to the spread "How do you spot overfitting?" (当てはまりすぎを、どう見抜くか), on AIC and cross-validation.
Sources & confidence (2)From population formula to probability
1838Verified in the original
The logistic function (S-shaped curve): in 1838 Verhulst publishes it as an equation for how a population grows; it gets its name in 1845
The logistic function is an equation that grows at nearly a doubling pace at first but slows as it nears an upper limit, drawing an S shape. In later logistic regression, this S shape is used to express a probability (a number between 0 and 1). In 1838, in a short paper in volume 10 of the journal , edited by Quetelet, Verhulst (P.-F. Verhulst) expressed the speed at which a population p grows by the equation dp/dt = mp − np². It is the doubling term mp minus a brake np² that acts more strongly the larger the population is. From the solution of this equation, the population approaches an upper limit m/n as time passes. He compared it with population data for France, Belgium, Essex and Russia, and wrote that this work was completed in 1833 and that he published it on Quetelet’s advice. He also noted that when the form of the brake was changed to np³, np⁴ or n log p, the error against the observations was almost the same, but the value of the upper limit changed greatly.
The popular story:People tend to think “logistic is a word for classification” The original source says:It was originally an equation for how populations grow. The 1838 paper did not yet have a name. Verhulst named it logistique in a 1845 paper of the Royal Academy of Belgium, and gave no reason (J. S. Cramer’s 2002 historical study). In 1920 Pearl and Reed, in a study of the US population, rediscovered the same curve without knowing Verhulst’s work.
What we still don’t know:Verhulst himself pointed out that even when the form of the brake was changed in various ways, the error against the observations stayed about the same, while the forecast of the upper limit changed greatly. He concluded that only the future can reveal the true form of the brake. The problem that a small amount of data cannot settle which equation is correct leads to the later discussion of model selection (AIC and so on).
Logistic regression — in 1958 D. R. Cox publishes a regression for 0-and-1 data; he himself said he did not know who did it first
Logistic regression is a method for data where the answer is 1 or 0, such as yes or no, and it predicts the probability of getting 1 from explanatory variables such as age. Ordinary linear regression can predict values below 0 or above 1, but with the S-shaped logistic curve the prediction always stays between 0 and 1, and the result can be read as an odds ratio (a ratio of how likely something is to happen). In his July 1958 paper , D. R. Cox (affiliation given in the paper: Birkbeck College, University of London) dealt with tests and estimation for the case where a sequence of 0s and 1s is observed and the probability of getting a 1 is thought to depend on one or more explanatory variables. The abstract says there is a fair amount of earlier work, which the paper also reviews. In a 1994 interview, Cox said his aim was to build a systematic framework, similar to least-squares regression, for analyzing 0-and-1 data, and that one thing that prompted it was the strong interest of experimental psychologists at Birkbeck College in analyzing 0-and-1 data.
The popular story:People tend to say “Cox invented logistic regression in 1958” The original source says:in a 1994 interview Cox himself said he did not know who first did logistic regression, and that Cornfield might have done it in the Framingham study a few years before 1958. In 1944, Berkson had used the logistic function in a bioassay (an experiment that measures the effect of a drug or similar by using living organisms) and coined the term (Cramer's 2002 historical study). A 2022 medical journal article writes that logistic regression was created by Berkson (1944) and polished by Cox (1958).
What we still don’t know:Who first used logistic regression is not clear even from the people involved. Cox named Cornfield, but also said he does not remember whether he was influenced by him.
What the figure shows:An S-shaped curve describing population growth later became a formula for predicting yes-or-no probabilities between 0 and 1.
This is the same mechanism as "Try it", written as Python code. Press ▶ Run and the result appears below. Change the numbers and run it again.
u09_logistic.py
051901, 1933 · Principal component analysis
Which line is closest to the points?
When both x and y have measurement error, how do you decide on a “fitting line”?
The least squares line assumes that only y has error and minimizes the vertical gaps. So you get a different line depending on whether you take x or y as the base. In 1901, Karl Pearson in London reasoned that it is enough to minimize the sum of the squares of the perpendicular lengths dropped from each point to the line. That line passes through the center of the points, and points along the longest axis of the ellipse that describes how the points scatter.
In the figure on the right, choose "Pearson's 10 points" to see the 10 points exactly as in the numerical example from Pearson's paper. Choose "Perpendiculars", rotate the line and look for the direction in which the average of the squared perpendiculars is smallest. The paper's answer is a slope of −0.546. Regression that predicts y from x gives −0.5396, and predicting x from y gives −0.5659, so the three lines are all different.
In 1933, Harold Hotelling of Columbia University named the method of taking out, one after another, the components that explain the most variation from the scores of many related tests .
The direction that minimizes Σ di² =the direction in which the points scatter most (the first principal component)
di
length of the perpendicular from the i-th point to the line
scatter ellipse
an ellipse that shows how the data scatter
What you think you knowThe popular story:Principal component analysis was devised by Hotelling in 1933 The original source says:That the line closest to the points is the longest axis of the scatter had already been shown by Pearson in 1901.
Why components and not factors
Hotelling's starting point was psychology, the study of pulling out the behind test scores. He wrote that he chose the name components so that it could be used outside psychology (for the properties of telephone poles, currency exchange rates and so on), and so that it would not collide with the mathematical word factor.
Sources & confidence (2)Which line is closest to the points?
November 1901Verified in the original
The prototype of principal component analysis: in 1901 Pearson shows that the line closest to a cloud of points runs along the longest axis of the spread of the points
Karl Pearson of University College, London dealt with the problem of finding the line or plane that best fits a cloud of points in a plane or in a space of three or more dimensions. In ordinary least squares (regression analysis, which assumes error only in y and makes the vertical gaps small), the line comes out different depending on whether x or y is taken as the base. Pearson said that when x and y both contain measurement error equally, one should minimize the sum of the squares of the lengths of the perpendiculars dropped from the points to the line. He then showed that this line passes through the center of gravity and runs along the longest axis of the ellipsoid that describes how the points are scattered. This corresponds to the direction of the first principal component in modern principal component analysis (a method of expressing many sets of numbers by a few numbers taken in order from the direction with the largest variation). The paper has a numerical example with 10 points: the slope of the closest line is −0.546, that of the regression line predicting y from x is −0.5396, and that predicting x from y is −0.5659, showing that the three lines differ.
The popular story:People tend to say “Hotelling devised principal component analysis in 1933 to summarize many variables in a few numbers” The original source says:In 1901 Pearson had already shown that the line closest to the points runs along the longest axis of their spread. But Pearson’s aim was to find the best-fitting line or plane when both x and y contain error, and his examples had only two and three variables. He also wrote that with four, five or more variables the calculation becomes troublesome (because one has to solve equations of the fourth, fifth or higher degree), but that it is not so hard to get a rough idea of the solution. The name principal components comes from Hotelling in 1933.
Principal component analysis: in 1933 Hotelling names , a way to express many numbers, such as test scores, by a few components taken in order of how much they vary
Harold Hotelling of Columbia University considered that behind many variables that are related to one another, such as children’s scores on arithmetic-speed and reading tests, there are fewer basic quantities that are unrelated to each other, and he showed a method for extracting them. First take the one component that explains the most of the variation in all the variables, then take, in turn, from the components unrelated to it, the one that best explains the remaining variation. Components with small effects can be thrown away. He named this , and explained that it is the same as aligning the coordinate axes with the principal axes of the ellipsoid of the distribution of points. What psychology called can also be used in fields outside psychology (such as the properties of utility poles and exchange rates), and he writes that he chose to call them components so as not to clash with the mathematical word factor.
The popular story:People tend to think “principal component analysis is a technique of the computer age for shrinking data to fewer dimensions” The original source says:Hotelling’s 1933 paper was psychology research on extracting the basic mental workings (in the terms of the time, factors) behind test scores, carried out under Thorndike’s committee (the Unitary Traits Committee) and the Carnegie Foundation. The name components was chosen so that it could be used outside psychology and not be confused with the mathematical factor. Note that Pearson had already shown in 1901 that the line closest to data points runs in the direction of greatest variance.
What we still don’t know:We have not checked on this site whether Hotelling knew Pearson’s 1901 paper (the copy we saved is machine-read text, and we have not read all the pages).
Whence Principal Components? (Lawrence Hubert · 2016 · p.3-6 Harold Hotelling and Truman Lee Kelley/Hotelling as a Quant… · Secondary summary)
What the figure shows:The line closest to the points, minimizing the sum of squared perpendicular distances, follows the longest direction of their spread. It differs from a regression line that measures only differences in y.
The direction of the line and the average squared perpendicular
This is the same mechanism as "Try it", written as Python code. Press ▶ Run and the result appears below. Change the numbers and run it again.
u09_pca.py
061950 · Japan · Quantification
Deciding the scores from the data
What score should we give to answers like “fits in well, ordinary, rebellious”?
Around the end of Japan's 1948 survey of literacy, Chikio Hayashi of the Institute of Statistical Mathematics was approached by a scholar of criminal law who wanted statistics to make the decision of whether to grant parole more reliable. Until then, people had decided scores for each answer, such as −2 to 2, and added them up.
Hayashi went to Yokohama Prison and read a trend: both model prisoners who fit in well and people who were openly rebellious were likely to reoffend, and people in the middle were least likely to reoffend. So he arrived at a method in which the numbers given to the answers themselves are decided from the data, so that the group that did not reoffend and the group that did are separated as clearly as possible (an English-language paper in 1950).
In the figure on the right, compare scores decided by a person with scores decided from the data (the headcounts are an example made for this site).
Choose the scores s so that the correlation ratio η² is as large as possible
s
the number given to an answer (category)
η²
the share of the overall variation that is explained by differences between groups (the form of today's Quantification Type II)
What you think you knowThe popular story:Quantification Types I to IV were created and named by Hayashi in that order The original source says:The first to be completed was the parole prediction method (today's Type II). The numbers I to IV were names given in 1964 by the social psychologist Hiroshi Akuto, and have nothing to do with the order in which they were made.
The consultation and the paper
The consultation was about deciding whether to grant parole to inmates who had served a third of their sentence. Hayashi heard that model prisoners who fit well into prison life and work hard come right back after release. Scores do not always line up in one straight direction. The 1950 English paper writes this as a problem of maximizing the hit rate of the prediction.
Quantification afterwards
After that Hayashi built methods for each problem in the field: a literacy survey (Type I), a survey of human relations among Japanese National Railways track maintenance workers (Type IV), and a design survey for labels on canned goods for export (Type III). Type III (1956) is deeply connected mathematically to Jean-Paul Benzécri's correspondence analysis (1973), and was published earlier.
How much Hayashi's numbers actually changed real parole decisions has not yet been confirmed.
Sources & confidence (1)Deciding the scores from the data
1950–1956Verified in the originalJapan
Quantification theory: in 1950 Chikio Hayashi writes up a method of giving numbers to qualitative data in order to predict parole
Chikio Hayashi of the Institute of Statistical Mathematics, as the 1948 survey of Japanese literacy was nearing its end, was consulted by a criminal law scholar at the Ministry of Justice training institute for corrections and parole. The request was to make the judgment of whether prisoners who had served a third of their sentence could be paroled statistically more reliable. The earlier method assigned scores such as −2, −1, 0, 1, 2 to the answers on each item by human judgment and added them up. Hayashi visited Yokohama Prison, where he heard that model prisoners who adapted well to prison life and worked hard would come back soon after release. He also read that very rebellious people were likely to reoffend, and that people in the middle were the least likely to. He saw that evaluation scores do not always line up in one simple direction. So he got the idea of deciding from the data the very numbers to assign to each answer (qualitative data), so that the group that did not reoffend and the group that did are separated most clearly. In his English paper of 1950 he set this up as a problem of maximizing the success rate of prediction. Hayashi then developed methods for field problems one by one, such as the literacy survey, a survey of human relations among track maintenance workers of Japanese National Railways, and a survey on the label design of canned goods for export, and these came to be called quantification theory together.
The popular story:People tend to think “Quantification Type I, II, III and IV are methods Chikio Hayashi made and named in that order” The original source says:The first to be made was the method for parole prediction (today’s Type II), while Type I came from the literacy survey, Type IV from the survey of track maintenance workers of Japanese National Railways, and Type III from the survey on label design for canned goods for export. The numbers Type I to Type IV are names given in 1964 by the social psychologist Hiroshi Akuto, and there is no direct relation between the order of the numbers and the order in which they were developed (Morimoto 2005). Morimoto 2005 also cites correspondence analysis by Benzécri of France (1973) as work mathematically closely related to Type III. Hayashi’s Type III (1956) is earlier than that.
What we still don’t know:We have not confirmed on this site how much the numbers (weights) Hayashi made actually changed parole decisions (whether they were used in the field).
What the figure shows:Letting data determine the scores assigned to answers, instead of choosing them by hand, captures relationships that do not follow a straight line, such as people in the middle range being least likely to reoffend.
This is the same mechanism as "Try it", written as Python code. Press ▶ Run and the result appears below. Change the numbers and run it again.
u09_quantification.py
071951, 1967 · k-nearest neighbors
Just copy the closest example
If you only copy the closest example you remember, how far can you get?
The k-nearest-neighbor method finds the k closest data points to a new point you want to classify, and assigns it to the class most of them belong to. It assumes nothing in advance about the shape of the distribution of the data. In February 1951, Fix and Hodges considered this kind of rule in a report from the U.S. Air Force School of Aviation Medicine. They showed that if k grows as the number of data points n grows (while k/n goes to 0), the rule approaches the best possible classification.
In 1967, Cover and Peter Hart of Stanford proved that even the nearest-neighbor rule (k=1), which copies just the single closest point, has with unlimited data an error rate that does not exceed twice the theoretical minimum (the Bayes error rate). In this sense, the paper says, half of the classification information in infinite data is contained in the one nearest point.
Place points in the figure on the right and change k. With k=1 every training point is classified correctly, but small islands form around the noisy points.
R* ≤ R ≤ R*(2 − M R* ÷(M − 1))≤ 2R*
R
error rate of the nearest-neighbor rule (with unlimited data)
R*
Bayes error rate (the error that remains even if the distribution is known perfectly)
M
number of classes
What you think you knowThe popular story:The k-nearest-neighbor method was devised by Cover and Hart in 1967 The original source says:The two themselves wrote that Fix and Hodges seem to have been the first to formulate the neighbor rule and study its properties. The 1951 report was never turned into a journal paper.
The condition for the guarantee
This guarantee of "at most twice" is about the case of unlimited data, and does not guarantee performance when data are few.
This site has not yet checked the 1951 report itself directly (it relies on the 1967 paper and on later summaries).
Sources & confidence (2)Just copy the closest example
February 1951Confirmed in summaries
The beginning of k-nearest neighbors: in 1951 Fix and Hodges devise it in a US Air Force report (it was never made into a journal paper)
The k-nearest neighbors method (k-NN) finds the k labeled data points closest to a new point you want to classify, and assigns the point to the group that is the majority among them. Because it does not fix the distribution of the data in advance, it is called a nonparametric method. In February 1951, E. Fix and J. L. Hodges, Jr. devised this kind of rule in a report of the USAF School of Aviation Medicine (Randolph Field, Texas), . They showed that if k is increased as the number of data points n increases, and k/n is made to approach 0, the rule is consistent (the more data there is, the closer it comes to the optimal classification). This is how the 1967 paper of Cover and Hart introduces it.
The popular story:People tend to say “k-nearest neighbors was devised by Cover and Hart in 1967” The original source says:Cover and Hart themselves write in their 1967 paper that it seems Fix and Hodges and others were the first to formulate the nearest neighbor rule and study its properties. That 1951 report was never published as a journal paper (Scholarpedia’s explanation).
What we still don’t know:We have not yet opened and checked the original 1951 report on this site. The content is based on the introductions in the 1967 paper and in an explanatory article.
A theoretical guarantee for the nearest neighbor method — in 1967, Cover and Hart prove that even a classifier that just matches the single nearest point has an error rate no more than twice the theoretical minimum
The nearest neighbor method (the k-nearest-neighbor method with k=1) is a technique that puts a new data point into the same class (group) as the nearest already-classified point. In their January 1967 paper , T. M. Cover (Stanford University) and P. E. Hart (Stanford Research Institute) showed that when there are enough classified points, the error rate R of this method cannot be smaller than R*, the minimum error rate when the distribution is perfectly known (the Bayes error rate), but does not exceed twice R* (with M the number of classes, R* ≤ R ≤ R*(2 − MR*/(M − 1))). The paper writes that, in this sense, half of the classification information in infinite data is contained in the single nearest point.
What we still don’t know:This theoretical guarantee is about the case where there is plenty of data (in the paper's words, a large sample), and it says nothing about performance when data is scarce.
What the figure shows:Simply copying the nearest example can give the right answer. Setting k to 1 gets every training point right, but creates small islands around noisy points.
R* and the upper bound on R (for two classes)
This is the same mechanism as "Try it", written as Python code. Press ▶ Run and the result appears below. Change the numbers and run it again.
u09_knn.py
081957, 1967 · k-means
Split, then move to the average
Without giving the right answers, can we split data into k groups?
k-means is a procedure that alternates two steps: (1) assign each point to the nearest representative, and (2) move each representative to the average of its group. It keeps shrinking the spread within the groups. With each repetition, the total spread either decreases or stays the same and stops.
Many people arrived at the same kind of procedure independently. In 1957, Lloyd of Bell Labs, working on sending a signal's voltage by replacing it with a fixed number of representative values (PCM), wrote up a procedure for finding the boundaries and representatives that give the smallest error, in a manuscript circulated inside the lab (it appeared in a journal about 25 years later). In 1967, James MacQueen of the University of California, Los Angeles gave the name k-means to a procedure that puts points one at a time into the nearest group and immediately updates the averages.
In the figure on the right, press "Assign" and "Move" in turn. If you change where the first representatives are placed, it may stop at a different answer.
W=Σgroup Σpoint x in the group ‖x − μgroup‖²
x, μ
a point, and the representative (average) of its group
‖ ‖²
squared distance
W
total spread within the groups
What you think you knowThe popular story:k-means was invented by MacQueen in 1967 The original source says:In the same paper he himself lists earlier examples: Sebestyen (1962), Forgy and Jennrich, and Wallace and Ball. Steinhaus in 1956 and Lloyd in 1957 had also put forward the same idea.
A paper 25 years later
Lloyd's manuscript appeared in a journal in 1982. He himself called it .
MacQueen's k-means and today's k-means
What MacQueen named was the sequential procedure that adds points one at a time. The two-stage loop you often see in textbooks today (reassign all the points, then recompute the averages) is introduced in the same paper as a method Forgy and Jennrich taught him around 1963 (as summarized by Bock, 2008).
Sources & confidence (2)Split, then move to the average
31 July 1957–1982Verified in the original
The seed of k-means — in 1957, Lloyd at Bell Laboratories writes down a procedure for finding the best way to represent a signal with a small set of representative values (it appeared in a journal only in 1982)
Lloyd of Bell Laboratories thought about PCM (pulse-code modulation: a way of sending a signal by replacing its voltage with one of a fixed number of representative values). The problem is to find the boundaries and representative values that make the mean of the squared error from the original signal (quantization noise) as small as possible. The procedure he called Method I alternates two steps: (1) fix the boundaries and take the centroid (average) of each interval as its representative value, (2) make the exact midpoint between neighboring representative values the new boundary. This has the same form as one-dimensional k-means (a procedure that repeats two steps: sort points into the group of the nearest representative, then move each representative to the average of its group). The manuscript was circulated inside Bell Laboratories in 1957, and part of it was presented at a meeting of a statistics society that September, but it appeared in a journal only about 25 years later, in 1982.
The popular story:People tend to say “k-means is a method MacQueen devised in 1967” The original source says:a procedure of the same form, which repeats sorting points to the nearest representative and moving the representative to the average, is already in the manuscript Lloyd circulated inside Bell Laboratories in 1957 (for quantizing one-dimensional signals). In computer science and pattern recognition it is often called Lloyd's algorithm (Bock 2008). However, the manuscript went unpublished in journals for 25 years, and in a note added to his 1982 paper Lloyd himself writes that it is an , a printing of the 1957 draft almost unchanged.
What we still don’t know:The 1957 manuscript itself (an internal Bell Laboratories record) has not been made public, and this site has only checked it through the note added to the 1982 paper. How much was presented at the 1957 meeting is not known.
k-means — in 1967, MacQueen names a procedure for splitting points into k groups k-means (k means)
J. MacQueen of the University of California, Los Angeles gave the name k-means to a procedure for dividing data (points) into k groups so that the scatter within the groups (the sum of squared distances from the mean) is small. In his procedure, the first k points each become a group of their own. After that, new points are added one at a time to the group whose mean is nearest, and that group's mean is updated right away. Because at every stage the k means are the means of the groups, he called it k-means. In the same paper MacQueen also describes the form that is often presented today as the k-means procedure (repeat two steps: reassign all points to the group of the nearest mean → recompute each group's mean). He says this is a method that Forgy and Jennrich each noticed separately and told him about around 1963.
The popular story:People tend to say “k-means was invented by MacQueen in 1967” The original source says:in the same paper, MacQueen himself writes as follows. As a clear use of the same idea, Sebestyen's 1962 book is the earliest he knows of. He heard of the two-step repeating method from Forgy and Jennrich. In a note added at the proof stage, he adds that Wallace and Ball had also used it separately. According to Bock's (2008) summary, what MacQueen named k-means was the sequential procedure of adding points one at a time, not the two-step repeating form commonly seen in textbooks today. Earlier still, in 1956 Steinhaus presented the same idea as a mechanics problem of dividing an object into k parts, and in 1957 Lloyd had written a procedure of the same form for quantizing signals.
What we still don’t know:It is known that several researchers independently arrived at the same procedure, but this site has not confirmed how much information passed between them (for example, who had read Lloyd's 1957 manuscript).
What the figure shows:Alternating between assigning points to their nearest representative and moving each representative to the average divides data into groups without being taught answers. Different starting points can lead to different final results.
This is the same mechanism as "Try it", written as Python code. Press ▶ Run and the result appears below. Change the numbers and run it again.
u09_kmeans.py
091958, 1963 · Hierarchical clustering
Link the most similar into a tree
How can you find clusters in data without deciding the number of groups?
In 1958, the entomologists Sokal and Michener of the University of Kansas compared 97 species of bees using 122 characters per species. They first grouped the species that resemble each other, then kept linking the groups they had made, and summarized the whole as a tree-shaped diagram (a dendrogram).
In 1963, Ward of an Air Force base laboratory turned this way of summarizing into a general procedure. Check every way of choosing two of the existing groups, and merge the two that make a number measuring how good the grouping is (the objective function) best. Repeat until everything is one group.
The figure on the right is exactly the example in Ward's paper. If the scores of 10 people, 2, 6, 5, 6, 2, 2, 2, 0, 0, 0, are put in one group, the sum of squared deviations from the average is 50.5. Split into four groups, {0,0,0} {2,2,2,2} {5} {6,6}, it is 0. Press "Merge two" and watch the tree grow.
The two groups to merge=the two groups whose merger increases ESS the least (the objective function in this figure)
ESS
the sum, over the groups, of the sum of squared deviations from each group's average
What you think you knowThe popular story:Ward's method merges the pair of groups with the smallest increase in the total squared deviation The original source says:The 1963 paper wrote it as a general procedure in which the objective function can be any function the researcher chooses to fit the purpose. The sum of squared deviations was just one of the examples.
The people who computed
The 1958 calculations used IBM punched-card machines and people calculating on desktop calculators, and the correlation coefficients alone took 160 person-hours. The paper uses the word computer for those people. The two authors also wrote that a digital computer of the day would take less than an hour.
People tend to think clustering was born in AI or information science, but one of its early forms was made by entomologists to make the classification of bees less dependent on human taste.
Sources & confidence (2)Link the most similar into a tree
20 March 1958Verified in the original
Hierarchical clustering — in 1958, Sokal and Michener group 97 species of bees into a tree diagram by similarity
Sokal and Michener of the Department of Entomology at the University of Kansas classified 97 species of solitary bees of the leafcutter bee family. They compared 122 characters per species (traits such as whether the abdomen is red or black) and expressed how alike species were with a correlation coefficient (a number showing how closely two lists of numbers rise and fall together). They then grouped similar species and joined the resulting groups step by step, summarizing everything in a tree diagram. This is the idea behind today's hierarchical clustering (a method that merges the closest items first and builds a tree of nested groups). Earlier classification gave weights to a few characters by human judgment, but the two men treated many characters with equal weight. For the calculations they used IBM punch-card machines and people operating desk calculators (the paper calls these human calculators “computer”), and the correlation coefficients alone took 160 person-hours. The two also wrote that a digital computer of the time could do it in under an hour.
The popular story:People tend to think “clustering (a way of grouping data into similar items) was born in AI and information science” The original source says:one of the early forms of hierarchical clustering was devised in 1958 by entomologists, so that the classification of bees would depend less on human preference. The calculations were done with punch-card machines and by people using desk calculators.
What we still don’t know:The idea of merging similar items step by step to build a tree is said to have appeared in other fields around the same time, but this site has not opened those original papers and has not confirmed it. Which came first has not been decided.
Ward's method — in 1963, Ward presents hierarchical clustering that picks an objective function and merges two groups at a time
Joe H. Ward Jr. of the Aerospace Medical Division at Lackland Air Force Base presented a procedure for dividing people or objects into similar groups (clusters) in large surveys of more than 100. Starting with n groups (at first, one item each), it examines all n(n−1)/2 ways of joining two of them, and merges the two for which the objective function (a number showing how good the grouping is; the analyst chooses it to suit the purpose) comes out best, leaving n−1 groups. Repeating this until only one group remains shows the whole hierarchical structure and the amount of information lost at each stage. In the paper's example, if the scores of 10 people, 2, 6, 5, 6, 2, 2, 2, 0, 0, 0, are represented as one group by their mean, the error sum of squares (the total of the squared deviations from the mean) is 50.5, but if they are divided into the four groups {0,0,0} {2,2,2,2} {5} {6,6}, it is 0.
The popular story:People tend to be taught that “Ward's method merges the pair of groups for which the increase in within-group scatter (the error sum of squares) is smallest” The original source says:Ward's 1963 paper presents a general procedure in which the objective function can be any function the analyst chooses to suit the purpose. The error sum of squares is only one example used in the explanation. At the paper's Personnel Research Laboratory, the same computer program was used with a different objective function for each of four purposes: grouping people who are alike, grouping jobs so that retraining time for reassignment is lower, reducing the number of job descriptions, and reducing the number of prediction equations.
What we still don’t know:This site checked the first two pages of the paper (p.236-237). It has not read the details of the flow chart and numerical example (p.238 onward).
What the figure shows:Joining the most similar items first brings everything into 1 branching tree, or dendrogram, without choosing the number of groups in advance.
This is the same mechanism as "Try it", written as Python code. Press ▶ Run and the result appears below. Change the numbers and run it again.
u09_ward.py
101966, 1975 · Hidden Markov models
Guessing hidden states from symbols
How do you recognize a signal like a voice, which wobbles every time even for the same word?
In 1966, Baum and Petrie of the Institute for Defense Analyses in the United States built a mathematical model in which states that cannot be seen from outside change with certain probabilities, and each state outputs a symbol with certain probabilities. From the sequence of output symbols alone, you can estimate the probabilities of the state changes and of the symbols being produced. This is today's hidden Markov model (HMM).
In 1975, James Baker used this mechanism in the speech understanding system DRAGON, and in 1976 Frederick Jelinek of IBM brought it into continuous speech recognition. It became widely known in 1989 through a tutorial by Rabiner of AT&T Bell Laboratories.
In the figure on the right, hidden states (voice, silence) change from one to another, and all we can observe is the loudness of the sound. From this we estimate the most probable sequence of states (the Viterbi algorithm).
Look for the sequence of states that maximizes P(states | symbols)
states and symbols
the states cannot be seen directly (voice, silence), only the symbols can be observed (loudness)
Viterbi
a dynamic programming method that keeps only the most probable route so far as it moves forward (the same idea as DP matching in Chapter 5)
What you think you knowThe popular story:With the arrival of hidden Markov models, the performance of speech recognition jumped The original source says:In Rabiner's (1989) experiments on words spoken one digit at a time, the error rates of the HMM and of the earlier pattern matching (DTW) were about the same (1.3% versus 2.0%, and 1.8% versus 1.1%). The strengths lay in its power to represent how sounds are produced, and in its ability to extend to recognition with large vocabularies.
A journal that speech engineers did not read
Rabiner (1989) wrote that the theory had appeared in specialist mathematics journals that speech engineers did not read, which delayed its spread into speech recognition. From the time it began to be used for speech in the 1970s, it took more than ten years to become firmly established.
The title of the 1966 paper did not say "hidden"; it was about probabilistic functions of Markov chains.
Sources & confidence (2)Guessing hidden states from symbols
December 1966–1970Verified in the original
The hidden Markov model — in 1966, Baum and Petrie work out the mathematics of a model in which unseen states change and symbols come out of them
Leonard E. Baum and Petrie of the Institute for Defense Analyses in Princeton considered a certain probabilistic model. In it, states that cannot be seen from outside move as a Markov chain (a sequence of changes in which the next state is determined probabilistically by the current state alone), and each state emits symbols probabilistically. They called this a probabilistic function of a Markov chain. They then proved that from the observed symbol sequence Y1, …, YT alone, one can estimate the table A of state transition probabilities and the table B of symbol emission probabilities, and that as data increases, maximum likelihood estimation (choosing the values that make the data most likely) approaches the correct values. This is one of the mathematical beginnings of today's hidden Markov model (HMM). In 1970 the same group (Baum with Petrie, Soules and Weiss) presented an iterative calculation that raises the likelihood step by step. Dempster and colleagues wrote in 1977 that this was a special case of the EM algorithm.
The popular story:People tend to think “the hidden Markov model (HMM) was a method made for speech recognition” The original source says:it began as a 1966 statistics paper by two researchers at the Institute for Defense Analyses. Its title spoke of a probabilistic function of a Markov chain, not of “hidden” anything. Rabiner (1989) writes that because this theory appeared in mathematics journals that speech engineers did not read, its application to speech recognition was delayed. According to a footnote in Rabiner 1989, the way of explaining the HMM through three basic problems (computing probabilities, the most likely sequence of states, and learning the parameters) was also presented in lectures and writings by Jack Ferguson of the IDA (written there as Institute for Defense Analysis).
What we still don’t know:This site has not confirmed who first used the name hidden Markov model, or in what year. The sources also do not say for what purpose this research was carried out at the Institute for Defense Analyses.
February 1975–1989Verified in the originalSources disagree
Speech recognition with hidden Markov models — Baker began using them in 1975 and Jelinek of IBM in 1976, and Rabiner's 1989 tutorial made them widely known
Even for the same word, speech varies with the person and the speaking speed. In 1975, J. Baker, in the speech understanding system DRAGON, expressed various kinds of linguistic knowledge in a single form, a probabilistic function of a Markov process (today's hidden Markov model: a model in which hidden states change with probabilities, and each state outputs features of sound with probabilities). In 1976, F. Jelinek of the IBM Thomas J. Watson Research Center published a paper on recognizing continuous speech by statistical methods. However, it was only in the last few years, as of 1988 when Rabiner wrote his manuscript, that this method came to be widely understood and used in speech research. In 1989, Lawrence R. Rabiner of AT&T Bell Laboratories pulled together the theory and its application to speech recognition in a single tutorial. In the experiments in that tutorial, for speaker-independent recognition of digits spoken one word at a time, the error rate of the hidden Markov model was about the same as that of conventional pattern matching (a method called DTW, which lines up differences in speaking speed with dynamic programming and then compares with learned sound patterns). This idea of pattern matching had also been published in Japan in 1971 by Hiroaki Sakoe and Seibi Chiba of NEC as DP matching (another event on this site).
The popular story:People tend to think “when hidden Markov models came in, speech recognition suddenly got much better” The original source says:in the table of Rabiner's own 1989 experiments, for speaker-independent recognition of digits spoken one word at a time, the error rates of the hidden Markov model and conventional pattern matching (DTW) were about the same (1.3% versus 2.0% on a new group of 100 speakers, and 1.8% versus 1.1% on another group of 100). The strengths Rabiner pointed to were the ability to express how speech is produced and the fact that it can be extended to recognition with a large vocabulary. Also, it took more than 10 years from when it began to be applied to speech recognition in the 1970s until it came into wide use.
What we still don’t know:Rabiner 1989 also points out that the assumption that sound observations are independent of each other and the Markov assumption do not fit real speech in some ways. Why it worked well in spite of that is not clearly explained in the tutorial.
3 sources and how they were checked
Where sources disagree:Baker's affiliation: Rabiner 1989 writes Baker at CMU, while the IEEE bibliographic record (for the 1975 paper) writes IBM Thomas J. Watson Research Center (the affiliation at the time of publication may differ from where the work was done, but we have not checked this).
What the figure shows:Hidden states, such as speaking or silence, change according to probabilities, while only sound volume is visible. Using connections over time identifies the sequence of states more accurately than guessing each time step separately.
This is the same mechanism as "Try it", written as Python code. Press ▶ Run and the result appears below. Change the numbers and run it again.
u09_hmm.py
111971–2010 · Japan · Choosing a model
How do you spot overfitting?
The more numbers you adjust, the better the fit to the data. Where should you stop?
The more numbers you can adjust (parameters), the more closely a formula fits the data in hand. But it tends to miss on new data (overfitting). In 1971, Hirotugu Akaike of the Institute of Statistical Mathematics published AIC, a criterion that adds a penalty for the number of parameters to the fit to the data. Smaller is better.
In 1974, Stone in London wrote up cross-validation, which splits the data into one part for learning and one part for checking, and chooses a model by the average score as the split is varied. In 1977 he showed that leave-one-out cross-validation and AIC behave the same way when there is a lot of data. In 2010, Sumio Watanabe of the Tokyo Institute of Technology created WAIC, which can be used even for neural networks and other models where AIC's assumptions do not hold, and proved that it approaches the same value as leave-one-out cross-validation.
Raise the degree of the polynomial in the figure on the right. The training error keeps falling, but the cross-validation error and AIC start to rise partway.
AIC= −2 log L̂ +2k
L̂
maximum likelihood (how likely the data are at the values that make them most likely)
k
number of parameters
What you think you knowThe popular story:The A in AIC is the A of Akaike The original source says:In his 1974 paper he wrote that he chose A so that similar criteria such as BIC and DIC could follow, and in a 1981 recollection he wrote that it stands for "an information criterion". The W in WAIC also stands for widely applicable in Watanabe's paper.
The first announcement and how it was received
According to his 1981 recollection, Akaike suddenly saw the connection on the morning of March 16, 1971, when he sat down on his commuter train. In his 2006 Kyoto Prize lecture, Akaike said (translated), and that the fields of application responded eagerly, but (translated).
The idea of cross-validation is old: Stone himself lists earlier examples such as Larson in 1931 and Mosteller and Tukey in 1968.
Sources & confidence (3)How do you spot overfitting?
1973–1974Verified in the originalSources disagreeJapan
The Akaike Information Criterion (AIC) — presented by Hirotugu Akaike in 1971 and put into a paper in 1973; it combines fit to the data and a penalty on the number of parameters in one
AIC is a criterion for deciding which to use among several candidate statistical models. It is calculated as AIC = (−2) × (maximum log-likelihood) + 2 × (number of parameters), and the model with the smaller value is judged better. The first term shows how well the model fits the data (the better the fit, the smaller the value), and the second term acts as a penalty for increasing the number of adjustable parameters too much. In 1968, Hirotugu Akaike of the Institute of Statistical Mathematics, working on a time-series model of the behavior of a cement rotary kiln, faced the problem of how many past values to use (the order), and came up with the expected prediction error (FPE). According to his 1981 recollection, on the morning of March 16, 1971, as he sat down on his commuter train, he suddenly noticed the link between maximum likelihood and the Kullback–Leibler information. In his 2006 Kyoto Prize commemorative lecture he said, (translated). He went on to recall that he presented it in detail the same year at the Second International Symposium on Information Theory held in Armenia, that it came out in the proceedings in 1973, and that in 1974 he published an explanation in a US IEEE journal on automatic control.
The popular story:People tend to say “the A in AIC is the A of Akaike” The original source says:Akaike himself wrote in his 1974 paper that IC stands for information criterion, and that the A was added so that similar criteria such as BIC and DIC could follow later. In his 1981 recollection, too, he states that AIC stands for “an information criterion.” On the other hand, the 2006 Kyoto Prize citation writes Akaike Information Criterion.
What we still don’t know:AIC was not accepted from the start. In his 2006 lecture, Akaike said that there was an eager response from applied fields, but (translated). In his 1981 recollection he also writes that it was still not fully accepted by professional statisticians. Also, AIC assumes a regular model, in which parameters and distributions correspond one to one. So how it can be used for learning machines such as neural networks led to later research (such as Sumio Watanabe's WAIC).
5 sources and how they were checked
Where sources disagree:The year it was proposed: Akaike himself writes that the first announcement was in 1971 (at the Japan Statistical Society and at an international symposium the same year) and that the proceedings were published in 1973. Genshiro Kitagawa's 2006 text says it was proposed in 1973, and the Kyoto Prize citation says the early 1970s.
Cross-validation — in 1974, Stone's paper on choosing and assessing prediction methods by splitting the data and checking. The idea itself dates from the 1930s
Cross-validation is a method in which you split the data you have into data for learning and data for checking, repeat the check while changing the split, and choose a prediction method from the average of the results. In the leave-one-out version, from n data points you remove one, build a model from the remaining n − 1, and test it on the removed one, repeating this in n ways. M. Stone (University College London) presented at the Royal Statistical Society on December 5, 1973, and in the paper , published in the Journal of the Royal Statistical Society, Series B in 1974, proposed using this criterion in a form that combines the choice of a prediction method with its assessment. He then showed it with examples of estimation from a single sample, multiple regression and analysis of variance.
The popular story:People tend to say “cross-validation was devised by Stone in 1974” The original source says:Stone himself notes at the start of the paper that this way of thinking about assessment existed before. He gave as earlier examples Larson's 1931 work that split the data into two, a 1951 symposium on the need for and methods of cross-validation, and the 1968 explanation by Mosteller and Tukey, who said to build a model from all but one, test on the one left out, and repeat this over all the data. In a note, Stone also wrote that validation (recognizing something as correct) sounds overconfident, so he prefers assessment.
What we still don’t know:In 1977, Stone showed that leave-one-out cross-validation and Akaike's AIC behave in the same way when there is a lot of data (they are asymptotically equivalent) (as introduced in Watanabe's 2010 paper). However, this is for regular models in which parameters and probability distributions correspond one to one. What happens with singular models such as neural networks was not known until Watanabe's 2010 paper.
WAIC — in 2010, Sumio Watanabe proves that an information criterion usable for neural networks and more, and leave-one-out cross-validation, approach the same value
WAIC (widely applicable information criterion) is a criterion devised by Sumio Watanabe (affiliation in the paper: Precision and Intelligence Laboratory, Tokyo Institute of Technology) for estimating the generalization loss (the prediction error on unseen data) of a Bayesian model. AIC assumes regular models, in which parameters and probability distributions correspond one to one. But many learning machines, such as neural networks, mixtures of normal distributions, hidden Markov models and Boltzmann machines, are singular models for which this assumption does not hold. In a December 2010 paper in the Journal of Machine Learning Research, Watanabe proved that even for such singular models, the loss of Bayesian leave-one-out cross-validation and WAIC approach the same value when the amount of data is large (they are asymptotically equivalent). For regular models, it was known, for example from Stone's 1977 paper, that leave-one-out cross-validation is asymptotically equivalent to AIC.
The popular story:People tend to say “the W in WAIC stands for Watanabe” (there is a name, Watanabe-Akaike information criterion) The original source says:In Watanabe's paper the name is widely applicable information criterion. It was Gelman and colleagues' 2013 overview paper that called it Watanabe-Akaike, and even that paper notes that Watanabe himself calls it widely applicable. Just as the A in AIC was not meant as the initial of Akaike, the W in WAIC was not a person's name for the one who named it.
What the figure shows:Adding more terms to a formula fits the available data better, but makes errors on new data more likely. AIC and cross-validation help find the right balance.
Training error, cross-validation error and AIC by degree
This is the same mechanism as "Try it", written as Python code. Press ▶ Run and the result appears below. Change the numbers and run it again.
u09_aic_cv.py
121977 · The EM algorithm
Fill in the unseen part and recalculate
When part of the data is hidden, how do you find the most plausible values?
In 1977, Arthur Dempster, Nan Laird and Donald Rubin of Harvard University and elsewhere presented a general method for finding plausible parameters when part of the data cannot be seen. It alternates an E step (expectation), which computes the expected value of the unseen part using the current values, with an M step (maximization), which recomputes the values from the filled-in data. With each repetition, the likelihood of the data does not go down.
The first example is data in which 197 animals split into 4 types by a genetic mechanism, as 125, 18, 20 and 34. In the figure on the right, choose "197-animal example" and press "E: Fill in" and "M: Recalculate" in turn. π settles from 0.5 to 0.6268214… in 8 steps, the same numbers as in Table 1 of the paper. In "Two peaks", mixed points are separated into their peaks.
the value we want (a probability inside the genetic model)
x₂
the expected number, among the 125 animals, that belong to the unseen second group
What you think you knowThe popular story:The EM algorithm was invented by Dempster and colleagues in 1977 The original source says:The three themselves wrote that it had been proposed many times before for special cases (Hartley 1958, Baum and colleagues, and others). Their contribution was to put it into a general form, give it a name, and prove that the likelihood does not go down.
Deliberately creating the unseen
In the 197-animal example, the first type (125 animals) is thought of as two groups, one with probability 1/2 and one with π/4. How many animals are in each is not visible, but splitting it up on purpose makes the calculation simple.
Sources & confidence (1)Fill in the unseen part and recalculate
8 December 1976–1977Verified in the original
The EM algorithm: in 1977, Dempster, Laird, and Rubin gathered a method for learning from data with hidden parts into one general form and gave it a name
A. P. Dempster, N. M. Laird, and D. B. Rubin, of Harvard University and the Educational Testing Service, presented a general method for estimating the most likely parameters (the numbers inside a model) when part of the data cannot be seen (is missing). In the E step (expectation), the current parameters are used to calculate the expected values of the unseen parts. In the M step (maximization), the filled-in data are used to recalculate the parameters to their most likely values. Because the two steps are repeated in turn, they named it the EM algorithm. They proved that with every repetition the likelihood of the data (how likely it is that these parameters would produce this data) never goes down. They also showed that many examples fit this form: data with missing values, distributions made by mixing several peaks (mixture distributions), factor analysis, and more. Their first example was data on 197 animals that fall into 4 types through a mechanism of heredity (an example from Rao's 1965 book).
The popular story:The EM algorithm is often said to have been invented in 1977 by Dempster, Laird, and Rubin The original source says:The three authors themselves wrote that the EM algorithm had been proposed many times before for special cases. They cited Hartley's example (1958) and the papers of Baum and colleagues (1967–1972) as separate derivations of the same method (Baum and colleagues' work was on a model where an unseen state changes as a Markov chain, which is today's hidden Markov model). They also wrote that Baum and colleagues had shown the same idea for a special case but seem not to have noticed how widely it applies. The three authors' contribution was to gather the scattered examples into one general form, name it EM, and prove in general that the likelihood never goes down.
What the figure shows:Alternating between filling in hidden parts using the current values (E) and recalculating values from the filled-in data (M) moves toward the most likely values.
π at each repetition
This is the same mechanism as "Try it", written as Python code. Press ▶ Run and the result appears below. Change the numbers and run it again.
u09_em.py
131979 · Japan · Otsu's binarization
Setting the white–black border by itself
When you split a grayscale image into white and black, how do you decide the brightness of the border?
To separate characters or objects from the background, you classify pixels into two groups at a certain brightness (the threshold). In 1979, Nobuyuki Otsu of the Electrotechnical Laboratory published, in an English-language paper, a method for deciding this threshold automatically. When the frequency distribution of brightness (the histogram) is split in two at a threshold k, choose the k that maximizes the variation between the groups (between-class variance). This is the same as minimizing the variation inside each group.
The calculation uses only the share of pixels up to that brightness and the sum of their brightness, so it is very simple. It is an unsupervised method that decides without using examples of correct answers. The paper gave as examples the letter A typed with a new ribbon and an old ribbon, a cloth pattern, and an image of cells.
In the figure on the right, move the threshold and look for the top of the peak in the graph (the image is an example made for this site).
σB²(k)=ω₀ ω₁(μ₀ − μ₁)² → choose the k with the largest σB²
ω₀, ω₁
share of pixels in the dark group and in the bright group
μ₀, μ₁
average brightness of each group
What you think you knowThe popular story:Otsu's binarization is a small trick for turning images into black and white The original source says:Otsu derived it from discriminant analysis and treated it as one of the basic problems of pattern recognition. The paper's conclusion also says that, as long as you have a histogram, it can be used for other unsupervised classification too.
The English paper came first
The English paper (1979) was published before the Japanese paper (1980). The abstract of the Japanese paper says it can easily be extended to classification into three or more groups (a cell image was split into three: background, cytoplasm and nucleus).
The English paper notes that a rigorous proof that the evaluation function always has only one peak has not yet been obtained.
Sources & confidence (1)Setting the white–black border by itself
January 1979–1980Verified in the originalJapan
Otsu's binarization: in 1979, Nobuyuki Otsu presented a method that automatically sets the dividing line between white and black in a grayscale image at the point where the two groups are best separated
Nobuyuki Otsu (his affiliation in the 1979 English paper was the Electrotechnical Laboratory, Tokyo) showed a method for automatically deciding the brightness boundary (the threshold) used to separate characters or objects from the background in a black-and-white grayscale image. When the histogram of brightness is split into two groups at a brightness k, the method picks the k at which the spread between the groups (the between-class variance) is largest. This comes from the idea of discriminant analysis (a statistical method for finding a yardstick that separates two groups well), and it is equivalent to making the spread within the groups (the within-class variance) smallest. The calculation uses only the zeroth- and first-order cumulative moments of the histogram (the fraction of pixels up to that brightness, and the sum of their brightness), so it is very easy to compute. It is an unsupervised method: it decides without being given the correct answers. The abstract of the 1980 Japanese paper says the result is also optimal in the least-squares sense when the original grayscale image is approximated by a few values, and that it extends easily to splitting into three or more groups. The paper's examples were the letter A typed with a new ribbon and with an old ribbon, a cloth pattern, and a cell image (split into three: background, cytoplasm, and nucleus).
The popular story:Otsu's binarization is often thought of as a small image-processing trick for turning pictures black and white The original source says:Otsu derived the choice of threshold from discriminant analysis and treated it as one of the basic problems of pattern recognition, not just image processing (the summary of the Japanese paper). The conclusion of the English paper also says it can be used for other unsupervised classification as long as there is a histogram of the feature that separates the objects. In other words, it was proposed as an unsupervised learning method that splits data into two groups without being given correct answers. Note that the English paper (1979) came out before the Japanese paper (1980).
What we still don’t know:The English paper says that a rigorous proof that the evaluation function always has a single peak (is unimodal) has not yet been obtained (in experiments it was always smooth and unimodal).
What the figure shows:The brightness threshold for separating black and white can be chosen automatically, without labeled examples, by finding where the 2 groups are most separated, measured by between-class variance.
Threshold k and the between-class variance σB²
This is the same mechanism as "Try it", written as Python code. Press ▶ Run and the result appears below. Change the numbers and run it again.
u09_otsu.py
141985 · Bayesian networks
Evidence spreads through the net
Can you reason from many clues without bending the rules of probability?
Many expert systems of the early 1980s (Chapter 6) used ad hoc calculation rules for uncertain reasoning. In 1982, Judea Pearl of the University of California wrote that many AI researchers had accepted the myth that the rules of probability must be bent, at least in part. He devised a graph in which matters are points, direct dependencies are arrows, and their strength is expressed by conditional probabilities. He showed that with this, points that exchange short messages only with their neighbors can update the probabilities of the whole correctly. In 1985 he called this form a Bayesian network.
The figure on the right is a small network in which "Rain" and "Sprinkler" are causes of "Grass is wet" (the probabilities are an example made for this site). Enter "The grass was wet" and both rain and the sprinkler become suspect, but if you then also enter "It rained", the probability of the sprinkler goes down. That is because one cause already explains it.
P(rain | wet)=P(wet | rain) × P(rain) ÷ P(wet)
P(A | B)
probability of A when B is known
arrow
direction from a direct cause to its effect
What you think you knowThe popular story:Bayesian networks were born as a statistical tool for learning from data The original source says:The title of Pearl's 1985 report was about a model of a self-activated memory that reasons from evidence, and it was presented at a cognitive science conference. The 1982 paper was written as an answer to the ad hoc calculations of expert systems.
The origin of the name and its limits
Pearl wrote that he called it a Bayes Network to stress that the probability values come from human judgment (he also called the same thing an Influence Network).
Pearl himself states in the abstract that correct updating of the whole by exchanges between neighbors alone works only for networks where there is just one path between any two points.
Sources & confidence (1)Evidence spreads through the net
April 1985Verified in the original
Bayesian networks: in 1985, Pearl named a method that shows links between causes and effects as an arrow graph with conditional probabilities (the idea dates from 1982)
Judea Pearl of the University of California, Los Angeles (UCLA) called Bayesian networks the graphs in which many propositions or variables are points (nodes), direct dependencies are arrows, and strengths are given by conditional probabilities (the probability that B happens when A holds). When new evidence (for example, a test result) arrives, each node exchanges short messages only with its neighbors, and the probabilities of the whole network (how strongly we believe things) are updated correctly. This makes it possible to reason from many clues while keeping to the rules of probability. In a 1982 AAAI paper, Pearl wrote that many AI researchers had accepted the myth that computational models of reasoning under uncertainty must bend the rules of probability, at least in part, and were using ad hoc rules. He then showed that for a tree-shaped network, the whole network can be updated according to the rules of probability with a single round of message passing. The 1985 report presented this mechanism at a cognitive science conference as a model of human memory and reasoning.
The popular story:Bayesian networks are often said to have been born as a statistical tool for learning from data The original source says:The title of Pearl's 1985 report was about a model of autonomous memory for reasoning from evidence, and it was presented at a cognitive science conference. The 1982 paper was written as an answer to the fact that AI systems such as expert systems used ad hoc calculations that bent the rules of probability (the references include Shortliffe and Buchanan 1975, a model of uncertain reasoning in medicine, and PROSPECTOR, an expert system for mineral exploration). About the name, too, he writes that he calls it a Bayes Network not so much because it uses Bayes' theorem but to stress that the probability numbers come from human judgment (he also calls the same thing an Influence Network).
What we still don’t know:In the abstract, Pearl himself writes that updating the whole network correctly by local message passing alone works only for networks in which there is just one path between any two nodes. For networks with multiple paths, this report only discussed the approach.
What the figure shows:Connecting causes and effects with arrows lets changes in likelihood spread through the network when evidence arrives. Once one cause explains the evidence, the other becomes less suspect.
This is the same mechanism as "Try it", written as Python code. Press ▶ Run and the result appears below. Change the numbers and run it again.
u09_bayesnet.py
151986 · Decision trees
Which question do you ask first?
How do you build a "tree" that gives an answer by lining up questions?
A decision tree lines up questions such as "What is the weather?" and "Is the humidity high?" in a branching shape, follows the branches according to the answers, and gives its answer at the leaf it reaches. In 1986, Quinlan of Australia wrote up a way of building one called ID3. It picks first the question with the largest amount of information (information gain, measured in bits) obtained by splitting the data, and repeats the same thing after each split.
The 14 cards in the figure on the right each represent a row of Table 1 of the paper. They show the row number and the answer P or N. Open "Details: 14 source rows, full tree, amount of information" to also read the weather, temperature, humidity and wind on a Saturday morning. Splitting by weather gives 0.246 bits of information, the largest, so the root of the tree is the weather. Press the questions and grow the tree yourself.
Information gain=I(whole) − Σ(share of each split group × I(that group)) I= −p log₂ p − n log₂ n
p, n
shares of P and N in the group
I
the amount of information needed to learn the answer (in bits)
What you think you knowThe popular story:Decision trees were born from AI (machine learning) research The original source says:The statistics side has its own separate lines: AID by Morgan and Sonquist in 1963 and CART by Breiman and colleagues in 1984. CART does not appear in the references of the 1986 paper, and the two lines grew separately. The paper has an example of the same form as the textbook "weather and tennis" example, but the word tennis does not appear in it.
From chess endgames
According to the paper, the ancestor of this family (the ID3 relatives) is CLS by Hunt and colleagues (a 1966 book). ID3 (1979) was built for a task posed by Donald Michie: recognizing chess endgame positions (king and rook against king and knight) from features of their shape alone.
Sources & confidence (1)Which question do you ask first?
March 1986Verified in the originalSources disagree
Decision trees: in 1986, Quinlan wrote up ID3 in a paper. Its sources are two streams, Hunt and colleagues' CLS (1966) and statistical trees dating from 1963
A decision tree is a classification method that arranges questions, such as what the weather is and whether the humidity is high, in a branching form. You follow the branches according to the answers, and when you reach a leaf, it outputs the result (the class). ID3 builds the tree by choosing which question to ask first as follows. It chooses the question that gives the greatest amount of information (information gain, measured in bits) when the data is split by it, and repeats the same process on each resulting branch. J. R. Quinlan (New South Wales Institute of Technology) worked through a concrete example in his 1986 paper . Using the weather conditions of 14 Saturday mornings (outlook, temperature, humidity and wind, and two classes, P and N), he showed that splitting by outlook gives 0.246 bits of information, the largest, so the root of the tree is outlook. According to the paper, the roots of this family of methods are Hunt and colleagues' CLS (a 1966 book). ID3 (1979) was developed while tackling a task posed by Donald Michie: identifying positions in a chess endgame (king and rook against king and knight) from features of their shape alone.
The popular story:It is often said that decision trees came out of AI (machine learning) research The original source says:There is another stream on the statistics side. The AID of Morgan and Sonquist in 1963 is regarded as the first regression tree to appear in the literature, and in 1984 Breiman and colleagues' CART regained interest with a method that grows a large tree and then prunes it by cross-validation (explanation by Loh, 2014). The reference list of Quinlan's 1986 paper does not include CART, and the two streams grew separately. Also, this paper has data of the same form as the weather-and-tennis example often seen in textbooks, but its classes are only P and N, and the word tennis does not appear.
3 sources and how they were checked
Where sources disagree:Publication year of the CART book: Loh (2014) cites Breiman et al. (1984), and Friedman's 1999 manuscript cites Breiman, Friedman, Olshen and Stone 1983 (this site has not opened the book itself).
What the figure shows:A decision tree is built by first asking the question that gives the most information about the answer, then repeating the same process at each branch.
This is the same mechanism as "Try it", written as Python code. Press ▶ Run and the result appears below. Change the numbers and run it again.
u09_id3.py
161988, 1989, 1994 · Reinforcement learning
Learning from rewards
Can a machine learn when it is never told the right answer and only now and then gets a "reward"?
Reinforcement learning is machine learning that picks up how to act by trial and error, using rewards as its clue. In 1988, Richard Sutton of GTE Laboratories put into a theory the TD method, which corrects today's prediction by its gap from tomorrow's prediction, without waiting for the final result. In 1989, Christopher Watkins presented in his doctoral thesis Q-learning, which corrects the value Q(x, a) of taking action a in state x a little at a time with each trial. He also proved that if every action is tried many times, the values approach the best ones.
In 1994, Kazuteru Miyazaki, Masayuki Yamamura and Shigenobu Kobayashi of the Tokyo Institute of Technology studied the profit sharing method, which, when a reward comes, strengthens the whole chain of actions that came before it. They proved the conditions under which it does not learn actions that only make a detour, and won a paper award from the Japanese Society for Artificial Intelligence.
In the maze in the figure on the right, train the route to the treasure (the reward) with Q-learning and with profit sharing.
Q(x, a) ← Q(x, a)+α(r+γ maxb Q(x′, b) − Q(x, a))
r, x′
the reward received and the next state
α, γ
the learning rate and the discount rate for future rewards (0–1)
What you think you knowThe popular story:TD learning was invented by Sutton in 1988 The original source says:Sutton himself wrote that the TD method had already been used, for example in Samuel's checkers program (1959, Chapter 3), but was not well understood. Q-learning, too, was born in the 1989 doctoral thesis, not in the 1992 paper.
From a question about animals
The starting point of Watkins' thesis was not games or robots but a question from behavioral ecology: how animals learn optimal behavior from experience. He wrote that the idea is very simple, and that, since the field has been studied for so long, he could not believe nobody had thought of it.
Miyazaki and colleagues' paper compares them: Q-learning needs an enormous number of trials to learn the environment correctly, whereas profit sharing needs far fewer. At the time it was not a situation where only one method was the right answer.
Sources & confidence (3)Learning from rewards
1988Verified in the original
TD learning: in 1988, Sutton puts together as a theory a way to learn from the error against the next prediction without waiting for the final result
Richard Sutton of GTE Laboratories put together, as a theory, the TD method (temporal-difference learning), a way to learn predictions of the future from experience. Ordinary supervised learning corrects a prediction using the error between the prediction and the final outcome. The TD method, on the other hand, corrects today's prediction using the error between “today's prediction” and “tomorrow's prediction.” The paper gives a weather forecaster as an example. Suppose you predict every day whether it will rain on Saturday, and you predicted 50% on Monday and 75% on Tuesday. The TD method raises the prediction for a day like Monday without waiting for Saturday's result. Sutton argued that this way, the computation can be done a little at a time each day, less memory is needed, and in many cases the predictions are more accurate. He then proved convergence and optimality for special cases. As a test example, he had the method learn, in a random-walk problem that starts at the center D, moves one step at random left or right, and ends on reaching A at the left end or G at the right end, the probability of reaching G from each point.
The popular story:It is often said that TD learning was invented by Sutton in 1988 The original source says:In the abstract and introduction, Sutton himself says that the TD method had already been used in Samuel's checkers program (1959), Holland's bucket brigade and Sutton's own Adaptive Heuristic Critic, but was not well understood. He gives Samuel's checkers as the earliest and best-known use, and writes that it corrected the evaluation of the earlier position by the difference between the evaluations of two successive positions. What was new in the 1988 paper was that it pulled out the TD method as a way of learning predictions, proved convergence and optimality, and compared it with supervised learning.
What we still don’t know:Sutton argued that most of the problems on which supervised learning was being used at the time were in fact prediction problems where the TD method can be used to advantage (his outlook in 1988). This site has not checked how far this outlook turned out to be true.
Q-learning: in 1989, Watkins presents it in a doctoral thesis on how animals learn optimal behavior (the detailed proof came in 1992)
Christopher Watkins, in his King's College doctoral thesis Learning from Delayed Rewards (May 1989), posed the question of how an animal can learn optimal behavior from experience. In behavioral ecology, the optimal behavior of animals is derived by dynamic programming (a method that looks ahead through all later results to compute the best action), but it is hard to think that animals do such a calculation. So Watkins presented one-step Q-learning, which gradually updates the action value Q(x, a) of taking action a in state x, trial by trial, using the reward r obtained and the largest Q value in the next state. In an appendix he proved that, even without a model of the environment, if every action is tried many times in every state, the values approach those of the optimal behavior. In 1992, Watkins and Dayan published in a journal a short paper that wrote out this proof in detail.
The popular story:It is often said that Q-learning was born in the 1992 paper by Watkins and Dayan The original source says:Q-learning itself is in Watkins's doctoral thesis of May 1989, and the 1992 paper was a short paper (Technical Note) that wrote out in detail the convergence proof outlined in the thesis. Also, the starting point of the thesis was not games or robots but the question from behavioral ecology of how animals learn optimal behavior from experience. Watkins himself writes that the idea is very simple and, as far as he knows, had not been proposed before, but that since it is a long-studied field he cannot believe nobody thought of it.
What we still don’t know:According to the preface of the thesis, Sutton pointed out a difficult point of Q-learning. This site has not checked what that was.
A theory of profit sharing: in 1994, Kazuteru Miyazaki, Masayuki Yamamura and Shigenobu Kobayashi of the Tokyo Institute of Technology proved the condition under which reinforcement learning that hands out rewards in a lump does not learn useless actions, and won the Japanese Society for Artificial Intelligence's paper award
Reinforcement learning (a method of learning how to act by trial and error, guided by rewards) includes TD learning and Q-learning, which update the value of an action at every step, and the profit sharing method, which, when a reward is obtained, reinforces the whole series of earlier actions (the episode) together. Kazuteru Miyazaki, Masayuki Yamamura and Shigenobu Kobayashi mathematically analyzed, for the profit sharing method, how much of the reward to hand out to the action taken how many steps before it (the reinforcement function). As a result, they proved the necessary and sufficient condition for not learning ineffective rules (actions) that cause detours (the ineffective rule suppression theorem), and also showed that under the same condition it can learn an action policy that keeps earning rewards. The paper compares this with Q-learning, which needs a huge number of trials to understand the environment accurately, and says that profit sharing needs far fewer trials to learn. This paper won the Japanese Society for Artificial Intelligence's paper award for fiscal 1994.
The popular story:“Research on reinforcement learning follows only the overseas line of TD learning and Q-learning,” people tend to think The original source says:In 1990s Japan, a theory of the profit sharing method, which reinforces a whole series of actions together when a reward comes, was being studied, and a 1994 paper won the Japanese Society for Artificial Intelligence's paper award. The paper itself compares them, saying that Q-learning needs many trials while profit sharing needs far fewer to learn, so it was not a situation in which only one method was the right answer.
What we still don’t know:This site has not confirmed how far this condition was later used in real robots and applications.
What the figure shows:Even without correct answers and with a reward only on reaching the treasure, repeated trials let the system learn a path to it.
Number of tries and steps taken to reach the treasure
This is the same mechanism as "Try it", written as Python code. Press ▶ Run and the result appears below. Change the numbers and run it again.
u09_q_learning.py
171990, 1995, 1999 · Boosting
Stacking weak learners on mistakes
Can you gather methods only slightly better than guessing and make a strong one?
It started as a theoretical question. Is weak learning, only slightly better than guessing, as powerful as strong learning, which can be made as accurate as you like? In 1990, Robert Schapire of MIT proved the two are equivalent and gave an algorithm that turns a weak method into a strong one. In 1995, Yoav Freund and Schapire published AdaBoost. It gives more weight to the data that was wrong in the previous round, has the next classifier learn from it, and finally combines the verdicts with a weighted majority vote.
In 1999, Friedman of Stanford generalized gradient boosting, in which the next small decision tree predicts the error (residual) between the current prediction and the right answer and the results are added up bit by bit, as an optimization that goes downhill.
In the figure on the right, press "Add one round" and the wrongly classified data grow larger, and the next boundary comes to classify them.
αt = ½ log((1 − εt) ÷ εt) weight of a misclassified example ← weight × eαt
εt
the weighted error rate of the t-th classifier
αt
that classifier's weight in the majority vote
What you think you knowThe popular story:Boosting does not overfit The original source says:In a 2013 commentary, Schapire himself showed an example where the test error keeps falling after the training error reaches 0 (8.4% at round 5 → 3.1% at round 1000), and yet wrote that boosting can overfit too.
A horse-racing analogy
The name AdaBoost comes from adaptively changing the weights to match the mistakes. The 1997 paper begins with an analogy: a gambler who keeps losing at the horse races redistributes his bets according to his friends' records.
Friedman himself wrote that the view of boosting as a downhill calculation had already been pointed out by Rätsch and colleagues (1998) and by Breiman (1999).
Sources & confidence (2)Stacking weak learners on mistakes
June 1990–1997Verified in the original
Boosting: in 1990, Schapire proves that weak learning can be turned into strong learning, and in 1995 Freund and Schapire announce AdaBoost
Boosting is a method that builds many weak classifiers that are only slightly better than guessing, trains each round with larger weights on the data that the previous round got wrong, and in the end combines them by a weighted majority vote into a strong classifier. Its origin was a question in computational theory. Kearns and Valiant left it as an open problem (the hypothesis boosting problem) whether weak learning, which is only slightly better than guessing, and strong learning, which can be made as accurate as you like, have the same power. R. E. Schapire (MIT Laboratory for Computer Science) proved in his 1990 paper that the two are the same, and showed a construction that turns a method of weak learning into a method of strong learning. Y. Freund and Schapire announced AdaBoost at a conference in Barcelona in March 1995 (the paper was published in an academic journal in 1997). They wrote in the paper that the name comes from the fact that, unlike earlier methods, it changes the weights adaptively according to the errors of the weak classifiers. The paper begins with an analogy of a person who keeps losing at horse racing and reallocates his bets according to the track records of his friends.
The popular story:Boosting is often thought to have come from a practical trick of gathering weak models together The original source says:The start was a question in computational theory. Kearns argued that simply repeating the same learning and taking a majority vote would not work. It is also often said that boosting does not overfit. But Schapire himself, in a 2013 explanatory article, shows an example in which the test error keeps falling after the training error reaches 0 (from 8.4% at round 5 to 3.1% at round 1000), while writing that boosting too can overfit.
What we still don’t know:Why does the test error keep falling even after the training error reaches 0? Schapire's 2013 explanation writes that the explanation using margins is the only theory that can account for the cases where overfitting does not occur. He also writes, though, that the numerical guarantees are very loose.
Gradient boosting: in 1999, Friedman generalized boosting, which fills in residuals little by little with small decision trees, as gradient descent
Gradient boosting is a method that improves predictions little by little by repeating the following many times: have the next small decision tree predict the error (the residual) between the current model's prediction and the correct answer, and add that tree a little at a time. Friedman (Jerome H. Friedman, Stanford University) wrote the paper dated February 24, 1999 (IMS 1999 Reitz Lecture, published in The Annals of Statistics in 2001). In it he regarded this piling up of additions as a computation that goes down a little at a time in the steepest direction, in function space, of the size of the error (the loss) (steepest descent). When the loss is squared error, the tree at each round simply learns the current residual. The same procedure can be applied to other losses, such as absolute error, the Huber loss, and the logistic likelihood for classification. He also writes that making the amount added at each round small with a learning rate ν (shrinkage) often gives better results.
The popular story:“Gradient boosting is a trick from recent machine-learning contests,” people tend to think The original source says:The basic paper is a manuscript dated February 1999, presented as the IMS 1999 Reitz Lecture and published in The Annals of Statistics in 2001. Friedman himself writes that the view of boosting as gradient descent had already been pointed out in earlier machine-learning work (Rätsch and colleagues in 1998, Breiman in 1999).
What the figure shows:Even classifiers only slightly better than guessing become strong when each new one learns with extra weight on previous mistakes and their answers are combined by voting.
Number of training errors by round
This is the same mechanism as "Try it", written as Python code. Press ▶ Run and the result appears below. Change the numbers and run it again.
u09_boosting.py
181992 · Support vector machines
Maximizing the margin from the border
You can draw many lines that split two groups. Which one is best?
In July 1992, Bernhard Boser, Isabelle Guyon and Vladimir Vapnik (working at AT&T Bell Laboratories) presented a new learning method at the conference COLT '92. It draws the boundary between two groups so that the distance to the nearest points (the margin) is as large as possible. The boundary is decided only by the few points closest to it, and does not change if other points are moved.
Rewrite it to use only a function that measures how similar points are to each other (a kernel), and the same calculation can maximize the margin even for a curved boundary. On handwritten digits from mail (7300 training images and 2000 test images), a fourth-degree polynomial kernel gave an error rate of 4.9%. That was slightly lower than the 5-layer neural network on the same data (5.1%).
Drag points in the figure on the right. Moving any point other than the ringed ones (the points that support the boundary) does not move the line unless it enters the margin band.
Maximize the margin ⇔ minimize ‖w‖ K(x, x′)=exp(−γ‖x − x′‖²)
w, b
direction and position of the boundary (every point satisfies yi(w·xi+b) ≥ 1, where y marks the group with+1 or −1)
K, γ
the kernel (similarity between points) and how narrow the range of influence is
What you think you knowThe popular story:The support vector machine was invented by Corinna Cortes and Vapnik in 1995 The original source says:The learning method that maximizes the margin and, with kernels, handles curved boundaries too had been presented by the three in 1992. The paper does not contain the term support vector; it calls the points that support the boundary supporting patterns. It also cites the potential function method of Aizerman and colleagues from 1964 as an example of the kernel idea.
Sources & confidence (1)Maximizing the margin from the border
July 1992Verified in the original
The support vector machine with kernels: in 1992 Boser, Guyon and Vapnik presented margin maximization in a form that also works for curved boundaries
A support vector machine (SVM) is a method that draws the boundary between two classes so that the distance (the margin) from the boundary to the nearest points is as large as possible. Boser (B. E. Boser, University of California, Berkeley; part of the work was done at AT&T Bell Laboratories), Guyon (I. M. Guyon) and Vapnik (V. N. Vapnik; both at AT&T Bell Laboratories), in their paper at the conference COLT '92 in July 1992, rewrote this problem so that it uses only a function that measures the similarity between points (the kernel) K(x, x') (the dual form). This lets the same computation maximize the margin even for curved boundaries such as polynomials or radial basis functions. The boundary is determined only by the few training examples closest to it (called supporting patterns in the paper). On handwritten digits taken from real mail (7,300 training and 2,000 test images), a 4th-degree polynomial kernel gave an error rate of 4.9%, slightly below the 5.1% reported for a 5-layer neural network on the same data.
The popular story:“The support vector machine was invented in 1995 by Cortes and Vapnik,” people tend to say The original source says:The learning method that maximizes the margin and uses kernels for curved boundaries had already been presented in 1992 at COLT by Boser, Guyon and Vapnik. That paper does not contain the term support vector; it calls the points that hold up the boundary supporting patterns. For the idea of the kernel itself, the paper cites as an example the potential function method of Aizerman, Braverman and Rozonoer from 1964.
What the figure shows:Chapter 8's widest road allowed only a straight line, but replacing the calculation with similarities between points, called a kernel, keeps the idea of maximizing the margin from the nearest points while allowing curved boundaries for groups a straight line cannot separate. Once again, only a few points nearest the boundary determine it.
This is the same mechanism as "Try it", written as Python code. Press ▶ Run and the result appears below. Change the numbers and run it again.
u09_svm.py
191995, 2001 · Random forests
A majority vote of randomly grown trees
One tree fits the training data too closely. What happens when you gather many?
If you grow a decision tree deep, it fits all the training data but memorizes even the noise. In 1995, Tin Kam Ho of AT&T Bell Laboratories showed that bundling many trees that each use a randomly chosen part of the features raises accuracy on handwritten digit recognition (random decision forests).
In 2001, Leo Breiman of the University of California, Berkeley, pulled random forests together. Trees are grown on randomly redrawn data, and at each split the best split is chosen from randomly selected features. As trees are added, the error rate settles at a certain value. With the roughly one-third of the data that was not chosen (out-of-bag), you can estimate the error rate without a separate test set.
In the figure on the right, as you increase the number of trees, the jagged border becomes smooth.
Answer of the forest=majority vote (answer of tree₁, answer of tree₂, …, answer of treeB)
B
number of trees
error of the forest
decided by the strength (accuracy) of each tree and the correlation (similarity) between trees
What you think you knowThe popular story:Random forests were thought up by Breiman alone The original source says:Random decision forests, with a very similar name, had been published by Ho in 1995. Breiman himself lists bagging (1996), Ho's method, and Amit and Geman (1997) as earlier work, and calls the combination random forests.
What Breiman wrote
He wrote that as trees are added the error rate settles at a certain value (by the strong law of large numbers), so overfitting is not a problem. Ho's trees improved accuracy on both the training data and data seen for the first time.
About one-third: if you pick n times at random from n items, allowing repeats, the share of items never picked is about (1 − 1
n)n ≒ 0.37. These are the out-of-bag data.
Sources & confidence (1)A majority vote of randomly grown trees
October 2001Verified in the original
Random forests: in 2001, Breiman combined the majority vote of many decision trees randomly varied a little. In 1995 Ho had already published random decision forests
A random forest is a method that grows many decision trees under slightly different conditions and decides the prediction by majority vote. Each tree is grown from examples chosen at random from the training data, and at each split the best way of dividing is chosen from a randomly chosen subset of the features. In his 2001 paper , Breiman (Leo Breiman, Statistics Department, University of California, Berkeley) wrote that as more trees are added, the error rate converges to a fixed value by the strong law of large numbers, so overfitting is not a problem. He also showed that the error rate of a forest is determined by the strength of the individual trees and the correlation between them. Using the roughly one third of examples that were not chosen at random (out-of-bag), the error rate can be estimated without separate test data or cross-validation. In 1995, Ho (Tin Kam Ho, AT&T Bell Laboratories) had shown in that bundling trees that each use a randomly chosen subset of the features improves accuracy in handwritten-digit recognition, both on the training data and on data seen for the first time.
The popular story:“Breiman thought up random forests all by himself,” people tend to say The original source says:Random decision forests, with a very similar name, were published by Ho in 1995. In his 2001 paper Breiman himself lists bagging (1996), methods that randomly choose the candidate splits, and Ho's random subspace method as earlier examples. He then writes that he was influenced by Amit and Geman's 1997 paper on handwritten characters, and calls these together random forests. Incidentally, the name Random Forests is also a trademark of Breiman and Cutler (Adele Cutler).
What the figure shows:A single deep tree memorizes noise too, but voting across many randomly grown trees smooths the boundary and makes the error rate more stable.
Number of trees and the error rate estimated out-of-bag
This is the same mechanism as "Try it", written as Python code. Press ▶ Run and the result appears below. Change the numbers and run it again.
u09_random_forest.py
201873, 1936 · Singular value decomposition
Any table can be broken into three
How can you represent a table of numbers (a matrix) well with few numbers?
Singular value decomposition (SVD) says that any matrix can be broken into three multiplications: rotate, stretch along each axis, and rotate once more. The stretch factors are called singular values. Eugenio Beltrami published it first in 1873, in a mathematics journal for university students in Italy. In 1874 Camille Jordan showed it independently.
If you keep only the large singular values, you can approximate the original table accurately with few numbers. In 1936, Carl Eckart and Gale Young of the University of Chicago, in connection with analyzing scores on psychological tests, solved the problem of finding the matrix that approximates best at a low rank. The answer is the one with the small singular values set to 0.
In "Circle (2×2)" in the figure on the right, a circle becomes an ellipse through "rotate, stretch, rotate". In "Picture (24×24)", a picture drawn from a 24×24 table of numbers is rebuilt keeping only k singular values (the picture is an example made for this site).
A=U Σ VT Ak =σ₁u₁v₁T + … +σkukvkT
U, V
rotations (they change direction)
Σ
stretching along each axis (a table with the singular values σ₁ ≥ σ₂ ≥ … along the diagonal)
Ak
the approximation that keeps only the largest k
What you think you knowThe popular story:The approximation theorem was found by Eckart and Young (1936) The original source says:According to Stewart (1992), Erhard Schmidt had proved it in 1907, and the two rediscovered it.
Computing it on a computer came much later
Beltrami writes at the start, in effect, that he hopes this will be a chance for young readers to become familiar with this kind of algebra. The decomposition itself is 19th-century mathematics. A method to compute it quickly and accurately on a computer appeared only in 1965 (Golub and Kahan) (Stewart 1992).
Sources & confidence (2)Any table can be broken into three
1873–1874Verified in the original
Singular value decomposition: Beltrami in 1873 and Jordan in 1874 show that any matrix can be split into three steps: rotate, stretch or shrink, rotate
Singular value decomposition (SVD) is the fact that any matrix (a table of numbers in rows and columns that expresses a transformation of vectors) can be split into three multiplications: change the direction (rotate) → stretch or shrink along each axis → rotate once more. The stretch or shrink factors are called singular values. If you keep only the large singular values, the original table can be approximated well with fewer numbers, so SVD is used for predicting tastes (recommendation features) and for computing principal component analysis. The first to publish this decomposition, in 1873, was Eugenio Beltrami, in a mathematics journal for Italian university students. At the start Beltrami wrote, in effect, that he hoped to give young readers a chance to become familiar with this kind of algebraic calculation. The next year, 1874, Camille Jordan separately presented the same decomposition in a more complete form.
The popular story:People tend to think “singular value decomposition is a new technique made for handling data with computers” The original source says:The decomposition itself is mathematical research by Beltrami in 1873 and Jordan in 1874. Beltrami’s paper appeared in a journal for university students as material for them to get used to algebra. A fast and accurate way to compute it on computers appeared much later, in 1965 (Golub and Kahan) (Stewart 1992).
What we still don’t know:According to Stewart 1992, the time when the term singular value settled into its present meaning was uncertain in the first half of the 20th century. We have not traced the original sources on this site.
Best approximation by a lower rank: in 1936 Eckart and Young solve it in connection with factor analysis (a rediscovery of Schmidt’s 1907 theorem)
Carl Eckart and Gale Young of the University of Chicago solved the problem of approximating a matrix by a matrix of lower rank (the number of mutually independent directions) so that the sum of squared errors is as small as possible. According to the abstract of the paper, the problem becomes easier to solve if the matrix is first brought to a canonical form (today’s singular value decomposition). A solution always exists and is usually unique. The solution is the singular value decomposition with the small singular values set to 0 (as explained by Stewart 1992). The abstract says the problem is deeply connected to a basic concept of factor analysis, which explains a table of test scores by a few factors. The paper appeared in volume 1 of Psychometrika, a journal of psychological measurement founded that year.
The popular story:People tend to say “Eckart and Young (1936) found the theorem that the best approximation by a lower rank is the matrix with the small singular values thrown away” (the name Eckart–Young theorem) The original source says:According to Stewart 1992, this approximation theorem had been proved in 1907 by Erhard Schmidt in his work on integral equations, and Eckart and Young rediscovered it. Stewart writes that the name is widely used but wrong.
What we still don’t know:Whether Eckart and Young knew Schmidt’s 1907 work has not been checked on this site (we have not opened the full text).
Whence Principal Components? (Lawrence Hubert · 2016 · p.1 The Journal of Educational Psychology as a Precursor to Psy… · Secondary summary)
What the figure shows:Any table can be split into 3 operations, rotate, stretch, and rotate; keeping only the largest stretches, called singular values, approximates the original table well using fewer numbers.
Size of the singular values (largest first)
This is the same mechanism as "Try it", written as Python code. Press ▶ Run and the result appears below. Change the numbers and run it again.
u09_svd.py
212006–2009 · The Netflix Prize
Learning a table full of blanks
Can you predict what score a person would give to a movie they have not seen?
In October 2006, Netflix started a contest. If you could predict the scores people give to movies 10% more accurately than the company's own recommendations, you won a $1 million prize. The data handed out were more than 100 million ratings, given by about 480,000 people to about 18,000 movies.
On December 11, 2006, a participant, Simon Funk, published his whole method on his blog. He represented people's tastes and movies' features with 40 numbers each, and predicted a rating from the sum of their products. Of the table's roughly 8.5 billion cells, 8.4 billion are blank, so a true singular value decomposition cannot be used. So he computed the error using only the ratings that are known, and corrected the values a little at a time. He himself called this a tweaked SVD.
In the figure on the right, press "Keep learning" and the blanks fill in (the table is an example made for this site).
prediction r̂ =pu · qm pu ← pu +η e qm (qm is updated the same way. e=r − r̂)
pu, qm
the tastes of person u and the features of movie m (sets of numbers)
r, η
the known rating and the learning rate
What you think you knowThe popular story:Netflix improved its recommendations with what was found in the contest The original source says:In the company's blog (2012), the winning method was not adopted because it was not worth the effort.
The details of the contest
Netflix was renting out DVDs of movies at the time. The rating data had names replaced with numbers. Funk's learning rate was 0.001. On September 21, 2009, BellKor's Pragmatic Chaos, which combined hundreds of models, won. The second-place team had the same score (0.8567), and the team that submitted about 20 minutes earlier won. Afterwards the center of the company moved to streaming, and became only one part of recommendation (the 2012 blog).
Even with names removed
It was thought that data with the names removed could not be traced to anyone. But researchers showed that individuals could be identified by matching against publicly posted movie ratings, and the second contest was cancelled.
Sources & confidence (2)Learning a table full of blanks
2 October 2006–2009Verified in the original
The $1 million contest to predict movie tastes, the “Netflix Prize” (2006–2009): the winner was decided by a 20-minute margin
In October 2006, Netflix, which then rented out movies on DVD, opened a contest. It offered a $1 million prize to anyone who could predict movie ratings 10% more accurately than its own recommendation system (Cinematch). It distributed more than 100 million ratings (1 to 5 stars) of about 18,000 movies by about 480,000 people, with names replaced by numbers. About three years later, on September 21, 2009, the team BellKor's Pragmatic Chaos won. The second-place team had the same score (error 0.8567, a 10.06% improvement), but the side that submitted about 20 minutes earlier became the winner (under the rule that a tie goes to whoever submitted first).
The popular story:“Netflix improved its recommendations with the system found in the $1 million contest,” people tend to think The original source says:According to Netflix's own blog (2012), it put two methods from the first year's progress prize into use, but the combination of hundreds of models that won the $1 million was not adopted, because the extra accuracy was not worth the effort. Meanwhile the company had moved from DVDs to streaming, and had itself become only one part of a recommendation system. Another point: people thought “Data with the names removed cannot be traced to anyone”, but researchers showed that individuals could be identified by matching against public movie ratings, and the second contest was cancelled (an FTC letter).
Factorizing the table of tastes: in 2006 Simon Funk published on his blog the method he used in the Netflix Prize, which learns from only the known entries of a table full of blanks
Simon Funk, who was taking part in the Netflix Prize, a contest to predict movie tastes, published his whole method (at that time tied for third) on his blog on December 11, 2006. The data provided was 100 million ratings by about 500,000 people of about 17,000 movies; as a user-by-movie matrix, 8.4 billion of its roughly 8.5 billion entries were missing. Funk represented each user's taste and each movie's features as 40-dimensional vectors and predicted a rating from their inner product (about 20 million parameters, 1/400 of the original matrix). This has the same form as an approximation by singular value decomposition (SVD) that keeps the top 40 singular values, but real SVD cannot be computed when values are missing. So he found the prediction error only for the known ratings, and repeatedly adjusted the parameters a little in the direction that made the error smaller (gradient descent; the learning rate was 0.001). He himself calls this a modified SVD.
The popular story:“The Netflix Prize used singular value decomposition (SVD),” it is said The original source says:What the participants called SVD was a matrix factorization that approximates a matrix with many missing values as the product of two long, thin matrices fitted to the known entries, and it is not the strict mathematical singular value decomposition itself. Funk himself explains that real SVD cannot be used with 8.4 billion blanks, and that his method computes the gradient of the error and learns from only the known ratings, and he calls his method a modified SVD. The Netflix blog also writes of matrix factorization (which participants generally called SVD).
What we still don’t know:This site has not confirmed who, before Funk, had devised a method that likewise factorizes a matrix with missing values to predict tastes.
2 sources and how they were checked
Netflix Update: Try This at Home (Simon Funk · 2006 · Body of the article (explanation of the problem, the 40 fea… · Original source)
What the figure shows:Representing people's preferences and films' features with a few numbers, then gradually adjusting them using only known ratings, makes it possible to predict blanks in a mostly empty table.
Number of repetitions and the error on the validation data
This is the same mechanism as "Try it", written as Python code. Press ▶ Run and the result appears below. Change the numbers and run it again.
u09_funk_mf.py
221770, 2005, 2009 · A person inside
Who put the correct labels on?
When "AI looked at photos and learned by itself", where did the right answers come from?
In November 2005, Amazon started Mechanical Turk, a service in which a program asks people to do small tasks and receives the results. The announcement called having computers get large numbers of people to do tasks and then using the results .
In 2009, Li and colleagues at Princeton University presented ImageNet (3.2 million images at that point), photos gathered for each cluster of nouns in WordNet, which organizes word meanings in a hierarchy. The photos were gathered with search engines, and whether the object was really in each photo was checked, one photo at a time, by people around the world who were paid for the work on Mechanical Turk, and decided by the votes of several people (page 4 of the paper). "A person inside" is not a story only about the 1770 chess machine (called the Turk).
The figure on the right shows how the voting works (the accuracy of each person is a provisional number made for this site).
Correctness of the majority vote=Σk>n/2 C(n, k) pk(1 − p)n−k
n
number of people who check one image (odd)
p
accuracy per person
C(n, k)
the number of ways to choose k people out of n
What you think you knowThe popular story:Amazon explains that it took the name from the 1770 chess machine (the Turk) The original source says:Amazon's 2005 announcement and its current top page, which this site checked, do not say where the name comes from.
What we still don't know
In 2009, ImageNet was 3.2 million images in 5247 clusters. The paper does not say how many people did the labeling, in which countries, or for how much pay.
Sources & confidence (2)Who put the correct labels on?
2 November 2005Verified in the original
Amazon launches , which lets computers ask people to do small jobs ()
Amazon announced , a system in which a program asks people to do small tasks and receives the results. The announcement calls a computer having many people do work and using the results . Later, in 2009, it was used for ImageNet to check whether a photo really showed a given object. The ImageNet paper describes it as a place where workers are paid.
The popular story:“Amazon explains that it took the name from the chess machine of 1770 with a person inside (the Turk),” people tend to say The original source says:Neither the Amazon announcement of 2005 opened for this site nor its current top page gives the origin of the name (the story of its origin is often told, but we have not yet found Amazon's original words). The words the announcement used were
Fei-Fei Li and colleagues at Princeton University present ImageNet, millions of named photos; the people who checked them were workers hired over the internet
Li and colleagues presented ImageNet, built on WordNet, a dictionary that organizes concepts in a hierarchy, and collecting an average of 500–1,000 photos for each of about 80,000 groups of nouns (the paper's example is ). At the time of publication it had reached 12 areas, 5,247 categories and 3.2 million images. The images were collected with search engines, and workers around the world, paid through Amazon Mechanical Turk, checked one by one whether the target object was really in each picture, deciding by the votes of several people.
The popular story:People tend to say “AI learned by looking at photos on its own” The original source says:The correct labels used for training were checked and attached one image at a time by large numbers of people paid for their work over the internet (paper, page 4). A person inside the machine is not only a story about the chess automaton of 1770.
What we still don’t know:The paper does not state the number of workers, the countries they lived in, or the amount they were paid (unconfirmed).
What the figure shows:When AI learned from photos, the correct labels came from paid workers around the world checking each image and deciding by a vote among several people.
Number of people and the accuracy of the majority vote
232009 · Follow the data
Check on separate training and test sets
An elaborate mechanism, or a huge amount of data: which wins?
In 2009, Alon Halevy, Peter Norvig and Fernando Pereira of Google wrote the paper . It said that , and urged:
But learned performance cannot be known unless you evaluate it on data not used for learning. Change the amount of training data in the figure on the right. When there is little, it gets 100% right on the training data but the accuracy on the test data is low (overfitting).
Correctness on the test=number of hits among the examples not used for learning ÷ number of those examples
practice (training data)
data used to train the model
test (test data)
data used only to check performance after learning
The paper's reasoning
The most successful machine learning on language, the paper pointed out, was statistical speech recognition and machine translation, because a large training set of existed in the world just as it was. It also mentions the 1 trillion words of text data Google released in 2006.
The size of data (open "Details: Size of data" in the figure)
The data that appear in this chapter are 14 cases (the 1986 decision tree table), 97 species (the 1958 bees), 405 abstracts (the 1961 abstracts), 7300 images (the 1992 handwritten digits), 100 million ratings (the 2006 Netflix), 1 trillion words (the 2006 Google), 3.2 million images (the 2009 ImageNet) and 10 million videos (the 2012 cat).
Sources & confidence (1)Check on separate training and test sets
March 2009Verified in the original
Three Google researchers present : stop aiming for a very neat theory, and make data your ally
Google's Alon Halevy, Peter Norvig and Fernando Pereira wrote that the most successful machine learning on language was statistical speech recognition and statistical machine translation. Translation (and transcribing speech is the same) is done routinely every day for real human needs, so, they said, there is naturally out in the world. They also mentioned that in 2006 Google released a data set of 1 trillion words of text (counts of word sequences up to 5 words long). And they stated that based on less data, and in the closing paragraph they wrote
What the figure shows:A learned system's performance must be checked on data it did not train on. With little data, it may score 100% in practice and still fail on new cases.
Size of the data that appear in this chapter (each tick is 10 times more)
This is the same mechanism as "Try it", written as Python code. Press ▶ Run and the result appears below. Change the numbers and run it again.
u09_digits_train_test.py
242006, 2009 · GPUs
The same calculation, many at once
Why did a part for drawing game graphics make AI learning faster?
The GPU, a part for drawing game graphics, can run the same calculation at the same time on many small cores (up to 128 in the 2006 GeForce 8800 GTX). In November 2006, NVIDIA announced CUDA, a system that lets this part be used in ordinary C for calculations other than graphics.
Most of neural network learning is the repetition of "multiply and add". Each one is the same kind of calculation, so it can be divided up and done at the same time. In 2009, Rajat Raina, Anand Madhavan and Andrew Ng of Stanford trained a deep network that learns without correct answers on a GeForce GTX 280 and made it up to 70 times faster than a 2-core CPU. Training a 4-layer network with 100 million parameters dropped from to .
The figure on the right is a race to fill in the cells of a multiplication table with a changing number of cores (a schematic made for this site). The more cores, the more cells are filled in the same time.
Time taken ≒ number of calculations ÷ (number of cores × speed of one core)
assumptions of the figure
every cell takes the same amount of work and the effort of dividing it up is taken as 0 (in reality, things like the time to move data are needed)
What you think you knowThe popular story:Deep learning spread thanks to GPUs The original source says:Google's big 2012 attempt (the cat paper and DistBelief) used thousands of ordinary computers, not GPUs. The paper wrote that GPUs of the time had difficulty speeding up large models that do not fit in memory (under 6 GB). There was more than one road.
Google Brain
A team inside Google that tried computations on huge neural networks. It began its work in 2011, and in a 2012 paper announced DistBelief, a system for training one network on thousands of computers. In April 2023 it merged with DeepMind to become Google DeepMind.
Sources & confidence (3)The same calculation, many at once
8 November 2006Verified in the originalSources disagree
NVIDIA announces CUDA, a system that lets the GPUs it makes for drawing images be used for other computing too
A GPU (a semiconductor chip for graphics) built to draw game images can do the same calculation at the same time on many small computing cores (up to 128 on the GeForce 8800 GTX). Along with the GeForce 8800, NVIDIA announced , a system that lets a GPU be used for computing other than graphics in an ordinary programming language (C). This later became a tool for training neural networks quickly (a 2009 Stanford paper, and the winning team of the 2012 image-recognition contest, also used GPUs).
1 sources and how they were checked
Where sources disagree:The year differs between sources: it was announced in NVIDIA's technical document of November 2006. In the list of candidate events made from memory when this site was first being built, it is 2007. 2007 may be the year the tools were made generally available, but we have not checked the original sources.
Andrew Ng and colleagues at Stanford show that GPUs can speed up neural network training by up to 70 times
Rajat Raina, Anand Madhavan and Andrew Ng trained deep networks that learn without labels (unsupervised learning) on a gaming GPU (NVIDIA GeForce GTX 280). It was up to 70 times faster than an ordinary CPU (2 cores), and cut the training of a 4-layer network with 100 million numbers from . The paper pointed out the problem that, although learning from large amounts of data is important, slow computation kept experiments small.
At Google, “Google Brain”, a team for trying out large neural networks, gets started
According to the official page of X, a division of Google (the footer reads and the title is ), the Google Brain team began work in 2011 and moved from X to Google at the end of 2012. In December 2011, researchers at Google and Stanford University released the first version of a paper on learning from unlabeled YouTube images (later known as Google's cat paper). In a 2012 paper they presented DistBelief, a system for training one large neural network on thousands of computers; they reported training a network 30 times larger than previously reported, and also speeding up the training of a commercial speech recognizer. In April 2023, Google announced that it would merge this team with DeepMind into Google DeepMind.
The popular story:It is often said that “deep learning spread thanks to GPUs (semiconductor chips for graphics)” The original source says:Google's big 2012 efforts (the cat paper and DistBelief) used not GPUs but thousands of ordinary computers (for example 16,000 CPU cores). The paper writes that GPUs of the time were hard to make fast for large models that do not fit in their memory (under 6 GB). The winning team of the image-recognition competition that same year, 2012, used GPUs, so there was more than one road
What the figure shows:Neural network training contains large amounts of the same kinds of calculations, so sharing them among a GPU's many small cores makes it much faster.
252009–2012 · Voices
Voices came before images
Was the first big success of AI using deep learning in images?
From around 2009, deep learning began to beat the statistical method that had been used for a long time in speech recognition (GMM). In 2012, four groups, the University of Toronto, Microsoft Research, Google and IBM, wrote a paper that , showing results such as cutting the word error rate from 27.4% to 18.5% on a conversation transcription test. By August 2012, Google had adopted it for speech recognition on Android. All of this came before the October 2012 image contest.
Around the same time, apps you could ask by voice appeared one after another (VoiceTra from NICT in August 2010, Siri from Apple in October 2011, and Shabette Concier from NTT Docomo on March 1, 2012; Japanese Siri came a week later).
The figure on the right shows how speech recognition errors are counted (word error rate).
Error rate= (substitutions+deletions+insertions) ÷ number of words in the correct text
how to count
line up two sentences with the fewest corrections and count them (dynamic programming, the same idea as DP matching in Chapter 5)
What you think you knowThe popular story:The first big success of deep learning was the image contest in the fall of 2012 The original source says:In speech recognition, results began to appear from 2009, and in August 2012 Google had adopted it on Android. Google's blog names, as the first attempt to use neural networks for speech, a paper from the late 1980s (co-authored by two researchers at ATR in Osaka).
Google's blog (August 6, 2012)
The new version of Android was introduced this way: when you speak to your device, .
The quiz-show Watson (2011)
In February 2011, IBM's Watson beat two all-time champions on the quiz show Jeopardy! ($77,147 against $24,000 and $21,600). It gathered evidence for candidate answers using machine learning, and pressed the buzzer only when its confidence was high enough. The development team itself wrote in 2010 that the Jeopardy! . In the last question of the first game it answered and was wrong.
Sources & confidence (5)Voices came before images
2009–2012Verified in the original
Deep-learning AI had a success in before images (2009–2012)
Deep learning (an AI method that learns with neural networks stacked in many layers) began, around 2009, to beat the statistical method that had long been used in speech recognition (GMM). In 2012, four research groups, at the University of Toronto, Microsoft Research, Google and IBM, wrote up a paper that , and showed results such as cutting the share of wrong words on a speech-transcription test from 27.4% to 18.5% (by about one third). In a blog post on August 6, 2012, Google wrote that with the new version of Android, when you speak to your phone, .
The popular story:People tend to say “deep learning's first big success was the image competition of autumn 2012 (ILSVRC 2012)” The original source says:In the original sources, results began to appear in speech recognition from 2009; in 2011 errors in transcribing conversations were cut by about one third (the four-group paper); and in August 2012 Google was using it for speech recognition on Android (Google's blog). Both came before the results of the image competition in October 2012. Google's blog also says the first attempt to use neural networks for speech was a paper from the late 1980s (co-written by two researchers at ATR in Osaka)
What we still don’t know:Looking back, the four-group paper writes that most of the gains came from the ability to (a frame is a short slice of sound) and from (predicting sound states that are split according to the sounds before and after, with similar ones grouped together). At the same time, it also writes that there is have already been found (PDF page 13). Why deep neural networks work better was still an unfinished explanation at the time.
5 August 2010Verified in the originalSources disagreeJapan
Japan's national research institute NICT releases “VoiceTra” worldwide, a smartphone app that translates what you say and answers in a voice
In the summer of 2010, NICT (National Institute of Information and Communications Technology), a national research institute of Japan, released the speech translation app VoiceTra and the text translation app TexTra free of charge worldwide on the App Store. With VoiceTra, you speak to your iPhone and it translates into another language and answers in a voice. The translation is computed on a server, not on the iPhone. The languages that worked by voice were Japanese, Mandarin, English, Indonesian and Vietnamese (Malay was in preparation). It was aimed at travel conversation, and NICT described the speech translation ability as (translated). The purpose of the release was to evaluate performance and to study ease of use. According to NICT's 2015 announcement, the technology was also used in real services such as Narita International Airport's app. It then became the central app of the Ministry of Internal Affairs and Communications' (April 2014; aiming for a society by 2020 in which foreign visitors to Japan do not have trouble with language).
2 sources and how they were checked
Where sources disagree:Statements about the release date: the August 5, 2010 announcement says “we have released it” (translated), while the 2015 announcement says “released in July 2010” (translated). The release date on the App Store and the date of the press announcement may have differed, but this site does not settle the question.
Apple unveils the iPhone 4S with “Siri”, which helps with tasks when you talk to it (at first )
On October 4, 2011, together with the iPhone 4S, Apple announced Siri, that helps with tasks when you ask by voice. In the example in the announcement, when you ask , it understands that you are looking for a weather forecast. At first it was (a trial version), and the only languages supported were English (for the U.S., the U.K. and Australia), French and German. In Japan the iPhone 4S went on sale on October 14, but Siri with Japanese support came only on March 8, 2012. Siri's origin is CALO, research that SRI International carried out under PAL, a program of DARPA, the U.S. Department of Defense's research agency. It was spun off as a company in 2007, and in April 2010 Apple acquired that company (according to SRI's and DARPA's explanations).
The popular story:People tend to think “Apple built Siri from scratch” The original source says:In the explanations of the research institute SRI and of DARPA themselves, it started as CALO, research under PAL, a program to support military decision-making (5 years, $150 million, more than 300 people). It was spun off in 2007 as the company Siri Inc., and Apple acquired that company in 2010. At the announcement it was also still (a trial version)
NTT Docomo launches the voice dialogue service “Shabette Concier” (March 2012, one week before Japanese Siri)
On March 1, 2012, NTT Docomo launched Shabette Concier, a “voice agent function” (translated) that lets you look up information and operate email, alarms and more just by talking to a smartphone, as a free app for Docomo's Android smartphones. It reads the speaker's “intention” (translated) from the combination of words and chooses a suitable service or answer. Complex processing such as speech recognition was done not on the device but on Docomo's network side (the cloud). Docomo described it as following the “interpreter phone” (translated) cloud service, whose trial offering began in November 2011. The launch of Japanese Siri came one week later, on March 8.
IBM's question-answering AI “Watson” beats two all-time champions on the U.S. quiz show Jeopardy!
In February 2011, Watson, a question-answering system (a program that reads a question in natural language and returns an answer) developed by an IBM research team, played two games against two all-time champions on the U.S. TV quiz show Jeopardy!. The final scores were $77,147 for Watson, $24,000 for Ken Jennings and $21,600 for Brad Rutter (IBM's history page). Watson was not connected to the internet; it produced candidate answers with many methods, gathered evidence, and pressed the buzzer only when the that machine learning (a way of learning from data) had put together was high enough. On the final question of the first game, a question about U.S. cities, it answered and was wrong.
The popular story:People tend to take it as “beating humans at a quiz means Watson understood language the way humans do” The original source says:The team that built it wrote in a 2010 paper that the Jeopardy! . IBM's history page also has the leader's words . In the real show, too, it got a low-confidence answer (Toronto) wrong
What the figure shows:Deep learning's first major success was in speech, before images. Speech recognition errors are counted by finding the fewest edits needed to make 2 sentences match.
262010–2012 · Deep networks
Keeping the signal alive in deep stacks
What happens when you stack 30 layers?
In deep neural networks, the signal that carries the error backward weakens or grows too large at every layer. In 2010–2012, three tricks to prevent this came together.
One is the activation function ReLU (if the input is positive it passes through unchanged, if negative it becomes 0; on the positive side the slope is always 1). Another is a starting range for the weights that is matched to the number of neurons (Glorot and Bengio, 2010). The third is dropout (2012), which at each training step makes half of the neurons in a hidden layer rest at random.
In the figure on the right, change the activation function and the starting values and look at the size of the signal at each layer.
ReLU(x)=max(0, x) weight range ±√6 ÷ √(nin +nout), standard deviation √(2 ÷ nin)
nin, nout
number of input and output neurons of the layer (the range is Glorot and Bengio's rule, the standard deviation is He and colleagues' rule)
What you think you knowThe popular story:ReLU was invented by Hinton and colleagues in 2010 The original source says:The 2010 paper itself cites an attempt by Jarrett and colleagues in 2009. An activation function of the same shape is in Kunihiko Fukushima's network of 1968–1969 (Chapter 4), but the three papers of 2010–2012 do not cite Fukushima.
Details of the three tricks
The S-shaped sigmoid has a slope of almost 0 when the input is large. ReLU was named in 2010 by Vinod Nair and Hinton of the University of Toronto, and the winning team of the 2012 image contest also used it. Glorot and Bengio's rule is meant to keep the spread of the signal about the same in every layer, and in 2015 Kaiming He and colleagues published a version for ReLU. Dropout (July 2012) prevents neurons from depending on one another and reduces overfitting. At test time all neurons are used and the output weights are halved. The detailed 2014 paper says that the idea was inspired by a theory of the role of sex in evolution (Livnat and colleagues, 2010).
The name Xavier
The 2010 paper calls this rule normalized initialization. The 2015 paper by He and colleagues writes that Xavier is its name in the tool Caffe (Xavier is the same as the first name of the lead author).
Dropout is often thought of as a trick of AlexNet, but the AlexNet paper uses it by citing a paper from the same lab in July 2012 as a .
Sources & confidence (3)Keeping the signal alive in deep stacks
2010–2012Verified in the original
How ReLU (0 if negative, unchanged if positive) became the standard in deep learning: named in 2010, widespread in 2011–2012. The same form appears in Kunihiko Fukushima's 1968 model
ReLU (rectified linear unit) is a unit (neuron) that outputs its input x as it is if x is greater than 0, and outputs 0 if x is 0 or less; as a formula, max(0, x). Sigmoid (S-shaped function) and tanh have a slope that becomes almost 0 (they saturate) when the input is large, but ReLU has a slope of always 1 for positive values, so the signal sent back through backpropagation is less likely to fade. In 2010, Nair and Hinton of the University of Toronto showed that if the 0-or-1 hidden units of a restricted Boltzmann machine are replaced by many copies of a unit with the same weights, shifted slightly one from another, their behavior can be approximated by noisy rectified linear units. They also reported better results in recognizing photos of objects and faces. In 2011, Glorot and colleagues at the Université de Montréal wrote that this unit is also a better model of real nerve cells and can form sparse representations, with many values exactly 0. They also wrote that, with plenty of labeled data, deep networks can learn well without prior unsupervised learning. The 2012 AlexNet paper, following Nair and Hinton, called it ReLU and wrote that it trained several times faster than the same network using tanh.
The popular story:People tend to say “ReLU was invented by Hinton and colleagues in 2010” and “it started with AlexNet in 2012” The original source says:The 2010 paper by Nair and Hinton itself says that Jarrett and colleagues had tried max(0,x) in 2009 under the name positive part. Moreover, a unit of the same form (output the input as it is if positive, 0 if negative) was used in networks modeled on the visual system in a 1968 Japanese-language explanation and a 1969 English paper by Kunihiko Fukushima of NHK. However, none of the three papers of 2010, 2011 and 2012 cites Fukushima's paper. What was new in 2010 was deriving this unit from the probabilistic model of the restricted Boltzmann machine and naming it; what was new in 2011 was showing that prior unsupervised learning is not needed.
What we still don’t know:I have not checked whether the 2010 paper was the first to use the name “rectified linear” (the 2010 paper writes that Jarrett and colleagues in 2009 called it positive part). Fukushima's 1969 English paper (IEEE) is behind a paywall and I have not opened its text. I have not opened the 2009 paper by Jarrett and colleagues either.
How to set initial weights: matching the spread to the number of neurons in the layers (Glorot and Bengio 2010), and a version for ReLU (He et al. 2015)
Before a deep network is trained, what values are put into its weights (initialization) greatly changes whether it can be trained at all. If the signal grows or shrinks by a fixed factor each time it passes through a layer, then after many layers the backpropagated gradient vanishes or blows up. In 2010, Glorot and Bengio of the Université de Montréal showed that with the commonly used method, choosing from the range ±1/√n where n is the number of neurons in the previous layer, the spread (variance) of the backpropagated signal gets smaller from layer to layer. They proposed normalized initialization, which uses the numbers of neurons on the input side and the output side of a layer, nj and nj+1, and chooses uniformly from the range ±√6/√(nj+nj+1). The aim is to keep the spread of the forward-propagated values and of the backpropagated signal about the same in every layer. In 2015, Kaiming He and colleagues at Microsoft Research pointed out that this derivation rests on the assumption that neurons are linear (output proportional to input), which does not fit ReLU. For ReLU they proposed choosing from a normal distribution with mean 0 and standard deviation √(2/n). They write that a 30-layer network could be trained this way, while training stalled with Glorot and Bengio's method.
The popular story:People tend to think the name was given by Glorot and Bengio's paper The original source says:The 2010 paper itself calls it normalized initialization. He and colleagues' 2015 paper writes that the name Xavier is what the tool Caffe (Jia et al., 2014) calls it (Xavier is the same as the first name of the lead author). It is also often said that “In 2006, pre-training made it possible for the first time to train deep networks”, but the 2010 paper writes that every success since 2006 came from either a new initialization or a new training method, and it investigated, from the side of the initial weights and the activation function (the function that decides a neuron's output), why standard gradient descent had not worked well.
What we still don’t know:I have not confirmed when, or by whom, the name Xavier was given inside Caffe (I have not opened the Caffe paper or the program's records). I have not confirmed who started the ±1/√n method, which the 2010 paper calls a commonly used heuristic, or when.
Dropout: randomly resting half of the hidden-layer neurons at each training step to prevent overfitting (2012)
Dropout is a method in which, each time one piece of data is shown during training, hidden-layer neurons are temporarily removed at random with a fixed probability (often 0.5), that is, their outputs are set to 0, while the network learns. It prevents specific neurons from depending on one another (co-adaptation) and suppresses overfitting, in which a network fits only the training data too closely and makes mistakes on new data. At use time, all neurons are used and the outgoing weights are halved (in the 2014 paper, the weights are multiplied by the retention probability p). The result is close to averaging the predictions of many thin networks that share weights. The July 2012 arXiv paper (Hinton, Srivastava, Krizhevsky, Sutskever and Salakhutdinov of the University of Toronto) tried it on handwritten digits (MNIST) and others, and wrote that it set records in speech and object recognition. The AlexNet paper of the same year also used it in the two fully connected layers, and wrote that without it overfitting became serious.
The popular story:People tend to think “dropout was a trick of AlexNet in 2012” The original source says:The AlexNet paper itself uses dropout as a , citing the July 2012 arXiv paper from the same lab. The detailed 2014 paper explains that the idea was inspired by a theory of the role of sexual reproduction in evolution (Livnat et al., 2010). It also writes that the idea of adding noise to neurons itself had already been used in denoising autoencoders, which add noise to the input and reconstruct it (Vincent et al., 2008 and 2010). What the paper claims as its own novelty is that it can be used in hidden layers too, and that it can be regarded as an average of a large number of networks.
What we still don’t know:I have not written down when dropout was conceived, or the several accounts Hinton himself later gave of how it came about, because I have not looked for the original talks.
What the figure shows:In a deep stack of layers, the backward signal weakens at each layer until it disappears. ReLU and carefully chosen starting weights allowed signals to travel through even 30 layers.
Size of the backpropagated signal at each layer (each tick is 10 times more)
This is the same mechanism as "Try it", written as Python code. Press ▶ Run and the result appears below. Change the numbers and run it again.
u09_relu_init.py
272012 · Three in one year
Traffic signs, and Google's cat
In 2012, how many results did deep learning produce?
In February, Dan Cireşan, Ueli Meier and Jürgen Schmidhuber of the Swiss lab IDSIA trained, on GPUs, a method that lines up several deep convolutional networks and averages their results. On German traffic sign recognition (12,569 test images) it reached an error rate of 0.54%. According to the paper, the human error rate was twice that. But there were only 43 answers, so it is a different level of difficulty from photos with 1000 categories.
In June, Quoc Le, Andrew Ng and colleagues of Google and Stanford University took images, one frame at a time from 10 million YouTube videos, and had a 9-layer network with 1 billion connections learn from them with no correct labels (3 days on 1000 computers and 16,000 cores). Even though nobody taught it, neurons that respond to human faces appeared, and so did ones that respond to cat faces and human bodies. Google's blog stressed that .
Then in October came the image recognition contest (the next spread). In the figure on the right, compare the cat neuron's score with "answering at random".
Table 1 of the paper: human face neuron 81.7%, cat face neuron 74.8%, answering at random 64.8% (the cat test). Read a score by comparing it with the points you would get without doing anything.
What you think you knowThe popular story:Google's AI learned to recognize cats on its own The original source says:The paper's main result is the human face neuron (81.7%), and the cat face neuron is 74.8%. It was a test in which random answers would score 64.8% (Table 1 of the paper).
What is going on inside
The paper examines from outside what a neuron responds to, by lining up the images that excite it most and by synthesizing images that make it respond more strongly. Tricks for looking at internal behavior had already started in this era. This question continues through Chapter 12.
Sources & confidence (2)Traffic signs, and Google's cat
13 February 2012Verified in the original
Switzerland's IDSIA: a deep convolutional neural network trained on GPUs gets half the human error rate at recognizing traffic signs (8 months before the 2012 competition)
Dan Cireşan, Ueli Meier and Jürgen Schmidhuber of the Swiss institute IDSIA trained, on gaming GPUs, a method that lines up several deep convolutional neural networks and averages their outputs (a multi-column deep network). On recognition of German traffic signs (12,569 test images) it recorded an error rate of 0.54%, and according to the paper the human error rate was twice that. On the handwritten-digit set MNIST it also reached 0.23%, which they write is the first time a machine came close to the human figure of about 0.2%.
The popular story:People tend to say “2015 was the first time AI surpassed humans at telling images apart” and “the first success of deep convolutional neural networks on GPUs was the University of Toronto in autumn 2012” The original source says:IDSIA's February 2012 report writes that a deep convolutional neural network trained on GPUs had half the human error rate at recognizing traffic signs (note, though, that the task was traffic-sign recognition and the paper's system has 43 outputs, so the number of classes differs from ImageNet, which has 1,000)
26 June 2012Verified in the originalSources disagree
Google's cat experiment: a network that learned only from unlabeled YouTube images developed a neuron that responds strongly to cat faces
Quoc Le, Andrew Ng and colleagues at Google and Stanford University showed a 9-layer neural network with 1 billion connections still images (200×200 pixels), one taken from each of 10 million YouTube videos, with no labels, and let it learn. It took 3 days on 1,000 computers (16,000 cores). Without being taught, it developed neurons that respond to human faces, and also ones that respond to cat faces and human bodies. In its blog, Google described how , and stressed that .
The popular story:People tend to say “Google's AI learned to recognize cats by itself” The original source says:The paper's main result is the neuron for human faces (81.7%); the cat-face neuron scored 74.8% on a test of telling cat faces from non-cat images, a test built so that random answers would score 64.8% (the paper's Table 1). The paper also claims not that it “saw many cats” but that it found, on its own, groups of features in unlabeled images. Google's blog wrote that the network “discovered” a cat's appearance from unlabeled images alone (the blog itself puts quotation marks around the word discovered), in its words
What we still don’t know:The paper examines what the trained neurons respond to, by lining up the images that excite them most and by constructing images that make them respond more strongly. Methods for looking inside the network were already being devised at that time.
2 sources and how they were checked
Where sources disagree:Training period: the paper says “3 days”, while Google's blog uses the figure of speech of showing it YouTube for a week. The cat figure: the upper row of the paper's Table 1 gives 74.6%, while the lower row and the main text give 74.8%
What the figure shows:Deep learning successes followed one another in 2012, and Google's cat neuron scored 74.8% on a test where random answers would score 64.8%. Results should be read against the score possible without doing anything useful.
282012 · The image contest
A huge win, with Tokyo in second
What was the “turning point” in the image contest of October 2012?
In ILSVRC, a contest to recognize 1000 kinds of objects from photos, SuperVision, the team of Alex Krizhevsky, Ilya Sutskever and Hinton of the University of Toronto, reached an error rate of 15.3% in the rule that allows 5 guesses (16.4% using only the data the contest provided). The team in second place, from the University of Tokyo (ISI), had 26.2%, a gap of about 11 points. They trained an 8-layer convolutional neural network with 60 million parameters on two gaming GPUs. The contest organizers wrote that for research on recognizing many kinds of objects.
Convolution is a calculation that slides a small window (a filter) over an image and responds strongly at places that have the same shape as the window. The structure continues from Kunihiko Fukushima's Neocognitron of 1979 (Chapter 5) and the network of Yann LeCun and colleagues in 1989. In the figure on the right, choose a window shape and slide it over the image. The graph below shows eight years of the contest's first-place error rate.
output(i, j)=Σa Σb window(a, b) × image(i+a, j+b)
window
a 3×3 table of weights (after learning, it takes shapes that respond to vertical lines, horizontal lines and so on)
What you think you knowThe popular story:Deep learning was suddenly born in 2012 The original source says:The winning mechanism was the network of LeCun and colleagues from 1989 (itself based on Fukushima's Neocognitron of 1979), made bigger and trained on GPUs. The paper itself writes that .
What is inside "humans: 5.1%"
The first-place error rate fell from 28.2% in 2010 to 2.3% in 2017. The organizers measured humans too. One person who practiced on 500 images got 5.1% on 1500 images, and another who practiced on only 100 images got 12.0% on 258 images. The Microsoft researchers themselves, who wrote in 2015 that machines had passed humans, also wrote that this .
Sources & confidence (2)A huge win, with Tokyo in second
13 October 2012Verified in the originalSources disagree
In the image-recognition competition ILSVRC 2012, the University of Toronto's deep convolutional neural network (SuperVision) wins by a wide margin. The runner-up is a University of Tokyo team
In a competition to identify 1,000 kinds of objects from photos, the team “SuperVision” of Krizhevsky, Sutskever and Hinton held its error rate, under a rule allowing 5 guesses, to 15.3% (16.4% if only the data supplied by the competition is used). The runner-up, a University of Tokyo team (ISI), scored 26.2%, a gap of about 11 points. They trained an 8-layer neural network with 60 million parameters on 2 gaming GPUs (graphics chips, made by NVIDIA). In a later paper, the competition's organizers wrote that for research on recognizing a wide variety of objects.
The popular story:People tend to say “deep learning was suddenly born in 2012” The original source says:The winning system was a large-scale version of the of LeCun and colleagues in 1989 (whose origin is Fukushima's Neocognitron of 1979), trained on GPUs. The paper itself writes that (prefaced by saying that every experiment points that way). People also tend to think “2012 was a Western story”, but the runner-up was a University of Tokyo team
4 sources and how they were checked
Where sources disagree:Training time: the team introduction on the competition page says “about 1 week”, while the paper says “5 to 6 days”. The first-place figure: there are two, 15.3% (with extra data) and 16.4% (supplied data only), and the table in the organizers' paper uses 16.4
In the image-recognition competition, the top team's fell from 28% to 2% in 8 years. The “human: 5.1%” figure is one person's measurement
In a competition to identify 1,000 kinds of objects from images (with a rule allowing up to 5 guesses), the winning error rate fell: 28.2% in 2010, 25.8% in 2011, 15.3% in 2012 (16.4% with only the supplied data), 11.7% in 2013 (11.2% when outside data was also used), 6.7% in 2014, 3.6% in 2015, 3.0% in 2016 and 2.3% in 2017. The organizers also measured the human error rate. It was a test with two people (A1 and A2): A1, who had practiced on 500 images, scored 5.1% on 1,500 images, and A2, who had practiced on only 100 images, scored 12.0% on 258 images.
The popular story:People tend to say “in 2015 AI surpassed humans (error rate 5.1%) at image recognition” The original source says:The 5.1% error rate is the result of one (A1) of the two skilled judges tested in the organizers' paper, who identified 1,500 images after practicing on 500. The other, who practiced on only 100 images, scored 12.0%. The researchers who themselves reported surpassing humans (Microsoft, 2015) also wrote that this result in general.
4 sources and how they were checked
Where sources disagree:How the “first place” result is counted differs by year. In 2012 and 2013 there were entries that used outside data and entries that did not, and the table in the organizers' paper uses the figures for the supplied data only (16.4 and 11.7). The 2010 and 2011 websites list the error rate for the same up-to-5-guesses method under the item name “flat cost”.
What the figure shows:At the 2012 image competition, a convolutional network beat the runner-up by about 11 percentage points. Inside was the same small-window sliding calculation developed from 1979 and 1989 onward.
First-place error rate of the contest (2010–2017)
This is the same mechanism as "Try it", written as Python code. Press ▶ Run and the result appears below. Change the numbers and run it again.
u09_convolution.py
292005–2012 · Japan · Shogi and checkers
Learning from records, voting on moves
Is the strength of a shogi program decided by how strong at shogi its maker is?
Most shogi programs until then had their makers decide by hand how to score whether a position was good or bad. Bonanza, made by Kunihito Hoki, a chemistry researcher (released in 2005), had the machine learn how to score so that it would play the same moves as professionals' game records. In May 2006 it won the World Computer Shogi Championship in its first appearance, and in March 2007 it played a public game against Ryuo Akira Watanabe (Watanabe won in 112 moves).
In October 2010, “Akara 2010”, from the Information Processing Society of Japan, beat Ichiyo Shimizu, the women's Osho title holder. It was a system that made the top four Japanese programs think about the same position and played the move chosen most often, by majority vote. At the first Shogi Den'o-sen in January 2012, Bonkras, based on Bonanza, beat Kunio Yonenaga, a retired professional and then president of the Japan Shogi Association, a lifetime title holder of Kisei, in 113 moves.
The figure on the right shows a learning method that updates weights so that the professional's move gets the highest score (a schematic made for this site).
w ← w+η(f(professional's move) − f(m))
m
a move that came to score higher than the professional's move
w, η
the weights for scoring and the learning rate
f
features of the position after making a move (material gain, king safety, how well the pieces work)
What you think you knowThe popular story:To build a strong shogi program, its maker has to be strong at shogi too The original source says:The author of Bonanza was a chemistry researcher who hardly knew shogi, yet by having the machine learn how to score from game records he won the world championship in his first entry.
The method was made public
How Bonanza works was made public (and in 2009 even the program's source code), and shogi programs based on Bonanza, such as Bonkras, were born.
Checkers was solved (2007)
Jonathan Schaeffer and colleagues at the University of Alberta kept computers running from 1989 and proved that in checkers The board has about 5×10²⁰ positions. One of the central tools was the search method (df-pn) from Ayumu Nagai's 2002 doctoral thesis at the University of Tokyo, and Akihiro Kishimoto is among the authors.
The first game between a professional and a shogi program is given as 2005, 2007 or 2013 depending on the source. Depending on whether you say in public, in an official setting, or officially, any of them can be called the first.
Sources & confidence (6)Learning from records, voting on moves
5 May 2006–2007Verified in the originalJapan
The shogi program Bonanza, which learns from professional game records, wins a world championship and the next year plays Ryuo Akira Watanabe in a public game
The shogi program Bonanza was developed by Kunihito Hoki, a chemistry researcher, and released free of charge in June 2005. Unlike many earlier shogi programs, it used “full-width search” (translated), which reads out every move allowed by the rules. It also adopted (translated), training the way positions are evaluated so as to play the same moves as professionals' game records. The author has said that he applied the idea from his research on controlling chemical reactions to learning position evaluation. In May 2006, at its first entry in the 16th World Computer Shogi Championship, it won the round robin of the 8 finalists with 6 wins and 1 loss. On March 21, 2007, it played a public game against Ryuo (Dragon King title holder) Akira Watanabe, with 2 hours of thinking time each. Bonanza was also said to have the upper hand in the middlegame, but the Ryuo won in 112 moves. Later its method (and in 2009 the inside of the program) was made public, and shogi programs based on Bonanza, such as Bonkras, were born.
The popular story:“To make strong shogi software, the maker must also be strong at shogi” (according to an interview article, until then developers built the evaluation by hand, so they were often strong amateur players) The original source says:Bonanza's author was a chemistry researcher who hardly knew shogi. He adopted a method in which the machine learns how to evaluate from professional game records, and won the world championship at his first entry
11 October 2010Verified in the originalSources disagreeJapan
The Information Processing Society of Japan's shogi software “Akara 2010” beats Women's Osho Ichiyo Shimizu
In April 2010 the Information Processing Society of Japan sent the Japan Shogi Association a challenge letter, (translated). The association's president, Kunio Yonenaga, accepted and chose as the first opponent Ichiyo Shimizu, the Women's Osho (a title of women's shogi) and the leading woman professional. The game was played on October 11 at the University of Tokyo's Hongo campus; Shimizu resigned at move 86 and Akara 2010 won. Akara 2010 had Japan's four top shogi programs (Gekisashi, GPS Shogi, Bonanza and YSS) think about the same position and played the move chosen most often, a (translated) system, and used 208 computers (658 cores in total).
The popular story:People tend to think “a shogi program's strength comes from one program” The original source says:Akara 2010 was a joint effort of Japan's computer shogi developers, a system in which four rival programs took a majority vote. The challenge letter described the society's shogi software as having (translated) (meaning strength on a par with the Meijin), but the first opponent was a woman professional chosen by the association
3 sources and how they were checked
Where sources disagree:Number of computers used: the August announcement says 169 computers and 676 cores; the addendum after the game says 208 computers and 658 cores (with a note that cores were counted as physical cores and that some were not used because of failures)
14 January 2012Verified in the originalSources disagreeJapan
1st Shogi Den-O-sen: Kunio Yonenaga, Eternal Kisei, loses in 113 moves to the shogi software “Bonkras”
On January 14, 2012, at the Shogi Kaikan (Shogi Hall) in Tokyo, Kunio Yonenaga, president of the Japan Shogi Association and a retired professional holding the lifetime Kisei title (Eternal Kisei), played the shogi software “Bonkras”, developed by Hidenori Ito of Fujitsu Laboratories, and Bonkras won in 113 moves. Bonkras is a version of the publicly available shogi software Bonanza, improved to search in parallel on several computers; that day it read about 18 million moves per second on 6 servers. The name is a coined word combining Bonanza and , which means a parallel computer (a system that calculates on several computers at the same time) (as explained in Fujitsu Laboratories' announcement).
The popular story:People tend to say “at the 1st Den-O-sen, a professional lost to software” The original source says:In the association's own summary, the 1st was a game against a retired player (then the association's president), and games against active professionals began with the 2nd, in 2013. The winning software was also developed by a single researcher from the publicly available Bonanza, which is another example of publicly shared knowledge building up
2 sources and how they were checked
Where sources disagree:Sources differ on “the first official game between a professional and software”: the Shogi Association's column (2017) calls the 2nd Den-O-sen in 2013 “officially for the first time”, but a note in Fujitsu's announcement (2012) says Ryuo Akira Watanabe played Bonanza in 2007. The difference seems to lie in what is called “official” (this site's estimate)
Checkers : computers prove it is always a draw if neither side makes a mistake
Jonathan Schaeffer and his colleagues at the University of Alberta kept dozens of computers running almost nonstop from 1989 and proved that Checkers has about 5×10²⁰ (500 quintillion) board positions, making it the hardest well-known game solved so far. One of the main tools in the proof was the search method (df-pn) from Ayumu Nagai's 2002 doctoral thesis at the University of Tokyo. Akihiro Kishimoto, who by the time of the paper had moved to a university in Hakodate, Hokkaido, is also among the co-authors.
The popular story:People tend to think “game AI is a Western story” The original source says:In the 2007 checkers proof, df-pn, the search method of Ayumu Nagai of Japan (University of Tokyo), was used as a core tool, and Akihiro Kishimoto of Japan was among the authors
3 sources and how they were checked
Checkers Is Solved (Science 317, 1518–1522) (Jonathan Schaeffer, Neil Burch, Yngvi Björnsson, Akihiro Kishimoto, Martin Müller, Robert Lake, Paul Lu, Steve Sutphen · 2007 · p.1518 abstract · Original source)
18 September 2005–2013Verified in the originalSources disagreeJapan
When did a professional first play a shogi program? The sources split between 2005, 2007 and 2013
For the (translated) game between a professional shogi player and a shogi program, the sources name three different years. (1) For the September 2005 game between TACOS and 5-dan Takayoshi Hashimoto, the newspaper of the time reported, as the words of its developer, Professor Hiroyuki Iida, (translated). (2) For the March 2007 game between Bonanza and Ryuo (Dragon King title holder) Akira Watanabe, a researcher's commentary wrote (translated), and a 2020 interview article wrote (translated). (3) For the 2nd Shogi Den-O-sen in 2013, a column by the Japan Shogi Association (2017) put it as (translated). Because of the difference in the words , and (translated), every one of them can be called the (translated).
The popular story:“The first game between a professional shogi player and a shogi program was in such-and-such a year,” people tend to say, pinning it to one year The original source says:The sources describe three years, 2005, 2007 and 2013, as the (translated). Because of the difference in the words (translated) (the 2005 TACOS game), (translated) (the 2007 Bonanza game) and (translated) (the 2013 Den-O-sen; the 1st in 2012 was a retired player), each of them can be said to be correct
4 sources and how they were checked
Where sources disagree:The year of the (translated): 2005 (The Sankei Shimbun, Professor Hiroyuki Iida's remarks; public setting, even game)/2007 (Hitoshi Matsubara's commentary, official setting; interview article, in public)/2013 (Japan Shogi Association column, officially). The conflicts section of the record on the 1st Shogi Den-O-sen in 2012 also has a mismatch between 2007 and 2013
The Japan Shogi Association decides that when professionals play software in public, they should (translated)
On September 18, 2005, in Komatsu, Ishikawa Prefecture, the shogi software (translated), developed by Professor Hiroyuki Iida and colleagues at the Japan Advanced Institute of Science and Technology, played a public even game (no handicap) against 5-dan Takayoshi Hashimoto. The professional won, but the software had the advantage until the endgame (The Sankei Shimbun). About a month later, on October 14, the Japan Shogi Association announced that, (translated), when players and women's professionals play shogi software in public, they should consult and confirm in advance with the organizer and the Association. The same announcement also says, (translated).
The popular story:“The Shogi Association banned professionals from playing computers,” people tend to say (a researcher's 2007 commentary also wrote that it (translated)) The original source says:The Association's own announcement says (translated), and it also writes (translated). In fact, a year and a half later, in March 2007, a public game between Ryuo (Dragon King title holder) Akira Watanabe and Bonanza took place (see the record → The shogi program Bonanza, which learns…)
What the figure shows:Instead of having developers choose shogi scores, training the machine to give the highest score to the move a professional would choose makes it stronger even if its developers are not shogi experts.
Number of training rounds and the share of professional moves chosen
This is the same mechanism as "Try it", written as Python code. Press ▶ Run and the result appears below. Change the numbers and run it again.
u09_bonanza_learn.py
302001, 2011 · Japan · Predictions
Deadline goals and market forecasts
What "by this year" targets were being put out in Japan in this era?
In November 2011, the National Institute of Informatics began (the Todai Robot Project). Its aims were to bring together AI research that had split into fine pieces since 1980, and to give young people a dream. As concrete markers it set goals with deadlines: (translated) and (translated) (the check against the results is in Chapter 10).
Ten years earlier, in 2001, a report by an industry association of Japan's machinery and robot industry forecast that the robot market would be about 3 trillion yen in 2010 (of which 1.5 trillion yen in the daily-life field) and about 8 trillion yen in 2025. In reality, though, the next-generation robot market outside manufacturing was 1.4–6.8 billion yen a year, not even reaching 10 billion yen, and in fact shrinking (analysis by Hirohiko Arai of the National Institute of Advanced Industrial Science and Technology, 2009). In the government's 2015 strategy, too, the whole robot market was written as (translated).
The figure on the right lines up the forecasts and the actual figures on a scale that goes up by tens.
forecast ÷ actual=how many times off
scale in tens
1, 10, 100 are spaced equally (the length between two points shows how many times off)
What you think you knowThe popular story: (translated) (the accepted view among robot researchers) The original source says:The industry association's own survey found that the next-generation robot market stayed at 1.4–6.8 billion yen a year and was in fact shrinking (−10.2% a year on average).
What the forecast was made of
These numbers were also used in explanations for the state to provide research funding. The forecast was computed in money terms from what share (5–50%) of the work people do now robots could take over. The person who compiled the report later wrote that it was (translated) and that (translated) (quoted in Arai's analysis).
In 2006–2008, the state agencies and the industry association itself lowered the forecast for 2010 to 50–75 billion yen.
Sources & confidence (2)Deadline goals and market forecasts
7 November 2011Verified in the originalSources disagreeJapan
The National Institute of Informatics launches (Todai Robot, Torobo-kun). The goals: a high score on the Center Test by 2016, passing the University of Tokyo entrance exam in 2021
In 2011 the National Institute of Informatics began the , which develops a computer program that solves university entrance exam questions. Its aims were to reintegrate AI research, which had split into specialized subfields since 1980, and to give young people a dream. As concrete targets it set (translated) and (translated).
What happened next:In the NII's announcement of November 14, 2016, the program scored 525 points (national average 454.8) on a mock Center Test (5 subject areas, 8 subjects), a standardized score (hensachi) of 57.1. That is a score with an 80% or higher chance of admission at 512 private universities (1,343 departments) and 23 national and public universities (30 departments). On the mock University of Tokyo entrance exam it reached a standardized score of 76.2 in mathematics (science track), fully solving 4 of the 6 problems without human help. For how the Center Test effort was treated at the briefing on the same day, see the 2016 Torobo-kun event (record → Torobo-kun (the Todai Robot): What Happ…). This site has not found an announcement that it passed the University of Tokyo entrance exam in 2021.
2 sources and how they were checked
Where sources disagree:The wording of the goal differs. The NII's 2011 announcement says (translated), while the 2012 paper says “by 2021, obtain a passing-level score on the University of Tokyo secondary exam questions” (translated).
May 2001–2015Verified in the originalSources disagreeJapan
A prediction: in 2010 Japan's robot market will be about 3 trillion yen (a 2001 forecast by industry groups). The home and hospital robot market actually did not reach 10 billion yen
In 2001, a report by the Japan Machinery Federation and the Japan Robot Association forecast the robot market in five fields. For 2010 it gave 1.5 trillion yen for daily life, 850 billion yen for manufacturing, 290 billion yen for public use, 260 billion yen for medical care and welfare, and 90 billion yen for biotechnology (about 3 trillion yen in total), and for 2025 about 8 trillion yen in total. According to a 2009 analysis by Hirohiko Arai of the National Institute of Advanced Industrial Science and Technology, this forecast was used as the basis for explaining research and development budgets by the Ministry of Economy, Trade and Industry and others, and it spread through news reports too. A 2005 Aichi Expo announcement also stated that (translated).
The popular story:“By expanding into non-manufacturing fields, the robot market will grow dramatically” (translated) (the words that Arai's analysis cites as a widely held belief among robot researchers; non-manufacturing fields means home, nursing care, public use and so on) The original source says:The industry association's own survey of actual conditions found that the next-generation robot market outside manufacturing was 1.4 to 6.8 billion yen a year over the period surveyed, never reaching 10 billion yen, and was in fact shrinking (an average of −10.2% a year, in Arai's analysis). The government's 2015 strategy also wrote that the whole robot market was (translated).
What happened next:According to Hirohiko Arai's 2009 analysis, the market for next-generation robots outside manufacturing stayed at 1.4 to 6.8 billion yen a year (under 10 billion yen) in surveys of actual conditions by the industry association and others, far from the 2010 forecast (about 2 trillion yen outside manufacturing). In 2006–2008, national agencies and the industry association itself revised the 2010 forecast downward to 50–75 billion yen. In 2015 the government's gave the size of the robot market as “600 billion yen at present” (translated) and set a new target of 2.4 trillion yen in 2020 (from the Ministry of Internal Affairs and Communications' 2015 white paper). This site has not checked the 2025 forecast.
4 sources and how they were checked
Where sources disagree:The same association's report of the year before (2000) had forecast the 2010 market outside manufacturing at 623.5 billion yen (Arai's analysis; we have not opened the original report). That is an increase of 3.3 times in one year.
What the figure shows:Japan in this period saw predictions and goals with target years. Robot market forecasts differed from reality by factors of several hundred or more.
Closing
Old tools, newly assembled
Most of the tools for learning from data are older than the words "artificial intelligence." Bayes (1763), least squares (1805), the S-shaped curve (1838), singular value decomposition (1873), principal components (1901, 1933). The tools that followed were also born, for the most part, in fields outside AI: bee classification, parole prediction, research on sending signals, psychological tests, questions about animal behavior. And many people arrived at the same methods separately. "Who was first" usually cannot be pinned on one person.
What was new in 2006–2012 was less the ideas than the conditions: fast computation with GPUs, huge amounts of data labeled by large numbers of people, and good starting weights and activation functions. In 2012, when these came together, deep learning produced result after result in speech, traffic signs, cats and the image contest.
There is still a lot we don't know. How many people labeled ImageNet, in which countries, for how much pay? What goes on inside a trained network (attempts to look from outside had only just begun)? Why do HMMs and naive Bayes, which treat things that are not independent as independent, still work well? Why, in boosting, does the test score keep improving after the training error reaches 0? Why is a deeper stack of layers better? None of these has a clear answer in this chapter's sources. In the next chapter, after 2012, AI beats professionals at Go and shogi, the Todai Robot's results are checked, and the forecast that "jobs will disappear" becomes a topic around the world.
Other events in this chapter (1; all with original sources and confidence grades)
2006Verified in the originalSources disagree
The two papers behind , which made Go programs strong
In Go, there are many possible moves in a single position, and evaluating a position is hard. So search methods like those used for chess did not make programs strong. In 2006, a paper by Coulom of France presented a new search method (Monte-Carlo Tree Search, MCTS). It plays out a game at random many times from a position to the end (Monte Carlo), compares moves by their win rates, and examines the more promising branches more deeply. With his 9×9 Go program Crazy Stone, Coulom won a computer Go tournament (KGS). In the same year, Kocsis and Szepesvári of Hungary proposed UCT, a method for deciding which branch to try and how much. Ten years later, the AlphaGo paper (2016) cited these two as the MCTS literature, and combined MCTS with two neural networks (one predicting the next move and one evaluating positions).
3 sources and how they were checked
Where sources disagree:Which edition of the conference the Kocsis and Szepesvári paper was at: the AlphaGo reference list says “15th European Conference on Machine Learning”, while Springer's bibliographic record says ECML 2006 was the 17th (we only confirmed this in search results and did not open the original). We have not settled which is right.