Fixed bugs in the Gibbs sampler. Fixed a bug which will set probF to 0 if the protoco...

[rsem.git] / Gibbs.cpp

diff --git a/Gibbs.cpp b/Gibbs.cpp
index 4b1e2d5e8857566fbf646af6b905ed3c999e365f..7727f37b73452910e7cb5ad0a5faa37e98a1a5bb 100644 (file)
--- a/Gibbs.cpp
+++ b/Gibbs.cpp
@@ -200,8 +200,7 @@ void Gibbs(char* imdName) {
                }
 
                if (ROUND > BURNIN) {
-                       if ((ROUND - BURNIN -1) % GAP == 0) writeCountVector(fo);
-                       writeCountVector(fo);
+                       if ((ROUND - BURNIN - 1) % GAP == 0) writeCountVector(fo);
                        for (int i = 0; i <= M; i++) { 
                          pme_c[i] += counts[i] - 1;
                          pme_theta[i] += counts[i] / totc;