fk_quanti:

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

Usage

fk_quanti(Xquanti, mean, sigma, Y, k = length(sigma), n = nrow(Xquanti), p = length(colnames(Xquanti)))

Arguments

Xquanti
mean
sigma
Y
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 (Xquanti, mean, sigma, Y, k = length(sigma), n = nrow(Xquanti), 
    p = length(colnames(Xquanti))) 
{
    K = length(unique(Y))
    fk_quanti_ik <- matrix(1, nrow = n, ncol = length(unique(Y)))
    fkx_q <- as.list(seq_len(length(Y)))
    for (i in 1:n) {
        for (k in 1:K) {
            fk_quanti_ik[i, k] <- (1/(((2 * pi)^(p/2)) * (det(as.matrix(sigma[[k]]))^(0.5)))) * 
                exp((-0.5) * as.matrix(Xquanti[i, ] - as.numeric(mean[k, 
                  -1])) %*% ginv(as.matrix(sigma[[k]])) %*% t(as.matrix(Xquanti[i, 
                  ] - as.numeric(mean[k, -1]))))
        }
        fkx_q[[i]] <- fk_quanti_ik[i, ]
    }
    return(fkx_q)
  }

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