My stuff
Learning From Noisy Singly-labeled Data

Supervised learning depends on annotated examples, which are taken to be the ground truth. But these labels often come from noisy crowdsourcing platforms, like Amazon Mechanical Turk. Practitioners typically collect multiple labels per example and aggregate the results to mitigate noise (the classic crowdsourcing problem). Given a fixed annotation budget and unlimited unlabeled data, redundant annotation comes at the expense of fewer labeled examples. This raises two fundamental questions: (1) How can we best learn from noisy workers? (2) How should we allocate our labeling budget to maximize the performance of a classi-fier? We propose a new algorithm for jointly modeling labels and worker quality from noisy crowd-sourced data. The alternating minimization proceeds in rounds, estimating worker quality from disagreement with the current model and then updating the model by optimizing a loss function that accounts for the current estimate of worker quality. Unlike previous approaches, even with only one annotation per example, our algorithm can estimate worker quality. We establish a generalization error bound for models learned with our algorithm and establish theoretically that it’s better to label many examples once (vs less multiply) when worker quality exceeds a threshold. Experiments conducted on both ImageNet (with simulated noisy workers) and MS-COCO (using the real crowdsourced labels) confirm our algorithm’s benefits.

Recent advances in supervised learning owe, in part, to the availability of large annotated datasets. For instance, the performance of modern image classifiers saturates only with millions of labeled examples. This poses an economic problem: Assembling such datasets typically requires the labor of human annotators. If we confined the labor pool to experts, this work might be prohibitively expensive. Therefore, most practitioners turn to crowdsourcing platforms such as Amazon Mechanical Turk (AMT), which connect employers with low-skilled workers who perform simple tasks, such as classifying images, at low cost.

Compared to experts, crowd-workers provide noisier annotations, possibly owing to high variation in worker skill; and a per-answer compensation structure that encourages rapid answers, even at the expense of accuracy. To address variation in worker skill, practitioners typically collect multiple independent labels for each training example from different workers. In practice, these labels are often aggregated by applying a simple majority vote. Academics have proposed many efficient algorithms for estimating the ground truth from noisy annotations. Research addressing the crowd-sourcing problem goes back to the early 1970s. Dawid & Skene (1979) proposed a probabilistic model to jointly estimate worker skills and ground truth labels and used expectation maximization (EM) to estimate the parameters. Whitehill et al. (2009); Welinder et al. (2010); Zhou et al. (2015) proposed generalizations of the Dawid-Skene model, e.g. by estimating the difficulty of each example.

Although the downstream goal of many crowdsourcing projects is to train supervised learning models, research in the two disciplines tends to proceed in isolation. Crowdsourcing research seldom accounts for the downstream utility of the produced annotations as training data in machine learning (ML) algorithms. And ML research seldom exploits the noisy labels collected from multiple human workers. A few recent papers use the original noisy labels and the corresponding worker identities together with the predictions of a supervised learning model trained on those same labels, to estimate the ground truth (Branson et al., 2017; Guan et al., 2017; Welinder et al., 2010). However, these papers do not realize the full potential of combining modeling and crowd-sourcing. In particular, they are unable to estimate worker qualities when there is only one label per training example.

This paper presents a new supervised learning algorithm that alternately models the labels and worker quality. The EM algorithm bootstraps itself in the following way: Given a trained model, the algorithm estimates worker qualities using the disagreement between workers and the current predictions of the learning algorithm. Given estimated worker qualities, our algorithm optimizes a suitably modified loss function. We show that accurate estimates of worker quality can be obtained even when only collecting one label per example provided that each worker labels sufficiently many examples. An accurate estimate of the worker qualities leads to learning a better model. This addresses a shortcoming of the prior work and overcomes a significant hurdle to achieving practical crowdsourcing without redundancy.

We give theoretical guarantees on the performance of our algorithm. We analyze the two alternating steps: (a) estimating worker qualities from disagreement with the model, (b) learning a model by optimizing the modified loss function. We obtain a bound on the accuracy of the estimated worker qualities and the generalization error of the model. Through the generalization error bound, we establish that it is better to label many examples once than to label less examples multiply when worker quality is above a threshold. Empirically, we verify our approach on several multi-class classification datasets: ImageNet and CIFAR10 (with simulated noisy workers), and MS-COCO (using the real noisy annotator labels). Our experiments validate that when the cost of obtaining unlabeled examples is negligible and the total annotation budget is fixed, it is best to collect a single label per training example for as many examples as possible. We emphasize that although this paper applies our approach to classification problems, the main ideas of the algorithm can be extended to other tasks in supervised learning.

The traditional crowdsourcing problem addresses the challenge of aggregating multiple noisy labels. A naive approach is to aggregate the labels based on majority voting. More sophisticated agreementbased algorithms jointly model worker skills and ground truth labels, estimating both using EM or similar techniques (Dawid & Skene, 1979; Jin & Ghahramani, 2003; Whitehill et al., 2009; Welinder et al., 2010; Zhou et al., 2012; Liu et al., 2012; Dalvi et al., 2013; Liu et al., 2012). Zhang et al. (2014) shows that the EM algorithm with spectral initialization achieves minimax optimal performance under the Dawid-Skene model. Karger et al. (2014) introduces a message-passing algorithm for estimating binary labels under the Dawid-Skene model, showing that it performs strictly better than majority voting when the number of labels per example exceeds some threshold. Similar observations are made by (Bragg et al., 2016). A primary criticism of EM-based approaches is that in practice, it’s rare to collect more than 3 to 5 labels per example; and with so little redundancy, the small gains achieved by EM over majority voting are not compelling to practitioners. In contrast, our algorithm performs well in the low-redundancy setting. Even with just one label per example, we can accurately estimate worker quality.

