A Unified Approach to Quantifying Algorithmic Unfairness: Measuring Individual & Group Unfairness via Inequality Indices

Till Speicher, Hoda Heidari, Nina Grgic-Hlaca, Krishna P. Gummadi, Adish Singla, Adrian Weller, Muhammad Bilal Zafar

Introduction

As algorithmic decision making systems are increasingly used in life-affecting scenarios such as criminal risk prediction (Angwin et al., 2016; Berk et al., 2017) and credit risk assessments (Furletti, 2002), concerns have risen about the potential unfairness of these decisions to certain social groups or individuals (Barocas and Selbst, 2016; Podesta et al., 2014; Romei and Ruggieri, 2014). In response, a number of recent works have proposed learning mechanisms for fair decision making by imposing additional constraints or conditions (Feldman et al., 2015; Zafar et al., 2017b, a; Hardt et al., 2016; Dwork et al., 2012; Joseph et al., 2016).

In this paper, we focus on a simple yet foundational question about unfairness of algorithms: Given two unfair algorithms, how should we determine which of the two is more unfair? Prior works on algorithmic fairness largely focus on formally defining conditions for fairness, but do not precisely define suitable measures for unfairness. That is, they can answer the binary question: is an algorithm fair or unfair?, but do not have a principled way to answer the nuanced question: if an algorithm is unfair, how unfair is it?

Figure 1 illustrates the questions we seek to answer through an example of two binary classifiers C1\mathcal{C}_{1} and C2\mathcal{C}_{2}, whose decisions affect 10 individuals belonging to 3 different groups. The figure shows that both C1\mathcal{C}_{1} and C2\mathcal{C}_{2} yield unequal false positive/negative rates across the 3 groups and are thus unfair at the level of groups—which set of unequal false positive/negative rates (C1\mathcal{C}_{1}’s or C2\mathcal{C}_{2}’s) are more unfair? Similarly, both C1\mathcal{C}_{1} and C2\mathcal{C}_{2} violate our individual-level fairness condition for “treating individuals deserving similar outcomes similarly”, but do so in different ways—whose violation of our individual fairness condition is more unfair?

We argue that how we address the unfairness measurement question has significant practical consequences. First, several studies have observed that satisfying multiple fairness conditions at the same time is infeasible (Kleinberg et al., 2017; Chouldechova, 2016; Corbett-Davies et al., 2017; Kearns et al., 2017). Hence in practice, designers often need to select the least unfair algorithm from a feasible set of unfair algorithms. Second, when training fair learning models, practitioners face a tradeoff between accuracy and fairness (Flores et al., 2016; Kleinberg et al., 2017). These tradeoffs rely on model-specific fairness measures (i.e., proxies chosen for computational tractability) that do not generalize across different models. Consequently, they cannot be used to compare accuracy-unfairness tradeoffs of models trained using different fair learning algorithms. Finally, designers of fair learning models make a number of ad hoc or implicit choices about fairness measures without explicit justification; for instance, it is unclear why in many previous works (Zafar et al., 2017b, a; Hardt et al., 2016; Dwork et al., 2012), the relative sizes of the groups in the population are not considered in estimating unfairness—even though these quantities matter when estimating accuracy.

In this paper, we propose to quantify unfairness using inequality indices that have been extensively studied in economics and social welfare (Atkinson, 1970; Cowell and Kuga, 1981; Kakwani, 1980). Traditionally, inequality indices such as Coefficient of Variation (Abdi, 2010), Gini (Gini, 1912; Bellù and Liberati, 2006), Atkinson (Atkinson, 1970), Hoover (Long and Nucci, 1997), and Theil (Theil, 1967), have been proposed to quantify how unequally incomes are distributed across individuals and groups in a population. Our interest in using these indices is rooted in the well-justified axiomatic basis for their designs. Specifically, we argue that many axioms satisfied by inequality indices such as anonymity, population invariance, progressive transfer preference, and subgroup decomposability are appealing properties for unfairness measures to satisfy. Thus, inequality indices are naturally well-suited as measures for algorithmic unfairness.

Our core idea is to use existing inequality indices in order to measure how unequally the outcomes of an algorithm benefit different individuals or groups in a population. This requires us to define a benefit function that maps the algorithmic output for each individual to a non-negative real number. By adapting the benefit function according to the desired fairness condition, we show that inequality indices can be applied generally to quantify unfairness across all the proposed fairness conditions shown in Figure 2. Since we quantify inequality of algorithmic outcomes, our measure is independent of the specifics of any learning model and can be used to compare unfairness of different algorithms.

We consider a family of inequality indices called generalized entropy indices, which includes Coefficient of Variation and Theil index as special cases. Generalized entropy indices have a useful property called subgroup decomposability. For any division of the population into a set of non-overlapping groups, the property guarantees that our unfairness measure over the entire population can be decomposed as the sum of a between-group unfairness component (computed imagining that all individuals in a group receive the group’s mean benefit) and a within-group unfairness component (computed as a weighted sum of inequality in benefits received by individuals within each group). Thus, inequality indices not only offer a unifying approach to quantifying unfairness at the levels of both individuals and groups, but they also reveal previously overlooked tradeoffs between individual-level and group-level fairness.

Further, the decomposition enables us to: (i) quantify how unfair an algorithm is along various sensitive attribute-based groups within a population (e.g., groups based on race, gender or age) and (ii) account for the “gerrymandered” unfairness affecting structured subgroups constructed from “intersecting” the sensitive attribute-groups (e.g., groups like young white women or old black men) (Kearns et al., 2017). Our empirical evaluations show that existing fair learning methods (Zafar et al., 2017a; Hardt et al., 2016), while successful in eliminating between-group unfairness, (a) may be targeting only a small fraction of the overall unfairness in the decision making algorithms and (b) can result in an increase in within-group unfairness, which paradoxically can lead to training algorithms whose overall unfairness is worse than those trained using traditional learning methods.

To summarize the contributions of this paper: (i) we propose inequality-indices based unfairness measures that offer a justified and generalizable framework to compare the fairness of a variety of algorithmic predictors against one another, (ii) we theoretically characterize and empirically illustrate the tradeoffs between individual fairness when measured using inequality indices and the prediction accuracy, and (iii) we study the relationship between individual- and group-level unfairness, showing that recently proposed learning mechanisms for mitigating (between-)group unfairness can lead to high within-group unfairness and consequently, high individual unfairness.

