Fixed a bug when find_start_line() does not locate exactly at the sta…

…rting CpG site, resulting in a false -nan.
smithlabcode · Oct 31, 2015 · f63954e · f63954e · mengzhou · Oct 31, 2015
1 parent 38d2b0b
commit f63954e
Showing 1 changed file with 7 additions and 0 deletions.
diff --git a/src/analysis/roimethstat.cpp b/src/analysis/roimethstat.cpp
@@ -266,6 +266,13 @@ get_cpg_stats(const bool METHPIPE_FORMAT,
   find_start_line(chrom, start_pos, cpg_in);
 
   GenomicRegion cpg;
+  // find_start_line not necessarily locate at the start site.
+  // in this case the file pointer needs to be move forward,
+  // a little bit hopefully.
+  while (load_cpg(METHPIPE_FORMAT, cpg_in, cpg) &&
+      (cpg.get_chrom() < chrom ||
+       (cpg.same_chrom(region) &&
+       cpg.get_end() < start_pos)));
   while (load_cpg(METHPIPE_FORMAT, cpg_in, cpg) &&
          (cpg.same_chrom(region) &&
           cpg.get_end() <= end_pos)) {