- .Call("EMGLLF", phiInit, rhoInit, piInit, gamInit,
- mini, maxi, gamma, lambda, X, Y, tau, PACKAGE="valse")
+ n = nrow(X) #nombre d'echantillons
+ p = ncol(X) #nombre de covariables
+ m = ncol(Y) #taille de Y (multivarié)
+ k = length(piInit) #nombre de composantes dans le mélange
+ .Call("EMGLLF",
+ phiInit, rhoInit, piInit, gamInit, mini, maxi, gamma, lambda, X, Y, tau,
+ phi=double(p*m*k), rho=double(m*m*k), pi=double(k), LLF=double(maxi),
+ S=double(p*m*k), affec=integer(n),
+ n, p, m, k,
+ PACKAGE="valse")