arXiv is now an independent nonprofit! Learn more
License: CC BY-NC-SA 4.0
arXiv:2406.00371v1 [cs.LG] 01 Jun 2024

Alternative Methods to SHAP Derived from Properties of Kernels: A Note on Theoretical Analysis

Kazuhiro Hiraki thanks: kazuhiro.hiraki86@gmail.com    Shinichi Ishihara thanks: ishihara5683@gmail.com    Junnosuke Shino thanks: Corresponding author: Waseda University, junnosuke.shino@waseda.jp
Abstract

This study first derives a general and analytical expression of AFA (Additive Feature Attribution) in terms of the kernel in LIME (Local Interpretable Model-agnostic Explanations). Then, we propose some new AFAs that have appropriate properties of kernels or that coincide with the LS prenucleolus in cooperative game theory. We also revisit existing AFAs such as SHAP (SHapley Additive exPlanations) and re-examine the properties of their kernels.
Keywords SHAP, LIME, Kernel, ES, FESP, LS prenucleolus, XAI, Alternatives to SHAP

1 Introduction

In the field of machine learning, Explainable Artificial Intelligence (XAI) refers to techniques and methods that make the decisions and predictions of machine learning models easier to understand. Among them, AFA (Additive Feature Attribution) is a method that decomposes a model’s prediction into the contributions of individual features. Notably, SHAP (SHapley Additive exPlanations), proposed by [5], which is based on the Shapley value [8] in cooperative game theory, is well-known in this context. Recently, research on SHAP has been rapidly expanding ([4]). To reduce the computational cost of SHAP, various methods such as Tree-SHAP[5] and Fast SHAP [3] have been proposed and applied to actual data (for example, [2]). As an alternative to SHAP, [1] considers ES (Equal Surplus) and FESP (Fair Efficient Symmetric Perturbation), both of which are based on solution concepts in cooperative game theory.

In this study, we investigate the relationship between AFA and the kernel in LIME (Local Interpretable Model-agnostic Explanations) as proposed by [6]. [5] characterizes SHAP in terms of the kernel (Kernel SHAP) and derive the expression of SHAP kernel explicitly. Intriguingly, the properties of the SHAP kernel seem different from those that the LIME kernel is expected to have. More specifically, in LIME, the kernel attaches a large weight as a perturbed sample gets closer to the instance being explained, which is different for that of SHAP. In this note, we first provide a general framework to relate an AFA with its associated kernel by deriving an analytical expression of an AFA in terms of its kernel. Then, we propose some new AFAs that have reasonable properties of kernels or that coincide with the LS prenucleolus in cooperative game theory. We also revisit existing AFAs such as SHAP (SHapley Additive exPlanations) and reexamine the properties of their kernels.

2 Preliminaries

Let t and nn be the number of the instances and the number of features, respectively. Suppose N={1,,n}N=\{1,...,n\}, T={1,,t}T=\{1,...,t\}. The feature input is a t×nt\times n matrix X=(X1,Xj,,Xn)X=(X_{1},...X_{j},...,X_{n}). The jjth feature vector is Xj=(x1,j,,xt,j)X_{j}=(x_{1,j},...,x_{t,j})^{\prime} and, for the τ\tauth instance of interest, the vector of features is xτ=(xτ,1,,xτ,j,,xτ,n)x_{\tau}=(x_{\tau,1},...,x_{\tau,j},...,x_{\tau,n}). Let ff be the original prediction model which takes xτx_{\tau} and produces a prediction.

For an element of the power set of NN, which is called a coalition in the cooperative game theory, S2NS\in 2^{N}, define xτ,S={xτ,j|jS}x_{\tau,S}=\{x_{\tau,j}|j\in S\}. xτ,Sx_{\tau,S} is a vector that consists of features in SS at τ\tauth instance. Similarly, for S2NS\in 2^{N}, define XS={Xj|jS}X_{S}=\{X_{j}|j\in S\}.

In cooperative game theory, a characteristic function form game is expressed as (N,v)(N,v) where N={1,,n}N=\{1,...,n\} is the set of players and vv is a real-valued function on the power set 2N2^{N}. For the τ\tauth instance and any coalition S2NS\in 2^{N}, when we define vτ:2Nv_{\tau}:2^{N}\longrightarrow\mathbb{R} as in (1), a characteristic function form game (N,vτ)(N,v_{\tau}) is specified for τ\tau:

