Minimization of Gini impurity via connections with the $k$-means problem

Abstract

The Gini impurity is one of the measures used to select attribute in Decision Trees/Random Forest construction. In this note we discuss connections between the problem of computing the partition with minimum Weighted Gini impurity and the $k$-means clustering problem. Based on these connections we show that the computation of the partition with minimum Weighted Gini is a NP-Complete problem and we also discuss how to obtain new algorithms with provable approximation for the Gini Minimization problem.

1 Introduction

Decision Trees and Random Forests are among the most popular methods for classification tasks. It is widely known that decision trees, specially small ones, are easy to interpret while random forest usually yield to more stable/accurate classifications.

A key decision during the construction of these structures is the selection of the attribute that is used for branching at each node. The standard approach for this selection is to evaluate the ability of each attribute to generate ’pure’ partitions, that is, partitions in which each branch is very homogeneous with respect to the class distribution of its examples. To measure how impure each branch is, impurity measures are often employed. An impurity measure maps a vector $\mathbf{u}=(u_{1},\ldots,u_{k})$, counting how many examples of each class we have in a node (branch), into a non-negative scalar . Arguably, two of the most classical impurity measures are the Gini impurity

$$i_{Gini}(\mathbf{u})=\sum_{i=1}^{k}\frac{u_{i}}{\|\mathbf{u}\|_{1}}\left(1-\frac{u_{i}}{\|\mathbf{u}\|_{1}}\right),$$

which is used in the CART package [3], and the Entropy impurity

$$i_{Entr}(\mathbf{u})=-\sum_{i=1}^{k}\frac{u_{i}}{\|\mathbf{u}\|_{1}}\log\left(\frac{u_{i}}{\|\mathbf{u}\|_{1}}\right),$$

that along with its variants is used in the C4.5 decision tree inducer [8].

Given an attribute, the goal is then to find a split for the attribute values that induces a partition of the set of examples with minimum weighted impurity, where the weights are given by the number of examples that lie into each of the branches.

Here, we discuss connections between the problem of computing the partition with minimum weighted Gini and the $k$-means clustering problem.

2 Connections between Gini minimization and $k$-means clustering

For a vector $\mathbf{v}$ where all components are non-negative the weighted Gini impurity $Gini(\mathbf{v})$ is defined as $Gini(\mathbf{v})=\|\mathbf{v}\|_{1}\cdot i_{Gini}(\mathbf{v})$. Let $A$ be a nominal attribute that may take $n$ possible values $a_{1},\ldots,a_{n}$. The $k$-ary Partition with Minimum Weighted Gini Problem ($k$-PMWGP) can be described abstractly as follows. We are given a collection of $n$ vectors $V\subset\mathbb{R}^{d}$, where the $i$th component of the $j$th vector counts the number of examples in class $i$ for which the attribute $A$ has value $a_{j}$. The goal is to find a partition ${\cal P}$ of $V$ into $k$ disjoint groups of vectors $V_{1},\ldots,V_{k}$ so as to minimize the sum of the weighted Gini impurities

$$Gini({\cal P})=\sum_{i=1}^{k}Gini\bigg(\sum_{\mathbf{v}\in V_{i}}\mathbf{v}\bigg).$$ (1)

Recently, we obtained simple constant approximation algorithms for this problem: an $O(n\log n+nd)$ time $2$-approximation for the case where $k=2$ [7] and a linear time $3$-approximation for arbitrary $k$ [5]. In fact these papers also handle a more general class of impurity measures that includes the Entropy impurity. The complexity of ($k$-PMWGP) remained open.

A problem that is equivalent to the above problem from the perspective of optimality but it is different from the perspective of approximation is the problem of finding the partition ${\cal P}$ of $V$ into $k$ groups that minimizes

$$Gini({\cal P})-\sum_{\mathbf{v}\in V}Gini(\mathbf{v}).$$ (2)

Using concavity properties of Gini one can prove that the above expression is always non-negative. An $\alpha$-approximation with respect to goal (2) implies an $\alpha$-approximation with respect to goal (1) but the converse is not necessarily true, so that approximations with respect to goal (2) are stronger.

In the geometric $k$-means problem we are given a set of vectors $V$ in $\mathbb{R}^{d}$ and the goal is to find a partition ${\cal P}$ of $V$ into $k$ groups $V_{1},\ldots,V_{k}$ and set of $k$ centers $\mathbf{c}_{1},\ldots,\mathbf{c}_{k}$ in $\mathbb{R}^{d}$ such that

