#' @title Regression BART (Bayesian Additive Regression Trees) Learner
#' @author ck37
#' @name mlr_learners_regr.bart
#'
#' @description
#' Bayesian Additive Regression Trees are similar to gradient boosting algorithms.
#' Calls [dbarts::bart()] from \CRANpkg{dbarts}.
#'
#' @template learner
#' @templateVar id regr.bart
#'
#' @section Custom mlr3 parameters:
#' * Parameter: offset
#' * The parameter is removed, because only `dbarts::bart2` allows an offset during training,
#' and therefore the offset parameter in `dbarts:::predict.bart` is irrelevant for
#' `dbarts::dbart`.
#' * Parameter: nchain, combineChains, combinechains
#' * The parameters are removed as parallelization of multiple models is handled by future.
#'
#' @section Initial parameter values:
#' * `keeptrees` is initialized to `TRUE` because it is required for prediction.
#'
#' @references
#' `r format_bib("sparapani2021nonparametric", "chipman2010bart")`
#'
#' @template seealso_learner
#' @template example
#' @export
LearnerRegrBart = R6Class("LearnerRegrBart",
inherit = LearnerRegr,
public = list(
#' @description
#' Creates a new instance of this [R6][R6::R6Class] class.
initialize = function() {
ps = ps(
ntree = p_int(default = 200L, lower = 1L, tags = "train"),
sigest = p_uty(default = NULL, tags = "train"),
sigdf = p_int(default = 3L, lower = 1L, tags = "train"),
sigquant = p_dbl(default = 0.90, lower = 0, upper = 1, tags = "train"),
k = p_dbl(default = 2.0, lower = 0, tags = "train"),
power = p_dbl(default = 2.0, lower = 0, tags = "train"),
base = p_dbl(default = 0.95, lower = 0, upper = 1, tags = "train"),
ndpost = p_int(default = 1000L, lower = 1L, tags = "train"),
nskip = p_int(default = 100L, lower = 0L, tags = "train"),
printevery = p_int(default = 100L, lower = 0L, tags = "train"),
keepevery = p_int(default = 1L, lower = 1L, tags = "train"),
keeptrainfits = p_lgl(default = TRUE, tags = "train"),
usequants = p_lgl(default = FALSE, tags = "train"),
numcut = p_int(default = 100L, lower = 1L, tags = "train"),
printcutoffs = p_int(default = 0, tags = "train"),
verbose = p_lgl(default = FALSE, tags = "train"),
nthread = p_int(default = 1L, tags = c("train", "threads")),
keeptrees = p_lgl(default = FALSE, tags = "train"),
keepcall = p_lgl(default = TRUE, tags = "train"),
sampleronly = p_lgl(default = FALSE, tags = "train"),
seed = p_int(default = NA_integer_, tags = "train", special_vals = list(NA_integer_)),
proposalprobs = p_uty(default = NULL, tags = "train"),
splitprobs = p_uty(default = NULL, tags = "train"),
keepsampler = p_lgl(default = NO_DEF, tags = "train")
)
ps$values = list(keeptrees = TRUE)
super$initialize(
id = "regr.bart",
packages = c("mlr3extralearners", "dbarts"),
feature_types = c("integer", "numeric", "factor", "ordered"),
# TODO: add "se" to the list of predict types.
predict_types = "response",
param_set = ps,
properties = c("weights"),
man = "mlr3extralearners::mlr_learners_regr.bart",
label = "Bayesian Additive Regression Trees"
)
}
),
private = list(
.train = function(task) {
pars = self$param_set$get_values(tags = "train")
# Extact just the features from the task data.
data = task$data(cols = task$feature_names)
# Convert from data.table to normal data.frame, just to be safe.
# Dbarts expects x.train to be a dataframe.
data.table::setDF(data)
# This will also extract a data.table
outcome = task$data(cols = task$target_names)
data.table::setDF(outcome)
# Outcome will now be a factor vector.
outcome = outcome[[1]]
if ("weights" %in% task$properties) {
pars$weights = task$weights$weight
}
# y.train should either be a binary factor or have values {0, 1}
invoke(dbarts::bart,
x.train = data, y.train = outcome,
.args = pars
)
},
.predict = function(task) {
pars = self$param_set$get_values(tags = "predict") # get parameters with tag "predict"
newdata = ordered_features(task, self)
# Other possible vars: offset.test, combineChains, ...
data.table::setDF(newdata)
# This will return a matrix of predictions, where each column is an observation
# and each row is a sample from the posterior.
p = invoke(predict, self$model, newdata = newdata, .args = pars)
# Transform predictions.
# TODO: confirm that this is the correct element name.
pred = colMeans(p)
list(response = pred)
}
)
)
.extralrns_dict$add("regr.bart", LearnerRegrBart)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.