Measuring Algorithmic Unfairness via Inequality Indices

We first formally describe the setup of a fairness-aware machine learning task; then proceed to show that by defining an appropriate benefit function, existing inequality indices can be applied across the board to quantify algorithmic unfairness. We describe important properties (axioms) which we suggest a reasonable measure of algorithmic unfairness must satisfy. We end this section by comparing our proposed approach with previous work.

We assume certain features (e.g., gender or race) are considered sensitive. Sensitive features specify an individual’s membership in socially salient groups (e.g., women or African-Americans). For simplicity of exposition, we assume there is just one sensitive feature. However, the discussion can be extended to account for multiple sensitive features. We denote the sensitive feature for each individual ii as zi∈Z={1,2,…,K}z_{i}\in\mathcal{Z}=\{1,2,\ldots,K\}. Note that ziz_{i} may or may not be part of the feature vector xi\mathbf{x}_{i}. One can define partitions of the dataset D\mathcal{D} based on the sensitive feature, that is, Dz={(xi,yi) ∣ zi=z}\mathcal{D}_{z}=\{(\mathbf{x}_{i},y_{i})\ |\ z_{i}=z\}. We refer to each partition Dz\mathcal{D}_{z} of the data as a sensitive feature group.

2. Unfairness as Inequality in Benefits

The core idea of our proposal is to quantify the unfairness of an algorithm by measuring how unequally the outcomes of the algorithm benefit different individuals or groups in a population. While intuitive, our proposal raises two key questions: (i) how should we map algorithmic predictions received by individuals or groups to benefits? and (ii) given a set of benefits received by individuals or groups, how should we quantify inequality in the benefit distribution? We now tackle the first question, related to defining a benefit function for an individual given an outcome. In Section 2.3, we propose inequality indices as the answer to the second question.

Our choice of the benefit function will be dictated by the type of fairness notion we wish to apply on the task at hand. Figure 2 summarizes the different fairness notions that have been defined in prior works and their corresponding benefit functions. We now explain the choice of our benefit functions for the different fairness notions in the context of binary classification. Formally, let yi∈Y={0,1}y_{i}\in\mathcal{Y}=\{0,1\} indicate the true label for individual ii. We assume that labels in the training data reflect ground truth, and thus, yiy_{i} is the label deserved by individual ii. Let y^i∈{0,1}\hat{y}_{i}\in\{0,1\} be the label the algorithm assigns to individual ii.

Intuitively, the algorithmic benefit an individual ii receives, bib_{i}, should capture the desirability of outcome y^i\hat{y}_{i} for the individual. The desirability of an individual’s outcome may be determined taking into account the individual’s own preferences or the broader societal good. For instance, consider the criminal risk prediction example, where the positive label (y^=1\hat{y}=1) indicates a low risk of criminal behavior and the negative label (y^=0\hat{y}=0) indicates a high risk of criminal behavior. An individual defendant would clearly prefer the former outcome over the latter. However, from a social good perspective, accurate outcomes (y^=y\hat{y}=y) would be more desirable than inaccurate outcomes. Furthermore, amongst the inaccurate outcomes, one might wish to distinguish between the desirability of false positives (where a high risk person is released) and false negatives (where an low risk person is withheld).

In our binary classification scenario, where all outcomes can be decomposed into true positives (y^=1,y=1\hat{y}=1,y=1), true negatives (y^=0,y=0\hat{y}=0,y=0), false positives (y^=1,y=0\hat{y}=1,y=0), and false negatives (y^=0,y=1\hat{y}=0,y=1), the choice of our benefit function crucially determines the relative desirability of these different types of outcomes and captures different notions of fairness. For instance, the notion of parity mistreatment considers accurate outcomes as more desirable than inaccurate ones – so we choose a benefit function that assigns higher value (bi=1b_{i}=1) to true positives and true negatives and a lower value (bi=0b_{i}=0) to false positives and false negatives. In contrast, the notion of parity impact considers a positive label outcome as more desirable than a negative label outcome – so we adapt the benefit function to assign higher value (bi=1b_{i}=1) to true positives and false positives and a lower value (bi=0b_{i}=0) to true negatives and false negatives. To capture group fairness, once we define the benefits for all individuals, b=(b1,⋯ ,bn)\mathbf{b}=(b_{1},\cdots,b_{n}), we can define the benefit for a subset/group gg of the population, denoted by μg\mu_{g}, as the mean value of the benefits received by individuals in the group: μg=1∣g∣∑i∈gbi\mu_{g}=\frac{1}{|g|}\sum_{i\in g}b_{i}.

To capture individual fairness, we propose defining the benefit function of an individual ii as the discrepancy between ii’s preference for the outcome ii truly deserves (i.e., yiy_{i}), and ii’s preference for the outcome the learning algorithm assigns (i.e., y^i\hat{y}_{i}). As an illustration, in this work we consider a benefit function that assigns the highest value (bi=2b_{i}=2) for false positives (i.e., individuals that receive the advantageous positive label undeservedly), moderate values (bi=1b_{i}=1) for true positives and true negatives (i.e., individuals that receive the labels they deserve) and lowest value (bi=0b_{i}=0) for false negatives (i.e., individuals that receive the disadvantageous negative label despite deserving the positive label). More precisely, we compute the benefit for individual ii as follows:

We make two observations about the values of the benefit functions for different types of outcomes. First, while different fairness notions specify a preference ordering for different types of outcomes (i.e., true positives, false positives, true negatives, and false negatives), the absolute benefit values could be specified differently. The choice of benefit values would depend on the context and task at hand and the difficulty of determining them may vary in practice. Second, as many existing measures of inequality in benefits are limited to handling non-negative values, we need to ensure that bi≥0b_{i}\geq 0 for i=1,⋯ ,ni=1,\cdots,n and that there exists j∈[n]j\in[n] such that bj>0b_{j}>0.

