Despite significant progress in post-hoc explanation methods for neural networks, many remain heuristic and lack provable guarantees. A key approach for obtaining explanations with provable guarantees is by identifying a cardinally-minimal subset of input features which by itself is provably sufficient to determine the prediction. However, for standard neural networks, this task is often computationally infeasible, as it demands a worst-case exponential number of verification queries in the number of input features, each of which is NP-hard. In this work, we show that for Neural Additive Models (NAMs), a recent and more interpretable neural network family, we can efficiently generate explanations with such guarantees. We present a new model-specific algorithm for NAMs that generates provably cardinally-minimal explanations using only a logarithmic number of verification queries in the number of input features, after a parallelized preprocessing step with logarithmic runtime in the required precision is applied to each small univariate NAM component. Our algorithm not only makes the task of obtaining cardinally-minimal explanations feasible, but even outperforms existing algorithms designed to find the relaxed variant of subset-minimal explanations - which may be larger and less informative but easier to compute - despite our algorithm solving a much more difficult task. Our experiments demonstrate that, compared to previous algorithms, our approach provides provably smaller explanations than existing works and substantially reduces the computation time. Moreover, we show that our generated provable explanations offer benefits that are unattainable by standard sampling-based techniques typically used to interpret NAMs.
Figures & tables
Figure 1: Comparison of a subset-minimal explanation and a cardinally-minimal explanation of a NAM. Both explanations (in yellow) in (b) and (c) are minimal, as perturbing any additional feature can lead the overall output to become negative.
Figure 2: Sufficient explanations in NAMs: (a) Users must examine the neighborhood of an input for proper interpretation; e.g, users might wrongly conclude that feature 1 from the FICO HELOC dataset alone determines a positive output, but small changes in features 15 or 18 can flip the classification. (b) Outputs of continuous neighborhoods can be misleading if not verified, since sampling may miss extrema. For example, users might wrongly believe based on sampling that all features but feature 1 only have positive outputs and thus fixing feature 1 will always lead to a positive output of the entire NAM; however, the output of feature 2 can become so negative to flip the classification.
Breast Cancer
CREDIT
FICO HELOC
Method
Size ( ↓ )
Time [s] ( ↓ )
Size ( ↓ )
Time [s] ( ↓ )
Size ( ↓ )
Time [s] ( ↓ )
Ours
4.00±
4.24
35.60±
1.34
3.76±
2.62
132.67±
36.76
5.59±
1.80
317.92±
222.07
Lexicographic
16.58±
5.44
634.92±
77.23
12.42±
6.45
473.63±
128.38
15.60±
7.53
146.16±
188.07
Sensitivity
16.27±
5.57
636.79±
87.44
3.82±
1.84
407.93±
126.63
9.45±
5.90
250.09±
148.44
Table 1: Comparison of average explanation size and computation time.
Figure 3: Explanation size over time for all datasets.
CREDIT
FICO HELOC
Method
Size ( ↓ )
Time [s] ( ↓ )
Sufficiency [%] ( ↑ )
Size ( ↓ )
Time [s] ( ↓ )
Sufficiency [%] ( ↑ )
Ours
3.76±
2.62
132.67±
36.76
100.00
5.59±
1.80
317.92±
222.07
100.00
Sampling
2.67±
3.10
5.10±
0.26
131.37
3.04±
3.05
6.89±
2.57
125.49
Table 2: Comparison against a purely sampling-based approach.
Appendix figures & tables7 assets
Supplementary material from the paper’s appendix.
Appendix
Figure 4: The computation of the binary decomposition defined in Equation 33 , with input x and the output vector o0,…,on−1 representing the binary decomposition.
Breast Cancer
CREDIT
FICO HELOC
ϵ
Size ( ↓ )
Time [s] ( ↓ )
Size ( ↓ )
Time [s] ( ↓ )
Size ( ↓ )
Time [s] ( ↓ )
0.01
4.00±
4.24
35.60±
1.34
—
—
—
—
0.1
4.45±
3.98
115.08±
73.48
1.79±
1.07
76.76±
88.21
2.77±
1.62
131.31±
141.90
0.2
6.29±
3.50
143.81±
121.27
2.80±
2.07
97.16±
34.45
3.88±
1.58
136.82±
77.97
0.5
9.76±
2.75
146.01±
64.76
3.76±
2.62
132.67±
36.76
5.59±
1.80
317.92±
222.07
Appendix
Table 3: Varying perturbation radius ϵ on all datasets.
Figure 5: Number of processed features over time.
Figure 6: A comparison of explanation generation with different number of processors ρ .
Figure 7: Frequencies of the required precision to sort the features according to their importance.
Figure 8: Impact of the required precision on the sorting time.
Bias Shift
Sorting Time [s]
1
0.08
0.1
7.58
0.01
12.71
0.001
16.24
0.0001
20.95
0.00001
28.37
Appendix
Table 4: Influence of near identical features on the sorting time on synthetic example.