- #selectedVariables: list where element j contains vector of selected variables in [1,m]
- selectedVariables = sapply(1:p, function(j) { ## je me suis permise de changer le type,
- ##une liste de liste ca devenait compliqué je trouve pour choper ce qui nous intéresse
- #from boolean matrix mxk of selected variables obtain the corresponding boolean m-vector,
- #and finally return the corresponding indices
- #seq_len(m)[ apply( abs(params$phi[j,,]) > thresh, 1, any ) ]
- c(seq_len(m)[ apply( abs(params$phi[j,,]) > thresh, 1, any ) ],
- rep(0, m-length(seq_len(m)[ apply( abs(params$phi[j,,]) > thresh, 1, any ) ] ) ))
- })
+ # selectedVariables: list where element j contains vector of selected variables
+ # in [1,m]
+ selectedVariables <- lapply(1:p, function(j) {
+ # from boolean matrix mxk of selected variables obtain the corresponding boolean
+ # m-vector, and finally return the corresponding indices
+ if (m>1) {
+ seq_len(m)[apply(abs(params$phi[j, , ]) > thresh, 1, any)]
+ } else {
+ if (any(params$phi[j, , ] > thresh))
+ 1
+ else
+ numeric(0)
+ }
+ })