Our proposal is to measure the overall individual-level unfairness of an algorithm by plugging bib_{i}’s (as defined above) into an existing inequality index (such as generalized entropy—to be defined shortly). Throughout the rest of the paper, we will use the terms “overall unfairness” and “individual unfairness” interchangeably, to refer to our proposed measure. Our approach can be further generalized to measuring (un)fairness beyond supervised learning tasks (e.g. for unsupervised tasks, such as clustering or ranking)—this only requires the specification a proper notion of benefit for individuals given their relative outcomes within the population. We leave a careful exploration of this direction for future, and focus on supervised learning tasks in the current work.

Next, we discuss how we can generally quantify the unfairness of an algorithm as the degree to which it distributes benefit unequally across individuals using inequality indices.

3. Axioms for Measuring Inequality

Borrowing insights from the rich body of work on the axiomatic characterization of inequality indices in economics and social science (Atkinson, 1970; Sen, 1973; Kolm, 1976a, b; Kakwani, 1980; Cowell and Kuga, 1981; Litchfield, 1999), we argue that many axioms satisfied by inequality indices are appealing properties for measures of algorithmic unfairness. Therefore, inequality indices are naturally well-suited as measures for algorithmic unfairness. In this section, we briefly overview these axioms.

Anonymity: The measure does not depend on any characteristics of the individuals other than their benefit, and is independent of who earns each level of benefit. Formally:

where (b(1),b(2),…,b(n))(b_{(1)},b_{(2)},\dots,b_{(n)}) is the benefit vector (b1,b2,⋯ ,bn)(b_{1},b_{2},\cdots,b_{n}) sorted in ascending order.

Transfer principle: Transferring benefit from a high-benefit to a low-benefit individual must decrease inequality. More precisely for any 1≤i<j≤n1\leq i<j\leq n and 0<δ<b(j)−b(i)20<\delta<\frac{b_{(j)}-b_{(i)}}{2},

Note that the transfer should not reverse the relative position of the two individuals ii and jj. The transfer principle is sometimes called the Pigou-Dalton principle (Pigou, 1912; Dalton, 1920).

See Figure 3 for an illustration of this property.

While not all inequality measures satisfy the decomposability property (e.g., the Gini Index does not), the property has been studied extensively in economics, as it allows economists to compare patterns and dynamics of inequality in different subpopulations (e.g., racial minorities (Conceição and Ferreira, 2000)).

Our measure of unfairness. For quantifying algorithmic unfairness, in this paper, we focus on a family of inequality indices called generalized entropy indices. For a constant α∉{0,1}\alpha\notin\{0,1\}, the generalized entropy of benefits b1,b2,⋯ ,bnb_{1},b_{2},\cdots,b_{n} with mean benefit μ\mu is defined as follows:

One can interpret generalized entropy as a measure of information theoretic redundancy in data. Generalized entropy satisfies the earlier properties of anonymity, population-invariance, the Pigou-Dalton transfer principle, and zero-normalization. Further it is subgroup decomposable (Cowell and Kuga, 1981), and also scale-invariant.A measure II is scale-invariant if for any constant c>0c>0, I(cb)=I(b)I(c\mathbf{b})=I(\mathbf{b}). In fact, Shorrocks (1980) show that generalized entropy is the only differentiable family of inequality indices that satisfies population- and scale-invariance. Our interest in this family of inequality indices is motivated by this result and by our aim of understanding the trade-offs between individual and group-level unfairness.

4. Comparison with Previous Work

Existing notions of algorithmic fairness can be divided into two distinct categories: group and individual fairness.

Group fairness. Group fairness notions require that given a classifier θ\theta, a certain group-conditional quality metric qz(θ)q_{z}(\theta) is the same for all sensitive feature groups. That is:

Different choices for qz(.)q_{z}(.) have led to different namings of the corresponding group fairness notions (see e.g., statistical parity (Kleinberg et al., 2017; Dwork et al., 2012; Corbett-Davies et al., 2017), disparate impact (Zafar et al., 2017b; Feldman et al., 2015), equality of opportunity (Hardt et al., 2016), calibration (Kleinberg et al., 2017), and disparate mistreatment (Zafar et al., 2017a)). Generally, these notions cannot guarantee fairness at the individual level, or when groups are further refined (see Kearns et al. (Kearns et al., 2017) for an illustrative example).

Existing group fairness notions are similar to the between-group component of fairness that we propose. However, these notions usually do not take into account the size of different groups, whereas our between-group measure considers the proportion of the groups relative to the total population as illustrated in Figure 3. For example consider a population divided into two groups AA and BB containing 70% and 30% of the population with the negative ground truth label respectively. Using generalized entropy with α=2\alpha=2, a classifier C1C_{1} achieving a false positive rate of 0.80.8 on AA and 0.60.6 on BB has a between-group inequality of 0.060.06, whereas a classifier C2C_{2} with false positive rates of 0.60.6 on AA and 0.80.8 on BB results in a lower between-group inequality of 0.040.04. However, when considering a group fairness measure based on differences in false positive rates between AA and BB, C1C_{1} and C2C_{2} would be equally fair.

Individual fairness. Dwork et al. (2012) first formalized the notion of individual fairness for classification tasks using Lipschitz conditions on the classifier outcomes. Their notion of individual fairness requires that two individuals who are similar with respect to the task at hand, receive similar classification outcomes. Dwork et al.’s definition is, therefore, formalized in terms of a similarity function between individuals. For instance, in practice given two individuals with feature values x\mathbf{x} and x′\mathbf{x}^{\prime}, and suitable distance functions DXD_{\mathcal{X}} and DYD_{\mathcal{Y}} (defined over X×X\mathcal{X}\times\mathcal{X} and Δ(Y)×Δ(Y)\Delta(\mathcal{Y})\times\Delta(\mathcal{Y}), respectively), Dwork et al.’s notion for individual fairness requires the following condition to hold:

Due to its dependence on the individual feature vectors x\mathbf{x}, Dwork et al.’s notion of individual fairness does not satisfy the anonymity principle.

