live.Rmd

---
title: "The live R package"
date: "`r format(Sys.time(), '%d-%m-%Y')`"
output:
  html_document:
    toc: true
    toc_float: true
---


This report aims to present the capabilities of the package `live`.

The document is a part of the paper "Landscape of R packages for eXplainable Machine Learning", S. Maksymiuk, A. Gosiewska, and P. Biecek. 
(https://arxiv.org/abs/2009.13248). It contains a real life use-case with a hand of [titanic_imputed](https://modeloriented.github.io/DALEX/reference/titanic.html) data set described in Section *Example gallery for XAI packages* of the article.

We did our best to show the entire range of the implemented explanations. Please note that the examples may be incomplete. If you think something is missing, feel free to make a pull request at the GitHub repository [MI2DataLab/XAI-tools](https://github.com/MI2DataLab/XAI-tools).

The list of use-cases for all packages included in the article is [here](http://xai-tools.drwhy.ai/).  

```{r setup, include=FALSE}
knitr::opts_chunk$set(echo = TRUE, message = FALSE, warning = FALSE)
```

Load [`titanic_imputed`](https://modeloriented.github.io/DALEX/reference/titanic.html) data set.

```{r}
data(titanic_imputed, package = "DALEX")

head(titanic_imputed)
```


```{r}
library(live)
library(ggplot2)
set.seed(123)
```

Fit a random forest to the titanic imputed data.

```{r}
rf_model <- ranger::ranger(survived~., data = titanic_imputed, classification = TRUE, probability = TRUE)
```





# Predict parts

## LIVE explanation


```{r}
similar <- sample_locally(data = titanic_imputed,
                          explained_instance = titanic_imputed[1, ],
                          explained_var = "survived",
                          size = 500)

pred_fun <- function(X.model, newdata) {
  predict(X.model, newdata)$predictions[,2]
}
similar1 <- add_predictions(to_explain = similar,
                            black_box_model = rf_model,
                            predict_fun = pred_fun)
explanation <- fit_explanation(live_object = similar1,
                             white_box = "regr.lm")
plot(explanation, type = "forest")
```




```{r}
similar <- sample_locally(data = titanic_imputed,
                          explained_instance = titanic_imputed[13, ],
                          explained_var = "survived",
                          size = 500)
pred_fun <- function(X.model, newdata) {
  predict(X.model, newdata)$predictions[,2]
}
similar1 <- add_predictions(to_explain = similar,
                            black_box_model = rf_model,
                            predict_fun = pred_fun)
explanation <- fit_explanation(live_object = similar1,
                             white_box = "regr.lm")

plot(explanation, type = "waterfall") + ylim(0.33, 0.47)
```

# Session info

```{r}
sessionInfo()
```