## Not run:
# x <- matrix(rnorm(100*100), nrow = 100, ncol = 100)
# y <- x[,1] * 2 + x[,2] * 2.5 + rnorm(100)
# sel <- bst.sel(x, y, q=10)
# library("hdi")
# fit.multi <- hdi(x, y, method = "multi.split",
# model.selector =bst.sel,
# args.model.selector=list(type="firstq", q=10))
# fit.multi
# fit.multi$pval[1:10] ## the first 10 p-values
# fit.multi <- hdi(x, y, method = "multi.split",
# model.selector =bst.sel,
# args.model.selector=list(type="cv"))
# fit.multi
# fit.multi$pval[1:10] ## the first 10 p-values
# ## End(Not run)
Run the code above in your browser using DataLab