bsnullinteract: Compute bootstrapped null interaction prediction rule...

View source: R/pre.R

bsnullinteractR Documentation

Compute bootstrapped null interaction prediction rule ensembles


bsnullinteract generates bootstrapped null interaction models, which can be used to derive a reference distribution of the test statistic calculated with interact.


  nsamp = 10,
  parallel = FALSE,
  penalty.par.val = "lambda.1se",
  verbose = FALSE,



object of class pre.


numeric. Number of bootstrapped null interaction models to be derived.


logical. Should parallel foreach be used to generate initial ensemble? Must register parallel beforehand, such as doMC or others.


character or numeric. Value of the penalty parameter λ to be employed for selecting the final ensemble. The default "lambda.min" employs the λ value within 1 standard error of the minimum cross-validated error. Alternatively, "lambda.min" may be specified, to employ the λ value with minimum cross-validated error, or a numeric value >0 may be specified, with higher values yielding a sparser ensemble. To evaluate the trade-off between accuracy and sparsity of the final ensemble, inspect pre_object$ and plot(pre_object$


logical. should progress be printed to the command line?


Further arguments to be passed to predict.pre.


Note that computation of bootstrapped null interaction models is computationally intensive. The default number of samples is set to 10, but for reliable results argument nsamp should be set to a higher value (e.g., ≥ 100).

See also section 8.3 of Friedman & Popescu (2008).


A list of length nsamp with null interaction models, to be used as input for interact.


Fokkema, M. (2020). Fitting prediction rule ensembles with R package pre. Journal of Statistical Software, 92(12), 1-30. doi: 10.18637/jss.v092.i12

Friedman, J. H., & Popescu, B. E. (2008). Predictive learning via rule ensembles. The Annals of Applied Statistics, 2(3), 916-954, doi: 10.1214/07-AOAS148.

See Also

pre, interact


airq.ens <- pre(Ozone ~ ., data=airquality[complete.cases(airquality),])
nullmods <- bsnullinteract(airq.ens)
interact(airq.ens, nullmods = nullmods, col = c("#7FBFF5", "#8CC876"))

pre documentation built on June 11, 2022, 1:10 a.m.