- if (! requireNamespace("capushe", quietly=TRUE))
- {
- warning("'capushe' not available: returning all models")
- return (models_list)
- }
+ if (!requireNamespace("capushe", quietly = TRUE)) {
+ warning("'capushe' not available: returning all models")
+ return(models_list)
+ }
+
+ # 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) print(tableauRecap)