Em_quali:

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

Usage

Em_quali(Xquali, alpha_ki = NULL, Xijh_i = NULL, pki = NULL, Y = NULL, methode = "init", e = 0.1, nbr_iteration = NULL, k)

Arguments

Xquali
alpha_ki
Xijh_i
pki
Y
methode
e
nbr_iteration
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 (Xquali, alpha_ki = NULL, Xijh_i = NULL, pki = NULL, 
    Y = NULL, methode = "init", e = 0.1, nbr_iteration = NULL, 
    k) 
{
    if (methode == "kmodes") {
        Y <- kmodes(Xquali, k)$cluster
    }
    else {
        Y <- init_class(iris, k)
    }
    if (is.null(alpha_ki)) {
        alpha_ki <- alphai(Y, Xquali)
    }
    if (is.null(Xijh_i)) {
        Xijh_i <- xijh(Xquali)
    }
    if (is.null(pki)) {
        pki <- as.numeric(table(Y)/length(Y))
    }
    Pki <- pki
    Alphaki <- alpha_ki
    p = length(colnames(Xquali))
    n <- nrow(Xquali)
    fkx <- fk_quali(Y, Xijh_i, alpha_ki, p = p, n = n)
    f <- fx(Pki, fkx, Y)
    tkx <- tk(fkx, Pki)
    li <- l(f)
    Qi <- Qt(tkx, fkx, Pki)
    BIC <- Bic(f, Y, Xquali = Xquali)
    ICL <- BIC - sum(tkx * log(tkx))
    i <- 1
    teta <- list()
    K <- k
    repeat {
        Y <- apply(tkx, MARGIN = 1, which.max)
        Nk <- nk(tkx)
        Pk <- pk(Nk, Xquali)
        Alphak <- alphak(tkx, Xijh_i, Y, Xquali)
        fkx <- fk_quali(Y, Xijh_i, Alphak, p = p, n = n)
        f <- fx(Pk, fkx, Y)
        lf <- l(f)
        tkx <- tk(fkx, Pk)
        Qf <- Qt(tkx, fkx, Pk)
        BIC <- Bic(f, Y, Xquali = Xquali)
        ICL <- BIC - sum(tkx * log(tkx))
        teta[[i]] <- list(iteration = i, Proba = Pk, Alphak = Alphak, 
            log_like = lf, Q = Qf, Bic = BIC, ICL = ICL, proba_individu = fkx, 
            part_MAP = tkx)
        if (!is.null(nbr_iteration)) {
            if (i >= nbr_iteration) 
                break
        }
        if (abs(li - lf) < e) {
            break
        }
        else {
            Pki <- Pk
            Alphaki <- Alphak
            li <- lf
            i <- i + 1
        }
    }
    return(teta)
  }

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