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

[rsem.git] / PairedEndModel.h

diff --git a/PairedEndModel.h b/PairedEndModel.h
index acd9f30a2b32ba0a3c30edf3104aa0e85f2a8c0e..bd046becb3fcd9cacbbc271863491800a6e67ee1 100644 (file)
--- a/PairedEndModel.h
+++ b/PairedEndModel.h
@@ -59,11 +59,11 @@ public:
                mw = NULL;
 
                if (isMaster) {
-                       ori = new Orientation(params.probF);
                        if (!estRSPD) rspd = new RSPD(estRSPD);
                        mld = new LenDist(params.mate_minL, params.mate_maxL);
                }
 
+               ori = new Orientation(params.probF);
                gld = new LenDist(params.minL, params.maxL);
                if (estRSPD) rspd = new RSPD(estRSPD, params.B);
                pro = new Profile(params.maxL);
@@ -321,7 +321,7 @@ void PairedEndModel::write(const char* outF) {
        npro->write(fo);
 
        if (mw != NULL) {
-         fprintf(fo, "%d\n", M);
+         fprintf(fo, "\n%d\n", M);
          for (int i = 0; i < M; i++) {
            fprintf(fo, "%.15g ", mw[i]);
          }