- .Call("EMGrank", Pi, Rho, mini, maxi, X, Y, tau, rank, PACKAGE="valse")
+ #TEMPORARY: use R version
+ return (EMGrank_R(Pi, Rho, mini, maxi, X, Y, tau, rank))
+
+ n = nrow(X) #nombre d'echantillons
+ p = ncol(X) #nombre de covariables
+ m = ncol(Y) #taille de Y (multivarié)
+ k = length(Pi) #nombre de composantes dans le mélange
+ .Call("EMGrank",
+ Pi, Rho, mini, maxi, X, Y, tau, rank,
+ phi=double(p*m*k), LLF=double(1),
+ n, p, m, k,
+ PACKAGE="valse")