Furthermore, Dwork et al.’s notion of individual fairness only provides a ‘yes/no’ answer to whether fairness conditions are satisfied, but does not provide a meaningful measure of algorithmic fairness when considered independent of prediction accuracy. We further illustrate this point with two examples: First, by this definition a model that assigns the same outcome to everyone is considered fair, regardless of people’s merit for different outcomes (e.g. awarding pretrial release to every defendant is considered fair, even though only some of them—those who appear for subsequent hearings and don’t commit a crimeFor making the pretrial release decisions, these two are the main criteria that the judges or the algorithms try to assess (Summers and Willis, 2010; Berk et al., 2017).—deserve to be awarded the pretrial release). Second, the definition does not take into account the difference in social desirability of various outcomes. For instance, if one flips the (binary) labels predicted by a fair classifier, the resulting classifier will be considered equally fair (e.g. a classifier that awards pretrial release to a defendant if and only if they go on to violate the release criteria is considered fair!). The measure we propose in equations 1 and 2 addresses these issues by offering a merit-based metric of fairness that seeks to equalize the benefit individuals receive as the result of being subject to algorithmic decision making. Finally, we remark that there has been interest in a similar axiomatic approach to methods for algorithmic interpretability (Sundararajan et al., 2017; Lundberg and Lee, 2017).

Theoretical Characterization

In this section, we characterize the conditions under which there is a tradeoff between accuracy and our notion of algorithmic fairness. Further, we shed light on the relationship between our notion of fairness and existing group measures, precisely connecting the two when the inequality index in use is additively decomposable (see (Subramanian, 2011) and the references therein). At a high level, we show that group unfairness is one piece of a larger puzzle: overall unfairness may be regarded as a combination of unfairness within- and between-groups. As the number of groupings increases, with each becoming smaller (eventually becoming single individuals), the between-group component grows to be an increasingly large part of the overall unfairness.

We begin by observing that the fairness optimal classifier is perfectly fair if and only if the accuracy optimal classifier is perfectly accurate. Given a classifier θ\theta and training data set D={(xi,yi)}i=1n\mathcal{D}=\{(\mathbf{x}_{i},y_{i})\}_{i=1}^{n}, let ID(θ)I_{\mathcal{D}}(\theta) specify the individual unfairness of θ\theta on D\mathcal{D}, that is, ID(θ)=I(b1θ,⋯ ,bnθ)I_{\mathcal{D}}(\theta)=I(b^{\theta}_{1},\cdots,b^{\theta}_{n}) where biθ=1+θ(xi)−yib^{\theta}_{i}=1+\theta(\mathbf{x}_{i})-y_{i}. LD(θ)L_{\mathcal{D}}(\theta) is the empirical loss of θ\theta on D\mathcal{D}.

Suppose I(.)I(.) is a zero-normalized inequality index and Θ\boldsymbol{\Theta} is closed under complements.In the context of a binary classification task with Y={0,1}\mathcal{Y}=\{0,1\}, Θ\boldsymbol{\Theta} is closed under complements if for any θ∈Θ\theta\in\boldsymbol{\Theta}, also 1−θ∈Θ1-\theta\in\boldsymbol{\Theta}. For any training data set D\mathcal{D}, there exists a classifier θ∈Θ\theta\in\boldsymbol{\Theta} for which ID(θ)=0I_{\mathcal{D}}(\theta)=0 if and only if there exists a classifier θ′\theta^{\prime} for which LD(θ′)=0L_{\mathcal{D}}(\theta^{\prime})=0. Proofs can be found in the appendix.

Proposition 3.1 may seem to suggest that our notion of fairness is entirely in harmony with prediction accuracy: by simply minimizing prediction error, unfairness will be automatically eliminated. While this is true in the special case of fully separable data (or when we have access to an oracle with 0 prediction error), it is not true in general. The following result shows that under broad conditions, the fairness optimal classifier may not coincide with the accuracy optimal classifier. For training data set D\mathcal{D}, let θDA\theta^{A}_{\mathcal{D}} be the accuracy optimal classifier, and θDF\theta^{F}_{\mathcal{D}} be the fairness optimal classifier:

Even though the fairness optimal and accuracy optimal classifiers do not necessarily coincide, one might wonder if the fairness optimal classifier always results in near-optimal accuracy. In fact, it does not. Example A.1 in the appendix shows that the accuracy of the fairness optimal classifier can be arbitrarily worse than that of the accuracy optimal classifier.

2. Individual vs. Group Fairness

Next, we focus on additive-decomposability and show how this property allows us to establish formally the existence of tradeoffs between individual- and group-level (un)fairness. Suppose we partition the population into ∣G∣|G| disjoint subgroups, where subgroup g∈Gg\in G consists of ngn_{g} individuals with the benefit vector bg=(b1g,⋯ ,bngg)\mathbf{b}^{g}=(b^{g}_{1},\cdots,b^{g}_{n_{g}}) and mean benefit μg\mu_{g}. Each partition could, for instance, correspond to a sensitive feature group (e.g., g=1g=1 consists of all African-American defendants and g=2g=2, all white defendants). One can re-write the Generalized Entropy as follows:

Note that imposing a constraint on a decomposable inequality measure, such as Eα(b)\mathcal{E}^{\alpha}(\mathbf{b}), guarantees both within-group and between-group inequality are bounded. Existing notions of group fairness, however, capture only the between-group component (when ∣G∣=2|G|=2, the between-group unfairness is minimized if and only if the two groups receive the same treatment on average). The problem with imposing a constraint on the between-group component (Eβα(b)\mathcal{E}^{\alpha}_{\beta}(\mathbf{b})) alone, is that it may drive up the within-group component, Eωα(b)\mathcal{E}^{\alpha}_{\omega}(\mathbf{b}). In fact, we show that if an individual-fairness optimal classifier is not group-fairness optimal, then optimizing for group fairness alone will certainly increase unfairness within groups sufficiently so as to raise the overall (individual) unfairness.

while minimizing only between-group unfairness corresponds to:

Let θ∗(δ)\theta^{*}(\delta) be an optimal solution for optimization (3)—if there are multiple optimal solutions, pick one with the lowest IβI_{\beta}. Let θβ∗(δ)\theta^{*}_{\beta}(\delta) be any optimal solution for optimization (4). The following holds.

