Consistent Individualized Feature Attribution for Tree Ensembles

Scott M. LundbergGabriel G. ErionSu-In Lee

article2018arXiv1,957 citations

Develops a fast, exact algorithm for computing SHAP values on tree ensembles, resolving fundamental inconsistency errors in traditional feature importance methods to provide mathematically sound, individualized explanations for models like XGBoost and Random Forests.

Listen

Tree-based machine learning models, such as gradient boosting machines and random forests, are widely deployed across high-stakes industries like healthcare and finance. However, explaining why these complex models make specific predictions remains a major challenge. The article addresses a critical flaw in standard feature attribution techniques, showing that widely used measures such as split count, gain, and existing individual tree explanation methods are inconsistent. Inconsistency means that when a model is modified to rely more heavily on a given feature, the importance score assigned to that feature can paradoxically decrease, making reliable comparisons between features impossible.

The article sets out to evaluate the mathematical consistency of existing tree interpretation methods and demonstrate a computationally efficient, exact approach based on game-theoretic principles for explaining individualized predictions in tree ensembles.

To achieve this, the authors applied SHapley Additive exPlanation (SHAP) values—the only additive attribution method mathematically proven to satisfy local accuracy, missingness, and consistency. Because standard exact calculations require evaluating an exponential number of feature combinations, the authors derived an exact, polynomial-time algorithm called Tree SHAP tailored to tree structures. They validated the method through mathematical proofs, runtime benchmarks on large ensemble models, human intuitive alignment studies via crowdsourced experiments, feature perturbation tests on sentiment data, and real-world demonstrations across census data and long-term mortality data from the NHANES I health study.

The findings establish that conventional gain, split count, and heuristic tree attribution methods frequently produce misleading feature rankings. In contrast, Tree SHAP provides guaranteed consistency and reduces computational complexity from exponential time to low-order polynomial time, allowing models with 1,000 deep trees and 100 features to be explained in just 0.08 seconds. In human evaluations, attributions calculated via SHAP matched human intuition significantly better than existing methods. Furthermore, the framework extends to pairwise interaction effects through SHAP interaction values and introduces powerful analytical tools, such as supervised clustering and rich dependence plots that separate main effects from complex feature interactions.

These findings have direct operational and governance implications. Organizations relying on flawed attribution metrics face hidden risks when making high-impact decisions, conducting feature selection, or ensuring regulatory compliance. By resolving inconsistency, Tree SHAP establishes a trustworthy, standardized baseline for model interpretability. The massive speed improvement also removes previous performance bottlenecks, enabling real-time, individualized explanations and full-dataset visual auditing at enterprise scale.

Organizations developing or deploying tree ensemble models should transition away from traditional split count and gain metrics toward consistent SHAP-based attributions for model debugging, reporting, and feature selection. Teams can readily adopt the open-source implementation integrated into major machine learning frameworks like XGBoost and LightGBM.

Confidence in the mathematical properties and computational speed of Tree SHAP is very high based on formal game-theoretic foundations and experimental validation. However, decision-makers should recognize that the algorithm explains the mechanics of the trained model itself, which may still inherit underlying data biases or spurious correlations present in the training set.

arXiv: 1802.03888slundberg/shap
  • Paper: A Unified Approach to Interpreting Model Predictions, Scott M. Lundberg et al. (2017). Introduces the unified SHAP framework and proves the axiomatic necessity of Shapley values for consistent additive feature attributions, which this paper adapts into an exact, polynomial-time algorithm for tree ensembles.
  • Paper: XGBoost: A Scalable Tree Boosting System, Tianqi Chen et al. (2016). Presents the foundational tree boosting architecture and objective formulation implemented by XGBoost, the primary tree ensemble system into which this paper embeds its exact SHAP algorithm.
  • Paper: LightGBM: A Highly Efficient Gradient Boosting Decision Tree, Guolin Ke et al. (2017). Describes the high-efficiency gradient boosted decision tree framework that serves as one of the key target ensemble systems integrated with the proposed Tree SHAP algorithm.
  • Paper: “Why Should I Trust You?”: Explaining the Predictions of Any Classifier, Marco Tulio Ribeiro et al. (2016). Establishes local surrogate models for individualized classifier explanation (LIME), motivating the need for consistent and locally accurate feature attribution methods.
  • Paper: Greedy function approximation: A gradient boosting machine, Jerome H. Friedman (2001). Formulates gradient boosting and tree-based additive expansions, providing the underlying predictive modeling paradigm that this paper aims to interpret consistently.
  • Paper: Random Forests, Leo Breiman (2001). Introduces random forests and classical impurity-based feature importance heuristics, establishing the baseline ensemble methods and inconsistent attribution metrics addressed in this work.
