diff --git a/src/rmn.cpp b/src/rmn.cpp
index 92a499a..a8e952b 100644
--- 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