X-Git-Url: https://git.auder.net/?p=valse.git;a=blobdiff_plain;f=pkg%2FR%2Fselectiontotale.R;h=bfb932564fe777eef9aadb39a51b1740e317e3da;hp=042c70b2ce6c5ea0c333d0a31aa5f1012851c03b;hb=25f737999f505619008eeaa0a35eac1e34a52d5b;hpb=610cfd0d7ee1061c884583c8f3f456b246eadbf8

diff --git a/pkg/R/selectiontotale.R b/pkg/R/selectiontotale.R
index 042c70b..bfb9325 100644
--- a/pkg/R/selectiontotale.R
+++ b/pkg/R/selectiontotale.R
@@ -40,7 +40,7 @@ selectiontotale = function(phiInit,rhoInit,piInit,gamInit,mini,maxi,gamma,glambd
       selectedVariables[[lambdaIndex]] = sapply(1:p, function(j) {
         #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(apply( abs(params$phi[j,,]) > thresh, 1, any ) )))
       })
       Rho[[lambdaIndex]] = params$Rho
       Pi[[lambdaIndex]] = params$Pi