$$Cost_{KM}({\cal P})=\sum_{i=1}^{k}\sum_{\mathbf{v}\in V_{i}}\|\mathbf{v}-\mathbf{c}_{i}\|_{2}^{2}$$

is minimized.

It is well known that if $U$ is a set of vectors then the vector $\mathbf{c}$ for which $\sum_{\mathbf{v}\in U}\|(\mathbf{v}-\mathbf{c})\|_{2}^{2}$ is minimum is the centroid of $U$, that is, $\mathbf{c}=(\sum_{\mathbf{v}\in U}\mathbf{v})/|U|$.

We argue that the following connections between $k$-PMWGP and $k$-means hold:

• Let $V$ be an instance of $k$-means where all vectors have the same $\ell_{1}$ norm. If ${\cal P}$ is an optimal partition for instance $V$ then ${\cal P}$ is also an optimal partition for instance $V$ of $k$-PMWGP

• there exists a pseudo-polynomial time reduction from $k$-PMWGP to the geometric $k$-means problem

The key observation for establishing C1 and C2 is the following lemma.

Lemma 1.

Let $X$ be a set of vectors, all of them with $\ell_{1}$ norm equal to $L$. Then,

$$Gini\bigg(\sum_{\mathbf{v}\in X}\mathbf{v}\bigg)-\sum_{\mathbf{v}\in X}Gini(\mathbf{v})=L\times\left(\sum_{\mathbf{v}\in X}\|\mathbf{v}-\mathbf{c}\|_{2}^{2}\right),$$

where $\mathbf{c}$ is the centroid of the set of vectors in $X$.

Proof.

Let $\mathbf{u}=\sum_{\mathbf{v}\in X}\mathbf{v}$. We have that

$$Gini(\mathbf{u})-\sum_{\mathbf{v}\in X}Gini(\mathbf{v})=\|\mathbf{u}\|_{1}\left(\sum_{i=1}^{d}\left(1-\frac{u_{i}}{\|\mathbf{u}\|_{1}}\right)\left(\frac{u_{i}}{\|\mathbf{u}\|_{1}}\right)\right)-\sum_{\mathbf{v}\in X}\sum_{i=1}^{d}\|\mathbf{v}\|_{1}\left(1-\frac{v_{i}}{\|\mathbf{v}\|_{1}}\right)\left(\frac{v_{i}}{\|\mathbf{v}\|_{1}}\right)$$

On the other hand,

$$\sum_{\mathbf{v}\in X}\|\mathbf{v}-\mathbf{c}\|_{2}^{2}=\sum_{i=1}^{d}\sum_{\mathbf{v}\in X}(v_{i}-c_{i})^{2}$$

Thus, it suffices to show that for any $i$

$$\|\mathbf{u}\|_{1}\left(1-\frac{u_{i}}{\|\mathbf{u}\|_{1}}\right)\left(\frac{u_{i}}{\|\mathbf{u}\|_{1}}\right)-\sum_{\mathbf{v}\in X}\|\mathbf{v}\|_{1}\left(1-\frac{v_{i}}{\|\mathbf{v}\|_{1}}\right)\left(\frac{v_{i}}{\|\mathbf{v}\|_{1}}\right)=L\left(\sum_{\mathbf{v}\in X}(v_{i}-c_{i})^{2}\right)$$

The left side is equal to

$$u_{i}-\frac{(u_{i})^{2}}{\|\mathbf{u}\|_{1}}-\left(\sum_{\mathbf{v}\in X}v_{i}-|X|\frac{\sum_{\mathbf{v}\in X}(v_{i})^{2}}{\|\mathbf{u}\|_{1}}\right)=\frac{|X|\sum_{\mathbf{v}\in X}(v_{i})^{2}}{\|\mathbf{u}\|_{1}}-\frac{(u_{i})^{2}}{\|\mathbf{u}\|_{1}}=\frac{\sum_{\mathbf{v}\in X}(v_{i})^{2}}{L}-\frac{(u_{i})^{2}}{L|X|}$$

Moreover, the righthand side is equal to