Several prior crowdsourcing papers incorporate the predictions of a supervised learning model, together with the noisy labels, to estimate the ground truth labels. Welinder et al. (2010) consider binary classification and frames the problem as a generative Bayesian model on the features of the examples and the labels. Branson et al. (2017) considers a generalization of the Dawid-Skene model and estimates its parameters using supervised learning in the loop. In particular, they consider a joint probability over observed image features, ground truth labels, and the worker labels and computes the maximum likelihood estimate of the true labels using alternating minimization. We also consider a joint probability model but it is significantly different from theirs as we assume that the optimal labeling function gives the ground truth labels. We maximize the joint likelihood using a variation of expectation maximization to learn the optimal labeling function and the true labels. Further, they train the supervised learning model using the intermediate predictions of the labels whereas we train the model by minimizing a weighted loss function where the weights are the intermediate posterior probability distribution of the labels. Moreover, with only one label per example, their algorithm fails and estimates all the workers to be equally good. They only consider binary classification, whereas we verify our algorithm on multi-class (ten classes) classification problem.

A rich body of work addresses human-in-loop annotation for computer vision tasks. However, these works assume that humans are experts, i.e., that they give noiseless annotations (Branson et al., 2010; Deng et al., 2013; Wah et al., 2011). We assume workers are unreliable and have varying skills. A recent work by Ratner et al. (2016) also proposes to use predictions of a supervised learning model to estimate the ground truth. However, their algorithm is significantly different than ours as it does not use iterative estimation technique, and their approach of incorporating worker quality parameters in the supervised learning model is different. Their theoretical results are limited to the linear classifiers.

Another line of work employs active learning, iteratively filtering out examples for which aggregated labels have high confidence and collect additional labels for the remaining examples (Whitehill et al., 2009; Welinder & Perona, 2010; Khetan & Oh, 2016). The underlying modeling assumption in these papers is that the questions have varying levels of difficulty. At each iteration, these approaches employ an EM-based algorithm to estimate the ground truth label of the remaining unclassified examples. For simplicity, our paper does not address example difficulties, but we could easily extend our model and algorithm to accommodate this complexity.

Several papers analyze whether repeated labeling is useful. Sheng et al. (2008) analyzed the effect of repeated labeling and showed that it depends upon the relative cost of getting an unlabeled example and the cost of labeling. Ipeirotis et al. (2014) shows that if worker quality is below a threshold then repeated labeling is useful, otherwise not. Lin et al. (2014a; 2016) argues that it also depends upon expressiveness of the classifier in addition to the factors considered by others. However, these works do not exploit predictions of the supervised learning algorithm to estimate the ground truth labels, and hence their findings do not extend to our methodology.

Another body of work that is relevant to our problem is learning with noisy labels where usual assumption is that all the labels are generated through the same noisy rate given their ground truth label. Recently Natarajan et al. (2013) proposed a generic unbiased loss function for binary clas-sification with noisy labels. They employed a modified loss function that can be expressed as a weighted sum of the original loss function, and gave theoretical bounds on the performance. However, their weights become unstably large when the noise rate is large, and hence the weights need to be tuned. Sukhbaatar et al. (2014); Jindal et al. (2016) learns noise rate as parameters of the model. A recent work by Guan et al. (2017) trains an individual softmax layer for each expert and then predicts their weighted sum where weights are also learned by the model. It is not scalable to crowdsourcing scenario where there are thousands of workers. There are works that aim to create noise-robust models (Joulin et al., 2016; Krause et al., 2016), but they are not relevant to our work.

Let D be the underlying true distribution generating pairs  (X, Y ) ∈ X × Kfrom which n i.i.d. samples  (X1, Y1), (X2, Y2), · · · , (Xn, Yn)are drawn, where K denotes the set of possible labels K := {1, 2, · · · , K}, and  X ⊆ Rddenotes the set of euclidean features. We denote the marginal distribution of Y by  {q1, q2, · · · , qK}, which is unknown to us. Consider a pool of m workers indexed by  1, 2, · · · , m. We use [m] to denote the set  {1, 2, · · · , m}. For each i-th sample  Xi, rworkers  {wij}j∈[r] ∈ [m]rare selected randomly, independent of the sample  Xi. Each selected worker provides a noisy label  Zijfor the sample  Xi, where the distribution of  Zijdepends on the selected worker and the true label  Yi. We call r the redundancy and, for simplicity, assume it to be the same for each sample. However, our algorithm can also be applied when redundancy varies across the samples. We use  Z(r)ito denote  {Zij}j∈[r], the set of r labels collected on the i-th example, and  w(r)ito denote  {wij}j∈[r].