Suppose I(.)I(.) is additively decomposable. For any δ∈\delta\in, if Iβ(θβ∗(δ))≠Iβ(θ∗(δ))I_{\beta}\left(\theta^{*}_{\beta}(\delta)\right)\neq I_{\beta}\left(\theta^{*}(\delta)\right), then Iω(θβ∗(δ))>Iω(θ∗(δ))I_{\omega}\left(\theta^{*}_{\beta}(\delta)\right)>I_{\omega}\left(\theta^{*}(\delta)\right) and I(θβ∗(δ))>I(θ∗(δ))I\left(\theta^{*}_{\beta}(\delta)\right)>I\left(\theta^{*}(\delta)\right).

Next, we show that the contribution of the between-group component to overall inequality, i.e. Iβ/II_{\beta}/{I}, depends—among other things—on the granularity of the groups. In particular, the between-group contribution increases with intersectionality: IβI_{\beta} is lower when computed over just race (African-Americans vs. Caucasians), and is higher when computed over the intersection of race and gender (female African-Americans, …, male Caucasians). More precisely, suppose G,G′G,G^{\prime} are two partitions of the population into disjoint groups. Let G×G′G\times G^{\prime} specify the Cartesian product of the two partitions: for g∈G,g′∈G′g\in G,g^{\prime}\in G^{\prime}, i∈(g,g′)i\in(g,g^{\prime}) if and only if i∈gi\in g and i∈g′i\in g^{\prime}. It is easy to show the following result.

Suppose I(.)I(.) is zero-normalized and additively decomposable. Suppose G,G′G,G^{\prime} are two different partitions of the population into disjoint groups. For any benefit distribution b\mathbf{b}, IβG(b)≤IβG×G′(b).I^{G}_{\beta}(\mathbf{b})\leq I^{G\times G^{\prime}}_{\beta}(\mathbf{b}).

If one continues refining the groups, eventually every individual will be in their own group and the between-group unfairness becomes equivalent to the overall individual unfairness. This offers a framework to interpolate between group and individual fairness. When the number of groups is small and people within each group receive highly unequal benefits, the contribution of the between-group component to overall unfairness is small, and narrowing down attention to reducing IβI_{\beta} alone may result in fairness gerrymandering (Kearns et al., 2017): while it is often easy to reduce group unfairness in this case, doing so will affect the overall unfairness in unpredictable ways—potentially making the within-group unfairness worse. On the other hand, as the number of groups increases or the treatment of people within a group becomes more uniform, the role that the between-group component plays in overall unfairness grows – but, as noted by Kearns et al. (Kearns et al., 2017), it also becomes computationally harder to control and limit the between-group unfairness.

Suppose I(.)I(.) is zero-normalized and additively decomposable. Suppose I(b)≠0I(\mathbf{b})\neq 0. For any partition GG of the population to disjoint subgroups, 0≤IβG(b)I(b)≤10\leq\frac{I^{G}_{\beta}(\mathbf{b})}{I(\mathbf{b})}\leq 1. Further, there exist benefit distributions b\mathbf{b} and b′\mathbf{b}^{\prime} such that IβG(b)I(b)=1\frac{I^{G}_{\beta}(\mathbf{b})}{I(\mathbf{b})}=1 and IβG(b′)I(b′)=0\frac{I^{G}_{\beta}(\mathbf{b}^{\prime})}{I(\mathbf{b}^{\prime})}=0.

We saw that the contribution of the group component to overall unfairness is a nuanced function of the granularity with which the groups are defined, as well as the unfairness within each group. If our goal is to reduce overall unfairness, note that existing fair learning models that exclusively focus on reducing between-group unfairness would help only when between-group accounts for a large part of the overall unfairness. Our measures of unfairness present a framework to examine this condition.

Empirical Analysis

In this section, we empirically validate our theoretical propositions from Section 3 on multiple real-world datasets. Specifically, in Section 4.1, our goal is to shed light on tradeoffs between overall individual-level unfairness and accuracy. In Section 4.2, we explore how the overall unfairness decomposes along the lines of sensitive attribute groups. We use the subgroup-decomposability of our proposed unfairness measures to study finer-grained fairness-accuracy tradeoffs at the levels of between-group and within-group unfairness. In Section 4.3, we empirically explore how methods to control between-group unfairness affect other unfairness components.

Setup and Datasets. We use the Generalized Entropy index (cf. Eq. 2) with α=2\alpha=2 (in other words, half the squared coefficient of variation) to measure unfairness:

As noted in Section 2, the Generalized Entropy index can be further decomposed into between-group and within-group unfairness as:

We will refer to the quantity E2\mathcal{E}^{2} as individual unfairness or overall unfairness interchangeably, Eβ2\mathcal{E}^{2}_{\beta} as between-group unfairness, and Eω2\mathcal{E}^{2}_{\omega} as within-group unfairness.

We experiment with two real-world datasets: (i) the Adult income dataset (Lichman, 2013), and (ii) the ProPublica COMPAS dataset (Larson et al., 2016). Both datasets have received previous attention (Zafar et al., 2017b, a; Corbett-Davies et al., 2017; Zemel et al., 2013; Feldman et al., 2015).

For the Adult income dataset, the task is to predict whether an individual earns more (positive class) or less (negative class) than 50,000 USD per year based on features like education level and occupation. We consider gender (female and male) and race (Black, White and Asian) as sensitive features. We filter out races (American Indian and Other) which constitute less than 1%1\% of the dataset. After the filtering, the dataset consists of 44,43444,434 subjects and 1111 features.

For the ProPublica COMPAS dataset, the task is to predict whether (negative class) or not (positive class) a criminal defendant would commit a crime within two years based on features like current charge degree or number of prior offenses. We use the same set of features as Zafar et al. (Zafar et al., 2017a). The sensitive features in this case are also gender (female and male) and race (Black, Hispanic, White). The dataset consists of 5,7865,786 subjects and 55 features.

For all experiments, we repeatedly split the data into 70%70\%-30%30\% train-test sets 10 times and report average statistics. All hyperparameters are validated using a further 70%70\%-30%30\% split of the train set into train and validation sets.

We begin by studying the tradeoff between the accuracy and the overall individual unfairness of a given classifier (E2(b)\mathcal{E}^{2}(\mathbf{b}) in Eq. 5). We use three standard classifier models: logistic regression, support vector machine with RBF kernel (SVM), and random forest classifier. Results in Section 4.1 and Section 4.2 are computed by optimizing these classifiers for accuracy.

