- 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))
+ LLH <- sapply( models, function(model) model$llh )
+ k == length(models[[1]]$pi)
+ # TODO: chuis pas sûr du tout des lignes suivantes...
+ # J'ai l'impression qu'il manque des infos
+ sumPen = sapply( models, function(model)
+ sum( model$pi^gamma * sapply(1:k, function(r) sum(abs(model$phi[,,r]))) ) )
+ data.frame(model=paste(i,".",seq_along(models),sep=""),
+ pen=sumPen/1000, complexity=sumPen, contrast=LLH)