lgb.importance.R 2.09 KB
Newer Older
1
#' Compute feature importance in a model
2
#'
3
#' Creates a \code{data.table} of feature importances in a model.
4
#'
5
6
#' @param model object of class \code{lgb.Booster}.
#' @param percentage whether to show importance in relative percentage.
7
#'
8
#' @return
9
#'
10
11
12
13
14
15
16
#' For a tree model, a \code{data.table} with the following columns:
#' \itemize{
#'   \item \code{Feature} Feature names in the model.
#'   \item \code{Gain} The total gain of this feature's splits.
#'   \item \code{Cover} The number of observation related to this feature.
#'   \item \code{Frequency} The number of times a feature splited in trees.
#' }
17
#'
18
#' @examples
19
20
#' library(lightgbm)
#' data(agaricus.train, package = "lightgbm")
21
22
23
#' train <- agaricus.train
#' dtrain <- lgb.Dataset(train$data, label = train$label)
#'
24
#' params <- list(objective = "binary",
25
26
27
#'                learning_rate = 0.01, num_leaves = 63, max_depth = -1,
#'                min_data_in_leaf = 1, min_sum_hessian_in_leaf = 1)
#' model <- lgb.train(params, dtrain, 10)
28
29
30
#'
#' tree_imp1 <- lgb.importance(model, percentage = TRUE)
#' tree_imp2 <- lgb.importance(model, percentage = FALSE)
31
#'
32
#' @importFrom magrittr %>% %T>% extract
33
34
35
#' @importFrom data.table :=
#' @export
lgb.importance <- function(model, percentage = TRUE) {
36

37
  # Check if model is a lightgbm model
38
  if (!inherits(model, "lgb.Booster")) {
39
40
    stop("'model' has to be an object of class lgb.Booster")
  }
41

42
  # Setup importance
43
  tree_dt <- lgb.model.dt.tree(model)
44

45
  # Extract elements
46
47
  tree_imp <- tree_dt %>%
    magrittr::extract(.,
48
                      i = ! is.na(split_index),
49
50
51
52
                      j = .(Gain = sum(split_gain), Cover = sum(internal_count), Frequency = .N),
                      by = "split_feature") %T>%
    data.table::setnames(., old = "split_feature", new = "Feature") %>%
    magrittr::extract(., i = order(Gain, decreasing = TRUE))
53

54
  # Check if relative values are requested
55
56
57
58
59
  if (percentage) {
    tree_imp[, ":="(Gain = Gain / sum(Gain),
                    Cover = Cover / sum(Cover),
                    Frequency = Frequency / sum(Frequency))]
  }
60

61
  # Return importance table
62
  return(tree_imp)
63

64
}