Cover for Consistent Individualized Feature Attribution for Tree Ensembles

Abstract

Interpreting predictions from tree ensemble methods such as gradient boosting machines and random forests is important, yet feature attribution for trees is often heuristic and not individualized for each prediction. Here we show that popular feature attribution methods are inconsistent, meaning they can lower a feature's assigned importance when the true impact of that feature actually increases. This is a fundamental problem that casts doubt on any comparison between features. To address it we turn to recent applications of game theory and develop fast exact tree solutions for SHAP (SHapley Additive exPlanation) values, which are the unique consistent and locally accurate attribution values. We then extend SHAP values to interaction effects and define SHAP interaction values. We propose a rich visualization of individualized feature attributions that improves over classic attribution summaries and partial dependence plots, and a unique "supervised" clustering (clustering based on feature attributions). We demonstrate better agreement with human intuition through a user study, exponential improvements in run time, improved clustering performance, and better identification of influential features. An implementation of our algorithm has also been merged into XGBoost and LightGBM, see this http URL for details.

Table of Contents

  • 1 Introduction
  • 2 Inconsistencies in current feature attribution methods
  • 2.1 SHAP values as the only consistent and locally accurate individualized feature attributions
  • 3 Tree SHAP: Fast SHAP value computation for trees
  • 3.1 Estimating SHAP values directly in O⁡(T​L​2M)O(TL2^{M}) time
  • 3.2 Estimating SHAP values in O⁡(T​L​D2)O(TLD^{2}) time
  • 4 SHAP Interaction Values
  • 5 Experiments and Applications
  • 5.1 Agreement with Human Intuition
  • 5.2 Computational Performance
  • 5.3 Supervised Clustering
  • 5.4 Identification of Influential Features
  • 5.5 SHAP Plots
  • 5.5.1 SHAP Summary Plots
  • 5.5.2 SHAP Dependence Plots
  • 6 Conclusion
  • References

