Fair Algorithms for Learning in Allocation Problems

2018·Arxiv

Abstract

Abstract

Settings such as lending and policing can be modeled by a centralized agent allocating a scarce resource (e.g. loans or police officers) amongst several groups, in order to maximize some objective (e.g. loans given that are repaid, or criminals that are apprehended). Often in such problems fairness is also a concern. One natural notion of fairness, based on general principles of equality of opportunity, asks that conditional on an individual being a candidate for the resource in question, the probability of actually receiving it is approximately independent of the individual’s group. For example, in lending this would mean that equally creditworthy individuals in different racial groups have roughly equal chances of receiving a loan. In policing it would mean that two individuals committing the same crime in different districts would have roughly equal chances of being arrested.

In this paper, we formalize this general notion of fairness for allocation problems and investigate its algorithmic consequences. Our main technical results include an efficient learning algorithm that converges to an optimal fair allocation even when the allocator does not know the frequency of candidates (i.e. creditworthy individuals or criminals) in each group. This algorithm operates in a censored feedback model in which only the number of candidates who received the resource in a given allocation can be observed, rather than the true number of candidates in each group. This models the fact that we do not learn the creditworthiness of individuals we do not give loans to and do not learn about crimes committed if the police presence in a district is low.

As an application of our framework and algorithm, we consider the predictive policing problem, in which the resource being allocated to each group is the number of police officers assigned to each district. The learning algorithm is trained on arrest data gathered from its own deployments on previous days, resulting in a potential feedback loop that our algorithm provably overcomes. In this case, the fairness constraint asks that the probability that an individual who has committed a crime is arrested should be independent of the district in which they live. We empirically investigate the performance of our learning algorithm on the Philadelphia Crime Incidents dataset.

1 Introduction

The bulk of the literature on algorithmic fairness has focused on classification and regression problems (see e.g. [3, 4, 6–8, 10, 14, 16, 17, 19, 20, 25–27] for a collection of recent work), but fairness concerns also arise naturally in many resource allocation settings. Informally, a resource allocation problem is one in which there is a limited supply of some resource to be distributed across multiple groups with differing needs. Resource allocation problems arise in financial applications (e.g. allocating loans), disaster response (allocating aid), and many other domains — but the primary example that we will focus on in this paper is policing. In the predictive policing problem, the resource to be distributed is police officers, which can be dispatched to different districts. Each district has a different crime distribution, and the goal (absent additional fairness constraints) might be to maximize the number of crimes caught.

Of course, fairness concerns abound in this setting, and recent work (see e.g. [11, 12, 21]) has highlighted the extent to which algorithmic allocation might exacerbate those concerns. For example, Lum and Isaac [21] show that if predictive policing algorithms such as PredPol are trained using past arrest data to predict future crime, then pernicious feedback loops can arise, which misestimate the true crime rates in certain districts, leading to an overallocation of police.Since the communities that Lum and Isaac [21] showed to be overpoliced on a relative basis were primarily poor and minority, this is especially concerning from a fairness perspective. In this work, we study algorithms that both avoid this kind of under-exploration and can incorporate quantitative fairness constraints.

In the predictive policing setting, Ensign et al. [11] implicitly consider an allocation to be fair if police are allocated across districts in direct proportion to the district’s crime rate; generally extended, this definition asks that units of a resource are allocated according to the group’s share of the total candidates for that resource. In our work, we study a different notion of allocative fairness that has a similar motivation to the notion of equality of opportunity proposed by Hardt et al. [14] in classification settings. Informally speaking, it asks that the probability that a candidate for a resource be allocated a resource should be independent of his group. In the predictive policing setting, it asks that conditional on committing a crime, the probability that an individual is apprehended should not depend on the district in which they commit the crime.

To illustrate that our notions of fairness do not depend on whether individuals would prefer to receive or not receive the good, we highlight another setting in which allocative fairness is a natural concern: hiring.Suppose a company wishes to recruit machine learning programmers by advertising on a social media platform. Many such platforms offer the ability to advertise to different demographics of users and charge by the number of times the advertisement is shown to different users (i.e., the number of impressions); a fixed advertising budget can then be viewed as a number of impressions to allocate. Depending on how well the platform can identify programmers within each demographic, the ad may be shown to a higher or lower number of programmers. In this setting, our notion of allocative fairness asks that the probability a programmer is exposed to the hiring ad (and thus, receives the opportunity to apply for a job) does not depend on the programmer’s demographic, and the allocation problem is to maximize the number of programmers reached via the choice of impressions across each demographic, subject to fairness constraints.

1.1 Our Results

To define the extent to which an allocation satisfies our fairness constraint, we must model the specific mechanism by which resources deployed to a particular group reach their intended targets. We study two such discovery models, and we view the explicit framing of this modeling step as one of the contributions of our work; the implications of a fairness constraint depend strongly on the details of the discovery model, and choosing one is an important step in making one’s assumptions transparent.

We study two discovery models which capture two extremes of targeting ability. In the random discovery model, however many units of the resource are allocated to a given group, all individuals within that group are equally likely to be assigned a unit, regardless of whether they are a candidate for the resource or not. In other words, the probability that a candidate receives a resource is equal to the ratio of the number of units of the resource assigned to his group to the size of his group (independent of the number of candidates in the group).

At the other extreme, in the precision discovery model, units of the resource are given only to actual candidates within a group, as long as there is sufficient supply of the resource. In other words, the probability that a candidate receives a resource is equal to the ratio of the number of units of the resource assigned to his group, to the number of candidates within his group.

In the policing setting, these models can be viewed as two extremes of police targeting ability for an intervention like stop-and-frisk. In the random model, police are viewed as stopping people uniformly at random. In the precision model, police have the omniscient ability to identify individuals with contraband, and stop only them. Of course, reality lies somewhere in between.

These different discovery models have different implications for fairness. In the random model, fairness constrains resources to be distributed in amounts proportional to group sizes, regardless of the distribution of candidates, and so is uninteresting from a learning perspective. On the other hand, the precision model yields an interesting fairness-constrained learning problem when the distribution of the number of candidates in each group must be learned via observation, and what counts as a ‘fair’ allocation depends greatly on these distributions.

