Wrapper function for applying bootstrap estimation using gradient boosting.

## Bootstrap boosting.
bboost(..., data, type = 1, cores = 1,
  n = 2, prob = 0.623, fmstop = NULL,
  trace = TRUE, drop = FALSE, replace = FALSE)

## Plotting function.
bboost_plot(object, col = NULL)

## Predict method.
# S3 method for bboost
predict(object, newdata, ..., cores = 1, pfun = NULL)



Arguments passed to bamlss and predict.bamlss.


The data frame to be used for modeling.


Type of algorithm, type = 1 uses all observations and samples with replacement, type = 2 uses only a fraction specified in prob and samples with replacement.


The number of cores to be used.


The number of bootstrap iterations.


The fraction that should be used to fit the model in each bootstrap iteration.


The function that should return the optimum stopping iteration. The function must have two arguments: (1) the model end (2) the data. The function must return a list with two named arguments: (1) "mstop" the optimum stopping iteration and (2) a vector of the objective criterion that should be evaluated by the hold out sample data during each bootstrap iteration. See the examples.


Prints out the current state of the bootstrap algorithm.


Should only the best set of parameters be saved?


Sampling with replacement, or sampling ceiling(nobs * prob) rows of the data for fitting the n models.


The "bboost" object used for prediction and plotting.


The color that should be used for plotting.


The data frame predictions should be made for.


The prediction function that should be used, for example predictn could be used, too. Note that this is experimental.


A list of bamlss objects.

See also


if (FALSE) ## Simulate data.
d <- GAMart()

## Estimate model.
f <- num ~ s(x1) + s(x2) + s(x3) + s(lon,lat)

## Function for evaluation of hold out sample
## criterion to find the optimum mstop.
fmstop <- function(model, data) {
  p <- predict(model, newdata = data, model = "mu")
  mse <- NULL
  for(i in 1:nrow(model$parameters))
    mse <- c(mse, mean((data$num - p[, i])^2))
  list("MSE" = mse, "mstop" = which.min(mse))

## Bootstrap boosted models.
b <- bboost(f, data = d, n = 50, cores = 3, fmstop = fmstop)

## Plot hold out sample MSE.

## Predict for each bootstrap sample.
nd <- data.frame("x2" = seq(0, 1, length = 100))
p <- predict(b, newdata = nd, model = "mu", term = "x2")
plot2d(p ~ x2, data = nd)
#> Error in h(simpleError(msg, call)): error in evaluating the argument 'x' in selecting a method for function 'as.matrix': object 'p' not found