Bic:

In SaraTouzani/EM: Algorithme EM ("Espérance-Maximisation")

Usage

Bic(f, Y = NULL, nbr_par = NULL, Xquali = NULL, Xquanti = NULL, K = length(unique(Y)))

Arguments

f
Y
nbr_par
Xquali
Xquanti
K

Examples

##---- Should be DIRECTLY executable !! ----
##-- ==>  Define data, use random,
##--	or do  help(data=index)  for the standard data sets.

## The function is currently defined as
function (f, Y = NULL, nbr_par = NULL, Xquali = NULL, Xquanti = NULL, 
    K = length(unique(Y))) 
{
    nbr_par <- 0
    if (!is.null(Xquali)) {
        p <- ncol(Xquali)
        for (j in 1:p) {
            nbr_par <- nbr_par + length(unique(Xquali[, j])) - 
                1
        }
        nbr_par <- nbr_par * K
    }
    if (!is.null(Xquanti)) {
        p <- ncol(Xquanti)
        nbr_par <- K * (p * (p + 1)/2 + p + 1) - 1 + nbr_par
    }
    return(-2 * sum(log(f)) + nbr_par * log(length(Y)))
  }

SaraTouzani/EM documentation built on Dec. 16, 2019, 12:31 a.m.