Following Dawid & Skene (1979), we assume the probability that the a-th worker labels an item in class  k ∈ Kas class  s ∈ Kis independent of any particular chosen item, that is, it is a constant over i ∈ [n]. Let us denote this constant by  πks; by definition, �s∈K πks = 1for all  k ∈ K, and we call π(a) ∈ [0, 1]K×Kthe confusion matrix of the a-th worker. In particular, the distribution of Z is:


The diagonal entries of the confusion matrix correspond to the probabilities of correctly labeling an example. The off-diagonal entries represent the probability of mislabeling. We use  πto denote the collection of confusion matrices  {π(a)}a∈[m].

We assume nr workers  w1,1, w1,2, · · · , wn,rare selected uniformly at random from a pool of m workers with replacement and a batch of r workers are assigned to each of the examples  X1, X2, · · · , Xn. The corrupted labels along with the worker information (X1, Z(r)1 , w(r)1 ), · · · , (Xn, Z(r)n , w(r)n )are what the learning algorithm sees.

Let F be the hypothesis class, and  f ∈ F, f : X → RK, denote a vector valued predictor function. Let  ℓ(f(X), Y )denote a loss function. For a predictor f, its  ℓ-risk under D is defined as


Given the observed samples  (X1, Z(r)1 , w(r)1 ), · · · , (Xn, Z(r)n , w(r)n ), we want to learn a good pre- dictor function �f ∈ Fsuch that its risk under the true distribution  D, Rℓ,D( �f)is minimal. Having access to only noisy labels  Z(r)by workers  w(r), we compute �fas the one which minimizes a suitably modified loss function  ℓ�π,�q(f(X), Z(r), w(r)). Where  �πdenote an estimate of confusion matrix  π, and  �qan estimate of q, the prior distribution on Y . We define  ℓ�π,�qin the following section.

Assume that there exists a function  f ∗ ∈ Fsuch that  f ∗(Xi) = Yifor all  i ∈ [n]. Under the Dawid-Skene model (described in previous section), the joint likelihood of true labeling function f ∗(Xi)and observed labels  {Zij}i∈[n],j∈[r]as a function of confusion matrices of workers  πcan be written as


qk’s are the marginal distribution of the true labels  Yi’s. We estimate the worker confusion matrices πand the true labeling function  f ∗by maximizing the likelihood function  L(π; f ∗(X), Z). Observe that the likelihood function  L(π; f ∗(X), Z)is different than the standard likelihood function of Dawid-Skene model in that we replace each true hidden labels  Yiby  f ∗(Xi). Like the EM algorithm introduced in (Dawid & Skene, 1979), we propose ‘Model Bootstrapped EM’ (MBEM) to estimate confusion matrices  πand the true labeling function  f ∗. EM converges to the true confusion matrices and the true labels, given an appropriate spectral initialization of worker confusion matrices (Zhang et al., 2014). We show in Section 4.4 that MBEM converges under mild conditions when the worker quality is above a threshold and the number of training examples is sufficiently large. In the following two subsections, we motivate and explain our iterative algorithm to estimate the true labeling function  f ∗, given a good estimate of worker confusion matrices  πand vice-versa.



To begin, we ask, what is the optimal approach to learn the predictor function �fwhen for each worker we have  �π, a good estimation of the true confusion matrix  π, and  �q, an estimate of the prior? A recent paper, Natarajan et al. (2013) proposes minimizing an unbiased loss function specifically, a weighted sum of the original loss over each possible ground truth label. They provide weights for binary classification where each example is labeled by only one worker. Consider a worker with confusion matrix  π, where  πy > 1/2and  π−y > 1/2represent her probability of correctly labeling the examples belonging to class y and  −yrespectively. Then their weights are  π−y/(πy + π−y − 1)for class y and  −(1 − πy)/(πy + π−y − 1)for class  −y. It is evident that their weights become unstably large when the probabilities of correct classification  πyand  π−yare close to 1/2, limiting the method’s usefulness in practice. As explained below, for the same scenario, our weights would be  πy/(1+πy −π−y)for class y and  (1−π−y)/(1+πy −π−y)for class  −y. Inspired by their idea, we propose weighing the loss function according to the posterior distribution of the true label given the  Z(r)observed labels and an estimate of the confusion matrices of the worker who provided those labels. In particular, we define  ℓ�π,�qto be


If the observed label is uniformly random, then all weights are equal and the loss is identical for all predictor functions f. Absent noise, we recover the original loss function. Under the Dawid-Skene model, given the observed noisy labels  Z(r), an estimate of confusion matrices  �π, and an estimate of prior  �q, the posterior distribution of the true labels can be computed as follows:


where I[.] is the indicator function which takes value one if the identity inside it is true, otherwise zero. We give guarantees on the performance of the proposed loss function in Theorem 4.1. In practice, it is robust to noise level and significantly outperforms the unbiased loss function. Given ℓ�π,�q, we learn the predictor function �fby minimizing the empirical risk



The next question is: how do we get a good estimate  �πof the true confusion matrix  πfor each worker. If redundancy r is sufficiently large, we can employ the EM algorithm. However, in practical applications, redundancy is typically three or five. With so little redundancy, the standard applications of EM are of limited use. In this paper we look to transcend this problem, posing the question: Can we estimate confusion matrices of workers even when there is only one label per example? While this isn’t possible in the standard approach, we can overcome this obstacle by incorporating a supervised learning model into the process of assessing worker quality.

