X-Git-Url: https://git.auder.net/?a=blobdiff_plain;f=R%2FinitSmallEM.R;fp=R%2FinitSmallEM.R;h=8cfb7e87c2fe65638e639c2b188a8bff0a00d4ab;hb=f2a9120810d7e1e423c7b5c2c4320f4e27221f50;hp=1fa2d9b9a5daf0b366b3314f35eba412200f2675;hpb=d1531659214edd6eaef0ac9ec835455614bba16c;p=valse.git

diff --git a/R/initSmallEM.R b/R/initSmallEM.R
index 1fa2d9b..8cfb7e8 100644
--- a/R/initSmallEM.R
+++ b/R/initSmallEM.R
@@ -22,12 +22,11 @@ initSmallEM = function(k,X,Y,tau)
   LLFinit1 = list()
   
   require(MASS) #Moore-Penrose generalized inverse of matrix
+  require(mclust) # K-means with selection of K
   for(repet in 1:20)
   {
-    clusters = hclust(dist(y)) #default distance : euclidean
-    #cutree retourne les indices (? quel cluster indiv_i appartient) d'un clustering hierarchique
-    clusterCut = cutree(clusters,k)
-    Zinit1[,repet] = clusterCut
+    clusters = Mclust(matrix(c(X,Y),nrow=n),k) #default distance : euclidean
+    Zinit1[,repet] = clusters$classification
     
     for(r in 1:k)
     {