vτ(S)=E[f^(xτ,S,XNS)].v_{\tau}(S)=E\left[\hat{f}(x_{\tau,S},X_{N\setminus S})\right]. (1)

vτ(S)v_{\tau}(S) is interpreted as the prediction that ff produces for the τ\tauth instance, when (i) features xτ,jx_{\tau,j} where jSj\in S are known but (ii) features xτ,kx_{\tau,k} where kNSk\in N\setminus S are unknown. Note that vτ(N)=E[f^(xτ,1,,xτ,n)]=f^(xτ,1,,xτ,n)v_{\tau}(N)=E\left[\hat{f}(x_{\tau,1},...,x_{\tau,n})\right]=\hat{f}(x_{\tau,1},...,x_{\tau,n}) and vτ()=E[f^(X1,,Xn)]=E[f^(X)]v_{\tau}(\emptyset)=E\left[\hat{f}(X_{1},...,X_{n})\right]=E\left[\hat{f}(X)\right], where the former is the prediction when all features at τ\tauth instances are known and the latter is the prediction when none of the features are known. It should be noted that, while standard cooperative game theory assumes that v()=0v(\emptyset)=0, this is not necessarily satisfied under this machine learning (ML) setting.

With this setup, Additive Feature attribution (AFA) is the method to decompose vτ(N)vτ()v_{\tau}(N)-v_{\tau}(\emptyset) into features at τ\tau, depending on their “contributions.” More precisely, for a characteristic function form game (N,vτ)(N,v_{\tau}) associated with the τ\tauth instance and for the feature (player) jj, define a real-valued function Ψτ(j):N\Psi_{\tau}(j):N\longrightarrow\mathbb{R}. We hereafter use Ψτ(j)\Psi_{\tau}(j) and Ψτ,j\Psi_{\tau,j} interchangeably and let Ψτ=(Ψτ,1,,Ψτ,n)\Psi_{\tau}=(\Psi_{\tau,1},...,\Psi_{\tau,n}). When Ψτ\Psi_{\tau} satisfies jNΨτ,j=vτ(N)vτ()\sum_{j\in N}\Psi_{\tau,j}=v_{\tau}(N)-v_{\tau}(\emptyset), then Ψτ\Psi_{\tau} is called Additive Feature Attribution (AFA), denoted by ΨτAFA\Psi_{\tau}^{AFA}.

3 A Brief Review on LIME and Kernel

Here we review [5] and [6], specifically the parts concerning the relationship between LIME and SHAP. In their notation, xx is the original representation of an instance being explained and zz is a perturbed sample from xx. They use a binary vector xx^{\prime} and a mapping x=hx(x)x=h_{x}(x^{\prime}), but in this study, just for simplicity, x=xx=x^{\prime} and z=zz=z^{\prime} i.e., the original instances are simplified ([5]), interpretable ([6]) or binary from the beginning.

[5] considers the following minimization problem (LIME, proposed by [6]).

ξ(x)=argmingGL(f,g,πx)+Ω(g)where\xi(x)=\mathop{\rm arg~min}\limits_{g\in G}L(f,g,\pi_{x})+\Omega(g)\ \ \ \ \ where
  • ff: the original prediction model.

  • gg: the explanation model defined as g(z)=ϕ0+i=1nϕizig(z)=\phi_{0}+\sum_{i=1}^{n}\phi_{i}z_{i}, where ϕi\phi_{i}\in\mathbb{R} and nn is the number of the features. Let GG be the set of all ggs and let ϕ=(ϕ1,,ϕn)n\phi=(\phi_{1},...,\phi_{n})\in\mathbb{R}^{n}.

  • xx: instance being explained.

  • zz: perturbed sample from xx. Let ZZ be the set of all zzs, including xx.

  • πx\pi_{x}: local kernel.

and LL is the loss function of the minimization problem and Ω(g)\Omega(g) is a measure of complexity of gg (a more complex gg is penalized). Note that, regarding the kernel πx\pi_{x}, [6] assumes it increases as the distance between xx and zz decreases, that is, as zz gets closer to xx, a larger weight is attached to zz.

Based on this setup, [5] assumes Ω(g)=0\Omega(g)=0 and L(f,g,πx)=zZ[f(z)g(z)]2πx(z)L(f,g,\pi_{x})=\sum_{z\in Z}\left[f(z)-g(z)\right]^{2}\pi_{x}(z). Therefore, the minimization problem of (2) is:

