- # Get summary "tableauRecap" from models ; TODO: jusqu'à ligne 114 à mon avis là c'est faux :/
- tableauRecap = t( sapply( models_list, function(models) {
- llh = do.call(rbind, lapply(models, function(model) model$llh))
- LLH = llh[-1,1]
- D = llh[-1,2]
- c(LLH, D, rep(k, length(model)), 1:length(model))
- } ))
- if (verbose)
- print('Model selection')
- tableauRecap = tableauRecap[rowSums(tableauRecap[, 2:4])!=0,]
- tableauRecap = tableauRecap[!is.infinite(tableauRecap[,1]),]
- data = cbind(1:dim(tableauRecap)[1], tableauRecap[,2], tableauRecap[,2], tableauRecap[,1])
+ # Get summary 'tableauRecap' from models
+ tableauRecap <- do.call(rbind, lapply(seq_along(models_list), function(i)
+ {
+ models <- models_list[[i]]
+ # For a collection of models (same k, several lambda):
+ LLH <- sapply(models, function(model) model$llh[1])
+ k <- length(models[[1]]$pi)
+ sumPen <- sapply(models, function(model) k * (dim(model$rho)[1] + sum(model$phi[,
+ , 1] != 0) + 1) - 1)
+ data.frame(model = paste(i, ".", seq_along(models), sep = ""), pen = sumPen/n,
+ complexity = sumPen, contrast = -LLH)
+ }))
+ tableauRecap <- tableauRecap[which(tableauRecap[, 4] != Inf), ]
+
+
+ if (verbose == TRUE)
+ print(tableauRecap)
+ modSel <- capushe::capushe(tableauRecap, n)
+ indModSel <- if (selecMod == "DDSE")
+ {
+ as.numeric(modSel@DDSE@model)
+ } else if (selecMod == "Djump")
+ {
+ as.numeric(modSel@Djump@model)
+ } else if (selecMod == "BIC")
+ {
+ modSel@BIC_capushe$model
+ } else if (selecMod == "AIC")
+ {
+ modSel@AIC_capushe$model
+ }
+
+ listMod <- as.integer(unlist(strsplit(as.character(indModSel), "[.]")))
+ modelSel <- models_list[[listMod[1]]][[listMod[2]]]
+ modelSel$tableau <- tableauRecap
+
+ if (plot)
+ print(plot_valse(X, Y, modelSel, n))