Merge pull request #66 from haizi-zh/master

🐛 Fix a critical bug in hic2mat/bed2mat
TaoYang-dev · Dec 13, 2020 · e485dfa · e485dfa
2 parents ce2c3a8 + a42551c
commit e485dfa
Showing 1 changed file with 4 additions and 4 deletions.
diff --git a/R/bed2mat.R b/R/bed2mat.R
@@ -28,11 +28,11 @@
 bed2mat <- function(bed, resol = "NONE"){ 
   # n the max number of bin
   if (resol == "NONE") {
-    N = max(bed[,c(1:2)])
+    N = max(bed[,c(1:2)]) + 1
   } else {
-    N = max(bed[,c(1:2)])/resol
-    bed[,1] = bed[,1]/resol
-    bed[,2] = bed[,2]/resol
+    N = max(bed[,c(1:2)])/resol + 1
+    bed[,1] = bed[,1]/resol + 1
+    bed[,2] = bed[,2]/resol + 1
   }
 
   mat = matrix(0, N, N)