R Under development (unstable) (2025-04-16 r88149 ucrt) -- "Unsuffered Consequences"
Copyright (C) 2025 The R Foundation for Statistical Computing
Platform: x86_64-w64-mingw32/x64

R is free software and comes with ABSOLUTELY NO WARRANTY.
You are welcome to redistribute it under certain conditions.
Type 'license()' or 'licence()' for distribution details.

R is a collaborative project with many contributors.
Type 'contributors()' for more information and
'citation()' on how to cite R or R packages in publications.

Type 'demo()' for some demos, 'help()' for on-line help, or
'help.start()' for an HTML browser interface to help.
Type 'q()' to quit R.

> #########################
> ## Stability selection ##
> #########################
> 
> library(hdi)
Loading required package: scalreg
Loading required package: lars
Loaded lars 1.3

> 
> set.seed(123)
> 
> x <- matrix(rnorm(100*1000), nrow = 100, ncol = 1000)
> y <- x[,1] + x[,2] + rnorm(100)
> 
> fit.stab <- stability(x, y, EV = 1)
> fit.tmp <- stability(x, y, EV = 1, verbose = TRUE) ## to check verbose
...Subsample 1 
...Subsample 2 
...Subsample 3 
...Subsample 4 
...Subsample 5 
...Subsample 6 
...Subsample 7 
...Subsample 8 
...Subsample 9 
...Subsample 10 
...Subsample 11 
...Subsample 12 
...Subsample 13 
...Subsample 14 
...Subsample 15 
...Subsample 16 
...Subsample 17 
...Subsample 18 
...Subsample 19 
...Subsample 20 
...Subsample 21 
...Subsample 22 
...Subsample 23 
...Subsample 24 
...Subsample 25 
...Subsample 26 
...Subsample 27 
...Subsample 28 
...Subsample 29 
...Subsample 30 
...Subsample 31 
...Subsample 32 
...Subsample 33 
...Subsample 34 
...Subsample 35 
...Subsample 36 
...Subsample 37 
...Subsample 38 
...Subsample 39 
...Subsample 40 
...Subsample 41 
...Subsample 42 
...Subsample 43 
...Subsample 44 
...Subsample 45 
...Subsample 46 
...Subsample 47 
...Subsample 48 
...Subsample 49 
...Subsample 50 
...Subsample 51 
...Subsample 52 
...Subsample 53 
...Subsample 54 
...Subsample 55 
...Subsample 56 
...Subsample 57 
...Subsample 58 
...Subsample 59 
...Subsample 60 
...Subsample 61 
...Subsample 62 
...Subsample 63 
...Subsample 64 
...Subsample 65 
...Subsample 66 
...Subsample 67 
...Subsample 68 
...Subsample 69 
...Subsample 70 
...Subsample 71 
...Subsample 72 
...Subsample 73 
...Subsample 74 
...Subsample 75 
...Subsample 76 
...Subsample 77 
...Subsample 78 
...Subsample 79 
...Subsample 80 
...Subsample 81 
...Subsample 82 
...Subsample 83 
...Subsample 84 
...Subsample 85 
...Subsample 86 
...Subsample 87 
...Subsample 88 
...Subsample 89 
...Subsample 90 
...Subsample 91 
...Subsample 92 
...Subsample 93 
...Subsample 94 
...Subsample 95 
...Subsample 96 
...Subsample 97 
...Subsample 98 
...Subsample 99 
...Subsample 100 
> 
> fit.stab
Selected predictors:
--------------------
[1] 1 2
--------------------
Expected number of false positives controlled at level 1 
> fit.stab$freq[1:10] ## selection frequency of the first 10 predictors
 [1] 0.99 1.00 0.00 0.00 0.01 0.01 0.00 0.01 0.13 0.01
> 
> ncores <- if(.Platform$OS.type == "windows") 1 else getOption("mc.cores", 2L)
> fit.stab2 <- stability(x, y, EV = 1, parallel = TRUE, ncores = ncores,
+                        verbose = TRUE)
...starting parallelization of bootstrap samples
> fit.stab2
Selected predictors:
--------------------
[1] 1 2
--------------------
Expected number of false positives controlled at level 1 
> fit.stab2$freq[1:10]
 [1] 0.99 1.00 0.00 0.00 0.02 0.04 0.02 0.00 0.13 0.03
>  
> stopifnot(all.equal(fit.stab$select, fit.stab2$select))
> 
> 
> 
> 
> 
> 
> 
> 
> proc.time()
   user  system elapsed 
   3.71    0.48    4.18