We study learning in a censored feedback setting: each round, the algorithm can choose a feasible deployment of resources across groups. Then the number of candidates for the current round in each group is drawn from a fixed, but unknown group-dependent distribution (which might be not be independent from the distributions in other groups). The algorithm does not observe the number of candidates present in each group, but only the number of candidates that received the resource. In the policing setting, this corresponds to the algorithm being able to observe the number of arrests, but not the actual number of crimes in each of the districts. Thus, the extent to which the algorithm can learn about the distribution in a particular group is limited by the number of resources it deploys there. The goal of the algorithm is to converge to an optimal fairness-constrained allocation, where here both the objective value of the solution, and the constraints imposed on it, depend on the unknown distributions.

One trivial solution to the learning problem is to sequentially deploy all of one’s resources to each group in turn for a sufficient amount of time to accurately learn the candidate distributions. This would reduce the learning problem to an offline constrained optimization problem, which we show can be efficiently solved by a greedy algorithm. But this algorithm is unreasonable: it has a large exploration phase in which it uses nonsensical deployments, vastly overallocating to some groups and underallocating to others. A much more realistic, natural approach is a greedy-style learning algorithm, which at each round simply uses its current best-guess estimate for the distribution in each group and deploys an optimal fairness-constrained allocation according to these estimates. Unfortunately, as we show, if one makes no assumptions on the underlying distributions, any algorithm that has a guarantee of converging to a fair allocation must behave like the trivial algorithm, deploying vast numbers of resources to each group in turn.

This impossibility result motivates us to consider the learning problem in which the unknown distributions are from a known parametric family. The natural greedy algorithm uses an optimal fair deployment at each round given the maximum likelihood estimates of candidate distributions given its (censored) observations so far; for concreteness, we analyze this algorithm in case of the Poisson distribution, and show that it converges to an optimal fair allocation, but our analysis generalizes for any single-parameter Lipschitz-continuous family of distributions.

Finally, we conduct an empirical evaluation of our algorithm on the Philadelphia Crime Incidents dataset, which records all crimes reported to the Philadelphia Police Department’s INCT system between 2006 and 2016. We verify that the crime distributions in each district are in fact well-approximated by Poisson distributions, and that our algorithm converges quickly to an optimal fair allocation (as measured according to the empirical crime distributions in the dataset). We also systematically evaluate the Price of Fairness, and plot the Pareto curves that trade off the number of crimes caught versus the slack allowed in our fairness constraint, for different sizes of police force, on this dataset. For the random discovery model, we prove worst-case bounds on the Price of Fairness.

1.2 Further Related Work

Our precision discovery model is inspired by and has technical connections to Ganchev et al. [13], which models the dark pool problem from quantitative finance, in which a trader wishes to execute a specified number of trades across a set of exchanges of unknown but independently distributed liquidity. In Ganchev et al. [13], the authors design an optimal allocation algorithm under the censored feedback of the precision model. It is straightforward to map their setting onto ours, but they assume independence between different exchanges, while the candidate distributions in our setting need not be independent. Regardless, we show that their allocation algorithm can be used to compute an optimal allocation (ignoring fairness) even when the independence assumption is relaxed (see Remark 1). Later, Agarwal et al. [1] extend the dark pool problem to an adversarial (rather than distributional) setting. This is quite closely related to the work of Ensign et al. [12] who also consider the precision model (under a different name) in an adversarial predictive policing setting. They provide no-regret algorithms for this setting by reducing the problem to learning in a partial monitoring environment. Since their setting is equivalent to that of Agarwal et al. [1], the algorithms in Agarwal et al. [1] can be directly applied to the problem studied by Ensign et al. [12].

Our desire to study the natural greedy algorithm rather than an algorithm which uses “unreasonable” allocations during an exploration phase is an instance of a general concern about exploration in fairnessrelated problems [5]. Recent works have studied the performance of greedy algorithms in different settings for this reason [2, 18, 24].

Lastly, the term fair allocation appears in the fair division literature (see e.g. [23] for a survey), but that body of work is technically quite distinct from the problem we study here.

2 Setting

We study an allocator who has V units of a resource and is tasked with distributing them across a population partitioned into G groups. Each group is divided into candidates, who are the individuals the allocator would like to receive the resource, and non-candidates, who are the remaining individuals. We let denote the total number of individuals in group i. The number of candidates in group i is a random variable drawn from a fixed but unknown distribution called the (marginal) candidate distribution. We do not make any assumptions about the relationship between the candidate distributions across different groups and in particular these distributions need not be independent. We use M to denote the total size of all groups (i.e., ). An allocation ) is a partitioning of these V units, where denotes the units of resources allocated to group i. Every allocation is bound by a feasibility constraint which requires that Σ.

A discovery model is a (possibly randomized) function disc() mapping the number of units allocated to group i and the number of candidates in group i to the number of candidates discovered in group i. In the learning setting, upon fixing an allocation v, the learner will get to observe (a realization of) disc() for the realized value of for each group i. Fixing an allocation v, a discovery model disc() and candidate distributions for all groups , we define the total expected number of discovered candidates, disc(as

where the expectation is taken over and any randomization in the discovery model disc(). When the discovery model and the candidate distributions are fixed, we will simply write ) for brevity. We also use the total expected number of discovered candidates and (expected) utility exchangeably. We refer to an allocation that maximizes the expected number of discovered candidates over all feasible allocations as an optimal allocation and denote it by .

2.1 Allocative Fairness

For the purposes of this paper, we say that an allocation is fair if it satisfies approximate equality of candidate discovery probability across groups. We call this discovery probability for brevity. This formalizes the intuition that it is unfair if candidates in one group have an inherently higher probability of receiving the resource than candidates in another. Formally, we define our notion of allocative fairness as follows.

Definition 1. Fix a discovery model disc() and the candidate distributions C. For an allocation v, let

for all pairs of groups i and j.

When it is clear from the context, for brevity, we write ) for the discovery probability in group i. We emphasize that this definition (1) depends crucially on the chosen discovery model, and (2) requires nothing about the treatment of non-candidates. We think of this as a minimal definition of fairness, in that one might want to further constrain the treatment of non-candidates — but we do not consider that extension.

