fk_quali:

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

Usage

fk_quali(Y, Xijh, alpha_ki, k = length(alpha_ki), n = nrow(Xquali), p = length(colnames(Xquali)))

Arguments

Y
Xijh
alpha_ki
k
n
p

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 (Y, Xijh, alpha_ki, k = length(alpha_ki), n = nrow(Xquali), 
    p = length(colnames(Xquali))) 
{
    if (p == 0) {
        p <- 1
        n <- length(Xquali)
    }
    fk_quali_ik <- matrix(1, nrow = n, ncol = length(unique(Y)))
    K = length(unique(Y))
    fkx_q <- as.list(seq_len(length(Y)))
    for (i in 1:n) {
        for (k in 1:K) {
            for (j in 1:p) {
                fk_quali_ik[i, k] <- fk_quali_ik[i, k] * (as.numeric(alpha_ki[[k]][[j]][["proba"]] %*% 
                  as.numeric(Xijh[[j]][i, -length(Xijh[[j]][i, 
                    ])])))
            }
        }
        fkx_q[[i]] <- fk_quali_ik[i, ]
    }
    return(fkx_q)
  }

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