meank:

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

Usage

meank(Nk, tkx, Xquanti, Y)

Arguments

Nk
tkx
Xquanti
Y

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 (Nk, tkx, Xquanti, Y) 
{
    K <- length(Nk)
    p <- ncol(Xquanti)
    moyenne <- matrix(0, ncol = p, nrow = K)
    moyenne[, 1] <- unique(Y)
    for (k in 1:K) {
        for (j in 1:p) {
            moyenne[k, j] <- sum((tkx[, k] * Xquanti)[, j])/Nk[k]
        }
    }
    moyenne <- cbind(sort(unique(Y)), moyenne)
    return(moyenne)
  }

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