Since discovery probabilities ) and ) are in [0, 1], the absolute value of their difference is in [0, 1]. By setting require exact fairness.

We refer to an allocation v that maximizes ) subject to -fairness and the feasibility constraint as an and denote it by . In general, ) is a non-increasing quantity in , since as diminishes, the utility maximization problem becomes more constrained.

Remark 1. We note that both the utility and discovery probabilities can be written solely in terms of the marginal candidate distributions in each of the groups, even when these distributions are not independent. This is because we have (implicitly) assumed that the number of candidates discovered in a group depends only on the number of candidates in the group and the allocation to that group, regardless of the allocations to and the number of candidates in other groups. This assumption together with the linearity of expectation allows us to write the expected utility as in the right hand side of Equation 1.

3 The Precision Discovery Model

We begin by describing the precision model of discovery. Allocating units to group i in the precision model results in the discovery of disc(min() candidates. This models the ability to perfectly discover and reach candidates in a group with resources deployed to that group, limited only by the number of deployed resources and the number of candidates present.

The precision model results in censored observations that have a particularly intuitive form. Recall that in general, a learning algorithm at each round gets to choose an allocation v and then observes disc() for each group i. In the precision model, this results in the following kind of observation: when is larger than , the allocator learns the number of candidates present on that day exactly. We refer to this kind of feedback as an uncensored observation. But when is smaller than , all the allocator learns is that the number of candidates is . We refer to this kind of feedback as a censored observation.

The rest of this section is organized as follows. In Sections 3.1 and 3.2 we characterize optimal and optimal fair allocations for the precision model when the candidate distributions are known. In Section 3.3 we focus on learning an optimal fair allocation when these distributions are unknown. We show that any learning algorithm that is guaranteed to find a fair allocation in the worst case over candidate distributions must have the undesirable property that at some point, it must allocate a vast number of its resources to each group individually. To bypass this hurdle, in Section 3.4 we show that when the candidate distributions have a parametric form, a natural greedy algorithm which always uses an optimal fair allocation for the current maximum likelihood estimates of the candidate distributions converges to an optimal fair allocation.

3.1 Optimal Allocation

We first describe how an optimal allocation (absent fairness constraints) can be computed efficiently when the candidate distributions are known. In Ganchev et al. [13], the authors provide an algorithm for computing an optimal allocation when the distributions over the number of shares present in each dark pool are known and the trader wishes to maximize the expected number of traded shares. While they assume that the distributions of shares across different dark pools are independent, our formulation does not require this assumption of independence. Regardless, we can use the same algorithm as in Ganchev et al. [13] to compute an optimal allocation in our setting; this is because, as stated in Remark 1, the utility in both settings can be written solely in terms of the (marginal) candidate distributions even when the candidate distributions are not independent across different groups. Here, we present the high level ideas of their algorithm in the language of our model, and provide full details for completeness in Appendix B.

Let ] denote the probability that there are at least c candidates in group i. We refer to ) as the . Recall that the value of the cumulative distribution function (CDF) of at c is defined to be

So ) can be written in terms of CDF values as 1). First, observe that the expected total number of candidates discovered by an allocation in the precision model can be written in terms of the tail probabilities of the candidate distributions i.e.

Since the objective function is concave (as ) is a non-increasing function in c for all i), a greedy algorithm which iteratively allocates the next unit of the resource to a group in

3.2 Optimal Fair Allocation

We next show how to compute an optimal -fair allocation in the precision model when the candidate distributions are known and do not need to be learned.

To build intuition for how the algorithm works, imagine that the group i has the highest discovery probability in , and the allocation to that group is somehow known to the algorithm ahead of time. The constraint of -fairness then implies that the discovery probability for each other group j in must satisfy ]. This in turn implies upper and lower bounds on the feasible allocations to group j. The algorithm is then simply a constrained greedy algorithm: subject to these implied constraints, it iteratively allocates units so as to maximize their marginal probability of reaching another candidate. Since the group i maximizing the discovery probability in and the corresponding allocation are not known ahead of time, the algorithm simply iterates through each possible choice of i.

Pseudocode is given in Algorithm 1. We prove that Algorithm 1 returns an optimal -fair allocation in Theorem 1. We defer the proof of Theorem 1 and all the other omitted proofs in the section to Appendix B.

Theorem 1. Algorithm 1 computes an optimal -fair allocation for the precision model in time O(GV(GV + M)).

3.3 Learning Fair Allocations Generally Requires Brute-Force Exploration

In Sections 3.1 and 3.2 we assumed the candidate distributions were known. When the candidate distributions are unknown, learning algorithms intending to converge to optimal -fair allocations must learn a sufficient amount about the distributions in question to certify the fairness of the allocation they finally output. Because learners must deal with feedback in the censored observation model, this places constraints on how

they can proceed. Unfortunately, as we show in this section, if candidate distributions are allowed to be worst-case, this will force a learner to engage in what we call “brute-force exploration” — the iterative deployment of a large fraction of the resources to each subgroup in turn. This is formalized in Theorem 2.

Theorem 2. Define to be the size of the largest group and assume 6 for all i and 2. Let [0)), (0, 1/2), and A be any learning algorithm for the precision model which runs for a finite number of rounds and outputs an allocation. Suppose that there is some group i for which A has not allocated at least 2 units for at least k ln(1) rounds upon termination, where k is an absolute constant. Then there exists a candidate distribution such that, with probability at least outputs an allocation that is not -fair.

Sketch of the Proof. Let i denote a group in which A has not allocated at least 2 units for at least k ln(1) rounds upon its termination and let v denote an arbitrary allocation. We will design two candidate distributions for group i which have true discovery probabilities that are at least 2apart given , but which are indistinguishable given the observations of the algorithm with probability at least . If the A cannot distinguish between and , it cannot distinguish between and , and thus cannot guarantee whether group i’s discovery probability is indeed within of every other group’s discovery probability.

To design these candidate distributions, consider distributions and which satisfy the following four conditions.

1. and agree on all values less than 2.

2. The total mass of both distributions below 2 is 1 .

3. The remaining 2mass of is on the value 1.

4. The remaining 2mass of is on the value .