argmingGzZ[f(z)g(z)]2πx(z)\displaystyle\mathop{\rm arg~min}\limits_{g\in G}\sum_{z\in Z}\left[f(z)-g(z)\right]^{2}\pi_{x}(z) =\displaystyle= argminϕnzZ[f(z){ϕ0+i=1nϕizi}]2πx(z)\displaystyle\mathop{\rm arg~min}\limits_{\phi\in\mathbb{R}^{n}}\sum_{z\in Z}\left[f(z)-\left\{\phi_{0}+\sum_{i=1}^{n}\phi_{i}z_{i}\right\}\right]^{2}\pi_{x}(z) (2)
=\displaystyle= argminϕnzZ[i=1nϕizi{f(z)ϕ0}]2πx(z).\displaystyle\mathop{\rm arg~min}\limits_{\phi\in\mathbb{R}^{n}}\sum_{z\in Z}\left[\sum_{i=1}^{n}\phi_{i}z_{i}-\left\{f(z)-\phi_{0}\right\}\right]^{2}\pi_{x}(z).

Now recall zz is perturbed sample from xx and xZx\in Z. Therefore, summation over ZZ in (2) coincides with the summation over 2N2^{N} under our notation, and the summation of i=1nϕizi\sum_{i=1}^{n}\phi_{i}z_{i} coincides with iSϕi\sum_{i\in S}\phi_{i}. Therefore, under our notation, (2) falls into the following:

argminϕnS2N[iSϕi{vτ(S)vτ()}]2πxτ(S).\mathop{\rm arg~min}\limits_{\phi\in\mathbb{R}^{n}}\sum_{S\in 2^{N}\setminus\emptyset}\left[\sum_{i\in S}\phi_{i}-\left\{v_{\tau}(S)-v_{\tau}(\emptyset)\right\}\right]^{2}\pi_{x_{\tau}}(S). (3)

Furthermore, [5] impose a local accuracy condition (or called efficiency) on this optimization problem: f(x)=g(x)=ϕ0+i=1nϕixif(x)=g(x)=\phi_{0}+\sum_{i=1}^{n}\phi_{i}x_{i} for xx. We also impose this condition, so (3) becomes:

ΨτAFA=argminϕnwithiNϕi=vτ(N)vτ()S2N[iSϕi{vτ(S)vτ()}]2πxτ(S).\Psi_{\tau}^{AFA}=\mathop{\rm arg~min}\limits_{\phi\in\mathbb{R}^{n}with\sum_{i\in N}\phi_{i}=v_{\tau}(N)-v_{\tau}(\emptyset)}\sum_{S\in 2^{N}\setminus\emptyset}\left[\sum_{i\in S}\phi_{i}-\left\{v_{\tau}(S)-v_{\tau}(\emptyset)\right\}\right]^{2}\pi_{x_{\tau}}(S). (4)

For the following analysis, we consider 4 to examine relationship between SHAP and other AFAs alternatives to SHAP and the associated kernel πxτ(S)\pi_{x_{\tau}}(S).

4 Results

4.1 Conditions on Kernel

Regarding the kernel πxτ(S)\pi_{x_{\tau}}(S) in (4), we impose the following two conditions:

SN{i,j}πxτ(S{i})=1(i,jNwithij)\displaystyle\sum_{S\subseteq N\setminus\{i,j\}}\pi_{x_{\tau}}(S\cup\{i\})=1\qquad\Big(\forall i,j\in N\ with\ i\neq j\Big) (5)
πxτ(S)=πxτ(T)(S,T2Nwith|S|=|T|)\displaystyle\pi_{x_{\tau}}(S)=\pi_{x_{\tau}}(T)\qquad\Big(\forall S,T\in 2^{N}\ with\ |S|=|T|\Big) (6)

(5) is only necessary for normalization and is not essential. (6) states that, in terms of the number of features, when SS and TT are equidistant from NN, the kernel must assign the same weight to SS and TT. This can be considered a form of symmetry, which is a naturally acceptable condition.

4.2 Analytical Solution to the Optimization Problem

Here, we derive ΨτAFA\Psi_{\tau}^{AFA} in (4) analytically. The Lagrangian of (4) is:

