bg_from_sequence_dna will now return a probability distribution

jhkorhonen · Aug 9, 2016 · 0a18cbb · 0a18cbb
1 parent ac1ae42
commit 0a18cbb
Showing 1 changed file with 3 additions and 1 deletion.
diff --git a/core/moods_tools.cpp b/core/moods_tools.cpp
@@ -58,8 +58,10 @@ vector<double> bg_from_sequence_dna(const std::string &seq, const double ps)
         counts[m[c]]++;
     }
 
+    unsigned int count_all = counts[0] + counts[1] + counts[2] + counts[3];
+
     for (unsigned int j = 0; j < alphabet_size; ++j){
-        bg[j] = (((double)counts[j] + ps)/ ((double)seq.size() + alphabet_size * ps));
+        bg[j] = (((double)counts[j] + ps)/ ((double)count_all + alphabet_size * ps));
     }
 
     return bg;