Authors: Emanuele Borgonovo, Francesco Cappelli, Xuefei Lu, Elmar Plischke, Cynthia Rudin
Organizations: Department of Decision Sciences and Bocconi Institute for Data Science and Analytics, Bocconi University, Milan, Italy · Department of Decision Science, Bocconi University, Milan, Italy · SKEMA Business School, Université Côte d’Azur, France · Clausthal University of Technology, Clausthal-Zellerfeld, Germany · Department of Computer Science, Duke University, Durham NC, 27708 USA
Feature importance methods using unrestricted permutations are flawed due to extrapolation errors; such errors appear in all non-trivial variable importance approaches. We propose three new approaches: conditional model reliance and Knockoffs with Gaussian transformation, and restricted ALE plot designs. Theoretical and numerical results show our strategies reduce/eliminate extrapolation.
Figures & tables
Figure 1: Original (black) and permuted data (red). (Left) with unrestricted permutations; (Middle) with residual imputation without Gaussian transformation; (Right) with residual imputation with Gaussian transformation (our proposed approach).
Figure 2: ALE design (left graph) and Jansen design (right graph) for a correlated case. Legend: ∙ original points; red ( ∘ ), new points after varying X1 ; blue ( ∘ ), new points after varying X2 .
Feature
X1 , X2
X3 , X4 , X5
X6
X7
X8
X9
X10
ρX1,X2=0.00
0.1667
0.1667
0
0.0417
0.1067
0.2400
0.3750
ρX1,X2=0.90
0.0332
0.1756
0
0.0438
0.1124
0.2528
0.3950
Table 1: Ground Truth values for νj=2τj in the Hooker test case.
Figure 3: Comparison between estimated feature importance and rankings using νj (the baseline approach), νjGCMR (ours) and νjGKnock (also ours) for the Hooker test case. The symbol ∘ indicates the theoretical rank of the features. Black is ground truth. Each of these lines is the average ranking over 50 replications.
Figure 4: Estimates (left subpanels) and ranking (right subpanels) of τjALE (left) and κjALE (right) for the Hooker test case.
Figure 5: Feature importance measures for the Linear Regression, νj , νGCMR and νGKnock .
Figure 6: Linear Regression: prediction densities with free permutations of X2 (left) and with permutations restricted by GCMR (right, our approach).
Figure 7: SHAP results for the Boston Housing case based on the trained linear regression. Panel (a) shows the SHAP values based on the linear regression model. Panel (b) shows the Predictions associated with calculating the SHAPs of TAX.
Figure 8: Permutation based importance measures νj , and νjGCMR (ours), νjGKnock (also ours) for the neural network.
Reliable estimation of feature contributions in machine learning models is essential for transparency, algorithmic fairness, and regulatory compliance. While permutation feature importance is widely used, classical implementations rely on repeated Monte Carlo shuffling, introducing significant computational overhead and stochastic instability. In this paper, we show that replacing B random permutations with a single, max-min rank-optimal deterministic permutation maintains or improves correlation with ground-truth importance while eliminating estimation variance and reducing complexity from O(B⋅n⋅p) to O(n⋅p). Under location-scale feature distributions, we formally prove exact recovery of scale-adjusted linear regression coefficients, alongside improved importance estimation under concave model sensitivity. We extend this deterministic framework along two complementary dimensions. First, Systemic Feature Importance (SFI) integrates empirical feature correlations to quantify indirect feature reliance through proxy variables. Second, Importance Direction extends scalar importance to a signed, directional representation by measuring concordance between covariate displacements and output shifts. Extensive empirical validation across nearly 200 simulation scenarios demonstrates superior bias-variance trade-offs in high-dimensional and low signal-to-noise regimes. Finally, two real-world credit risk case studies show how coupling SFI with Importance Direction enables practitioners and regulators to audit models for both the magnitude and net sign of hidden reliance on protected attributes, delivering a principled, transparent, and scalable framework for model governance.
Feature importance is central to interpretable machine learning, but the term "importance" encompasses several fundamentally different notions of relevance. We develop a unified perspective based on null importance: a population-level characterization of when a feature is irrelevant under a specified notion of relevance. We consider standard notions of null importance arising from marginal and conditional statistical relevance, predictive risk, functional invariance, and causal effects, and show how these notions answer different scientific questions. We illustrate the framework in two applications in which the distinction is particularly consequential: algorithmic fairness, where common fairness criteria correspond to different notions of null importance, and genomic perturbation modeling, where different notions of relevance lead to different conclusions about what a prediction model has learned. The framework connects three aspects of feature analysis: the scientific question defining relevance, the data and model assumptions that shape how different null notions relate, and the methods used to assess importance. We establish sufficient conditions under which null notions coincide and give counterexamples showing how they diverge when those conditions fail. We then characterize which nulls different method families target and when their zero-importance statistics identify those targets. Finally, simulations spanning feature dependence, redundancy, nonlinearity, hidden features and other standard phenomena, along with case studies on image and multiomics data, provide empirical evidence for these theoretical distinctions and their practical consequences. Taken together, these results provide a common statistical language for relating scientific questions, data-generating assumptions, and algorithms, and clarify the conclusions that feature-importance analyses can support.
Garvesh Raskutti, Kris Sankaran, Jiaxin Ye
Department of Statistics, University of Wisconsin–Madison
Explainable artificial intelligence promises to yield insights into relevant features, thereby enabling humans to examine and scrutinize machine learning models or even facilitating scientific discovery. Considering the widespread technique of Shapley values, we find that purely data-driven operationalization of multivariate feature importance is unsuitable for such purposes. Even for simple problems with two features, spurious associations due to collider bias and suppression arise from considering one feature only in the observational context of the other, which can lead to misinterpretations. Causal knowledge about the data-generating process is required to identify and correct such misleading feature attributions. We propose cc-Shapley (causal context Shapley), an interventional modification of conventional observational Shapley values leveraging knowledge of the data's causal structure, thereby analyzing the relevance of a feature in the causal context of the remaining features. We show theoretically that this eradicates spurious association induced by collider bias. We compare the behavior of Shapley and cc-Shapley values on various, synthetic, and real-world datasets. We observe nullification or reversal of associations compared to univariate feature importance when moving from observational to cc-Shapley.