(ϕ1,ϕn,λ)=S2N{}[iSϕi{vτ(S)vτ()}]2πxτ(S)λ[iNϕivτ(N)+vτ()].\mathcal{L}(\phi_{1},...\phi_{n},\lambda)=\sum_{S\in 2^{N}\setminus\{\emptyset\}}\left[\sum_{i\in S}\phi_{i}-\{v_{\tau}(S)-v_{\tau}(\emptyset)\}\right]^{2}\cdot\ \pi_{x_{\tau}}(S)-\lambda\left[\sum_{i\in N}\phi_{i}-v_{\tau}(N)+v_{\tau}(\emptyset)\right].

The F.O.C. on ϕj\phi_{j} is:

S2N{}:jS2(iSϕi{vτ(S)vτ()})πxτ(S)λ=0.\displaystyle\sum_{S\in 2^{N}\setminus\{\emptyset\}:j\in S}2\left(\sum_{i\in S}\phi_{i}-\{v_{\tau}(S)-v_{\tau}(\emptyset)\}\right)\cdot\ \pi_{x_{\tau}}(S)-\lambda=0.

Therefore, for any i,jNi,j\in N with iji\neq j, the following holds:

S2N{}:iS(kSϕk{vτ(S)vτ()})πxτ(S)=S2N{}:jS(kSϕk{vτ(S)vτ()})πxτ(S)\displaystyle\sum_{S\in 2^{N}\setminus\{\emptyset\}:i\in S}\left(\sum_{k\in S}\phi_{k}-\{v_{\tau}(S)-v_{\tau}(\emptyset)\}\right)\cdot\pi_{x_{\tau}}(S)=\sum_{S\in 2^{N}\setminus\{\emptyset\}:j\in S}\left(\sum_{k\in S}\phi_{k}-\{v_{\tau}(S)-v_{\tau}(\emptyset)\}\right)\cdot\pi_{x_{\tau}}(S)
\displaystyle\Longleftrightarrow SN{i,j}(kS{i}ϕk{vτ(S{i})vτ()})πxτ(S{i})\displaystyle\sum_{S\subseteq N\setminus\{i,j\}}\left(\sum_{k\in S\cup\{i\}}\phi_{k}-\{v_{\tau}(S\cup\{i\})-v_{\tau}(\emptyset)\}\right)\cdot\pi_{x_{\tau}}(S\cup\{i\})
=SN{i,j}(kS{j}ϕk{vτ(S{j})vτ()})πxτ(S{j})\displaystyle\qquad\qquad\qquad\qquad\qquad\qquad\qquad=\sum_{S\subseteq N\setminus\{i,j\}}\left(\sum_{k\in S\cup\{j\}}\phi_{k}-\{v_{\tau}(S\cup\{j\})-v_{\tau}(\emptyset)\}\right)\cdot\pi_{x_{\tau}}(S\cup\{j\})
\displaystyle\Longleftrightarrow SN{i,j}(πxτ(S{i})ϕiπxτ(S{j})ϕj)\displaystyle\sum_{S\subseteq N\setminus\{i,j\}}\bigg(\pi_{x_{\tau}}(S\cup\{i\})\cdot\phi_{i}-\pi_{x_{\tau}}(S\cup\{j\})\cdot\phi_{j}\bigg)
=SN{i,j}(πxτ(S{i}){vτ(S{i})πxτ(S{j}){vτ(S{j}))\displaystyle\qquad\qquad\qquad\qquad\qquad\qquad\qquad=\sum_{S\subseteq N\setminus\{i,j\}}\bigg(\pi_{x_{\tau}}(S\cup\{i\})\cdot\{v_{\tau}(S\cup\{i\})-\pi_{x_{\tau}}(S\cup\{j\})\cdot\{v_{\tau}(S\cup\{j\})\bigg)
\displaystyle\Longleftrightarrow ϕiϕj=SN{i,j}(πxτ(S{i}){vτ(S{i})πxτ(S{j}){vτ(S{j})),\displaystyle\phi_{i}-\phi_{j}=\sum_{S\subseteq N\setminus\{i,j\}}\bigg(\pi_{x_{\tau}}(S\cup\{i\})\cdot\{v_{\tau}(S\cup\{i\})-\pi_{x_{\tau}}(S\cup\{j\})\cdot\{v_{\tau}(S\cup\{j\})\bigg),

which implies:

ϕ1S:1SNπxτ(S)vτ(S)==ϕnS:nSNπxτ(S)vτ(S).\displaystyle\phi_{1}-\sum_{S:1\in S\neq N}\pi_{x_{\tau}}(S)\cdot v_{\tau}(S)=...=\displaystyle\phi_{n}-\sum_{S:n\in S\neq N}\pi_{x_{\tau}}(S)\cdot v_{\tau}(S). (7)

Therefore, ϕ=(ϕ1,,ϕj,,ϕn)\phi=(\phi_{1},...,\phi_{j},...,\phi_{n}) that satisfies (7) and jNϕj=vτ(N)vτ()\sum_{j\in N}\phi_{j}=v_{\tau}(N)-v_{\tau}(\emptyset) is derived as:

ϕj=S:jSNπxτ(S)vτ(S)+vτ(N)vτ()iN{S:iSNπxτ(S)vτ(S)}n.\displaystyle{\phi_{j}=\sum_{S:j\in S\neq N}\pi_{x_{\tau}}(S)\cdot v_{\tau}(S)+\frac{v_{\tau}(N)-v_{\tau}(\emptyset)-\sum_{i\in N}\left\{\sum_{S:i\in S\neq N}\pi_{x_{\tau}}(S)\cdot v_{\tau}(S)\right\}}{n}}. (8)

(8) expresses the AFA ϕ\phi as a function of the associated kernels πxτ(S)\pi_{x_{\tau}}(S) enabling the establishment of an AFA from any kernels. In the following sections, we examine several AFAs, some proposed by existing research, while others are newly proposed and generated by kernels having appropriate properties.

4.3 SHAP

In [5], the kernel of SHAP is specified as follows:

πxτ(S)=n1C|S|n|S|(n|S|).\pi_{x_{\tau}}(S)=\displaystyle{\frac{n-1}{{}_{n}C_{|S|}\cdot|S|\cdot(n-|S|)}}. (9)

Given that the solution of the minimization problem (4) remains unchanged even if a constant in the kernel is altered, we slightly modify the kernel as follows,

πxτ(S)=nC|S|n|S|(n|S|)\pi_{x_{\tau}}(S)=\displaystyle{\frac{n}{{}_{n}C_{|S|}\cdot|S|\cdot(n-|S|)}} (10)

so that (5) is satisfied. By substituting (10) into (8), we obtain the following:

Ψτ,jSHAP=ϕj=SNj|S|!(n|S|1)!n!(vτ(S{j})vτ(S)).\Psi_{\tau,j}^{SHAP}=\phi_{j}=\displaystyle\sum_{S\subseteq N\setminus j}\frac{|S|!(n-|S|-1)!}{n!}\left(v_{\tau}(S\cup\{j\})-v_{\tau}(S)\right).

That is, SHAP is derived as an AFA generated from the kernel expressed in (10). Therefore, it may be more appropriate to consider (10) rather than (9) as the kernel for SHAP. Additionally, it should be noted that the kernel of (9) or (10) reaches its maximum if |S|=0|S|=0 and |S|=n|S|=n, and it has a concave shape regarding |S||S|, which is different from [6] where the weight assigned by the kernel increases as a perturbed sample gets closer to the instance being explained.

4.4 ES and FESP in [1]

As alternative AFAs to SHAP, [1] proposes ES (Equal Surplus) and FESP (Fair Efficient Symmetric Perturbation), based on the solution concepts in cooperative game theory.

First, consider the following kernel:

πxτ(S)={1if|S|=10if2|S|n1ifS=N.\pi_{x_{\tau}}(S)=\left\{\begin{array}[]{ll}1&\mbox{if}\quad|S|=1\\ 0&\mbox{if}\quad 2\leq|S|\leq n-1\\ \infty&\mbox{if}\quad S=N.\\ \end{array}\right. (11)

Similarly to the previous case, by substituting (11) into (8), ϕi\phi_{i} becomes as follows:

Ψτ,jES=ϕj=vτ({j})+vτ(N)vτ()kNvτ({k})n.\Psi_{\tau,j}^{ES}=\phi_{j}=\displaystyle v_{\tau}(\{j\})+\frac{v_{\tau}(N)-v_{\tau}(\emptyset)-\sum_{k\in N}v_{\tau}(\{k\})}{n}.

That is, ϕi\phi_{i} coincides with ES.

Next, suppose the following kernel:

πxτ(S)={wτif|S|=10if2|S|n21wτif|S|=n1ifS=N.\pi_{x_{\tau}}(S)=\left\{\begin{array}[]{ll}w_{\tau}&\mbox{if}\quad|S|=1\\ 0&\mbox{if}\quad 2\leq|S|\leq n-2\\ 1-w_{\tau}&\mbox{if}\quad|S|=n-1\\ \infty&\mbox{if}\quad S=N.\\ \end{array}\right. (12)

Then, (8) follows that the associated solution of the minimization problem is FESP:

Ψτ,jFESP=ϕj=wτ(vτ({j})vτ())+(1wτ)(vτ()vτ(N\{j})).\Psi_{\tau,j}^{FESP}=\phi_{j}=\displaystyle w_{\tau}\Big(v_{\tau}(\{j\})-v_{\tau}(\emptyset)\Big)+(1-w_{\tau})\Big(v_{\tau}(\emptyset)-v_{\tau}(N\backslash\{j\})\Big).

Note that the kernel of (11) and (12) also have a concave shape regarding |S||S|.

4.5 AFA based on LS preucleolus

Consider the following kernel:

πxτ(S)={12n2ifSNifS=N.\pi_{x_{\tau}}(S)=\left\{\begin{array}[]{ll}\displaystyle{\frac{1}{2^{n-2}}}&\mbox{if}\quad S\neq N\\ \infty&\mbox{if}\quad S=N.\end{array}\right. (13)

Note that the shape of this kernel is not concave with respect to |S||S|, although it is still different from [6] in that the shape is flat unless S=NS=N. By substituting (13) into (8), the resulting ϕi\phi_{i} is:

Ψτ,jPNucl=ϕj=2(12n1S:jSvτ(S))+vτ(N)vτ()iN{2(12n1S:iSvτ(S))}n\Psi_{\tau,j}^{PNucl}=\phi_{j}=\displaystyle 2\left(\frac{1}{2^{n-1}}\displaystyle\sum_{S:j\in S}v_{\tau}(S)\right)+\frac{v_{\tau}(N)-v_{\tau}(\emptyset)-\sum_{i\in N}\left\{2\left(\frac{1}{2^{n-1}}\sum_{S:i\in S}v_{\tau}(S)\right)\right\}}{n}

Intriguingly, this solution is identical to that in the following minimization problem in which a kernel does not appear, thus coinciding with the LS prenucleolus proposed by [7]:

argminϕRn:iNϕi=vτ(N)vτ()S2N[iSϕi{vτ(S)vτ()}]2.\displaystyle\mathop{\rm arg~min}\limits_{\phi\in R^{n}:\sum_{i\in N}\phi_{i}=v_{\tau}(N)-v_{\tau}(\emptyset)}\sum_{S\in 2^{N}\setminus\emptyset}\left[\sum_{i\in S}\phi_{i}-\{v_{\tau}(S)-v_{\tau}(\emptyset)\}\right]^{2}.

4.6 AFA with a reasonable kernel (I)

The next kernel we consider is as follows:

πxτ(S)={|S|n2n3ifSNifS=N.\pi_{x_{\tau}}(S)=\left\{\begin{array}[]{ll}\displaystyle{\frac{|S|}{n\cdot 2^{n-3}}}&\mbox{if}\quad S\neq N\\ \infty&\mbox{if}\quad S=N.\\ \end{array}\right. (14)

This kernel satisfies the conditions of (5) and (6). Furthermore, this is increasing in |S||S| and thus consistent with the condition on the kernel in [6]. By substituting (14) into (8), we have:

Ψτ,jLnK=ϕj=S:jSN|S|n2n3vτ(S)+vτ(N)vτ()iN{S:iSN|S|n2n3vτ(S)}n,\Psi_{\tau,j}^{LnK}=\displaystyle{\phi_{j}=\sum_{S:j\in S\neq N}\frac{|S|}{n\cdot 2^{n-3}}\cdot v_{\tau}(S)+\frac{v_{\tau}(N)-v_{\tau}(\emptyset)-\sum_{i\in N}\left\{\sum_{S:i\in S\neq N}\frac{|S|}{n\cdot 2^{n-3}}\cdot v_{\tau}(S)\right\}}{n}},

which is the first AFA we propose as an alternative to SHAP. The superscript LnKLnK stands for linealy increasing kernel.

4.7 AFA with a reasonable kernel (II)

Lastly, assume the following kernel:

πxτ(S)={2|S|13n2ifSNifS=N.\pi_{x_{\tau}}(S)=\left\{\begin{array}[]{ll}\displaystyle{\frac{2^{|S|-1}}{3^{n-2}}}&\mbox{if}\quad S\neq N\\ \infty&\mbox{if}\quad S=N.\\ \end{array}\right.

This kernel also satisfies (5) and (6) and is increasing in |S||S|. More specifically, each time |S||S| increases by 11, the value of the kernel doubles. Then, we get the following expression, which is our second proposed AFA alternative to SHAP.

Ψτ,jExK=ϕj=S:jSN2|S|13n2vτ(S)+vτ(N)vτ()iN{S:iSN2|S|13n2vτ(S)}n\Psi_{\tau,j}^{ExK}=\displaystyle{\phi_{j}=\sum_{S:j\in S\neq N}\frac{2^{|S|-1}}{3^{n-2}}\cdot v_{\tau}(S)+\frac{v_{\tau}(N)-v_{\tau}(\emptyset)-\sum_{i\in N}\left\{\sum_{S:i\in S\neq N}\frac{2^{|S|-1}}{3^{n-2}}\cdot v_{\tau}(S)\right\}}{n}}

The superscript ExKExK stands for exponentially increasing kernel.

5 Conclusion

In this study, we first derive an analytical and general expression of an AFA as a function of its associated kernel. Next, we compute several AFAs based on representations of several different specific kernels. Among the existing AFAs, we show that for SHAP, by slightly modifying the kernel into an appropriate form, the generated AFA coincides with SHAP. Additionally, for ES and FESP, we derive the representations of the corresponding kernels. The last three kernels and the AFAs generated from them are proposed for the first time in this study. ΨτPNucl\Psi_{\tau}^{PNucl} has a kernel that is not concave and coincides with the notion of the LS prenucleolus in the cooperative game theory. ΨτLnK\Psi_{\tau}^{LnK} and ΨτExK\Psi_{\tau}^{ExK} are generated as AFAs from kernels that have desirable properties and consistent with the idea from [6] that the kernel assigns a large weight as a perturbed sample gets closer to the instance being explained.

The extent to which these AFAs show different decomposition patterns in experiments using actual data is an empirical question of great importance and one that should be addressed promptly. Another important theme is how the newly presented Ψτ,jLnK\Psi_{\tau,j}^{LnK} and Ψτ,jExK\Psi_{\tau,j}^{ExK} in this study can be characterized from the perspective of cooperative game theory, for example, whether they can be axiomatized, is also worth investigating.

References

  • [1] C. Condevaux, S. Harispe, and S. Mussard. 2023. Fair and Efficient Alternatives to Shapley-based Attribution Methods. Joint European Conference on Machine Learning and Knowledge Discovery in Databases.
  • [2] S. B. Jabeur, S. M Wali, and J-L. Viviani. 2024. Forecasting golod price with the XGBoost algorithm and SHAP interaction values. Annals of Operational Research 334: 679-699.
  • [3] N. Jethani, M. Sudarshan, I.C. Covert, S.-I. Lee, and R. Ranganath. 2021 Fastshap: real-time Shapley value estimation. In International Conference on Learning Representations. S. B. Jabeur, S. M Wali, and J-L. Viviani. 2024. Forecasting golod price with the XGBoost algorithm and SHAP interaction values. Annals of Operational Research 334: 679-699.
  • [4] M. Li, H. Sun, Y. Huang and H. Chen. 2024. Shapley value: from cooperative game to explainable artificial intelligence. Autonomous Intelligent Systems 4, No.2.
  • [5] S. M. Lundberg and S-I Lee. 2016. A unified approach to interpreting model predictions. Advances in Neural Information Processing Systems 30.
  • [6] Ribeiro, Marco Tulio, Singh, Sameer, and Guestrin, Carlos. 2016. Why Should I Trust You?: Explaining the Predictions of Any Classifier. Proceedings of the 22Nd ACM SIGKDD International Conference on Knowledge Discovery and Data Mining: 1135–1144. New York, NY, USA. ACM. ISBN 978-1- 4503-4232-2. doi: 10.1145/2939672.2939778.
  • [7] L. M Ruiz, F. Valenciano, and J. M. Zarzuelo. 1996. The Least Square Prenucleolus and the least Square Nucleolus. Two values for TU Games Based on the Excess Vector. International Journal of Game Theory 25: 113–134.
  • [8] L. S. Shapley. 1953. A Value for n-Person Games. Annals of Mathematics Studies 28: 307–318.