Under the Dawid-Skene model, the EM algorithm estimates the ground truth labels and the confusion matrices in the following way: It alternately fixes the ground truth labels and the confusion matrices by their estimates and and updates its estimate of the other by maximizing the likelihood of the observed labels. The alternating maximization begins by initializing the ground truth labels with a majority vote. With only 1 label per example, EM estimates that all the workers are perfect.

We propose using model predictions as estimates of the ground truth labels. Our model is initially trained on the majority vote of the labels. In particular, if the model prediction is  {ti}i∈[n], where ti ∈ K, then the maximum likelihood estimate of confusion matrices and the prior distribution is given below. For the a-th worker,  �π(a)ksfor  k, s ∈ K, and  �qkfor  k ∈ K, we have,


The estimate is effective when the hypothesis class F is expressive enough and the learner is robust to noise. Thus the model should, in general, have small training error on correctly labeled examples and large training error on wrongly labeled examples. Consider the case when there is only one label per example. The model will be trained on the raw noisy labels given by the workers. For simplicity, assume that each worker is either a hammer (always correct) or a spammer (chooses labels uniformly random). By comparing model predictions with the training labels, we can identify which workers are hammers and which are spammers, as long as each worker labels sufficiently many examples. We expect a hammer to agree with the model more often than a spammer.


Building upon the previous two ideas, we present ‘Model Bootstrapped EM’, an iterative algorithm for efficient learning from noisy labels with small redundancy. MBEM takes data, noisy labels, and the corresponding worker IDs, and returns the best predictor function �fin the hypothesis class F. In the first round, we compute the weights of the modified loss function  ℓ�π,�qby using the weighted majority vote. Then we obtain an estimate of the worker confusion matrices  �πusing the maximum likelihood estimator by taking the model predictions as the ground truth labels. In the second round, weights of the loss function are computed as the posterior probability distribution of the ground truth labels conditioned on the noisy labels and the estimate of the confusion matrices obtained in the previous round. In our experiments, only two rounds are required to achieve substantial improvements over baselines.




The following result gives guarantee on the excess risk for the learned predictor function �fin terms of the VC dimension of the hypothesis class F. Recall that risk of a function f w.r.t. loss function  ℓis defined to be  Rℓ,D(f) := E(X,Y )∼D [ℓ(f(X), Y )], Equation (2). We assume that the classification problem is binary, and the distribution q, prior on ground truth labels Y , is uniform and is known to us. We give guarantees on the excess risk of the predictor function �f, and accuracy of  �πestimated in the second round. For the purpose of analysis, we assume that fresh samples are used in each round for computing function �fand estimating  �π. In other words, we assume that �fand  �πare each computed using n/4 fresh samples in the first two rounds. We define  αand  βϵto capture the average worker quality. Here, we give their concise bound for a special case when all the workers are identical, and their confusion matrix is represented by a single parameter,  0 ≤ ρ < 1/2. Where πkk = 1 − ρ, and  πks = ρfor  k ̸= s. Each worker makes a mistake with probability  ρ. βϵ ≤(ρ + ϵ)r �ru=0�ru�(τ u + τ r−u)−1, where  τ := (ρ + ϵ)/(1 − ρ − ϵ). αfor this special case is  ρ. A general definition of  αand  βϵfor any confusion matrices  πis provided in the Appendix.

Theorem 4.1. Define N := nr to be the number of total annotations collected on n training examples with redundancy r. Suppose  minf∈F Rℓ,D(f) ≤ 1/4. For any hypothesis class F with a finite VC dimension V , and any  δ < 1, there exists a universal constant C such that if N is large enough and satisfies


then for binary classification with 0-1 loss function  ℓ, �fand  �πreturned by Algorithm 1 after T = 2 iterations satisfies


and  ∥�π(a) − π(a)∥∞ ≤ ϵ1for all  a ∈ [m], with probability at least  1 − δ. Where  ϵ := 24γ +28�m log(26mδ)/N, and  γ := minf∈F Rℓ,D(f) + C(√V +�log(1/δ))/((1 − 2α)�N/r). ϵ1is defined to be  ϵwith  αin it replaced by  βϵ.


The price we pay in generalization error bound on �fis  (1 − 2βϵ). Note that, when n is large,  ϵgoes to zero, and  βϵ ≤ 2ρ(1 − ρ), for r = 1.

If  minf∈F Rℓ,D(f)is sufficiently small, VC dimension is finite, and  ρis bounded away from 1/2 then for n = O(m log(m)/r), we get  ϵ1to be sufficiently small. Therefore, for any redundancy r, error in confusion matrix estimation is small when the number of training examples is sufficiently large. Hence, for N large enough, using Equation (9) and the bound on  βϵ, we get that for fixed total annotation budget, the optimal choice of redundancy r is 1 when the worker quality  (1 − ρ)is above a threshold. In particular, if  (1−ρ) ≥ 0.825then label once is the optimal strategy. However, in experiments we observe that with our algorithm the choice of r = 1 is optimal even for much smaller values of worker quality.

