X-Git-Url: https://git.auder.net/?p=valse.git;a=blobdiff_plain;f=pkg%2FR%2FselectVariables.R;h=869e7bfc99967f5f4e79b131200c625de39e0a4e;hp=b4fc0abbcc02e8634671cadfb7591777904d8926;hb=4cc632c9a1e1d93e9a43a402d1361f23afc50e5e;hpb=086ca318ed5580e961ceda3f1e122a2da58e4427 diff --git a/pkg/R/selectVariables.R b/pkg/R/selectVariables.R index b4fc0ab..869e7bf 100644 --- a/pkg/R/selectVariables.R +++ b/pkg/R/selectVariables.R @@ -14,6 +14,7 @@ #' @param Y matrix of responses #' @param thres threshold to consider a coefficient to be equal to 0 #' @param tau threshold to say that EM algorithm has converged +#' @param ncores Number or cores for parallel execution (1 to disable) #' #' @return a list of outputs, for each lambda in grid: selected,Rho,Pi #' @@ -22,7 +23,7 @@ #' @export #' selectVariables = function(phiInit,rhoInit,piInit,gamInit,mini,maxi,gamma,glambda, - X,Y,thresh,tau, ncores=1) #ncores==1 ==> no // + X,Y,thresh,tau, ncores=3) { if (ncores > 1) { @@ -54,7 +55,8 @@ selectVariables = function(phiInit,rhoInit,piInit,gamInit,mini,maxi,gamma,glambd out <- if (ncores > 1) parLapply(cl, glambda, computeCoefs) - else lapply(glambda, computeCoefs) + else + lapply(glambda, computeCoefs) if (ncores > 1) parallel::stopCluster(cl)