X-Git-Url: https://git.auder.net/?p=valse.git;a=blobdiff_plain;f=pkg%2FR%2FselectVariables.R;h=2d1c9b7b71c16d5e7c15453be7d37d65fae501f9;hp=a4bc0f403960428f0934efd181ab7e97a22276ff;hb=fb3557f39487d9631ffde30f20b70938d2a6ab0c;hpb=6775f5b98ffc7eae7ce9d4081b23b39ce66d3c0b diff --git a/pkg/R/selectVariables.R b/pkg/R/selectVariables.R index a4bc0f4..2d1c9b7 100644 --- a/pkg/R/selectVariables.R +++ b/pkg/R/selectVariables.R @@ -1,6 +1,6 @@ #' selectVariables #' -#' It is a function which constructs, for a given lambda, the sets for each cluster of relevant variables. +#' For a given lambda, construct the sets of relevant variables for each cluster. #' #' @param phiInit an initial estimator for phi (size: p*m*k) #' @param rhoInit an initial estimator for rho (size: m*m*k) @@ -15,13 +15,11 @@ #' @param thresh real, threshold to say a variable is relevant, by default = 1e-8 #' @param eps threshold to say that EM algorithm has converged #' @param ncores Number or cores for parallel execution (1 to disable) +#' @param fast boolean to enable or not the C function call #' #' @return a list of outputs, for each lambda in grid: selected,Rho,Pi #' -#' @examples TODO -#' #' @export -#' selectVariables <- function(phiInit, rhoInit, piInit, gamInit, mini, maxi, gamma, glambda, X, Y, thresh = 1e-08, eps, ncores = 3, fast) { @@ -68,8 +66,7 @@ selectVariables <- function(phiInit, rhoInit, piInit, gamInit, mini, maxi, gamma if (ncores > 1) parallel::stopCluster(cl) - print(out) - # Suppress models which are computed twice + # Suppress models which are computed twice # sha1_array <- lapply(out, digest::sha1) out[ duplicated(sha1_array) ] selec <- lapply(out, function(model) model$selected) ind_dup <- duplicated(selec)