- # Get summary "tableauRecap" from models
- tableauRecap = do.call( rbind, lapply( models_list, function(models) {
- #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) #TODO: LLF? harmonize between EMGLLF and EMGrank?
- LLH = llh[-1,1]
- D = llh[-1,2]
- k = length(models[[1]]$pi)
- cbind(LLH, D, rep(k, length(models)), 1:length(models))
- } ) )
- 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')
- 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
- models_list[[tableauRecap[indModSel,3]]][[tableauRecap[indModSel,4]]]