Distinguishing between and requires at least one uncensored observation beyond 2. However, conditioned on allocating at least 2 units, the probability of observing an uncensored observation is at most 2. So to distinguish between and with confidence 1 , and therefore to guarantee an -fair allocation, a learning algorithm must allocate at least 2 units to group i for k ln(1) rounds.

Recall that we used to denote the size of the largest group. When , then Theorem 2 implies that no algorithm can guarantee -fairness for sufficiently small . Moreover, even when 2V, Theorem 2 shows that in general, if we want algorithms that have provable guarantees for arbitrary candidate distributions, it is impossible to avoid something akin to brute-force search (recall that there is a trivial algorithm which simply allocates all resources to each group in turn, for sufficiently many rounds to approximately learn the CDF of the candidate distribution, and then solves the offline problem). In the next section, we circumvent this by giving an algorithm with provable guarantees, assuming that the candidate distributions have a known parametric form.

3.4 Poisson Distributions and Convergence of the MLE

In this section, we assume that all the candidate distributions have a particular and known parametric form but that the parameters of the these distributions are not known to the allocator. Concretely, we assume that the candidate distribution for each group is Poisson(denoted by )) and write ) for the true underlying parameters of the candidate distributions; this choice appears justified, at least in the predictive policing application, as the candidate distributions in the Philadelphia Crime Incidents dataset are well-approximated by Poisson distributions (see Section 4 for further discussion). This assumption allows an algorithm to learn the tails of these distributions without needing to rely on brute-force search, thus circumventing the limitation given in Theorem 2. Indeed, we show that (a small variant of) the natural greedy algorithm incorporating these distributional assumptions converges to an optimal fair allocation.

For simplicity, we assume a parametric form on the marginal candidate distribution in each of the groups. We could have equivalently assumed that the candidates across groups are drawn from a multivariate Poisson distribution to highlight the (potential) correlation between candidates distributions. However, since for a given multivariate Poisson distribution the marginal distribution on each group is itself a Poisson distribution [15], we made our parametric assumption directly on these marginal distributions.

At a high level, in each round, our algorithm uses Algorithm 1 to calculate an optimal fair allocation with respect to the current maximum likelihood estimates of the group distributions; then, it uses the new observations it obtains from this allocation to refine these estimates for the next round. This is summarized in Algorithm 2. The algorithm differs from this pure greedy strategy in one respect, to overcome the following subtlety: there is a possibility that Algorithm 1, when operating on a preliminary estimate for the candidate distributions, will suggest sending zero units to some group, even when the optimal allocation for the true distributions sends some units to every group. Such a deployment would result in the algorithm receiving no feedback for the zero-allocated group that round. If this suggestion is followed and a lack of feedback is allowed to persist indefinitely, the algorithm’s parameter estimate for the zero-allocated group will also stop updating — potentially at an incorrect value. In order to avoid this problem and continue making progress in learning, our algorithm chooses another allocation in this case. As we show, any allocation that allocates positive resources to all groups will suffice; in particular, our algorithm makes the natural choice of simply repeating the allocation from the previous round.

Notice that Algorithm 2 chooses an allocation at every round which is fair with respect to its estimates of the parameters of the candidate distributions; hence, asymptotic convergence of its output to an optimal allocation follows directly from the convergence of the estimates to true parameters. However, we seek a stronger, finite sample guarantee, as stated in Theorem 3.

Theorem 3. Let 0. Suppose that the candidate distributions are Poisson distributions with unknown parameters in the vector , where lies in the known interval [. Suppose we run Algorithm 2 for t > ˜O(ln(rounds, where ) is some distribution specific functionto get an allocation ˆv and estimated parameters ˆfor all groups i. Then with probability at least 1

1. For all i in [G], .

2. Let ))) where denotes the total variation distance between two dis- tributions. Then ˆv

• has utility at most 4DGV smaller than the utility of an optimal ()-fair allocation i.e. .

Remark 2. Theorem 3 implies that in the limit, the allocation from Algorithm 2 will converge to an optimal -fair allocation. As , ˆfor all i, meaning 0 and more importantly, ˆv will be -fair and optimal.

The rest of this section is dedicated to the proof of Theorem 3. First, we introduce notation. Since we assumed the candidate distribution for each group is Poisson, the probability mass function (PMF) and the CDF of the candidate distribution for group i can be written as

Given an allocation of units of the resource to group i we use to denote the (possibly censored) observation received by Algorithm 2. So while the candidates in group i are generated according to ), the observations of Algorithm 2 follow a censored Poisson distribution which we abbreviate by ). We can write the PMF of this distribution as

where 1; ) is the CDF value of ) at 1. Since Algorithm 2 operates in rounds, we use the superscript t throughout to denote the round. For each round t, denote the history of the units allocated to group i and observations received (candidates discovered)

in rounds up to t by ). We use ) to denote the history for all groups. All the probabilities and expectations in this section are over the randomness of the observations drawn from the censored Poisson distributions unless otherwise noted; we suppress related notation for brevity. Finally, an allocation function A in round t is a mapping from the history of all groups to the number of units to be allocated to each group i.e. . For convenience, we use to denote the allocation at round t. We are now ready to define the likelihood functions.

Definition 2. Let ] denote the (censored) likelihood of discovering candidates given an allocation to group i assuming the candidate distribution follows ). We write ) as log ). So, given any history , the empirical log-likelihood function for group i is

The expected log-likelihood function given the history of allocations but over the randomness of the candidacy distribution can be written as

where the expectation is over the randomness of drawn from ).

Proof of Theorem 3

To prove Theorem 3, we first show that any sequence of allocations selected by Algorithm 2 will eventually recover the true parameters. There are two conceptual difficulties here: the first is that standard convergence results typically leverage the assumption of independence, which does not hold in this case as Algorithm 2 computes adaptive allocations which depend on the allocations in previous rounds; the second is the censoring of the observations. Despite these difficulties, we give quantifiable rates with which the estimates converge to the true parameters. Next, we show that computing an optimal -fair allocation using the estimated parameters will result in an allocation that is (+ 4D)-fair with respect to the true candidate distributions where D denotes the maximum total variation distance between the true and estimated Poisson distributions across all groups. Finally, we show that this allocation also achieves a utility that is comparable to the utility of an optimal ()-fair allocation. We note that while Theorem 3 is only stated for Poisson distributions, our results can be generalized to any single parameter Lipschitz-continuous family of distributions (see Remark 3).

