print(tableauRecap)
tableauRecap = tableauRecap[which(tableauRecap[,4]!= Inf),]
+
modSel = capushe::capushe(tableauRecap, n)
indModSel <-
if (selecMod == 'DDSE')
modSel@BIC_capushe$model
else if (selecMod == 'AIC')
modSel@AIC_capushe$model
-
+
mod = as.character(tableauRecap[indModSel,1])
listMod = as.integer(unlist(strsplit(mod, "[.]")))
modelSel = models_list[[listMod[1]]][[listMod[2]]]
-
+
##Affectations
Gam = matrix(0, ncol = length(modelSel$pi), nrow = n)
for (i in 1:n){
Gam = Gam/rowSums(Gam)
modelSel$affec = apply(Gam, 1,which.max)
modelSel$proba = Gam
-
+
if (plot){
print(plot_valse(X,Y,modelSel,n))
}
-
+
return(modelSel)
}