Knowls

  1. Knowl 1 — Tree SHAP Algorithm for Fast Exact Tree Ensemble Feature Attribution

    algorithm

    Tree SHAP computes exact SHapley Additive exPlanation (SHAP) values for a decision tree or an ensemble of decision trees in O(TLD2)O(T L D^2) time and O(D2+M)O(D^2 + M) memory, where TT is the number of trees, LL is the maximum number of leaves in any tree, DD is the maximum tree depth (D=log⁡LD = \log L for balanced trees), and MM is the number of input features. For an ensemble, tree attributions sum linearly.

    The algorithm operates by recursively maintaining subset weights during tree traversal rather than computing exponential subset evaluations. A path mm records four attributes for each split: feature index dd, zero-path proportion zz, one-path proportion oo, and cardinality weights ww.

    Input: Input instance xx, decision tree tree={v,a,b,t,r,d}tree = \{v, a, b, t, r, d\}
    Output: Vector ϕ∈RM\phi \in \mathbb{R}^M of feature attribution values
    procedure TreeSHAP(xx, treetree)
        ϕ←\phi \leftarrow array of MM zeros
        
        procedure Extend(mm, pzp_z, pop_o, pip_i)
            l←len(m)l \leftarrow \text{len}(m)
            m′←copy(m)m' \leftarrow \text{copy}(m)
            ml+1′.(d,z,o,w)←(pi,pz,po,if l=0 then 1 else 0)m'_{l+1}.(d, z, o, w) \leftarrow (p_i, p_z, p_o, \text{if } l=0 \text{ then } 1 \text{ else } 0)
            for i←l−1i \leftarrow l - 1 down to 11 do
                mi+1′.w←mi+1′.w+po⋅mi′.w⋅(i/l)m'_{i+1}.w \leftarrow m'_{i+1}.w + p_o \cdot m'_i.w \cdot (i / l)
                mi′.w←pz⋅mi′.w⋅((l−i)/l)m'_i.w \leftarrow p_z \cdot m'_i.w \cdot ((l - i) / l)
            end for
            return m′m'
        end procedure
        procedure Unwind(mm, ii)
            l←len(m)l \leftarrow \text{len}(m)
            n←ml.wn \leftarrow m_l.w
            m′←copy(m1…l−1)m' \leftarrow \text{copy}(m_{1 \dots l-1})
            for j←l−1j \leftarrow l - 1 down to 11 do
                if mi.o≠0m_i.o \neq 0 then
                    t←mj′.wt \leftarrow m'_j.w
                    mj′.w←n⋅l/(j⋅mi.o)m'_j.w \leftarrow n \cdot l / (j \cdot m_i.o)
                    n←t−mj′.w⋅mi.z⋅((l−j)/l)n \leftarrow t - m'_j.w \cdot m_i.z \cdot ((l - j) / l)
                else
                    mj′.w←(mj′.w⋅l)/(mi.z⋅(l−j))m'_j.w \leftarrow (m'_j.w \cdot l) / (m_i.z \cdot (l - j))
                end if
            end for
            for j←ij \leftarrow i to l−1l - 1 do
                mj′.(d,z,o)←mj+1′.(d,z,o)m'_j.(d, z, o) \leftarrow m'_{j+1}.(d, z, o)
            end for
            return m′m'
        end procedure
        procedure Recurse(jj, mm, pzp_z, pop_o, pip_i)
            m←Extend(m,pz,po,pi)m \leftarrow \text{Extend}(m, p_z, p_o, p_i)
            if vj≠internalv_j \neq \text{internal} then
                for i←2i \leftarrow 2 to len(m)\text{len}(m) do
                    w←∑Unwind(m,i).ww \leftarrow \sum \text{Unwind}(m, i).w
                    ϕmi.d←ϕmi.d+w⋅(mi.o−mi.z)⋅vj\phi_{m_i.d} \leftarrow \phi_{m_i.d} + w \cdot (m_i.o - m_i.z) \cdot v_j
                end for
            else
                if xdj≤tjx_{d_j} \le t_j then
                    h←ajh \leftarrow a_j
                    c←bjc \leftarrow b_j
                else
                    h←bjh \leftarrow b_j
                    c←ajc \leftarrow a_j
                end if
                iz←1i_z \leftarrow 1
                io←1i_o \leftarrow 1
                k←FindFirst(m.d,dj)k \leftarrow \text{FindFirst}(m.d, d_j)
                if k≠nothingk \neq \text{nothing} then
                    iz←mk.zi_z \leftarrow m_k.z
                    io←mk.oi_o \leftarrow m_k.o
                    m←Unwind(m,k)m \leftarrow \text{Unwind}(m, k)
                end if
                Recurse(hh, mm, iz⋅rh/rji_z \cdot r_h / r_j, ioi_o, djd_j)
                Recurse(cc, mm, iz⋅rc/rji_z \cdot r_c / r_j, 00, djd_j)
            end if
        end procedure
        Recurse(11, [ ][\,], 11, 11, 00)
        return ϕ\phi
    end procedure
  2. Knowl 2 — SHAP Interaction Values for Tree Models

    model/method

    Pairwise interaction effects for individualized predictions are captured using SHAP interaction values based on the Shapley interaction index. For a model prediction f(x)f(x) with feature set NN of size MM, the interaction attribution between distinct features ii and jj (i≠ji \neq j) is defined as:

    Φi,j=∑S⊆N∖{i,j}∣S∣!(M−∣S∣−2)!2(M−1)!∇ij(S)\Phi_{i,j} = \sum_{S \subseteq N \setminus \{i,j\}} \frac{|S|!(M - |S| - 2)!}{2(M - 1)!} \nabla_{ij}(S)

    where the discrete second-order difference is:

    ∇ij(S)=fx(S∪{i,j})−fx(S∪{i})−fx(S∪{j})+fx(S)\nabla_{ij}(S) = f_x(S \cup \{i,j\}) - f_x(S \cup \{i\}) - f_x(S \cup \{j\}) + f_x(S)

    Here fx(S)=E[f(x)∣xS]f_x(S) = \mathbb{E}[f(x) \mid x_S]. The interaction effect is split equally so that Φi,j=Φj,i\Phi_{i,j} = \Phi_{j,i}, and the total interaction between features ii and jj is Φi,j+Φj,i\Phi_{i,j} + \Phi_{j,i}. The main effect Φi,i\Phi_{i,i} represents the impact of feature ii with all pairwise interactions removed:

    Φi,i=ϕi−∑j≠iΦi,j\Phi_{i,i} = \phi_i - \sum_{j \neq i} \Phi_{i,j}

    Because ∇ij(S)=[fx(S∪{i,j})−fx(S∪{j})]−[fx(S∪{i})−fx(S)]\nabla_{ij}(S) = [f_x(S \cup \{i, j\}) - f_x(S \cup \{j\})] - [f_x(S \cup \{i\}) - f_x(S)], Φi,j\Phi_{i,j} is the difference between the SHAP value of feature ii computed with feature jj present and with feature jj absent. Using Tree SHAP twice per feature pair computes the full interaction matrix Φ∈RM×M\Phi \in \mathbb{R}^{M \times M} in O(TMLD2)O(T M L D^2) time.

  3. Knowl 3 — Inconsistency of Heuristic Tree Feature Attribution Methods

    theoretical result

    A feature attribution method is defined as consistent if, whenever a model is modified such that a feature's true marginal impact on the model's output increases or stays the same, that feature's assigned attribution value never decreases.

    Three widely used tree ensemble feature attribution methods fail consistency:

    1. Gain (Gini/Impurity decrease): Computes the total loss or impurity reduction across all splits on a given feature. Inconsistent because a feature with higher impact can receive lower attribution depending on the greedy split hierarchy chosen by the tree.
    2. Split Count: Counts the frequency with which a feature is split upon. Inconsistent because a critical feature may require only a single split near the root while less important features split repeatedly in deeper subtrees.
    3. Saabas Method: Computes individualized attributions by taking the difference between expected node values along the single decision path taken by an input. Inconsistent because attribution depends entirely on the split order along that single path rather than averaging over all possible feature subsets.

    Only SHAP values and permutation-based methods guarantee consistency.

  4. Knowl 4 — Axiomatic Characterization of SHAP Values for Individualized Tree Explanations

    definition

    Additive feature attribution methods explain a model prediction f(x)f(x) using a linear explanation model gg over binary coalition vectors z′∈{0,1}Mz' \in \{0, 1\}^M:

    g(z′)=ϕ0+∑i=1Mϕizi′g(z') = \phi_0 + \sum_{i=1}^M \phi_i z'_i

    where MM is the number of input features, zi′=1z'_i = 1 indicates feature ii is observed, zi′=0z'_i = 0 indicates feature ii is missing, and ϕi∈R\phi_i \in \mathbb{R} is the feature attribution.

    SHAP values are the unique attribution values satisfying three properties:

    1. Local Accuracy: ∑i=0Mϕi=f(x)\sum_{i=0}^M \phi_i = f(x), where ϕ0=E[f(x)]\phi_0 = \mathbb{E}[f(x)].
    2. Missingness: If zi′=0z'_i = 0 in the explanation domain, then ϕi=0\phi_i = 0.
    3. Consistency: If a model ff changes to f′f' such that for all feature subsets S⊆N∖{i}S \subseteq N \setminus \{i\}, fx′(S∪{i})−fx′(S)≥fx(S∪{i})−fx(S)f'_x(S \cup \{i\}) - f'_x(S) \ge f_x(S \cup \{i\}) - f_x(S), then ϕi(f′,x)≥ϕi(f,x)\phi_i(f', x) \ge \phi_i(f, x).

    Using conditional expectations to model missingness (fx(S)=E[f(x)∣xS]f_x(S) = \mathbb{E}[f(x) \mid x_S]), the unique consistent attribution values are:

    ϕi=∑S⊆N∖{i}∣S∣!(M−∣S∣−1)!M![E[f(x)∣xS∪{i}]−E[f(x)∣xS]]\phi_i = \sum_{S \subseteq N \setminus \{i\}} \frac{|S|!(M - |S| - 1)!}{M!} [\mathbb{E}[f(x) \mid x_{S \cup \{i\}}] - \mathbb{E}[f(x) \mid x_S]]

  5. Knowl 5 — Supervised Clustering Using Individualized SHAP Attributions

    model/method

    Supervised clustering applies clustering algorithms (such as hierarchical agglomerative clustering) to vectors of individualized feature attribution values ϕ(x)=[ϕ1(x),…,ϕM(x)]∈RM\phi(x) = [\phi_1(x), \dots, \phi_M(x)] \in \mathbb{R}^M rather than raw feature vectors x∈RMx \in \mathbb{R}^M.

    Because every ϕi(x)\phi_i(x) is denominated in the units of the model prediction f(x)f(x) (such as log-odds or risk):

    1. Attributions across heterogeneous feature types (e.g., dollars, age, biometric metrics) share a common scale, eliminating the need to define arbitrary feature scaling or distance metrics.
    2. Feature variations that have no impact on the model's prediction do not influence the distance between instances.
    3. Instances are clustered according to shared predictive explanations, grouping samples that arrive at similar predictions for the same underlying reasons.
  6. Knowl 6 — Recursive Conditional Expectation Estimation in Decision Trees

    algorithm

    For a single decision tree tree={v,a,b,t,r,d}tree = \{v, a, b, t, r, d\} (where vv is node values, aa and bb are left and right child indices, tt is thresholds, rr is node covers/sample counts, and dd is split feature indices), the conditional expectation E[f(x)∣xS]\mathbb{E}[f(x) \mid x_S] for a subset S⊆{1,…,M}S \subseteq \{1, \dots, M\} is computed recursively in O(L)O(L) time, where LL is the number of leaves.

    Input: Input instance xx, feature subset SS, decision tree tree={v,a,b,t,r,d}tree = \{v, a, b, t, r, d\}
    Output: Expected value E[f(x)∣xS]\mathbb{E}[f(x) \mid x_S]
    procedure ExpValue(xx, SS, treetree)
        procedure G(jj, ww)
            if vj≠internalv_j \neq \text{internal} then
                return w⋅vjw \cdot v_j
            else
                if dj∈Sd_j \in S then
                    if xdj≤tjx_{d_j} \le t_j then
                        return G(aja_j, ww)
                    else
                        return G(bjb_j, ww)
                    end if
                else
                    return G(aja_j, w⋅raj/rjw \cdot r_{a_j} / r_j) + G(bjb_j, w⋅rbj/rjw \cdot r_{b_j} / r_j)
                end if
            end if
        end procedure
        return G(1, 1)
    end procedure

    Evaluating this for all 2M2^M subsets SS yields an exact SHAP computation baseline in O(TL2M)O(T L 2^M) time for an ensemble of TT trees.

  7. Knowl 7 — SHAP Summary Plots and SHAP Dependence Plots

    model/method

    Two visualization formats utilize individualized SHAP attributions across datasets:

    1. SHAP Summary Plot: Replaces standard global feature importance bar charts. Features are ordered vertically by global mean absolute attribution ∑j=1N∣ϕi(j)∣\sum_{j=1}^N |\phi_i^{(j)}|. For each feature, each dataset sample is displayed as a point plotted horizontally at its attribution value ϕi(j)\phi_i^{(j)} and colored by its feature value (e.g., low blue to high red). Points pile up vertically to indicate density, visualizing value distribution, impact sign, magnitude, and asymmetry.

    2. SHAP Dependence Plot: Replaces 1D partial dependence plots. For a chosen feature ii, each sample is plotted with its feature value xix_i on the x-axis and its attribution ϕi\phi_i on the y-axis. Vertical dispersion at a single feature value indicates interaction effects, which are illuminated by coloring points according to the value of an interacting feature xkx_k. Plotting the main effect Φi,i\Phi_{i,i} eliminates vertical dispersion, while plotting Φi,k\Phi_{i,k} isolates the interaction component.

  8. Knowl 8 — Influential Feature Identification via Feature Perturbation

    empirical result

    Attribution methods were evaluated on an XGBoost sentiment analysis model (50 trees, maximum depth 30) trained on 11,712 airline tweets with 1,686 bag-of-words features, tested across 2,928 samples. For each method, the most influential negative feature for each sample was identified and replaced with that feature's value from a random training tweet (mimicking missingness), after which the updated prediction was computed.

    Tracking the cumulative increase in model prediction across all 2,928 test tweets demonstrated that SHAP values achieved the highest total increase in model output (approaching ∼550\sim 550 total increase), followed by Saabas (∼200\sim 200), Gain, Permutation, and Split Count (sim0\\sim 0). This confirmed that SHAP values most accurately identify the features exerting the strongest local influence on model predictions.

  9. Knowl 9 — Supervised Clustering Variance Explanation on Census and Gene Expression Datasets

    empirical result

    The fidelity of individualized attributions was quantified by measuring the proportion of model output variance explained (R2R^2) during agglomerative hierarchical clustering as sample clusters were iteratively merged from NN clusters (R2=1.0R^2 = 1.0) to 1 cluster (R2=0.0R^2 = 0.0).

    1. UCI Census Dataset: Using an XGBoost model (500 trees, max depth 6, shrinkage η=0.005\eta = 0.005) predicting annual income ≥$50K\ge \$50\text{K}, SHAP-based clustering achieved an Area Under the R2R^2 Curve (AUC) of 0.98 compared to 0.97 for Saabas attributions.
    2. Alzheimer's Cognitive Score Dataset: Using an XGBoost model predicting CERAD cognitive scores from 200 gene expression module levels, SHAP-based clustering achieved an AUC of 0.96 compared to 0.88 for Saabas attributions, maintaining higher R2R^2 as the number of clusters was reduced.
  10. Knowl 10 — Human Credit Allocation Intuition User Study

    empirical result

    A user study with 34 Amazon Mechanical Turk participants evaluated intuitive credit assignment on a two-feature disease risk decision tree (Model A: disease risk equals 20 at baseline, and 80 when both Fever = Yes and Cough = Yes). Participants allocated the 60-point increase in risk between Fever and Cough without constraints.

    Results:

    • Equal split (30 to Fever, 30 to Cough), matching SHAP value attribution, was chosen by the largest group of participants (12 out of 34).
    • Allocation giving greater weight to the root feature Fever (e.g., 40/20) was the second most common.
    • Allocation matching the Saabas heuristic (which splits credit along the decision path, allocating 20 to Fever and 40 to Cough) had fewer adherents than the equal split.

