Plurality-voting ensembles aggregate the predictions of multiple base classifiers and output the most voted class, with ties typically broken at random. Despite the simple aggregation rule, attributing the impact of each base classifier to the ensemble's predictions is far from straightforward. That is, a single base classifier's impact depends on how it interacts with the votes of all other classifiers in the ensemble, for example, in resolving ties, and not only on its own predictions. We address the attribution problem using established tools from cooperative game theory, studying the Banzhaf and Shapley power indices for new games modeling plurality voting and random tie-breaking, for binary and multiclass classification settings. Computing power indices is, in general, intractable. Surprisingly, our methods, collectively named ENPOWER, enable the exact computation of Banzhaf and Shapley values in linear time for the binary classification setting and in low-degree polynomial time for the multiclass setting. The algorithms in ENPOWER rely on new combinatorial formulations of the Banzhaf and Shapley indices, yielding closed-form solutions for the binary case and ordinary generating function representations for the multiclass case. We validate empirically the methods in ENPOWER demonstrating their high applicability, providing attribution values not captured by existing methods. We apply the attribution of ENPOWER to: ensemble pruning and prompt-based ensembling of vision-language models.
Figures & tables
Figure 1: Runtime (mean and standard deviation) on our largest datasets, binary classification.
Figure 2: For each pair (dataset, baseline) encoded by a pair (shape, color), we display the Kendall’s τ with respect to the exact rankings obtained by enpower - ϕ and enpower - β on the x -axis, and y -axis respectively, averaged over five splits. Each plot corresponds to a different ensemble.
Figure 3: Generalization metrics for sub-ensemble selection, averaged over five random splits.
Figure 4: Accuracy and F1 -score test metrics for the prompt ensembling of CLIP on Dtest (see Sec. 5.4 ), for varying k . We report mean and standard deviation (over multiple splits).
Appendix figures & tables16 assets
Supplementary material from the paper’s appendix.
Appendix
Dataset (Label)
∣Dval∣+∣Dtest∣
Classes ( ℓ )
y=1
y=0
Blood Transfusion ( BTSC )
187
2
143
44
Credit-G ( CR )
250
2
179
71
CNAE-9 ( C9 )
270
9
38
232
Titanic ( TT )
328
2
210
118
Bank Marketing ( BM )
11 303
2
9 994
1 309
Covertype ( CT )
145 253
7
70 628
74 625
Appendix
Table 1: Datasets used for binary and multiclass classification settings in Sec. 5 . Classes ( ℓ ), denotes the original number of labels, which are also retained for the multiclass experiments. The last two columns display the number of points in Dval∪Dtest with label y=1 and y=0 after binarization.
Dataset
Total images
Classes ( ℓ )
∣Dval∣
∣Dtest∣
∣U∣
CIFAR (Cifar-100)
10 000
100
5 000
5 000
675.2±6.6
PETS (Oxford-IIIT Pets)
3 669
37
1 834
1 835
123.6±5.4
IMAGENET-A (ImageNet-A)
7 500
200
3 750
3 750
861.0±11.6
Appendix
Table 2: Dataset statistics for the zero-shot multiclass classification experiment. For each dataset we display its total number of images, and the number of possible classes. We also show ∣Dval∣ , used to compute the ranking metrics, and ∣Dtest∣ used to evaluate the generalization of model selection. Column ∣U∣ reports the number of unique vote profiles observed on Dval (the ensemble has n=22 base predictors), averaged over five splits. Other columns do not depend on the split.
Dataset
Classes ( ℓ )
min
max
IR
CV (%)
Hnorm
CIFAR (Cifar-100)
100
100
100
1.00
0.0
1.000
PETS (Oxford-IIIT Pets)
37
88
100
1.14
2.6
1.000
IMAGENET-A (ImageNet-A)
200
3
100
33.33
78.2
0.946
Appendix
Table 3: Class-imbalance statistics for the multiclass datasets of Sec. 5.4 (computed on Dval∪Dtest , hence equal on every split). ℓ is the number of classes. “ min / max ” are the smallest/largest number of points (images) with the same class; IR=max/min is the imbalance ratio; CV is the coefficient of variation computed over the cardinality of all classes; and Hnorm∈(0,1] is the normalized Shannon entropy of the class distribution ( 1 corresponds to a uniform distribution).
Figure 5: Sample images from the three datasets considered in Sec. 5.4 . CIFAR consists of tiny, low resolution (32 × 32 pixel) images. PETS consists of high quality images of different pets with fine-grained classes, and IMAGENET-A is a notorious challenging dataset for image classification.
Dataset
Ensemble
Index
εabs
εrel
Kendall’s ( τ )
Spearman ( ρ )
BTSC
RF
ϕ
0.0045±0.0007
0.42±0.15
0.720±0.036
0.889±0.027
β
0.0003±0.0000
3.34±4.51
0.967±0.004
0.997±0.001
AB
ϕ
0.0013±0.0002
0.26±0.18
0.819±0.030
0.930±0.017
β
0.0007±0.0001
0.07±0.02
0.794±0.036
0.917±0.023
HET
ϕ
0.0111±0.0021
0.51±0.10
0.515±0.094
0.676±0.105
β
0.0005±0.0002
13.72±25.07
0.886±0.073
0.964±0.034
Appendix
Table 4: Approximation error of the MC baseline. For each dataset and power index we report the maximum absolute (εabs) and the maximum relative approximation error (εrel) , over all models from C . We further report the Kendall’s ( τ ) and Spearman ( ρ ) ranking coefficients, with the exact scores computed by enpower . The results are the average and standard deviation over five random splits.
Ensemble
Method
C9
CT
RF ( n=50 )
enpower - β
(3.7±0.4)⋅10−3
0.10±0.01
enpower - ϕ
0.14±0.02
2.6±0.2
enpower MC-ϕ ++
(2.8±0.4)⋅10−2
0.48±0.04
AB ( n=150 )
enpower - β
(5.0±0.9)⋅10−2
1.9±0.4
enpower - ϕ
7.8±1.4
265±50
enpower MC-ϕ ++
2.7±0.5
77±14
Appendix
Table 5: Runtime (seconds) of enpower on Dval in the multiclass setting. The results are the average and standard deviation over five random splits. enpower MC-ϕ ++ corresponds to the refined algorithm of App. H .
Figure 6: Runtime comparison (mean and standard deviation) on the missing datasets from Table 1 not reported in Sec. 5 due to space constraints.
Ensemble
BTSC
CR
C9
TT
BM
CT
RF
0.697±0.033
0.648±0.035
0.979±0.014
0.884±0.032
0.868±0.003
0.886±0.005
AB
0.510±0.251
0.505±0.176
0.534±0.376
0.535±0.156
0.545±0.311
0.519±0.025
HET
0.745±0.032
0.700±0.036
0.940±0.062
0.834±0.121
0.895±0.008
0.768±0.068
Appendix
Table 6: Accuracy on Dval of the base classifiers of each ensemble, binary classification. Each entry is the mean ± the standard deviation over the base classifiers of the ensemble, both averaged over five random splits.
Figure 7: Generalization metrics for sub-ensemble selection, missing results from Sec. 5.3 .
Figure 8: Accuracy and F1 -score test metrics for the prompt ensembling of CLIP on Dtest (see Sec. 5.4 ), as a function of k . Missing results from Sec. 5.4 .
Dataset
Ranking
best acc.
k
avg. acc.
best F1
k
avg. F1
CIFAR100
Acc
0.648±0.005
14
0.645±0.004
0.641±0.004
16
0.637±0.003
enpower - β
0.649±0.004
10
0.646±0.004
0.641±0.003
16
0.638±0.003
enpower - ϕ
0.649±0.004
8
0.646±0.004
0.641±0.004
17
0.638±0.003
Tr-PS
0.647±0.005
15
0.643±0.006
0.639±0.004
15
0.635±0.004
Tr-MLE
0.648±0.004
14
0.643±0.005
0.640±0.003
14
0.636±0.003
IMAGENET-A
Acc
0.308±0.004
4
0.302±0.004
0.288±0.005
1
0.279±0.003
Appendix
Table 7: Prompt ensembling of CLIP, metrics on Dtest of the sub-ensemble for varying k (see Sec. 5.4 ). Columns report, best: the maximum over k of the average over the five random splits, and the k attaining it. Avg.: the average over k∈[n] . Bold marks the best value for each configuration.
Method
CIFAR100
IMAGENET-A
PETS
enpower - β
(1.17±0.01)⋅10−2
(1.37±0.04)⋅10−2
(1.64±0.09)⋅10−3
enpower - ϕ
(3.61±0.07)⋅10−1
(5.14±0.18)⋅10−1
(3.45±0.20)⋅10−2
Tr-PS
(1.33±0.13)⋅100
(8.64±0.44)⋅10−1
(3.77±0.04)⋅10−1
Tr-MLE
(2.37±0.05)⋅10−2
(1.84±0.08)⋅10−2
(1.26±0.02)⋅10−2
Appendix
Table 8: Runtime (seconds) of the ranking methods on Dval for the prompt ensembling of CLIP, average ± standard deviation over five random splits.
Prompt template
Accuracy
Banzhaf
Shapley
Δβ
Δϕ
Value
Rank
Value
Rank
Value
Rank
A pixelated photo of a {}.
0.6338
1
0.00296
1
0.03285
1
−
−
A heavily compressed JPEG photo of a {}.
0.6308
2
0.00262
2
0.03226
2
−
−
A high-resolution, crisp photo of a {}.
0.6248
3
0.00113
5
0.03047
4
↓ 2
↓ 1
A cropped photo of a {}.
0.6242
4
-0.00003
15
0.02941
10
↓ 11
↓ 6
A photo of a {}, cluttered background.
0.6238
5
0.00161
3
0.03077
3
↑ 2
↑ 2
Appendix
Table 9: Prompt ensemble evaluation for LLM zero-shot multiclass image classification (CIFAR100). Each prompt corresponds to a different model of C . All values and rankings are computed on the validation split Dval . Columns Δβ and Δϕ display the difference in the validation rankings induced by the power indices β(GMC) and ϕ(GMC) with respect to validation accuracy.
Prompt template
Accuracy
Banzhaf
Shapley
Δβ
Δϕ
Value
Rank
Value
Rank
Value
Rank
A poor quality, grainy photo of a {}.
0.8648
1
0.00434
1
0.04529
1
−
−
A blurry photo of a {}.
0.8539
2
0.00328
2
0.04335
2
−
−
A close-up photo of a {}.
0.8522
3
0.00327
3
0.04312
3
−
−
A cropped photo of a {}.
0.8506
4
0.00188
4
0.04175
4
−
−
A heavily compressed JPEG photo of a {}.
0.8462
5
0.00095
9
0.04059
7
↓ 4
↓ 2
Appendix
Table 10: Prompt ensemble evaluation for LLM zero-shot multiclass image classification (PETS). Each prompt corresponds to a different model of C . All values and rankings are computed on the validation split Dval . Columns Δβ and Δϕ display the difference in the validation rankings induced by the power indices β(GMC) and ϕ(GMC) with respect to validation accuracy.
Prompt template
Accuracy
Banzhaf
Shapley
Δβ
Δϕ
Value
Rank
Value
Rank
Value
Rank
A poor quality, grainy photo of a {}.
0.3136
1
0.00359
1
0.01936
1
−
−
A pixelated photo of a {}.
0.3029
2
0.00242
3
0.01738
3
↓ 1
↓ 1
A photo of a {}, with a watermark.
0.2989
3
0.00096
8
0.01581
6
↓ 5
↓ 3
A blurry photo of a {}.
0.2989
4
0.00184
4
0.01656
4
−
−
An unusual, hard to recognize photo of a {}.
0.2981
5
0.00295
2
0.01750
2
↑ 3
↑ 3
Appendix
Table 11: Prompt ensemble evaluation for LLM zero-shot multiclass image classification (IMAGENET-A). Each prompt corresponds to a different model of C . All values and rankings are computed on the validation split Dval . Columns Δβ and Δϕ display the difference in the validation rankings induced by the power indices β(GMC) and ϕ(GMC) with respect to validation accuracy.
Figure 9: Runtime to process a unique vote profile on synthetic instances, where p=0.1 (top), p=0.5 (middle) and p=0.9 (bottom). The y-axis reports the measured time to compute the power indices under a vote profile p∈U with n players and ℓ labels. See Sec. G.6 for a detailed description of each panel.
Shapley value and its priority-aware extensions are widely used for valuation in machine learning, but existing methods require pairwise priority to be binary and acyclic, a restriction spectacularly violated in real-data examples such as aggregated human preferences and multi-criterion comparisons. We introduce the generalized priority-aware Shapley value (GPASV), a random order value defined on arbitrary directed weighted priority graphs, in which pairwise edges penalize rather than forbid order violations. GPASV covers a range of classical models as boundary cases. We establish GPASV through an axiomatic characterization, develop the associated computational methods, and introduce a priority sweeping diagnostic extending PASV's. We apply GPASV to LLM ensemble valuation on the cyclic Chatbot Arena preference graph, illustrating that priority-aware valuation is not a one-button operation: different balances of pairwise graph priority versus individual soft priority produce substantively different valuations of the same data.
Kiljae Lee, Ziqi Liu, Weijing Tang +1
The Ohio State University · Carnegie Mellon University
Originating in game theory, Shapley values are widely used for explaining a machine learning model's prediction by quantifying the contribution of each feature's value to the prediction. This requires a scalar prediction as in binary classification, whereas a multiclass probabilistic prediction is a discrete probability distribution, living on a multidimensional simplex. In such a multiclass setting the Shapley values are typically computed separately on each class in a one-vs-rest manner, ignoring the compositional nature of the output distribution. In this paper, we introduce Shapley compositions as a well-founded way to properly explain a multiclass probabilistic prediction, using the Aitchison geometry from compositional data analysis. We prove that the Shapley composition is the unique quantity satisfying linearity, symmetry and efficiency on the Aitchison simplex, extending the corresponding axiomatic properties of the standard Shapley value. We demonstrate this proper multiclass treatment in a range of scenarios.
In this article, we provide an axiomatic characterization of feature attribution for multi-output predictors within the Shapley framework. While SHAP explanations are routinely computed independently for each output coordinate, the theoretical necessity of this practice has remained unclear. By extending the classical Shapley axioms to vector-valued cooperative games, we establish a rigidity theorem showing that any attribution rule satisfying efficiency, symmetry, dummy player, and additivity must necessarily decompose component-wise across outputs. Consequently, any joint-output attribution rule must relax at least one of the classical Shapley axioms. This result identifies a previously unformalized structural constraint in Shapley-based interpretability, clarifying the precise scope of fairness-consistent explanations in multi-output learning. Numerical experiments on a biomedical benchmark illustrate that multi-output models can yield computational savings in training and deployment, while producing SHAP explanations that remain fully consistent with the component-wise structure imposed by the Shapley axioms.
Umberto Biccari, Alain Ibáñez de Opakua, José María Mato +3