X-Git-Url: https://git.auder.net/?p=valse.git;a=blobdiff_plain;f=pkg%2FR%2FselectVariables.R;h=b8ea1a07b0ba39aa42b461c3ecfa92c1a51ca729;hp=a4bc0f403960428f0934efd181ab7e97a22276ff;hb=6af1d4897dbab92a7be05068e0e15823378965d9;hpb=04845e3300b5450629bf1a2c3344d2f9419e91a6 diff --git a/pkg/R/selectVariables.R b/pkg/R/selectVariables.R index a4bc0f4..b8ea1a0 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,12 @@ #' @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 +#' @return a list, varying lambda in a grid, with selected (the indices of variables that are selected), +#' Rho (the covariance parameter, reparametrized), Pi (the proportion parameter) #' #' @export -#' selectVariables <- function(phiInit, rhoInit, piInit, gamInit, mini, maxi, gamma, glambda, X, Y, thresh = 1e-08, eps, ncores = 3, fast) { @@ -68,8 +67,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)