Each of the above models computes the likelihood of belonging to the positive class for every instance. We denote this likelihood by pip_{i} for an individual ii. We compare the fairness and accuracy of these classifiers with that of an “oracle” that can perfectly predict the label for every instance (and assigns pi∈{0,1}p_{i}\in\{0,1\}). To predict a label in {0,1}\{0,1\} for individuals, we first rank all instances (in increasing order) according to their pip_{i} values with ties broken randomly; then we designate a decision ranking threshold 0≤τ≤10\leq\tau\leq 1 and output a label of 11 for individual ii if and only if rank(i)≥nτrank(i)\geq n\tau, where rank(i)rank(i) denotes the rank of individual ii in the sorted list. In other words, an increasing value of τ\tau corresponds to the classifier rejecting more people from the sorted list (in order of their positive class likelihood pip_{i}). As we vary τ\tau, we expect both accuracy as well as unfairness of the resulting predictions to change as discussed below and shown in Figure 5.

For the oracle, as expected from Proposition 3.1, a perfect accuracy corresponds to zero unfairness: with an increasing τ\tau, the accuracy increases while the unfairness decreases. After a certain optimal value of threshold τ\tau (close to 0.750.75 in the Adult data and 0.450.45 in the COMPAS data), the trend reverses. We note that 0.750.75 and 0.450.45 represent the fraction of instances in the negative class in the respective datasets. Hence, at these optimal thresholds, all of the oracle’s predictions are accurate (since the points are ranked based on their positive class likelihood) resulting in unfairness.

However, for all other (non-oracle) classifiers, as expected via Proposition 3.2, the trend is very different: the optimal threshold for (imperfect) accuracy is far from the optimal threshold for unfairness. Moreover, with increasing τ\tau, while unfairness continually increases, for accuracy we initially see an increase followed by a drop. We note that the overall unfairness is not always a monotone function of the decision ranking threshold as illustrated in Example A.1.

2. Fairness Decomposability

As Figure 3 and Eq. 5 show, the individual unfairness of a predictor can be decomposed into between-group and within-group unfairness. In this part, we study the between-group unfairness component (Eβ2(b)\mathcal{E}^{2}_{\beta}(\mathbf{b}) in Eq. 5) as we change τ\tau. To this end, we consider two sensitive features: gender and race. We split each of the datasets into all possible disjoint groups based on these sensitive features (e.g., White women, Hispanic men, Black women).

Figure 7 shows between-group unfairness along with overall unfairness for different values of τ\tau. We notice that for the Adult dataset, the between-group unfairness follows a multi-modal trend: it starts from a non-zero value at τ=0\tau=0, falls to almost for most classifiers (except for the oracle) at around τ=0.2\tau=0.2, reaches a local peak again around τ=0.5\tau=0.5, and finally completes another cycle to fall and then reach its maximum value at τ=1.0\tau=1.0. The COMPAS dataset also shows a similar trend, albeit to a lesser extent.

Between-group unfairness and overall unfairness. Comparing the between-group unfairness and overall unfairness in Figure 7 reveals a very interesting insight: for the same value of τ\tau, the between-group unfairness (solid lines) is a very small fraction of the overall unfairness (dotted lines). For example, considering the performance of the logistic regression classifier on the COMPAS dataset, the maximum value of the overall unfairness is close to 0.60.6 whereas the maximum value of the between-group unfairness is merely 0.010.01. We hypothesize that since the number of sensitive feature-based groups is much smaller than the number of all individuals in the dataset, the individual unfairness value dominates the between-group unfairness.

To test this hypothesis, we experiment with the following setup: We take the three sensitive features present in the COMPAS dataset, namely gender, race and age, and form sensitive feature groups based on all possible combinations of these features. For example, groups formed based on gender would be men and women; groups formed based on race would be Black, White and Hispanic; whereas groups formed based on gender as well as race would be Black men, Black women, White men and so on. For each of these sensitive feature combinations we plot in Figure 8 the percentage of contribution that the between-group unfairness has towards the overall unfairness. Figure 8 shows that as the number of sensitive feature groups increases, the between-group unfairness contributes more and more towards the overall unfairness. This result is also in line with the implications of Proposition 3.4.

Figure 8 also shows the following interesting insight: Even though the overall unfairness and accuracy of all the classifiers is very similar (cf. Figure 5), the random forest classifier leads of significantly smaller contribution of between-group unfairness as compared to other classifier (e.g., gender, gender+age). In other words even for similar levels of accuracy and overall unfairness, classifiers have very different between-group unfairness across different feature sets.

Accuracy and between-group unfairness. We also study the between group unfairness in Figure 7 and corresponding accuracy in Figure 5. We notice that there is no definitive correlation between the accuracy and the between-group unfairness: In the Adult dataset, the highest level of accuracy (around τ=0.8\tau=0.8) corresponds to one of the lowest values of between-group unfairness for all classifiers. However, this doesn’t hold in the case of COMPAS dataset.

3. Interaction Between Different Types of Unfairness

In this section, we revisit the literature on fairness-aware machine learning and investigate how methods proposed to control between-group unfairness (which is what most existing methods focus on (Zafar et al., 2017b, a; Feldman et al., 2015; Kamiran and Calders, 2009; Kamishima et al., 2013; Zemel et al., 2013; Hardt et al., 2016)) can affect the overall/individual and the within-group unfairness. Specifically, we study how the overall unfairness (E2(b)\mathcal{E}^{2}(\mathbf{b}) in Eq. 5) and the within-group unfairness (Eω2(b)\mathcal{E}^{2}_{\omega}(\mathbf{b}) in Eq. 5) would change when training a constrained classifier to minimize the between-group unfairness (Eβ2(b)\mathcal{E}^{2}_{\beta}(\mathbf{b}) in Eq. 5). Our study is motivated by the fact that while several methods focus on designing constraints to remove the between-group unfairness (e.g., see (Zafar et al., 2017a; Hardt et al., 2016; Zemel et al., 2013)), to the best of our knowledge, no prior work in fairness-aware machine learning has studied the effect of these constraints on the overall and the within-group unfairness.