We experimentally investigate our algorithm, MBEM, on multiple large datasets. On CIFAR-10 (Krizhevsky & Hinton, 2009) and ImageNet (Deng et al., 2009), we draw noisy labels from synthetic worker models. We confirm our results on multiple worker models. On the MS-COCO dataset (Lin et al., 2014b), we accessed the real raw data that was used to produce this annotation. We compare MBEM against the following baselines:

MV: First aggregate labels by performing a majority vote, then train the model. weighted-MV: Model learned using weighted loss function with weights set by majority vote. EM: First aggregate labels using EM. Then train model in the standard fashion. (Dawid & Skene, 1979)

weighted-EM: Model learned using weighted loss function with weights set by standard EM. oracle weighted EM: This model is learned by minimizing  ℓπ, using the true confusion matrices. oracle correctly labeled: This baseline is trained using the standard loss function  ℓbut only using those training examples for which at least one of the r workers has given the true label.

Note that oracle models cannot be deployed in practice. We show them to build understanding only. In the plots, the dashed lines correspond to MV and EM algorithm. The black dashed-dotted line shows generalization error if the model is trained using ground truth labels on all the training examples. For experiments with synthetic noisy workers, we consider two models of worker skill:

hammer-spammer: Each worker is either a hammer (always correct) with probability  γor a spammer (chooses labels uniformly at random).

class-wise hammer-spammer: Each worker can be a hammer for some subset of classes and a spammer for the others. The confusion matrix in this case has two types of rows: (a) hammer class: row with all off-diagonal elements being 0. (b) spammer class: row with all elements being 1/|K|. A worker is a hammer for any class  k ∈ Kwith probability  γ.

We sample m confusion matrices  {π(a)}a∈[m]according to the given worker skill distribution for a given  γ. We assign r workers uniformly at random to each example. Given the ground truth labels, we generate noisy labels according to the probabilities given in a worker’s confusion matrix, using Equation (1). While our synthetic workers are sampled from these specific worker skill models, our algorithms do not use this information to estimate the confusion matrices. A Python implementation of the MBEM algorithm is available for download at https://github.com/khetan2/MBEM.

CIFAR-10 This dataset has a total of 60K images belonging to 10 different classes where each class is represented by an equal number of images. We use 50K images for training the model and 10K images for testing. We use the ground truth labels to generate noisy labels from synthetic workers. We choose m = 100, and for each worker, sample confusion matrix of size  10 × 10according to the worker skill distribution. All our experiments are carried out with a 20-layer ResNet


Figure 1: Plots for CIFAR-10. Line colors- black: oracle correctly labeled, red: oracle weighted EM, blue: MBEM, green: weighted EM, yellow: weighted MV.

which achieves an accuracy of 91.5%. With the larger ResNet-200, we can obtain a higher accuracy of 93.5% but to save training time we restrict our attention to ResNet-20. We run MBEM 1 for T = 2 rounds. We assume that the prior distribution  �qis uniform. We report mean accuracy of 5 runs and its standard error for all the experiments.

Figure 1 shows plots for CIFAR-10 dataset under various settings. The three plots in the first row correspond to “hammer-spammer” worker skill distribution and the plots in the second row correspond to “class-wise hammer-spammer” distribution. In the first plot, we fix redundancy r = 1, and plot generalization error of the model for varying hammer probability  γ. MBEM significantly outperforms all baselines and closely matches the Oracle weighted EM. This implies MBEM recovers worker confusion matrices accurately even when we have only one label per example. When there is only one label per example, MV, weighted-MV, EM, and weighted-EM all reduce learning with the standard loss function  ℓ.

In the second plot, we fix hammer probability  γ = 0.2, and vary redundancy r. This plot shows that weighted-MV and weighted-EM perform significantly better than MV and EM and confirms that our approach of weighing the loss function with posterior probability is effective. MBEM performs much better than weighted-EM at small redundancy, demonstrating the effect of our bootstrapping idea. However, when redundancy is large, EM works as good as MBEM.

In the third plot, we show that when the total annotation budget is fixed, it is optimal to collect one label per example for as many examples as possible. We fixed hammer probability  γ = 0.2. Here, when redundancy is increased from 1 to 2, the number of of available training examples is reduced by 50%, and so on. Performance of weighted-EM improves when redundancy is increased from 1 to 5, showing that with the standard EM algorithm it might be better to collect redundant annotations for fewer example (as it leads to better estimation of worker qualities) than to singly annotate more examples. However, MBEM always performs better than the standard EM algorithm, achieving lowest generalization error with many singly annotated examples. Unlike standard EM, MBEM can estimate worker qualities even with singly annotated examples by comparing them with model predictions. This corroborates our theoretical result that label-once is the optimal strategy when worker quality is above a threshold. The plots corresponding to class-wise hammer-spammer workers follow the same trend. Estimation of confusion matrices in this setting is difficult and hence the gap between MBEM and the baselines is less pronounced.

ImageNet The ImageNet-1K dataset contains 1.2M training examples and 50K validation examples. We divide test set in two parts: 10K for validation and 40K for test. Each example belongs to one of the possible 1000 classes. We implement our algorithms using a ResNet-18 that achieves top-


