Ceteris Paribus Plots

Please, note that there are two methods that create Ceteris Paribus Plots. The older solution is to use what_if() function. The new solution is to use ceteris_paribus() function. These functions have slightly different interfaces.

The old interface (what_if)

Here is the one with what_if().

##      m2.price construction.year surface floor no.rooms    district
## 1001     4644              1976     131     3        5 Srodmiescie
wi_rf <- what_if(explainer_rf, observation = new_apartment)
wi_rf
##       y_hat new_x             vname   x_quant quant relative_quant
## 0% 4255.354  1920 construction.year 0.6268889  0.00     -0.6268889
## 1% 4300.702  1921 construction.year 0.6268889  0.01     -0.6168889
## 2% 4301.926  1922 construction.year 0.6268889  0.02     -0.6068889
## 3% 4305.352  1923 construction.year 0.6268889  0.03     -0.5968889
## 4% 4305.352  1923 construction.year 0.6268889  0.04     -0.5868889
## 5% 4267.723  1924 construction.year 0.6268889  0.05     -0.5768889
##           label
## 0% randomForest
## 1% randomForest
## 2% randomForest
## 3% randomForest
## 4% randomForest
## 5% randomForest

wi_lm <- what_if(explainer_lm, observation = new_apartment)
wi_lm
##           y_hat new_x             vname   x_quant quant relative_quant
## 1001   4832.833  1920 construction.year 0.6268889  0.00     -0.6268889
## 1001.1 4832.604  1921 construction.year 0.6268889  0.01     -0.6168889
## 1001.2 4832.375  1922 construction.year 0.6268889  0.02     -0.6068889
## 1001.3 4832.146  1923 construction.year 0.6268889  0.03     -0.5968889
## 1001.4 4832.146  1923 construction.year 0.6268889  0.04     -0.5868889
## 1001.5 4831.917  1924 construction.year 0.6268889  0.05     -0.5768889
##        label
## 1001      lm
## 1001.1    lm
## 1001.2    lm
## 1001.3    lm
## 1001.4    lm
## 1001.5    lm

The new interface (ceteris_paribus)

Here is the one with ceteris_paribus().

cp_rf <- ceteris_paribus(explainer_rf, observation = new_apartment)
cp_rf
## Top profiles    : 
##        m2.price construction.year surface floor no.rooms    district
## 1001       4644              1920     131     3        5 Srodmiescie
## 1001.1     4644              1921     131     3        5 Srodmiescie
## 1001.2     4644              1922     131     3        5 Srodmiescie
## 1001.3     4644              1923     131     3        5 Srodmiescie
## 1001.4     4644              1923     131     3        5 Srodmiescie
## 1001.5     4644              1924     131     3        5 Srodmiescie
##          _yhat_           _vname_ _ids_      _label_
## 1001   4255.354 construction.year  1001 randomForest
## 1001.1 4300.702 construction.year  1001 randomForest
## 1001.2 4301.926 construction.year  1001 randomForest
## 1001.3 4305.352 construction.year  1001 randomForest
## 1001.4 4305.352 construction.year  1001 randomForest
## 1001.5 4267.723 construction.year  1001 randomForest
## 
## 
## Top observations:
##      m2.price construction.year surface floor no.rooms    district  _yhat_
## 1001     4644              1976     131     3        5 Srodmiescie 4160.84
##           _label_
## 1001 randomForest

cp_lm <- ceteris_paribus(explainer_lm, observation = new_apartment)
cp_lm
## Top profiles    : 
##        m2.price construction.year surface floor no.rooms    district
## 1001       4644              1920     131     3        5 Srodmiescie
## 1001.1     4644              1921     131     3        5 Srodmiescie
## 1001.2     4644              1922     131     3        5 Srodmiescie
## 1001.3     4644              1923     131     3        5 Srodmiescie
## 1001.4     4644              1923     131     3        5 Srodmiescie
## 1001.5     4644              1924     131     3        5 Srodmiescie
##          _yhat_           _vname_ _ids_ _label_
## 1001   4832.833 construction.year  1001      lm
## 1001.1 4832.604 construction.year  1001      lm
## 1001.2 4832.375 construction.year  1001      lm
## 1001.3 4832.146 construction.year  1001      lm
## 1001.4 4832.146 construction.year  1001      lm
## 1001.5 4831.917 construction.year  1001      lm
## 
## 
## Top observations:
##      m2.price construction.year surface floor no.rooms    district
## 1001     4644              1976     131     3        5 Srodmiescie
##        _yhat_ _label_
## 1001 4820.009      lm