Closeness of the Estimated Parameters Our argument can be stated at a high level as follows: for any group i and any history , the empirical log-likelihood converges to the expected log-likelihood for any sequence of allocations made by Algorithm 2 as formalized in Lemma 1. We then show in Lemma 2 that the closeness of the empirical and expected log-likelihoods implies that the maximizers of these quantities (corresponding to the estimated and true parameters) will also become close. Since in our analysis we consider the groups separately, we fix a group i throughout the rest of this section and drop the subscript i for convenience.

We start by studying the rate of convergence of the empirical log-likelihood to the expected log-likelihood.

Lemma 1. With probability at least 1 , for any t and any observed by Algorithm 2

The true and estimated parameters for each group correspond to the maximizers of the expected and empirical log-likelihoods, respectively (see Corollary 1 in Appendix B). We next show that closeness of the empirical and expected log-likelihoods implies that the true and estimated parameters are also close.

Lemma 2. Let ˆdenote the estimate of the Algorithm 2 after rounds. Then with probability at least 1 .

Proof. Since Corollary 1(in Appendix B) gives that ) has a unique maximizer at and Corollary 3(in Appendix B) gives that there exists some ) so that for any such that ), we must have that . We denote ) by for brevity. We define the empirical maximizer ˆto be the maximizer of ˆ) i.e.

Applying Lemma 1 implies that for any with and ], with probability at least 1 ,

In particular, we must have that

Since ˆis a maximizer of Equation 2, we have that

where the last inequality is by Equation 3. This implies that

where the last inequality is by Equation 3. So ˆ, and thus Corollary 3 gives that .

Combining Lemma 2 with a union bound over all groups show that, with probability 1 , if Algorithm 2 is run for rounds, then , for all i. Note that as , the maximum total variation distance D between the estimated and the true distribution will converge in probability to 0.

Fairness of the Allocation In this section, we show that the fairness violation (i.e. the maximum difference in discovery probabilities over all pairs of groups) is linear in terms of D. Therefore, as the running time of the Algorithm 2 increases and hence, 0, the fairness violation of ˆv approaches . This is stated formally as follows.

Lemma 3. Let ˆv denote the allocation returned by Algorithm 2 after rounds. Then with probability at least 1 + 4

Proof. For any ] we have that

The first inequality follows from the triangle inequality. In the second inequality, the second term can be bounded because Algorithm 1 returns an -fair allocation with respect to its input distribution. The first and third term in the second inequality can be bounded by Lemma 12 (in Appendix B). Lemma 12 shows that for any fixed allocation the difference between the discovery probability with respect to the true and estimated candidate distributions in group i is proportional to the total variation distance between the true and estimated distributions.

Utility of the Allocation In this section we analyze the utility of the allocation returned by Algorithm 2. Once again, note that as 0, which happens as the running time of Algorithm 2 increases, ˆv will become optimal and -fair.

Lemma 4. Let ˆv denote the allocation returned by Algorithm 2 after rounds. Then with probability at least 1

Proof. Consider the following optimization problem, ).

We can think of the above optimization problem as the case where the underlying candidate distributions used for the objective value and the fairness constraints are different. Let us write ) to denote an optimal allocation in the above optimization problem, ). So an optimal fair allocation and the allocation returned by Algorithm 2 can be written as ) and ), respectively.

Note that for any fixed allocation v

where ˆis the tail probability of C(ˆ). This is because ()). In other words, even when the underlying candidate distribution changes for the objective value, an allocation value can change by at most 2GVD.

Now observe that

The inequalities in the first and third lines are by Equation 4, which shows how the utility deteriorates when the underlying distribution for the objective function changes. The inequality in the second line follows from Lemma 3, as any () fair allocation is a feasible allocation to P(V, {C(ˆ(ˆ), and ) is an optimal solution to this problem.

Remark 3. Although we assumed Poisson distributions in this section, all our results hold for any single-parameter Lipschitz-continuous distribution whose parameter is drawn from a compact set. However, the convergence rate of Theorem 3 depends on the quantity ) which depends on the family of distributions used to model the candidate distributions.

4 Experiments

In this section, we apply our allocation and learning algorithms for the precision model to the Philadelphia Crime Incidents dataset, and complement the theoretical convergence guarantee of Algorithm 2 to an optimal fair allocation with empirical evidence suggesting fast convergence in practice. We also study the empirical trade-off between fairness and utility in the dataset.

4.1 Experimental Design

The Philadelphia Crime Incidents datasetcontains all the crimes reported to the Police Department’s INCT system between 2006 and 2016. The crimes are divided into two types. Type I crimes include violent offenses such as aggravated assault, rape, and arson among others. Type II crimes include simple assault, prostitution, gambling and fraud. For simplicity, we aggregate all crime of both types, but in practice, an actual police department would of course treat different categories of crime differently. We note as a caveat that these incidents are reported and may not represent the entirety of committed crimes.

Figure 1: Frequencies of the number of reported crimes in each district in the Philadelphia Crime Incidents dataset. The red curves display the best Poisson fit to the data.

To create daily crime frequencies in Figure 1, we first calculate the daily counts of criminal incidents in each of the 21 geographical police districts in Philadelphia by grouping together all the crime reports with the same date; we then normalize these counts to get frequencies.Each subfigure in Figure 1 represents a police district. The horizontal axis of the subfigure corresponds to the number of reported incidents in a day and the vertical axis represents the frequency of each number on the horizontal axis. These frequencies approximate the true (marginal) distributions of the number of reported crimes in each of the districts in Philadelphia. Therefore, throughout this section we take these frequencies as the ground truth candidate distributions for the number of reported incidents in each of the districts.

Figure 1 shows that crime distributions in different districts can be quite different; e.g., the average number of daily reported incidents in District 15 is 43.5, which is much higher than the average of 11.35 in District 1 (see Table 1 in Appendix C for more details). Despite these differences, each of the crime distributions can be approximated well by a Poisson distribution. The red curves overlayed in each subfigure correspond to the Poisson distribution obtained via maximum likelihood estimation on data from that district. Throughout, we refer to such distributions as the best Poisson fit to the data (see Table 2 in Appendix C for details about the goodness of fit).

In our experiments, we take the police officers assigned to the districts as the resource to be distributed, the ground truth crime frequencies as candidate distributions, and aim to maximize the sum of the number of crimes discovered under the precision model of discovery.

4.2 Results

We can quantify the extent to which fairness degrades utility in the dataset through a notion we call Price of Fairness (PoF henceforth). In particular, given the ground truth crime distributions and the precision model of discovery, for a fairness level , we define PoF(). The PoF is simply the ratio of the expected number of crimes discovered by an optimal allocation to the expected number of crimes discovered by an optimal -fair allocation. Since ) for all , the PoF is at least one. Furthermore, the PoF is monotonically non-increasing in . We can apply the algorithms given in Sections 3.1 and 3.2 respectively for computing optimal unconstrained, and optimal fair allocations with the with ground truth distributions as input and numerically compute the PoF. This is illustrated in Figure 2. The x axis corresponds to different values and the y axis displays 1/PoF(). Each curve corresponds to a different number of total police officers denoted by V. Because feasible allocations must be integral, there can sometimes be no feasible -fair allocation for small . Since the PoF in these cases is infinite we instead opt to display the inverse, 1/PoF, which is always bounded in [0, 1]. Higher values of inverse PoF are more desirable.

