is mu correct?

chung · Jan 20, 2016 · 1fae975 · 1fae975
1 parent a16e00f
commit 1fae975
Showing 1 changed file with 2 additions and 1 deletion.
diff --git a/src/rmn.cpp b/src/rmn.cpp
@@ -57,7 +57,7 @@ void estimateZ(
   for (i = 1; i < L + 1; ++i, ++F_it, ++B_it){
 
     if (i >= w){
-      std::cout << i << '\t' << *forlnEmPN_it << std::endl;
+
       *F_it = logSumExp(0, *forlnEmPN_it++ + lnTPN + *(F_it - w) - *(F_it - 1));
       *B_it = logSumExp(0, *revlnEmPN_it++ + lnTPN + *(B_it - w) - *(B_it - 1));
 
@@ -239,6 +239,7 @@ void estimate(
         calcSumU(meanU[iContig].slice(iExperiment), s[iContig], sumU);
       }
       mu.slice(iExperiment) = sumU / sumZ;
+      std::cout << mu.slice(iExperiment) << std::endl;
     }
     // calculate lnEmPN
     // re-estimate z