X-Git-Url: https://git.auder.net/?p=valse.git;a=blobdiff_plain;f=pkg%2FR%2Fmain.R;h=93f8e3f3b444cf3d2597ab8bb2e87071f83a6622;hp=bff2ec5b70e971cae6f7bfa64cd7a728949144c8;hb=b9b0b72a2c8f7f0d1a3216528aefcec0a92c6c99;hpb=08f4604c778da8af7e26b52b1d433a6be82c3139

diff --git a/pkg/R/main.R b/pkg/R/main.R
index bff2ec5..93f8e3f 100644
--- a/pkg/R/main.R
+++ b/pkg/R/main.R
@@ -106,7 +106,7 @@ valse = function(X, Y, procedure='LassoMLE', selecMod='DDSE', gamma=1, mini=10,
 		#Pour un groupe de modeles (même k, différents lambda):
 		llh = matrix(ncol = 2)
 		for (l in seq_along(models))
-			llh = rbind(llh, models[[l]]$llh)
+			llh = rbind(llh, models[[l]]$llh) #TODO: LLF? harmonize between EMGLLF and EMGrank?
 		LLH = llh[-1,1]
 		D = llh[-1,2]
 		k = length(models[[1]]$pi)
@@ -115,7 +115,7 @@ valse = function(X, Y, procedure='LassoMLE', selecMod='DDSE', gamma=1, mini=10,
 	tableauRecap = tableauRecap[rowSums(tableauRecap[, 2:4])!=0,]
   tableauRecap = tableauRecap[(tableauRecap[,1])!=Inf,]
   data = cbind(1:dim(tableauRecap)[1], tableauRecap[,2], tableauRecap[,2], tableauRecap[,1])
-
+browser()
   modSel = capushe::capushe(data, n)
   indModSel <-
 		if (selecMod == 'DDSE')