#'
#' export
constructionModelesLassoRank = function(pi, rho, mini, maxi, X, Y, tau, A1, rangmin,
- rangmax, ncores, verbose=FALSE)
+ rangmax, ncores, fast=TRUE, verbose=FALSE)
{
n = dim(X)[1]
p = dim(X)[2]
if (ncores > 1)
{
- cl = parallel::makeCluster(ncores)
+ cl = parallel::makeCluster(ncores, outfile='')
parallel::clusterExport( cl, envir=environment(),
varlist=c("A1","Size","Pi","Rho","mini","maxi","X","Y","tau",
"Rank","m","phi","ncores","verbose") )
for (j in 1:Size)
{
res = EMGrank(Pi[,lambdaIndex], Rho[,,,lambdaIndex], mini, maxi,
- X[,active], Y, tau, Rank[j,])
+ X[,active], Y, tau, Rank[j,], fast)
llh = rbind(llh,
c( res$LLF, sum(Rank[j,] * (length(active)- Rank[j,] + m)) ) )
phi[active,,,] = rbind(phi[active,,,], res$phi)