############################### Feature Interaction Constraints ############################### The decision tree is a powerful tool to discover interaction among independent variables (features). Variables that appear together in a traversal path are interacting with one another, since the condition of a child node is predicated on the condition of the parent node. For example, the highlighted red path in the diagram below contains three variables: :math:`x_1`, :math:`x_7`, and :math:`x_{10}`, so the highlighted prediction (at the highlighted leaf node) is the product of interaction between :math:`x_1`, :math:`x_7`, and :math:`x_{10}`. .. plot:: :nofigs: from graphviz import Source source = r""" digraph feature_interaction_illustration1 { graph [fontname = "helvetica"]; node [fontname = "helvetica"]; edge [fontname = "helvetica"]; 0 [label=10 < -1.5 ?>, shape=box, color=red, fontcolor=red]; 1 [label=2 < 2 ?>, shape=box]; 2 [label=7 < 0.3 ?>, shape=box, color=red, fontcolor=red]; 3 [label="...", shape=none]; 4 [label="...", shape=none]; 5 [label=1 < 0.5 ?>, shape=box, color=red, fontcolor=red]; 6 [label="...", shape=none]; 7 [label="...", shape=none]; 8 [label="Predict +1.3", color=red, fontcolor=red]; 0 -> 1 [labeldistance=2.0, labelangle=45, headlabel="Yes/Missing "]; 0 -> 2 [labeldistance=2.0, labelangle=-45, headlabel="No", color=red, fontcolor=red]; 1 -> 3 [labeldistance=2.0, labelangle=45, headlabel="Yes"]; 1 -> 4 [labeldistance=2.0, labelangle=-45, headlabel=" No/Missing"]; 2 -> 5 [labeldistance=2.0, labelangle=-45, headlabel="Yes", color=red, fontcolor=red]; 2 -> 6 [labeldistance=2.0, labelangle=-45, headlabel=" No/Missing"]; 5 -> 7; 5 -> 8 [color=red]; } """ Source(source, format='png').render('../_static/feature_interaction_illustration1', view=False) Source(source, format='svg').render('../_static/feature_interaction_illustration1', view=False) .. figure:: ../_static/feature_interaction_illustration1.svg :align: center :figwidth: 80 % When the tree depth is larger than one, many variables interact on the sole basis of minimizing training loss, and the resulting decision tree may capture a spurious relationship (noise) rather than a legitimate relationship that generalizes across different datasets. **Feature interaction constraints** allow users to decide which variables are allowed to interact and which are not. Potential benefits include: * Better predictive performance from focusing on interactions that work -- whether through domain specific knowledge or algorithms that rank interactions * Less noise in predictions; better generalization * More control to the user on what the model can fit. For example, the user may want to exclude some interactions even if they perform well due to regulatory constraints. **************** A Simple Example **************** Feature interaction constraints are expressed in terms of groups of variables that are allowed to interact. For example, the constraint ``[0, 1]`` indicates that variables :math:`x_0` and :math:`x_1` are allowed to interact with each other but with no other variable. Similarly, ``[2, 3, 4]`` indicates that :math:`x_2`, :math:`x_3`, and :math:`x_4` are allowed to interact with one another but with no other variable. A set of feature interaction constraints is expressed as a nested list, e.g. ``[[0, 1], [2, 3, 4]]``, where each inner list is a group of indices of features that are allowed to interact with each other. In the following diagram, the left decision tree is in violation of the first constraint (``[0, 1]``), whereas the right decision tree complies with both the first and second constraints (``[0, 1]``, ``[2, 3, 4]``). .. plot:: :nofigs: from graphviz import Source source = r""" digraph feature_interaction_illustration2 { graph [fontname = "helvetica"]; node [fontname = "helvetica"]; edge [fontname = "helvetica"]; 0 [label=0 < 5.0 ?>, shape=box]; 1 [label=2 < -3.0 ?>, shape=box]; 2 [label="+0.6"]; 3 [label="-0.4"]; 4 [label="+1.2"]; 0 -> 1 [labeldistance=2.0, labelangle=45, headlabel="Yes/Missing "]; 0 -> 2 [labeldistance=2.0, labelangle=-45, headlabel="No"]; 1 -> 3 [labeldistance=2.0, labelangle=45, headlabel="Yes"]; 1 -> 4 [labeldistance=2.0, labelangle=-45, headlabel=" No/Missing"]; } """ Source(source, format='png').render('../_static/feature_interaction_illustration2', view=False) Source(source, format='svg').render('../_static/feature_interaction_illustration2', view=False) .. plot:: :nofigs: from graphviz import Source source = r""" digraph feature_interaction_illustration3 { graph [fontname = "helvetica"]; node [fontname = "helvetica"]; edge [fontname = "helvetica"]; 0 [label=3 < 2.5 ?>, shape=box]; 1 [label="+1.6"]; 2 [label=2 < -1.2 ?>, shape=box]; 3 [label="+0.1"]; 4 [label="-0.3"]; 0 -> 1 [labeldistance=2.0, labelangle=45, headlabel="Yes"]; 0 -> 2 [labeldistance=2.0, labelangle=-45, headlabel=" No/Missing"]; 2 -> 3 [labeldistance=2.0, labelangle=45, headlabel="Yes/Missing "]; 2 -> 4 [labeldistance=2.0, labelangle=-45, headlabel="No"]; } """ Source(source, format='png').render('../_static/feature_interaction_illustration3', view=False) Source(source, format='svg').render('../_static/feature_interaction_illustration3', view=False) .. |fig1| image:: ../_static/feature_interaction_illustration2.svg :scale: 7% :align: middle .. |fig2| image:: ../_static/feature_interaction_illustration3.svg :scale: 7% :align: middle +-----------+---------+ | |fig1| | |fig2| | +-----------+---------+ | forbidden | allowed | +-----------+---------+ **************************************************** Enforcing Feature Interaction Constraints in XGBoost **************************************************** It is very simple to enforce feature interaction constraints in XGBoost. Here we will give an example using Python, but the same general idea generalizes to other platforms. Suppose the following code fits your model without feature interaction constraints: .. code-block:: python model_no_constraints = xgb.train(params, dtrain, num_boost_round = 1000, evals = evallist, early_stopping_rounds = 10) Then fitting with feature interaction constraints only requires adding a single parameter: .. code-block:: python params_constrained = params.copy() # Use nested list to define feature interaction constraints params_constrained['interaction_constraints'] = '[[0, 2], [1, 3, 4], [5, 6]]' # Features 0 and 2 are allowed to interact with each other but with no other feature # Features 1, 3, 4 are allowed to interact with one another but with no other feature # Features 5 and 6 are allowed to interact with each other but with no other feature model_with_constraints = xgb.train(params_constrained, dtrain, num_boost_round = 1000, evals = evallist, early_stopping_rounds = 10) ************************** Using feature name instead ************************** XGBoost's Python package supports using feature names instead of feature index for specifying the constraints. Given a data frame with columns ``["f0", "f1", "f2"]``, the feature interaction constraint can be specified as ``[["f0", "f2"]]``. ************** Advanced topic ************** The intuition behind interaction constraints is simple. Users may have prior knowledge about relations between different features, and encode it as constraints during model construction. But there are also some subtleties around specifying constraints. Take the constraint ``[[1, 2], [2, 3, 4]]`` as an example. The second feature appears in two different interaction sets, ``[1, 2]`` and ``[2, 3, 4]``. So the union set of features allowed to interact with ``2`` is ``{1, 3, 4}``. In the following diagram, the root splits at feature ``2``. Because all its descendants should be able to interact with it, all 4 features are legitimate split candidates at the second layer. At first sight, this might look like disregarding the specified constraint sets, but it is not. .. plot:: :nofigs: from graphviz import Source source = r""" digraph feature_interaction_illustration4 { graph [fontname = "helvetica"]; node [fontname = "helvetica"]; edge [fontname = "helvetica"]; 0 [label=2>, shape=box, color=black, fontcolor=black]; 1 [label={1, 2, 3, 4}>, shape=box]; 2 [label={1, 2, 3, 4}>, shape=box, color=black, fontcolor=black]; 3 [label="...", shape=none]; 4 [label="...", shape=none]; 5 [label="...", shape=none]; 6 [label="...", shape=none]; 0 -> 1; 0 -> 2; 1 -> 3; 1 -> 4; 2 -> 5; 2 -> 6; } """ Source(source, format='png').render('../_static/feature_interaction_illustration4', view=False) Source(source, format='svg').render('../_static/feature_interaction_illustration5', view=False) .. figure:: ../_static/feature_interaction_illustration4.png :align: center :figwidth: 80 % ``{1, 2, 3, 4}`` represents the sets of legitimate split features. This has lead to some interesting implications of feature interaction constraints. Take ``[[0, 1], [0, 1, 2], [1, 2]]`` as another example. Assuming we have only 3 available features in our training datasets for presentation purpose, careful readers might have found out that the above constraint is the same as simply ``[[0, 1, 2]]``. Since no matter which feature is chosen for split in the root node, all its descendants are allowd to include every feature as legitimate split candidates without violating interaction constraints. For one last example, we use ``[[0, 1], [1, 3, 4]]`` and choose feature ``0`` as split for the root node. At the second layer of the built tree, ``1`` is the only legitimate split candidate except for ``0`` itself, since they belong to the same constraint set. Following the grow path of our example tree below, the node at the second layer splits at feature ``1``. But due to the fact that ``1`` also belongs to second constraint set ``[1, 3, 4]``, at the third layer, we are allowed to include all features as split candidates and still comply with the interaction constraints of its ascendants. .. plot:: :nofigs: from graphviz import Source source = r""" digraph feature_interaction_illustration5 { graph [fontname = "helvetica"]; node [fontname = "helvetica"]; edge [fontname = "helvetica"]; 0 [label=0>, shape=box, color=black, fontcolor=black]; 1 [label="...", shape=none]; 2 [label=1>, shape=box, color=black, fontcolor=black]; 3 [label={0, 1, 3, 4}>, shape=box, color=black, fontcolor=black]; 4 [label={0, 1, 3, 4}>, shape=box, color=black, fontcolor=black]; 5 [label="...", shape=none]; 6 [label="...", shape=none]; 7 [label="...", shape=none]; 8 [label="...", shape=none]; 0 -> 1; 0 -> 2; 2 -> 3; 2 -> 4; 3 -> 5; 3 -> 6; 4 -> 7; 4 -> 8; } """ Source(source, format='png').render('../_static/feature_interaction_illustration6', view=False) Source(source, format='svg').render('../_static/feature_interaction_illustration7', view=False) .. figure:: ../_static/feature_interaction_illustration6.png :align: center :figwidth: 80 % ``{0, 1, 3, 4}`` represents the sets of legitimate split features.