few updates

diff --git a/pkg/R/selectiontotale.R b/pkg/R/selectiontotale.R
index 042c70b..bfb9325 100644 (file)
--- 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