add some comments after meeting 17-11-2016

[epclust.git] / code / stage2 / src / 02_cluster_2009.r

diff --git a/code/stage2/src/02_cluster_2009.r b/code/stage2/src/02_cluster_2009.r
index f08d82b..3514ce7 100644 (file)
--- a/code/stage2/src/02_cluster_2009.r
+++ b/code/stage2/src/02_cluster_2009.r
@@ -24,7 +24,7 @@ sdcontrib <- apply(matcontrib0, 1, sd)
 lims      <- quantile(sdcontrib, probs = c(.005, .995)) # obtain 1%-extreme data
 is_normal <- which((sdcontrib > lims[1]) & (sdcontrib < lims[2]))
 
-matcontri_ext <- matcontrib0[-is_normal, ]""
+matcontri_ext <- matcontrib0[-is_normal, ]#""
 matcontrib    <- matcontrib0[is_normal, ]     # wipe out aberrant data
 
 matcontrib <- t(apply(matcontrib, 1, function(x) x / sum(x)))
@@ -39,6 +39,9 @@ selvar <- ci$selectv
 ## c. Clustering  ##########
 K <- 200
 system.time(
+
+#TODO: cette partie en C
+
   clfit <- clara(x        = tdata[, selvar], 
                  k        = K, 
                  sampsize = 4000,