$$\sum_{\mathbf{v}\in X}(v_{i})^{2}-\frac{(\sum_{\mathbf{v}\in X}v_{i})^{2}}{|X|}=\sum_{\mathbf{v}\in X}(v_{i})^{2}-\frac{(u_{i})^{2}}{|X|},$$

which established the lemma. ∎

We should note that the result presented in the previous lemma is mentioned in the appendix of [4] where the Gini index is discussed.

The connection C1 is a direct consequence of Lemma 1 since it implies that for all $k$-partitions ${\cal P}$ of $V$

$$Gini({\cal P})=L\cdot Cost_{KM}({\cal P})+\sum_{\mathbf{v}\in V}Gini(\mathbf{v})$$

From the connection C1 and the hardness of geometric $k$-means established in [2] we obtain:

Theorem 1.

The Partition with Minimum Weighted Gini Problem (PMWGP) is NP-Complete with respect to goal (1) and APX-Hard with respect to goal (2).

Proof.

The result follows from [2], where a polynomial time reduction from the vertex cover problem on triangle free graphs to the $k$-means problem is presented. In this reduction, given a graph $G=(V,E)$, every edge $e$ in $E$ is mapped into a vector $\mathbf{v}$ in $\mathbb{R}^{|V|}$ where the $i$-th component $v_{i}$ is $1$ if $i$ is incident on $e$ and it is 0, otherwise. It is proved that if the minimum vertex cover of $G$ has size $k$ then the optimum cost of the corresponding $k$-means problem is at most $|E|-k$ and if the minimum vertex cover has size at least $(1+\epsilon)k$ then the minimum cost is at least $|E|-(1-\Omega(\epsilon)k$.

Our result follows from Lemma 1 and from the fact that in the instance of $k$-means above described all vectors have $\ell_{1}$ norm equals 2. ∎

With regards to the connection C2, let $V$ be an input of $k$-PMWGP and let $V^{\prime}$ be an instance of $k$-means obtained from $V$ as follows: for each vector $\mathbf{v}\in V$ we add to the input set $V^{\prime}$ exactly $\|\mathbf{v}\|_{1}$ copies of vector $\mathbf{v}^{\prime}=\mathbf{v}/\|\mathbf{v}\|_{1}$. Using Lemma 1 and also the fact that in any optimal solution for $k$-means identical vectors are in the same partition, we conclude that the optimum value of $V$ and $V^{\prime}$ differ by exactly $\sum_{\mathbf{v}\in V}Gini(\mathbf{v})$. Note that instance $V^{\prime}$ is obtained from $V$ in pseudo-poytime.

From this reduction one we can obtain new algorithms for $k$-PMWGP with provable approximation. As an example, we discuss how to obtain a PTAS for $k$-PMWGP with respect to the objective function (2) when $k$ is fixed. First, in our reduction, we keep each distinct vector and its multiplicity rather than all the $\sum_{\mathbf{v}\in V}\|\mathbf{v}\|_{1}$ vectors of instance $V^{\prime}$. Thus, we can construct the instance $V^{\prime}$ from $V$ in polytime. Next, we run over instance $V^{\prime}$ an adapted version of the recursive PTAS for $k$-means proposed in [6]. This version efficiently handles copies of the same vector and it is explained refering to the presentation of the PTAS that is given in Figure 1 of [1].

Let $W=\sum_{\mathbf{v}\in V}\|\mathbf{v}\|_{1}$. The adapted version is as follows:

1. The set of vectors in $V^{\prime}$ is represented using the distinct vectors and its multiplicities.

2. In the step 6 of the algorithm described in Figure 1 of [1] a constant number of vectors is sampled from a set of at most $n$ vectors. In our adaptation we sample from a set of $W$ vectors, with at most $n$ of them being distinct. This incurs an extra $O(\log W)$ factor to the running time.

3. At Step 12 of the same Figure one needs to compute, from a set of vectors $R$, the $|R|/2$ closest vectors to a given set of centroids $C$. Since this computation can be performed with time complexity proportional to the number of distinct vectors in $R$, rather than in $O(|R|)$ time, we do not incur any additional cost.

This adapted version also incurs an extra factor of $O(\log(W)^{k})$ with respect to the original one (executed over $n$ vectors) due to the number of nodes in the recursion tree. Thus, it runs in polynomial time when $k$ is fixed. Without efficiently handling the copies we would have a pseudo-polynomial time approximation scheme.