Figure 2: Inverse PoF plots for the Philadelphia Crime Incidents dataset. Smaller values indicate greater sacrifice in utility to meet the fairness constraint.

Figure 2 shows a diverse set of utility/fairness trade-offs depending on the number of available police officers. It also illustrates that the cost of fairness is rather low in most regimes. For example, in the worst case, with only 50 police officers (the black curve) (which is much smaller than the average number of daily reported crimes: 563.88), the inverse PoF is 1 for 1, which corresponds to a 10% difference in the discovery probability across districts. When we increase the number of available police officers to 400 (the magenta curve), tolerating only a 4% difference in the discovery probability across districts is sufficient to guarantee no loss in the utility. Figure 2 also shows that for any fixed , the inverse PoF() tends to increase as the number of police increases (i.e. the cost of fairness decreases).This captures the intuition that fairness becomes a less costly constraint when resources are in greater supply. Finally, we observe a thresholding phenomenon in Figure 2; in each curve, increasing beyond a threshold will significantly increase the inverse PoF. This is due to discretization effects, since only integral allocations are feasible.

We next turn into analyzing the performance of Algorithm 2 in practice. We run the algorithm instantiated to fit Poisson distributions, but use observations from the ground truth distribution at each round. As we have shown in Figure 1 and Table 2, the ground truth is well approximated by a Poisson distribution.

We measure the performance of Algorithm 2 as follows. First, we fix a police budget V and unfairness budget and run Algorithm 2 for 2000 rounds using the dataset as the ground truth. That is, we simulate each round’s crime count realizations in each of the districts as being sampled from the ground truth distributions, and return censored observations under the precision model to Algorithm 2 according to the algorithm’s allocations and the drawn realizations. The algorithm returns an allocation after termination and we can measure the expected number of crimes discovered and fairness violation (the maximum differ-ence in discovery probabilities over all pairs of districts) of the returned allocation using the ground truth distributions. Varying while fixing V allows us to trace out the Pareto frontier of the utility/fairness trade-off for a fixed police budget. Similarly, for any fixed V and , we can run Algorithm 1 (the offline algorithm for computing an optimal fair allocation) with the ground truth distributions as input and trace out a Pareto curve by varying . We refer to these two Pareto curves by the learned and optimal Pareto curves, respectively.So to measure the performance of Algorithm 2, we can compare the learned and optimal Pareto curves.

Figure 3: Pareto frontier of expected crimes discovered versus fairness violation.

In Figure 3, each curve corresponds to a police budget. The x and y axes represent the expected number of crimes discovered and fairness violation for allocations on the Pareto frontier, respectively. In our simulations we varied between 0 and 0.15. For each police budget V, the ‘x’ s connected by the dashed lines show the learning Pareto frontier. Similarly, the circles connected by solid lines show the optimal Pareto frontier. We point out that while it is possible for the fairness violations in the learned Pareto curves to be higher than the level of set as an input to Algorithm 2, the fairness violations in the optimal Pareto curves are always bounded by .

The disparity between the optimal and learned Pareto curves are due to the fact that the learning algorithm has not yet fully converged. This can be attributed to the large number of censored observations received by Algorithm 2, which are significantly less informative than uncensored observations. Censoring happens frequently because the number of police used in every case plotted is less than the daily average of 563.88 crimes across all the districts in the dataset — so it is unavoidable that in any allocation, there will be significant censoring in at least some districts.

Figure 3 shows that while the learning curves are dominated by the optimal curves, the performance of the learning algorithm approaches the performance of the offline optimal allocation as V increases. Again, this is because increasing V generally has the effect of decreasing the frequency of censoring.

We study the V = 500 regime in more detail, to explore the empirical rate of convergence. In Figure 4, we study the round by round performance of the allocation computed by Algorithm 2 in a single run with the choice of 05.

Figure 4: The per round expected number of crimes discovered and fairness violation of Algorithm 2. V = 500 and 05.

In Figure 4, the x axis labels progression of rounds of the algorithm. The y axis measures the fairness violation (left) and expected number of crimes discovered (right) of the allocation deployed by the algorithm, as measured with respect to the ground truth distributions. The black curves represent Algorithm 2. For comparison we also show the same quantities for an offline optimal fair allocation as computed with respect to the ground truth (red line), and an offline optimal fair allocation as computed with respect to the best Poisson fit to the ground truth (blue line). Note that in the limit, the allocations chosen by Algorithm 2 are guaranteed to converge to the blue baselines — but not the red baseline, because the algorithm is itself learning a Poisson approximation to the ground truth. The disparity between the red and blue lines quantifies the degradation in performance due to using Poisson approximations, rather than due to non-convergence of the learning process.