Coverage note — None was omitted; all core theoretical definitions, algorithms (Tree SHAP and conditional expectation calculation), interaction formulations, visualizations, and experimental results were included.

References

  1. 1.Lidia Auret and Chris Aldrich. 2011. Empirical comparison of tree ensemble variable importance measures. Chemometrics and Intelligent Laboratory Systems 105, 2 (2011), 157–170.
  2. 2.David Baehrens, Timon Schroeter, Stefan Harmeling, Motoaki Kawanabe, Katja Hansen, and Klaus-Robert Müller. 2010. How to explain individual classification decisions. Journal of Machine Learning Research 11, Jun (2010), 1803–1831.
  3. 3.Leo Breiman, Jerome Friedman, Charles J Stone, and Richard A Olshen. 1984. Classification and regression trees. CRC press.
  4. 4.Safiye Celik, Benjamin Logsdon, and Su-In Lee. 2014. Efficient dimensionality reduction for high-dimensional network estimation. In International Conference on Machine Learning. 1953–1961.
  5. 5.S Chebrolu, A Abraham, and J Thomas. 2005. Feature deduction and ensemble design of intrusion detection systems. Computers & security 24, 4 (2005), 295–307.
  6. 6.Tianqi Chen and Carlos Guestrin. 2016. XGBoost: A scalable tree boosting system. In Proceedings of the 22Nd ACM SIGKDD International Conference on Knowledge Discovery and Data Mining. ACM, 785–794.
  7. 7.Christine S Cox, Jacob J Feldman, Cordell D Golden, Madelyn A Lane, Jennifer H Madans, Michael E Mussolino, and Sandra T Rothwell. 1997. Plan and operation of the NHANES I Epidemiologic Followup Study, 1992. (1997).
  8. 8.Crowdflower. 2015. Twitter US Airline Sentiment. https://www.kaggle.com/crowdflower/twitter-airline-sentiment. (2015). Accessed: 2018-02-06.
  9. 9.Anupam Datta, Shayak Sen, and Yair Zick. 2016. Algorithmic transparency via quantitative input influence: Theory and experiments with learning systems. In Security and Privacy (SP), 2016 IEEE Symposium on. IEEE, 598–617.
  10. 10.R Díaz-Uriarte and S De Andres. 2006. Gene selection and classification of microarray data using random forest. BMC bioinformatics 7, 1 (2006), 3.
  11. 11.Jerome Friedman, Trevor Hastie, and Robert Tibshirani. 2001. The elements of statistical learning. Vol. 1. Springer series in statistics Springer, Berlin.
  12. 12.Katsushige Fujimoto, Ivan Kojadinovic, and Jean-Luc Marichal. 2006. Axiomatic characterizations of probabilistic and cardinal-probabilistic interaction indices. Games and Economic Behavior 55, 1 (2006), 72–99.
  13. 13.A Irrthum, L Wehenkel, P Geurts, et al. 2010. Inferring regulatory networks from expression data using tree-based methods. PloS one 5, 9 (2010), e12776.
  14. 14.Hemant Ishwaran et al. 2007. Variable importance in binary regression trees and forests. Electronic Journal of Statistics 1 (2007), 519–537.
  15. 15.M. Lichman. 2013. UCI ML Repository. (2013). http://archive.ics.uci.edu/ml
  16. 16.Scott M Lundberg and Su-In Lee. 2017. A Unified Approach to Interpreting Model Predictions. In Advances in Neural Information Processing Systems 30. Curran Associates, Inc., 4768–4777. http://papers.nips.cc/paper/7062-a-unified-approach-to-interpreting-model-predictions.pdf
  17. 17.Scott M Lundberg, Bala Nair, Monica S Vavilala, Mayumi Horibe, Michael J Eisses, Trevor Adams, David E Liston, Daniel King-Wai Low, Shu-Fang Newman, Jerry Kim, et al. 2017. Explainable machine learning predictions to help anesthesiologists prevent hypoxemia during surgery. bioRxiv (2017), 206540.
  18. 18.Henry W Miller. 1973. Plan and operation of the health and nutrition examination survey, United States, 1971-1973. DHEW publication no.(PHS)-Dept. of Health, Education, and Welfare (USA) (1973).
  19. 19.Suzanne S Mirra, A Heyman, D McKeel, SM Sumi, Barbara J Crain, LM Brownlee, FS Vogel, JP Hughes, G Van Belle, L Berg, et al. 1991. The Consortium to Establish a Registry for Alzheimer's Disease (CERAD) Part II. Stand. of the neuropathologic assessmen== of Alzheimer's disease. Neurology 41, 4 (1991), 479–479.
  20. 20.F Pedregosa, G Varoquaux, A Gramfort, V Michel, B Thirion, O Grisel, M Blondel, P Prettenhofer, R Weiss, V Dubourg, et al. 2011. Scikit-learn: Machine learning in Python. JMLR 12, Oct (2011), 2825–2830.
  21. 21.Marco Tulio Ribeiro, Sameer Singh, and Carlos Guestrin. 2016. Why should i trust you?: Explaining the predictions of any classifier. In Proceedings of the 22nd ACM SIGKDD. ACM, 1135–1144.
  22. 22.Greg Ridgeway. 2010. Generalized boosted regression models. Documentation on the R Package 'gbm', version 1.6–3. (2010).
  23. 23.W Rodenburg, G Heidema, J Boer, I Bovee-Oudenhoven, E Feskens, E Mariman, and J Keijer. 2008. A framework to identify physiological responses in microarray-based gene expression studies: selection and interpretation of biologically relevant genes. Physiological genomics 33, 1 (2008), 78–90.
  24. 24.Ando Saabas. 2014. Interpreting random forests. http://blog.datadive.net/interpreting-random-forests/. (2014). Accessed: 2017-06-15.
  25. 25.Marco Sandri and Paola Zuccolotto. 2008. A bias correction algorithm for the Gini variable importance measure in classification trees. Journal of Computational and Graphical Statistics 17, 3 (2008), 611–628.
  26. 26.C Strobl, A Boulesteix, T Kneib, T Augustin, and A Zeileis. 2008. Conditional variable importance for random forests. BMC bioinformatics 9, 1 (2008), 307.
  27. 27.Erik Štrumbelj and Igor Kononenko. 2014. Explaining prediction models and individual predictions with feature contributions. Knowledge and information systems 41, 3 (2014), 647–665.

