Density of model residuals.

plot_residual_density(object, ..., split = FALSE, variable = NULL)

plotResidualDensity(object, ..., split = FALSE, variable = NULL)

Arguments

object

An object of class 'auditor_model_residual' created with model_residual function.

...

Other 'auditor_model_residual' objects to be plotted together.

split

Logical. Indicates whenever plot should be splitted by variable.

variable

Variable name o split. Optional.

Value

ggplot object

See also

Examples

dragons <- DALEX::dragons[1:100, ] # fit a model model_lm <- lm(life_length ~ ., data = dragons) # use DALEX package to wrap up a model into explainer exp_lm <- DALEX::explain(model_lm, data = dragons, y = dragons$life_length)
#> Preparation of a new explainer is initiated #> -> model label : lm (default) #> -> data : 100 rows 8 cols #> -> target variable : 100 values #> -> predict function : yhat.lm will be used (default) #> -> predicted values : numerical, min = 585.8311 , mean = 1347.787 , max = 2942.307 #> -> residual function : difference between y and yhat (default) #> -> residuals : numerical, min = -88.41755 , mean = -1.489291e-13 , max = 77.92805 #> A new explainer has been created!
# validate a model with auditor library(auditor) mr_lm <- model_residual(exp_lm) # plot results plot_residual_density(mr_lm)
plot(mr_lm, type = "residual_density")
library(randomForest) model_rf <- randomForest(life_length~., data = dragons) exp_rf <- DALEX::explain(model_rf, data = dragons, y = dragons$life_length)
#> Preparation of a new explainer is initiated #> -> model label : randomForest (default) #> -> data : 100 rows 8 cols #> -> target variable : 100 values #> -> predict function : yhat.randomForest will be used (default) #> -> predicted values : numerical, min = 758.7355 , mean = 1341.421 , max = 2466.698 #> -> residual function : difference between y and yhat (default) #> -> residuals : numerical, min = -187.6278 , mean = 6.365703 , max = 431.7228 #> A new explainer has been created!
mr_rf <- model_residual(exp_rf) plot_residual_density(mr_lm, mr_rf)
plot(mr_lm, mr_rf, type = "residual_density")