diff --git a/src/rmn.cpp b/src/rmn.cpp
index 3265b18..7e914e3 100644
--- a/src/rmn.cpp
+++ b/src/rmn.cpp
@@ -268,6 +268,7 @@ Rcpp::List estimate(
         calcSumU(meanU[iContig].slice(iExperiment), s[iContig], sumU);
       }
       mu.slice(iExperiment) = sym(sumU / sumZ);
+      std::cout << "sumU: " << sumU.row(146) << std::endl;
       std::cout << "mu: " << mu.slice(iExperiment).row(146) << std::endl;
     }
     // calculate lnEmPN