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

[rsem.git] / PairedEndQModel.h

diff --git a/PairedEndQModel.h b/PairedEndQModel.h
index e328779b6675a3579aa221c32f490b99648053e4..5437eb4c0051974b501829aa6ad68c6f69072ac3 100644 (file)
--- a/PairedEndQModel.h
+++ b/PairedEndQModel.h
@@ -61,12 +61,12 @@ public:
                mw = NULL;
 
                if (isMaster) {
-                       ori = new Orientation(params.probF);
                        if (!estRSPD) rspd = new RSPD(estRSPD);
                        qd = new QualDist();
                        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);
                qpro = new QProfile();
@@ -335,7 +335,7 @@ void PairedEndQModel::write(const char* outF) {
        nqpro->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]);
          }