#'
#' export
constructionModelesLassoRank = function(pi, rho, mini, maxi, X, Y, tau, A1, rangmin,
#'
#' export
constructionModelesLassoRank = function(pi, rho, mini, maxi, X, Y, tau, A1, rangmin,
parallel::clusterExport( cl, envir=environment(),
varlist=c("A1","Size","Pi","Rho","mini","maxi","X","Y","tau",
"Rank","m","phi","ncores","verbose") )
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,
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)
llh = rbind(llh,
c( res$LLF, sum(Rank[j,] * (length(active)- Rank[j,] + m)) ) )
phi[active,,,] = rbind(phi[active,,,], res$phi)