Figure 2: Plots for ImageNet. Solid lines represent top-5 error, dashed-lines represent top-1 error. Line colors- blue: MBEM, green: weighted majority vote, yellow: majority vote


Figure 3: Results on raw MS-COCO annotations.

1 accuracy of 69.5% and top-5 accuracy of 89% on ground truth labels. We use m = 1000 simulated workers. Although in general, a worker can mislabel an example to one of the 1000 possible classes, our simulated workers mislabel an example to only one of the 10 possible classes. This captures the intuition that even with a larger number of classes, perhaps only a small number are easily confused for each other. Therefore, each workers’ confusion matrix is of size  10 × 10. Note that without this assumption, there is little hope of estimating a  1000 × 1000confusion matrix for each worker by collecting only approximately 1200 noisy labels from a worker. The rest of the settings are the same as in our CIFAR-10 experiments. In Figure 2, we fix total annotation budget to be 1.2M and vary redundancy from 1 to 9. When redundancy is 9, we have only (1.2/9)M training examples, each labeled by 9 workers. MBEM outperforms baselines in each of the plots, achieving the minimum generalization error with many singly annotated training examples.

MS-COCO These experiments use the real raw annotations collected when MS-COCO was crowdsourced. Each image in the dataset has multiple objects (approximately 3 on average). For validation set images (out of 40K), labels were collected from 9 workers on average. Each worker marks which out of the 80 possible objects are present. However, on many examples workers disagree. These annotations were collected to label bounding boxes but we ask a different question: what is the best way to learn a model to perform multi-object classification, using these noisy annotations. We use 35K images for training the model and 1K for validation and 4K for testing. We use raw noisy annotations for training the model and the final MS-COCO annotations as the ground truth for the validation and test set. We use ResNet-98 deep learning model and train independent binary classifier for each of the 80 object classes. Table in Figure 3 shows generalization F1 score of four different algorithms: majority vote, EM, MBEM using all 9 noisy annotations on each of the training examples, and a model trained using the ground truth labels. MBEM performs significantly better than the standard majority vote and slightly improves over EM. In the plot, we fix the total annotation budget to 35K. We vary redundancy from 1 to 7, and accordingly reduce the number of training examples to keep the total number of annotations fixed. When redundancy is r < 9 we select uniformly at random r of the original 9 noisy annotations. Again, we find it best to singly annotate as many examples as possible when the total annotation budget is fixed. MBEM significantly outperforms majority voting and EM at small redundancy.

We introduced a new algorithm for learning from noisy crowd workers. We also presented a new theoretical and empirical demonstration of the insight that when examples are cheap and annotations expensive, it’s better to label many examples once than to label few multiply when worker quality is above a threshold. Many avenues seem ripe for future work. We are especially keen to incorporate our approach into active query schemes, choosing not only which examples to annotate, but which annotator to route them to based on our models current knowledge of both the data and the worker confusion matrices.

Jonathan Bragg, Daniel S Weld, et al. Optimal testing for crowd workers. In Proceedings of the 2016 International Conference on Autonomous Agents & Multiagent Systems, pp. 966–974. International Foundation for Autonomous Agents and Multiagent Systems, 2016.

Steve Branson, Catherine Wah, Florian Schroff, Boris Babenko, Peter Welinder, Pietro Perona, and Serge Belongie. Visual recognition with humans in the loop. Computer Vision–ECCV 2010, pp. 438–451, 2010.

Steve Branson, Grant Van Horn, and Pietro Perona. Lean crowdsourcing: Combining humans and machines in an online system. In Proceedings of the IEEE Conference on Computer Vision and Pattern Recognition, pp. 7474–7483, 2017.

Nilesh Dalvi, Anirban Dasgupta, Ravi Kumar, and Vibhor Rastogi. Aggregating crowdsourced binary ratings. In Proceedings of the 22nd international conference on World Wide Web, pp. 285–294. ACM, 2013.

Alexander Philip Dawid and Allan M Skene. Maximum likelihood estimation of observer error-rates using the em algorithm. Applied statistics, pp. 20–28, 1979.

J. Deng, W. Dong, R. Socher, L.-J. Li, K. Li, and L. Fei-Fei. ImageNet: A Large-Scale Hierarchical Image Database. In CVPR09, 2009.

Jia Deng, Jonathan Krause, and Li Fei-Fei. Fine-grained crowdsourcing for fine-grained recognition. In Proceedings of the IEEE conference on computer vision and pattern recognition, pp. 580–587, 2013.

Melody Y Guan, Varun Gulshan, Andrew M Dai, and Geoffrey E Hinton. Who said what: Modeling individual labelers improves classification. arXiv preprint arXiv:1703.08774, 2017.

Panagiotis G Ipeirotis, Foster Provost, Victor S Sheng, and Jing Wang. Repeated labeling using multiple noisy labelers. Data Mining and Knowledge Discovery, 28(2):402–441, 2014.

Rong Jin and Zoubin Ghahramani. Learning with multiple labels. In Advances in neural information processing systems, pp. 921–928, 2003.

Ishan Jindal, Matthew Nokleby, and Xuewen Chen. Learning deep networks from noisy labels with dropout regularization. In Data Mining (ICDM), 2016 IEEE 16th International Conference on, pp. 967–972. IEEE, 2016.

