filter by seF and seR

Biosoft-ru · Mar 15, 2017 · 555696c · 555696c
1 parent e07582a
commit 555696c
Showing 1 changed file with 2 additions and 0 deletions.
diff --git a/output.c b/output.c
@@ -36,6 +36,8 @@ static void output(struct MixtureComp* comps, int32_t nComp, char* chr) {
     if(!score_filter_disabled && c.score <= Opt.score_cutoff) continue;
     if(c.delta <= 50 || c.delta >= 300) continue;
     if(c.se >= 50 || c.se <= 0) continue; //!!! ==0 ???
+    if(c.seF >= 50 || c.seF <= 0) continue; //!!! ==0 ???
+    if(c.seR >= 50 || c.seR <= 0) continue; //!!! ==0 ???
     if(c.sigmaSqF > 22500 || c.sigmaSqR > 22500) continue;
     int32_t from = (int32_t) (c.mu - c.delta/2 - 3*c.seF);
     int32_t to = (int32_t) (c.mu + c.delta/2 + 3*c.seR);