Add files via upload

hanjunlee21 · Aug 22, 2022 · 3a9b034 · 3a9b034
1 parent e7d66da
commit 3a9b034
Show file tree

Hide file tree

Showing 2 changed files with 12 additions and 12 deletions.
diff --git a/bedGraph2Cluster b/bedGraph2Cluster
diff --git a/bedGraph2Cluster.m b/bedGraph2Cluster.m
@@ -242,24 +242,24 @@ function bedGraph2Cluster(bedGraphs_Peak, bedGraphs_Cluster, bedGraphs_Noncluste
 chrs = ["chr1";"chr2";"chr3";"chr4";"chr5";"chr6";"chr7";"chr8";"chr9";"chr10";"chr11";"chr12";"chr13";"chr14";"chr15";"chr16";"chr17";"chr18";"chr19";"chr20";"chr21";"chr22";"chrX";"chrY"];
 fileID = fopen(convertStringsToChars(strcat(outputdir,"/peaks.bed")), 'w');
 for i = 1:size(X.peak.chr,1)
-    fprintf('%s\t',chrs(X.peak.chr(i,1),1));
-    fprintf('%s\t',convertCharsToStrings(num2str(X.peak.st(i,1)+1)));
-    fprintf('%s\t',convertCharsToStrings(num2str(X.peak.en(i,1))));
-    fprintf('%s\t',convertCharsToStrings(num2str(i)));
-    fprintf('%s\t',"1000");
-    fprintf('%s\n',".");
+    fprintf(fileID, '%s\t',chrs(X.peak.chr(i,1),1));
+    fprintf(fileID, '%s\t',convertCharsToStrings(num2str(X.peak.st(i,1)+1)));
+    fprintf(fileID, '%s\t',convertCharsToStrings(num2str(X.peak.en(i,1))));
+    fprintf(fileID, '%s\t',convertCharsToStrings(num2str(i)));
+    fprintf(fileID, '%s\t',"1000");
+    fprintf(fileID, '%s\n',".");
 end
 fclose(fileID);
 for clusterID = 1:k
     fileID = fopen(convertStringsToChars(strcat(outputdir,"/peaks.clust",convertCharsToStrings(num2str(k)),".",convertCharsToStrings(num2str(clusterID)),".bed")), 'w');
     for i = 1:size(X.peak.chr,1)
         if X.peak.(['clust',num2str(k)])(i,1) == clusterID
-            fprintf('%s\t',chrs(X.peak.chr(i,1),1));
-            fprintf('%s\t',convertCharsToStrings(num2str(X.peak.st(i,1)+1)));
-            fprintf('%s\t',convertCharsToStrings(num2str(X.peak.en(i,1))));
-            fprintf('%s\t',convertCharsToStrings(num2str(i)));
-            fprintf('%s\t',"1000");
-            fprintf('%s\n',".");
+            fprintf(fileID, '%s\t',chrs(X.peak.chr(i,1),1));
+            fprintf(fileID, '%s\t',convertCharsToStrings(num2str(X.peak.st(i,1)+1)));
+            fprintf(fileID, '%s\t',convertCharsToStrings(num2str(X.peak.en(i,1))));
+            fprintf(fileID, '%s\t',convertCharsToStrings(num2str(i)));
+            fprintf(fileID, '%s\t',"1000");
+            fprintf(fileID, '%s\n',".");
         end
     end
     fclose(fileID);