Armand Joulin, Laurens van der Maaten, Allan Jabri, and Nicolas Vasilache. Learning visual fea- tures from large weakly supervised data. In European Conference on Computer Vision, pp. 67–84. Springer, 2016.

David R Karger, Sewoong Oh, and Devavrat Shah. Budget-optimal task allocation for reliable crowdsourcing systems. Operations Research, 62(1):1–24, 2014.

Ashish Khetan and Sewoong Oh. Achieving budget-optimality with adaptive schemes in crowd- sourcing. In Advances in Neural Information Processing Systems, pp. 4844–4852, 2016.

Jonathan Krause, Benjamin Sapp, Andrew Howard, Howard Zhou, Alexander Toshev, Tom Duerig, James Philbin, and Li Fei-Fei. The unreasonable effectiveness of noisy data for fine-grained recognition. In European Conference on Computer Vision, pp. 301–320. Springer, 2016.

Alex Krizhevsky and Geoffrey Hinton. Learning multiple layers of features from tiny images. 2009.

Christopher H Lin, Daniel S Weld, et al. To re (label), or not to re (label). In Second AAAI conference on human computation and crowdsourcing, 2014a.

Christopher H Lin, M Mausam, and Daniel S Weld. Re-active learning: Active learning with rela- beling. In AAAI, pp. 1845–1852, 2016.

Tsung-Yi Lin, Michael Maire, Serge Belongie, James Hays, Pietro Perona, Deva Ramanan, Piotr Doll´ar, and C Lawrence Zitnick. Microsoft coco: Common objects in context. In European conference on computer vision, pp. 740–755. Springer, 2014b.

Qiang Liu, Jian Peng, and Alexander T Ihler. Variational inference for crowdsourcing. In Advances in neural information processing systems, pp. 692–700, 2012.

Nagarajan Natarajan, Inderjit S Dhillon, Pradeep K Ravikumar, and Ambuj Tewari. Learning with noisy labels. In Advances in neural information processing systems, pp. 1196–1204, 2013.

Alexander J Ratner, Christopher M De Sa, Sen Wu, Daniel Selsam, and Christopher R´e. Data pro- gramming: Creating large training sets, quickly. In Advances in Neural Information Processing Systems, pp. 3567–3575, 2016.

Victor S Sheng, Foster Provost, and Panagiotis G Ipeirotis. Get another label? improving data quality and data mining using multiple, noisy labelers. In Proceedings of the 14th ACM SIGKDD international conference on Knowledge discovery and data mining, pp. 614–622. ACM, 2008.

Sainbayar Sukhbaatar, Joan Bruna, Manohar Paluri, Lubomir Bourdev, and Rob Fergus. Training convolutional networks with noisy labels. arXiv preprint arXiv:1406.2080, 2014.

Catherine Wah, Steve Branson, Pietro Perona, and Serge Belongie. Multiclass recognition and part localization with humans in the loop. In Computer Vision (ICCV), 2011 IEEE International Conference on, pp. 2524–2531. IEEE, 2011.

Peter Welinder and Pietro Perona. Online crowdsourcing: rating annotators and obtaining cost- effective labels. In Computer Vision and Pattern Recognition Workshops (CVPRW), 2010 IEEE Computer Society Conference on, pp. 25–32. IEEE, 2010.

Peter Welinder, Steve Branson, Pietro Perona, and Serge J Belongie. The multidimensional wisdom of crowds. In Advances in neural information processing systems, pp. 2424–2432, 2010.

Jacob Whitehill, Ting-fan Wu, Jacob Bergsma, Javier R Movellan, and Paul L Ruvolo. Whose vote should count more: Optimal integration of labels from labelers of unknown expertise. In Advances in neural information processing systems, pp. 2035–2043, 2009.

Yuchen Zhang, Xi Chen, Denny Zhou, and Michael I Jordan. Spectral methods meet em: A provably optimal algorithm for crowdsourcing. In Advances in neural information processing systems, pp. 1260–1268, 2014.

Dengyong Zhou, Qiang Liu, John C Platt, Christopher Meek, and Nihar B Shah. Regularized minimax conditional entropy for crowdsourcing. arXiv preprint arXiv:1503.07240, 2015.

Denny Zhou, Sumit Basu, Yi Mao, and John C Platt. Learning from the wisdom of crowds by minimax entropy. In Advances in Neural Information Processing Systems, pp. 2195–2203, 2012.


Assuming the prior on Y , distribution q, to be uniform, we change the notation for the modified loss function  ℓ�π,�qto  ℓ�π. Observe that for binary classification,  Z(r) ∈ {±1}r. Let  ρdenote the posterior distribution of Y , Equation (5), when q is uniform. Let  τdenote the probability of observing an instance of  Z(r)as a function of the latent true confusion matrices  π, conditioned on the ground truth label Y = y.


Let W denote the uniform distribution over a pool of m workers, from which nr workers are selected i.i.d. with replacement, and a batch of r workers are assigned to each example  Xi. We define the following quantities which play an important role in our analysis.


