rsem-calculate-expression error message

[Yvonne Zhao]

Hi 

I've been trying to follow the tutorial with the data provided. I managed to prepared reference sequence but am stuck with the rsem-calculate-expression step:

Below is my command in Mac terminal: 

rsem-calculate-expression -p 1 \

            --paired-end \

            --bowtie2 \

            --bowtie2-path /Users/yvonnezhao/miniconda/bin \

            --estimate-rspd \

             --append-names \

             --output-genome-bam \

/Users/yvonnezhao/RSEM_tutorial/data/SRR937564_1.fastq \

/Users/yvonnezhao/RSEM_tutorial/data/SRR937564_2.fastq \

/Users/yvonnezhao/RSEM_tutorial/ref/mouse_ref  \

exp/LPS_6h

It first proceeds but then stopped with the error message:

[E::hts_open_format] Failed to open file "b" : No such file or directory

Assertion failed: (in != 0), function BamWriter, file ./BamWriter.h, line 54.

"rsem-run-em /Users/yvonnezhao/RSEM_tutorial/ref/mouse_ref 3 exp/LPS_6h exp/LPS_6h.temp/LPS_6h exp/LPS_6h.stat/LPS_6h -p 1 -b b exp/LPS_6h.temp/LPS_6h.bam 0 --append-names" failed! Plase check if you provide correct parameters/options for the pipeline!

Could you help me understand which part went wrong and how to fix it? 

[Ruban Rex Peter]

Were you able to fix this error?

[Yvonne Zhao]

Yes, I did! It worked out when I replace the RSEM package included in the tutorial package with an updated RSEM package. 

--
RSEM website: http://deweylab.biostat.wisc.edu/rsem/
---
You received this message because you are subscribed to a topic in the Google Groups "RSEM Users" group.
To unsubscribe from this topic, visit https://groups.google.com/d/topic/rsem-users/kgv3Hp03M-0/unsubscribe.
To unsubscribe from this group and all its topics, send an email to rsem-users+...@googlegroups.com.
To view this discussion on the web visit https://groups.google.com/d/msgid/rsem-users/e6e0e2e9-f462-4b35-9a64-747f44eec524n%40googlegroups.com.