diff abims_hclustering.r @ 1:36fc0a87d7fb draft default tip

planemo upload
author lecorguille
date Fri, 13 Nov 2015 08:43:11 -0500
parents 2f7381ee5235
children
line wrap: on
line diff
--- a/abims_hclustering.r	Tue Jun 30 06:36:09 2015 -0400
+++ b/abims_hclustering.r	Fri Nov 13 08:43:11 2015 -0500
@@ -21,10 +21,23 @@
 
     # -- Normalization: logratio --
     if (normalization) {
-	#meandata = apply(data,1,mean, na.rm=T)
-	#data = log2(data/meandata)
-	data=t(scale(t(data)))
+	    #meandata = apply(data,1,mean, na.rm=T)
+	    #data = log2(data/meandata)
+	    data=t(scale(t(data)))
+        
+        #AMAP: Unable to compute Hierarchical Clustering: missing values in distance matrix
+        #Erreur dans hcluster(x, method = method, link = link) :
+        #  Missing values in distance Matrix
+        #Calls: do.call -> <Anonymous> -> hclust2treeview -> hcluster
+        #Exécution arrêtée
+        data[is.nan(data)] = 0
     }
+    
+    #Erreur dans `[.default`(xj, i) :
+    #  les indices négatifs ne peuvent être mélangés qu'à des 0
+    #Calls: do.call ... r2cdt -> [ -> [.data.frame -> [ -> [.factor -> NextMethod
+    #Exécution arrêtée
+    data = data[!apply(data,1,sum)==0,]
 
     # -- hclust / output files for TreeView --
     file="hclust.cdt"