To this end, we use the methodology proposed by Zafar et al. (Zafar et al., 2017a) to remove the between-group unfairness based on false negative rates between different races (Whites and non-Whites) in the COMPAS dataset. Zafar et al. propose to remove the between-group unfairness by bounding the covariance between misclassification distance from the decision boundary and the sensitive feature value. The method operates by bounding the covariance of the unconstrained classifier by successive multiplicative factors between 1 and 0. A covariance multiplicative factor of 11 means that no fairness constraints are applied while training the classifier, whereas a factor of means the tightest possible constraints are applied. As done by Zafar et al., we train several logistic regression classifiers to limit the between-group unfairness; each classifier is trained with a covariance multiplicative factor in the range [1.00,0.95,0.90,…,0.05,0.00][1.00,0.95,0.90,\ldots,0.05,0.00].

Figure 10 shows the between-group unfairness, within-group unfairness, and overall/individual unfairness as the fairness constraints of Zafar et al. (Zafar et al., 2017a) are tightened towards . The figure shows the following key insights: (i) Reducing the between-group unfairness can in fact increase the within-group unfairness: the within-group unfairness for Whites almost monotonically increases as the between-group unfairness is reduced. This observation also follows Proposition 3.3. (ii) Reducing the between-group unfairness can exacerbate overall/individual unfairness: As the between-group unfairness decreases between the covariance multiplicative factor of 0.80.8 to 0.60.6 (on the x-axis), the overall unfairness in fact goes up. These insights point to possible significant tensions between these different components of unfairness.

Summary of empirical analysis. Experiments on multiple real-world datasets performed in this section support the theoretical analysis of Section 3. The empirical (as well as the theoretical) analysis brings out the inherent tensions between fairness and accuracy, as well as between different (between- and within-group) components of fairness. These results point to potential for situations where optimizing for one type of fairness can exacerbate the other.

Conclusion

We proposed using inequality indices from economics as a principled way to compute the scalar degree of total unfairness of any algorithmic decision system. The approach is based on well-justified principles (axioms), and is general enough so that by varying the benefit function, we can capture all previous notions of algorithmic fairness conditions as special cases, while also admitting interesting generalizations. The resulting measures of total unfairness unify previous concepts of group and individual fairness, and allow us to study quantitatively the behavior of earlier methods to mitigate unfairness. These earlier methods typically worry only about between-group unfairness, which may be justified for legal reasons, or in order to redress particular social prejudices. However, we demonstrate that minimizing exclusively between-group unfairness may actually increase overall unfairness.

Acknowledgements

AW acknowledges support from the David MacKay Newton research fellowship at Darwin College, The Alan Turing Institute under EPSRC grant EP/N510129/1 & TU/B/000074, and the Leverhulme Trust via the CFI. HH acknowledges support from the Innosuisse grant 27248.1 PFES-ES.

References

Appendix A Appendix: Technical Material

If there exists a classifier θ′\theta^{\prime} with LD(θ′)=0L_{\mathcal{D}}(\theta^{\prime})=0, then that classifier minimizes II: For all i=1,⋯ ,ni=1,\cdots,n, the benefit ii receives under θ′\theta^{\prime}, denoted by biθ′b^{\theta^{\prime}}_{i}, is equal to 1+θ′(xi)−yi=11+\theta^{\prime}(\mathbf{x}_{i})-y_{i}=1. That is everyone gets the same benefit under θ′\theta^{\prime}, and as the result, I(bθ′)=I(1)=0I(\mathbf{b}^{\theta^{\prime}})=I(\mathbf{1})=0.

If there exists a classifier θ\theta with I(bθ)=0I(\mathbf{b}^{\theta})=0, θ\theta must assign the same benefit to everyone: there exists b∈{0,1,2}b\in\{0,1,2\} such that for all i=1,⋯ ,ni=1,\cdots,n, biθ=1+θ(xi)−yi=bb^{\theta}_{i}=1+\theta(\mathbf{x}_{i})-y_{i}=b. Now let θ′(x)=θ(x)+1−b\theta^{\prime}(\mathbf{x})=\theta(\mathbf{x})+1-b. It is easy to verify that θ′\theta^{\prime} has zero error (i.e. LD(θ′)=0L_{\mathcal{D}}(\theta^{\prime})=0). ∎

Proof of Proposition 3.2

Now consider a probabilistic classifier θq\theta^{q} that randomly assigns label 1 to each instance in D\mathcal{D} with probability q∈(0,1)q\in(0,1). The resulting benefit distribution bq\mathbf{b}^{q} of θq\theta^{q} is as follows: p(1−q)p(1-q) fraction of the population receive benefit ; pq+(1−p)(1−q)pq+(1-p)(1-q) faction receive benefit 11; and (1−p)q(1-p)q faction receive benefit 22.

We claim that for q=(1−p)q=(1-p), ID(θA)>ID(θq)I_{\mathcal{D}}(\theta^{A})>I_{\mathcal{D}}(\theta^{q}). To see this, note that bq\mathbf{b}^{q} can be constructed from bA\mathbf{b}^{A} via a series of inequality-reducing operations:

b′=2×bA\mathbf{b}^{\prime}=2\times\mathbf{b}^{A}, so that b′\mathbf{b}^{\prime} consists of pp fraction of the population receiving benefit and the other (1−p)(1-p) fraction receiving 22. Due to the scale invariance property of II, we know I(b′)=I(bA)I(\mathbf{b}^{\prime})=I(\mathbf{b}^{A}).

Perform the following progressive transfer on b′\mathbf{b}^{\prime} to obtain b′′\mathbf{b}^{\prime\prime}: Take one unit of benefit from p(1−p)p(1-p) fraction of the population whose benefit is 2, and give it to p(1−p)p(1-p) fraction with benefit 0. The resulting distribution, b′′\mathbf{b}^{\prime\prime}, consists of p2p^{2} faction with benefit ; 2p(1−p)2p(1-p) faction with benefit 11; and (1−p)2(1-p)^{2} faction with benefit 22. Because II satisfies the Dalton principle and p(1−p)>0p(1-p)>0, we have that I(b′′)<I(b′)I(\mathbf{b}^{\prime\prime})<I(\mathbf{b}^{\prime}).