Citation

MLA
Lundberg, S. M., et al. “Consistent Individualized Feature Attribution for Tree Ensembles”. arXiv, 2018, http://arxiv.org/abs/1802.03888v3.
APA
Lundberg, S. M., Erion, G. G., & Lee, S.-I. (2018). Consistent Individualized Feature Attribution for Tree Ensembles. arXiv. http://arxiv.org/abs/1802.03888v3
Chicago
Lundberg, S. M., G. G. Erion, and S.-I. Lee. 2018. “Consistent Individualized Feature Attribution for Tree Ensembles”. arXiv. http://arxiv.org/abs/1802.03888v3.
Harvard
Lundberg, S.M., Erion, G.G. and Lee, S.-I. (2018) “Consistent Individualized Feature Attribution for Tree Ensembles”, arXiv [Preprint]. Available at: http://arxiv.org/abs/1802.03888v3.
Vancouver
1. Lundberg SM, Erion GG, Lee S-I (2018) Consistent Individualized Feature Attribution for Tree Ensembles. arXiv

BibTeX

@article{lundberg2018consistent,
  title = {Consistent Individualized Feature Attribution for Tree Ensembles},
  author = {Lundberg, Scott M. and Erion, Gabriel G. and Lee, Su-In},
  year = {2018},
  journal = {arXiv},
  url = {http://arxiv.org/abs/1802.03888v3},
  eprint = {1802.03888}
}
Metadata:arXiv

Source Code

This paper has an official code repository available. Click below to access the source code.

View Repository

Access the Paper

This paper is available from its original source. Click below to access the PDF.

Open PDF

License: Authors