Made the BAM output of RSEM have all information contained in the input alignments...

[rsem.git] / Gibbs.cpp

diff --git a/Gibbs.cpp b/Gibbs.cpp
index c1f4fe276b2100ca6e66a151303cb062e0ec0d0e..63d35a829c059e12ae79bd8f3246fcec5b50438b 100644 (file)
--- a/Gibbs.cpp
+++ b/Gibbs.cpp
@@ -404,15 +404,17 @@ void writeEstimatedParameters(char* modelF, char* imdName) {
 
 int main(int argc, char* argv[]) {
        if (argc < 7) {
-               printf("Usage: rsem-run-gibbs-multi reference_name sample_name sampleToken BURNIN NSAMPLES GAP [-p #Threads] [--var] [-q]\n");
+               printf("Usage: rsem-run-gibbs reference_name imdName statName BURNIN NSAMPLES GAP [-p #Threads] [--var] [-q]\n");
                exit(-1);
        }
 
+       strcpy(imdName, argv[2]);
+       strcpy(statName, argv[3]);
+
        BURNIN = atoi(argv[4]);
        NSAMPLES = atoi(argv[5]);
        GAP = atoi(argv[6]);
-       sprintf(imdName, "%s.temp/%s", argv[2], argv[3]);
-       sprintf(statName, "%s.stat/%s", argv[2], argv[3]);
+
        load_data(argv[1], statName, imdName);
 
        nThreads = 1;