X-Git-Url: https://git.auder.net/?p=valse.git;a=blobdiff_plain;f=R%2FmodelSelection.R;fp=R%2FmodelSelection.R;h=0000000000000000000000000000000000000000;hp=86e2efd563829e9474b6769a03a76803b907be16;hb=f87ff0f5116c0c1c59c5608e46563ff0f79e5d43;hpb=53fa233d8fbeaf4d51a4874ba69d8472d01d04ba diff --git a/R/modelSelection.R b/R/modelSelection.R deleted file mode 100644 index 86e2efd..0000000 --- a/R/modelSelection.R +++ /dev/null @@ -1,40 +0,0 @@ -#' Among a collection of models, this function constructs a subcollection of models with -#' models having strictly different dimensions, keeping the model which minimizes -#' the likelihood if there were several with the same dimension -#' -#' @param LLF a matrix, the first column corresponds to likelihoods for several models -#' the second column corresponds to the dimensions of the corresponding models. -#' -#' @return a list with indices, a vector of indices selected models, -#' and D1, a vector of corresponding dimensions -#' @export -#' -modelSelection = function(LLF) -{ - D = LLF[,2] - D1 = unique(D) - - indices = rep(1, length(D1)) - #select argmax MLE - if (length(D1)>2) - { - for (i in 1:length(D1)) - { - A = c() - for (j in 1:length(D)) - { - if(D[[j]]==D1[[i]]) - a = c(a, LLF[j,1]) - } - b = max(a) - #indices[i] : first indices of the binary vector where u_i ==1 - indices[i] = which.max(LLF == b) - } - } - - return (list(indices=indices,D1=D1)) -} - -#TODO: -## Programme qui sélectionne un modèle -## proposer à l'utilisation différents critères (BIC, AIC, slope heuristic)