Directly test whether off-manifold behavior explains model-dependent PDP distortion
Determine whether the differing refitting gains of tree ensembles and smooth models under strong feature dependence are caused by the models' distinct on-manifold and off-manifold behavior in partial dependence estimation, specifically by comparing partial dependence estimates evaluated on marginally generated feature combinations with the models' behavior on combinations supported by the joint data distribution.
References
We offer the following explanation, and flag that we have not tested it directly. A partial dependence function averages $f$ over the marginal distribution of the remaining features, so at $\rho = 0.85$ it evaluates the model at feature combinations that the joint distribution almost never produces. How badly the resulting curve is distorted therefore depends on how the model behaves away from the data. A tree ensemble is piecewise constant on axis-aligned cells and its predictions off the data manifold are bounded by the leaf values it learned on it, so its partial dependence curves stay close to the right scale. An RBF support vector machine and a multilayer perceptron are smooth global functions that extrapolate, and off-manifold they can return systematically inflated or deflated values, which distorts the amplitude of their curves. A single coefficient per feature is exactly the correction such a distortion needs, which is why refitting buys them so much. Testing this would require comparing on-manifold and off-manifold partial dependence estimates directly, which we did not do.