diff --git a/src/rmn.cpp b/src/rmn.cpp
index b33a083..4693c38 100644
--- a/src/rmn.cpp
+++ b/src/rmn.cpp
@@ -83,7 +83,7 @@ void estimateZ(
       z(i - w) = std::exp(F[i - w] + B[i + 1] + *forlnEmPN_it++ + lnTPN - lnL);
     }
     if(i == w){
-      std::cout << F[i - w] << '\t' << B[i + 1] << std::endl;
+      std::cout << z(i - w) << '\t' << z(i - 1) << '\t' << F[i - w] << '\t' << B[i + 1] << std::endl;
 
     }
   }