For any given  �πwith  |�π(a)ks − π(a)ks | ≤ ϵ, for all  a ∈ [m], k, s ∈ K, we can compute  βϵfrom the definition of  β�πsuch that  β�π ≤ βϵ. For the special case described in Section 4.4, we have the following bound on  βϵ.


It can easily be checked that  βϵ ≤ (ρ + ϵ)r �⌈r/2⌉u=0 �ru�(1 − ρ − ϵ)u(ρ + ϵ)−u.

We present two lemma that analyze the two alternative steps of our algorithm. The following lemma gives a bound on the excess risk of function �flearnt by minimizing the modified loss function  ℓ�π.

Lemma A.1. Under the assumptions of Theorem 4.1, the excess risk of function �fin Equation (6), computed with posterior distribution  P�π(5) using n training examples is bounded by


with probability at least  1 − δ1, where C is a universal constant. When  P�πis computed using majority vote, while initializing the iterative Algorithm 1, the above bound holds with  β�πreplaced by  α.

The following lemma gives an  ℓ∞norm bound on confusion matrices  �πestimated using model prediction �f(X)as the ground truth labels. In the analysis, we assume fresh samples are used for estimating confusion matrices in step 3, Algorithm 1. Therefore the function �fis independent of the samples  Xi’s on which  �πis estimated. Let K = |K|.

Lemma A.2. Under the assumptions of Theorem 4.1,  ℓ∞error in estimated confusion matrices  �πas computed in Equation (7), using n samples and a predictor function �fwith risk  Rℓ,D ≤ δ, is bounded by


with probability at least  1 − δ1.

First we apply Lemma A.1 with  P�πcomputed using majority vote. We get a bound on the risk of function �fcomputed in the first round. With this �f, we apply Lemma A.2. When n is sufficiently large such that Equation (8) holds, the denominator in Equation (18),  1/K − δ −8�m log(4mK2/δ1)/(nr) ≥ 1/8. Therefore, in the first round, the error in confusion matrix estimation is bounded by  ϵ, which is defined in the Theorem.

For the second round: we apply Lemma A.1 with  P�πcomputed as the posterior distribution (5). Where  ℓ∞error in  �πis bounded by  ϵ. This gives the desired bound in (9). With this �f, we apply Lemma A.2 and obtain  ℓ∞error in  �πbounded by  ϵ1, which is defined in the Theorem.

For the given probability of error  δin the Theorem, we chose  δ1in both the lemma to be  δ/4such that with union bound we get the desired probability of  δ.


Let  f ∗ := arg minf∈F Rℓ,D(f). Let’s denote the distribution of  (X, Z(r), w(r))by  DW,π,r. For ease of notation, we denote  DW,π,rby  Dπ. Similar to  Rℓ,D, risk of decision function f with respect to the modified loss function  ℓ�πis characterized by the following quantities:

1. �π-risk under  Dπ: Rℓ�π,Dπ(f) := E(X,Z(r),w(r))∼Dπ�ℓ�π(f(X), Z(r), w(r))�.

2. Empirical  ℓ�π-risk on samples: �Rℓ�π,Dπ(f) := 1n�ni=1 ℓ�π(f(Xi), Z(r)i , w(r)i ).

With the above definitions, we have the following,


where (19) follows from Equation (24). (20) follows from the fact that �fis the minimizer of �Rℓ�π,Dπas computed in (6). (21) follows from the basic excess-risk bound. V is the VC dimension of hypothesis class F, and C is a universal constant.

Following shows the inequality used in Equation (19). For binary classification, we denote the two classes by  Y, −Y.


where (22) follows from Equation (26). (23) follows from the fact that for 0-1 loss function ℓ(f(X), Y ) + ℓ(f(X), −Y ) = 1. (24) follows from the definition of  β�πdefined in Equation (12). When  ℓ�πis computed using weighted majority vote of the workers then (24) holds with  β�πreplaced by  α. αis defined in (14).

Following shows the equality used in Equation (22). Using the notations  ρ�πand  τπ, in the following, for any function  f ∈ F, we compute the excess risk due to the unbiasedness of the modified loss function  ℓ�π.


where  β�π(Y )is defined in (11). Where (25) follows from the definition of  ℓ�πgiven in Equation (4). Observe that when  ℓ�πis computed using weighted majority vote of the workers then Equation (26) holds with  β�π(Y )replaced by  α(y). α(y)is defined in (13).


Recall that we have


Let |K| = K. Define, for fixed  a ∈ [m], and  k, s ∈ K,


Note that A, B, C, D, E depend upon  a ∈ [m], k, s ∈ K. However, for ease of notations, we have not included the subscripts. We have,


Now, we have,


We have,


Observe that C is a sum of nr i.i.d. Bernoulli random variables with mean  δ/m. Using Chernoff bound we get that


for all  a ∈ [m], and  k ∈ Kwith probability at least  1 − δ1. Similarly, D is a sum of nr i.i.d. Bernoulli random variables with mean  πks/(mk). Again, using Chernoff bound we get that


for all  a ∈ [m], k, s ∈ Kwith probability at least  1 − δ1. From the bound on  |D − ¯A|, it follows that


Collecting Equations (33)-(38), we have for all  a ∈ [m], k, s ∈ K


with probability at least  1 − 2δ1.

Designed for Accessibility and to further Open Science