Combining the above two, we have I(b′′)<I(bA)I(\mathbf{b}^{\prime\prime})<I(\mathbf{b}^{A}). It only remains to note that b′′\mathbf{b}^{\prime\prime} is precisely bq\mathbf{b}^{q} for q=(1−p)q=(1-p). Therefore, we conclude I(b′′)=I(b(1−p))<I(bA)I(\mathbf{b}^{\prime\prime})=I(\mathbf{b}^{(1-p)})<I(\mathbf{b}^{A}). This finishes the proof. ∎

The following example shows that the accuracy of the fairness optimal classifier can be arbitrarily bad compared to that of the accuracy optimal classifier.

Consider the example in the proof of Proposition 3.2. Let II be the generalized entropy with α=2\alpha=2. We claim that the fairness optimal classifier is one that assigns label 1 to every instance. To see this, recall that under θq\theta^{q}, p(1−q)p(1-q) fraction of the population receive benefit ; pq+(1−p)(1−q)pq+(1-p)(1-q) faction receive benefit 11; and (1−p)q(1-p)q faction receive benefit 22. So the mean benefit μ\mu is equal to 1−p+q1-p+q. Taking derivative with respect to qq, we have

The derivative is positive for q<q∗q<q^{*} and negative for q>q∗q>q^{*}. The minimum therefore happens at either q=0q=0 or q=1q=1. Given that for 0<p<10<p<1, 11−p>4−3p(2−p)2\frac{1}{1-p}>\frac{4-3p}{(2-p)^{2}}, we obtain that q=1q=1 minimizes II.

The fairness optimal classifier assigns label 1 to every instance resulting in accuracy pp, whereas the accuracy optimal classifier can achieve accuracy (1−p)(1-p). The ratio 1−pp\frac{1-p}{p} can be arbitrarily large if pp is taken to be sufficiently small.

Proof of Proposition 3.3

Note that because of the optimality of θβ∗\theta^{*}_{\beta} for (4), if Iβ(θβ∗)≠Iβ(θ∗)I_{\beta}(\theta^{*}_{\beta})\neq I_{\beta}(\theta^{*}), it must be the case that Iβ(θβ∗)<Iβ(θ∗)I_{\beta}(\theta^{*}_{\beta})<I_{\beta}(\theta^{*}). If I(θβ∗)≤I(θ∗)I(\theta^{*}_{\beta})\leq I(\theta^{*}), then θβ∗\theta^{*}_{\beta} is an optimal solution to (3), and Iβ(θβ∗)<Iβ(θ∗)I_{\beta}(\theta^{*}_{\beta})<I_{\beta}(\theta^{*}). This is a contradiction with the choice of θ∗\theta^{*}. If Iω(θβ∗)≤Iω(θ∗)I_{\omega}(\theta^{*}_{\beta})\leq I_{\omega}(\theta^{*}), then combined with the fact that Iβ(θβ∗)<Iβ(θ∗)I_{\beta}(\theta^{*}_{\beta})<I_{\beta}(\theta^{*}), we have that I(θβ∗)<I(θ∗)I(\theta^{*}_{\beta})<I(\theta^{*}), which is a contradiction with the optimality of θ∗\theta^{*} for (3). ∎

Proof of Proposition 3.4

Suppose ∣G∣=m|G|=m and ∣G′∣=m′.|G^{\prime}|=m^{\prime}. Let b=(b(g1,g1′),⋯ ,b(gm,gm′′))\mathbf{b}=\left(\mathbf{b}^{(g_{1},g^{\prime}_{1})},\cdots,\mathbf{b}^{(g_{m},g^{\prime}_{m^{\prime}})}\right) where b(gi,gj′)\mathbf{b}^{(g_{i},g^{\prime}_{j})} specifies the benefit distribution for individuals in group (gi,gj′)(g_{i},g^{\prime}_{j}) and μ(gi,gj′)\boldsymbol{\mu}^{(g_{i},g^{\prime}_{j})} specifies the distribution in which each individual in (gi,gj′)(g_{i},g^{\prime}_{j}) receives the group’s mean benefit. Note that IβG×G′(b)I^{G\times G^{\prime}}_{\beta}(\mathbf{b}) can be written as

where to obtain the second line, we used the additive decomposability property of II; for the third line we used the definition of the between-group component, and finally to obtain the conclusion, we used the zero-normalization property of II. ∎

Proof of Proposition 3.5

Recall that I(b)=IβG(b)+IωG(b)I(\mathbf{b})=I^{G}_{\beta}(\mathbf{b})+I^{G}_{\omega}(\mathbf{b}) and IβG(b),IωG(b)≥0I^{G}_{\beta}(\mathbf{b}),I^{G}_{\omega}(\mathbf{b})\geq 0. Therefore, we have 0≤IβG(b)I(b)≤10\leq\frac{I^{G}_{\beta}(\mathbf{b})}{I(\mathbf{b})}\leq 1.

Consider a benefit distribution b\mathbf{b} in which members of group g1∈Gg_{1}\in G receive benefit 1, and everyone else receives benefit . It is easy to see that for this distribution IβG(b)I(b)=1\frac{I^{G}_{\beta}(\mathbf{b})}{I(\mathbf{b})}=1. Similarly, consider a benefit distribution b′\mathbf{b}^{\prime} that assigns a benefit of 1 to half of the population in each group, and 0 to everyone else. It is easy to that IβG(b′)I(b′)=0\frac{I^{G}_{\beta}(\mathbf{b}^{\prime})}{I(\mathbf{b}^{\prime})}=0. ∎

The following example shows that an added feature may in fact worsen the unfairness of the accuracy optimal classifier.

In the resulting benefit distribution, r2−ϵ\frac{r}{2}-\epsilon fraction of the total population receives benefit , p+2ϵp+2\epsilon fraction of the total population receives benefit 11, and the other (1−p−r2−ϵ)(1-p-\frac{r}{2}-\epsilon) fraction receives benefit 22. The mean benefit is, therefore, (2−p−r)(2-p-r) and GE is equal to

Take p=0.9p=0.9, r=0.2r=0.2 and ϵ=0.001\epsilon=0.001, and we have:

The above shows the addition of a new feature can worsen the fairness of the accuracy optimal classifier.