Update file make_heatmap.R

bin/make_heatmap.R

@@ -23,6 +23,12 @@ setkey(trait_name, "ID")
 Rho = fread(correlation_file)
 setkey(Rho, "V1")
 
+D = as.matrix(1-Rho[,-1])
+D[is.na(D)] = 0
+h = hclust(as.dist(D))
+order_trait <- trait_name$ID[h$order]
+
+
 Pval = fread(p_file)
 setkey(Pval, "V1")