selective-inference / R-software

R software for selective inference
GNU General Public License v2.0
34 stars 18 forks source link

randomized lasso with empty inactive set #40

Open jonathan-taylor opened 6 years ago

jonathan-taylor commented 6 years ago

One thing that seems to be a problem is if there is no inactive set.

X = matrix(rnorm(500), 50, 10)

Y = rnorm(50)

library(selectiveInference)

randomizedLassoInf(X, Y, 1.e-5)

Error in [<-(*tmp*, inactive_set[i], i, value = 1) :

subscript out of bounds