Figure 4 shows that Algorithm 2 converges to the Poisson approximation baseline well before the termination time of 2000, and substantially before the convergence bound guaranteed by our theory. Examining the estimated Poisson parameters used internally by Algorithm 2 reveals that although the allocation has converged to an optimal fair allocation, the estimated parameters have not yet converged to the parameters of the best Poisson fit in any of the districts. In particular, Algorithm 2 underestimates the parameters in all of the districts — but the degree of the underestimation is systematic: the correlation coefficient between the true and estimated parameters is 0.9975.

We see also in Figure 4 that convergence to the optimum expected number of discovered crimes occurs more quickly than convergence to the target fairness violation level. This is also apparent in Figure 3 where the learning and optimal Pareto curves are generally similar in terms of the maximum number of crimes discovered, while the fairness violations are higher in the learning curves.

5 The Random Discovery Model

Finally, we consider the random model of discovery. In the random model, when units are allocated to a group with candidates, the number of discovered candidates is a random variable corresponding to the number of candidates that appear in a uniformly random sample of individuals from a group of size . Equivalently, when units are allocated to a group of size with candidates, the number of candidates discovered by disc() is a random variable disc(where is drawn from the hypergeometric distribution with parameters and . Furthermore, the expected number of candidates discovered when allocating units to group i is E[disc(.

For simplicity, throughout this section, we assume for all i. This assumption can be completely relaxed (see the discussion in Appendix D). Moreover, let denote the expected fraction of candidates in group i. Without loss of generality, for the rest of this section, we assume .

5.1 Optimal Allocation

In this section, we characterize optimal allocations. Note that the expected number of candidates discovered by the allocation choice in group i is simply . This suggests a simple algorithm to compute : allocating every unit of the resource to group 1. More generally, let denote the subset of groups with the highest expected number of candidates. An allocation is optimal if and only if it only allocates all resources to groups in .

5.2 Properties of Fair Allocations

We next discuss the properties of fair allocations in the random discovery model. First, we point out that the discovery probability can be simplified as

So an allocation is -fair in the random model if for all groups i and j. Therefore, fair allocations (roughly) distribute resources in proportion to the size of the groups, essentially ignoring the candidate distributions within each group. We defer the full characterization to Appendix D.

5.3 Price of Fairness

Recall that PoF quantifies the extent to which constraining the allocation to satisfy -fairness degrades utility. While in Section 4 we study the PoF on the Philadelphia Crime Incidents dataset, we can define a worst-case variant as follows.

Definition 3. Fix the random model of crime discovery and let [0, 1]. We define the PoF as

where C ranges over all possible candidate distributions.

We can fully characterize this worst-case PoF in the random discovery model. We defer the proof of Theorem 4 to Appendix D.

Theorem 4. The PoF in the random discovery model is

The PoF in the random model can be as high as in the worst case. If all groups are identically sized, this grows linearly with the number of groups.

6 Conclusion and Future Directions

Our presentation of allocative fairness provides a family of fairness definitions, modularly parameterized by a “discovery model”. What counts as “fair” depends a great deal on the choice of discovery model, which makes explicit what would otherwise be unstated assumptions about the process of tasks like policing. The random and precision models of discovery studied in this paper represent two extreme points of a spectrum. In the predictive policing setting, the random model of discovery assumes that officers have no advantage over random guessing when stopping individuals for further inspection. The precision model assumes they can oracularly determine offenders, and stop only them. An interesting direction for future work is to study discovery models that lie in between these two.

We have also made a number of simplifying assumptions that could be relaxed. For example, we assumed the candidate distributions are stationary — fixed independently of the actions of the algorithm. Of course, the deployment of police officers can change crime distributions. Modeling this kind of dynamics, and designing learning algorithms that perform well in such dynamic settings would be interesting. Finally, we have assumed that the same discovery model applies to all groups. One friction to fairness that one might reasonably conjecture is that the discovery model may differ between groups — being closer to the precision model for one group, and closer to the random model for another. We leave the study of these extensions to future work.

Acknowledgements

We thank Sorelle Friedler for giving a talk at Penn which initially inspired this work. We also thank Carlos Scheidegger, Kristian Lum, Sorelle Friedler, and Suresh Venkatasubramanian for helpful discussions at an early stage of this work. We thank Richard Berk and Greg Ridgeway for helpful discussions about predictive policing. Finally, we thank the anonymous reviewer for helpful comments regarding Figure 2.

References

[1] Alekh Agarwal, Peter Bartlett, and Max Dama. Optimal allocation strategies for the dark pool problem. In Proceedings of the 13th International Conference on Artificial Intelligence and Statistics, pages 9–16, 2010.

[2] Hamsa Bastani, Mohsen Bayati, and Khashayar Khosravi. Exploiting the natural exploration in con- textual bandits. CoRR, abs/1704.09011, 2017.

[3] Richard Berk, Hoda Heidari, Shahin Jabbari, Matthew Joseph, Michael Kearns, Jamie Morgenstern, Seth Neel, and Aaron Roth. A convex framework for fair regression. CoRR, abs/1706.02409, 2017.

[4] Richard Berk, Hoda Heidari, Shahin Jabbari, Michael Kearns, and Aaron Roth. Fairness in criminal justice risk assessments: The state of the art. Sociological Methods & Research, 2018.

[5] Sarah Bird, Solon Barocas, Kate Crawford, Fernando Diaz, and Hanna Wallach. Exploring or exploiting? social and ethical implications of autonomous experimentation in AI. 2016.

[6] Toon Calders, Asim Karim, Faisal Kamiran, Wasif Ali, and Xiangliang Zhang. Controlling attribute effect in linear regression. In Proceedings of 13th International Conference on Data Mining, pages 71–80, 2013.

[7] Flavio Chierichetti, Ravi Kumar, Silvio Lattanzi, and Sergei Vassilvitskii. Fair clustering through fairlets. In Proceedings of the 31th Annual Conference on Neural Information Processing Systems, pages 5029–5037, 2017.

[8] Sam Corbett-Davies, Emma Pierson, Avi Feller, Sharad Goel, and Aziz Huq. Algorithmic decision making and the cost of fairness. In Proceedings of the 23rd ACM SIGKDD International Conference on Knowledge Discovery and Data Mining, pages 797–806, 2017.

[9] Cynthia Dwork and Christina Ilvento. Fairness under composition. In Proceedings of the 10th Innovations in Theoretical Computer Science Conference, 2019.

[10] Cynthia Dwork, Moritz Hardt, Toniann Pitassi, Omer Reingold, and Richard Zemel. Fairness through awareness. In Proceedings of the 3rd Innovations in Theoretical Computer Science, pages 214–226, 2012.

[11] Danielle Ensign, Sorelle Friedler, Scott Neville, Carlos Scheidegger, and Suresh Venkatasubramanian. Runaway feedback loops in predictive policing. In Conference on Fairness, Accountability and Transparency, pages 160–171, 2018.

[12] Danielle Ensign, Sorelle Frielder, Scott Neville, Carlos Scheidegger, and Suresh Venkatasubramanian. Decision making with limited feedback. In Proceedings of the 29th Conference on Algorithmic Learning Theory, pages 359–367, 2018.

[13] Kuzman Ganchev, Michael Kearns, Yuriy Nevmyvaka, and Jennifer Wortman Vaughan. Censored exploration and the dark pool problem. In Proceedings of the 25th Conference on Uncertainty in Artificial Intelligence, pages 185–194, 2009.

[14] Moritz Hardt, Eric Price, and Nathan Srebro. Equality of opportunity in supervised learning. In Proceedings of the 30th Annual Conference on Neural Information Processing Systems, pages 3315– 3323, 2016.

[15] David Inouye, Eunho Yang, Genevera Allen, and Pradeep Ravikumar. A review of multivariate distri- butions for count data derived from the poisson distribution. Wiley Interdisciplinary Reviews: Computational Statistics, 9, 2017.

[16] Shahin Jabbari, Matthew Joseph, Michael Kearns, Jamie Morgenstern, and Aaron Roth. Fairness in reinforcement learning. In Proceedings of the 34th International Conference on Machine Learning, pages 1617–1626, 2017.

[17] Matthew Joseph, Michael Kearns, Jamie Morgenstern, and Aaron Roth. Fairness in learning: classic and contextual bandits. In Proceedings of the 30th Annual Conference on Neural Information Processing Systems, pages 325–333, 2016.

[18] Sampath Kannan, Jamie Morgenstern, Aaron Roth, Bo Waggoner, and Zhiwei Steven Wu. A smoothed analysis of the greedy algorithm for the linear contextual bandit problem. In Proceedings of the 32nd Annual Conference on Neural Information Processing Systems, 2018.

[19] Jon Kleinberg, Sendhil Mullainathan, and Manish Raghavan. Inherent trade-offs in the fair determi- nation of risk scores. In Proceedings of the 8th Conference on Innovations in Theoretical Computer Science, pages 43:1–43:23, 2017.

[20] Lydia Liu, Sarah Dean, Esther Rolf, Max Simchowitz, and Moritz Hardt. Delayed impact of fair machine learning. In Proceedings of the 35th International Conference on Machine Learning, pages 3156–3164, 2018.

[21] Kristian Lum and William Isaac. To predict and serve? Significance, pages 14–18, October 2016.

[22] David MacKay. Information Theory, Inference and Learning Algorithms. Cambridge University Press, 2003.

[23] Ariel Procaccia. Cake cutting: Not just child’s play. Communications of the ACM, 56(7):78–87, 2013.

[24] Manish Raghavan, Aleksandrs Slivkins, Jennifer Wortman Vaughan, and Zhiwei Steven Wu. The exter- nalities of exploration and how data diversity helps exploitation. In Proceedings of the 31st Conference On Learning Theory, pages 1724–1738, 2018.

[25] Blake Woodworth, Suriya Gunasekar, Mesrob Ohannessian, and Nathan Srebro. Learning nondiscriminatory predictors. In Proceedings of the 30th Conference on Learning Theory, pages 1920–1953, 2017.

[26] Muhammad Bilal Zafar, Isabel Valera, Manuel Gomez-Rodriguez, and Krishna Gummadi. Fairness beyond disparate treatment & disparate impact: Learning classification without disparate mistreatment. In Proceedings of the 26th International Conference on World Wide Web, pages 1171–1180, 2017.

[27] Rich Zemel, Yu Wu, Kevin Swersky, Toni Pitassi, and Cynthia Dwork. Learning fair representations. In Proceedings of the 30th International Conference on Machine Learning, pages 325–333, 2013.

A Feasibility in Expectation

In this section, we show how to compute for any arbitrary but known candidate distributions C and known discovery model disc() in a relaxation where the feasibility constraint is satisfied in expectation.

The first observation is that when disc() and C are both known, for a group i and allocation of j units of resource to that group, the expected number of discovered candidates

and the discovery probability

can both be computed exactly. The second observation is that when allowing the feasibility condition to be satisfied in expectation, instead of allocating integral units of resources to each group, we can allocate resources to a group using a distribution.

Let denote the probability that j units of resource is allocated to group i. We can compute by

writing the following linear program with s as variables.

The objective function maximizes the number of candidates discovered given the allocation. The first constraint guarantees that the allocation is feasible in expectation. The second constraint (which is linear in ) ensures that -fairness is satisfied by the allocation. The last two constraints guarantees that for any values define a valid probability distribution on all the possible allocations to group i.

B Omitted Details from Section 3

B.1 Omitted Details from Section 3.1

We first show how the expected number of discovered candidates in a group in the precision model can be written as a function of the tail probabilities of the group’s candidate distribution.

Lemma 5 (Ganchev et al. [13]). The expected number of discovered candidates in the precision model when allocating units of resource to group i can be written as [min ().

Proof.

Note that we can perform the telescoping in the 3rd and 4th lines by observing that Pr1] + 1).

We then show that a greedy algorithm would find an optimal allocation in the precision model when the candidate distributions are known.

Theorem 5 (Theorem 1 in Ganchev et al. [13]). The allocation returned by greedily allocating the next unit of resource to a group in

where is the current allocation to group i at the round maximizes the expected number of candidates discovered.

Proof. Since the tail probability functions ) are all non-increasing (that is, for , we have )), the greedy allocation returns an allocation v which maximizes

Using Lemma 5 we have that