Issue running RSEM
I am trying to run RSEM on my output from STAR. I have vreated the RSEM ref and have it in my working diretory, however when I run RSEM:
rsem-calculate-expression --paired-end --bam -p 16 "$f" RSEMref/ rsem_out/"$base"
I get this error:
rsem-parse-alignments RSEMref/ rsem_out/MR13_S13.temp/MR13_S13 rsem_out/MR13_S13.stat/MR13_S13 MR13_S13_mappedAligned.toTranscriptome.out.bam 3 -tag XM
Cannot open RSEMref/.grp! It may not exist.
However when I look in RSEMref/ there is a .grp file there, why can it not see it?
1 answer
I think it wants more than just the directory where the files are kept; it wants the prefix of the files too. So if your file is called something.grp, your command line should be
rsem-calculate-expression --paired-end --bam -p 16 "$f" RSEMref/something rsem_out/"$base"
Log in to answer this question.
More posts like this
-
I have encountered some bugs in transcriptome analysis using STAR and RSEM.
written by Pallondyle •
-
RSEM getting TPM from STAR alignment
written by Emily •
-
RSEM error while calculating expression
written by Jonah •
-
The two mates do not align to a same transcript -- RSEM error
written by abhilashkannan2987 •
-
Error in running RSEM
written by bioinformatics.queries
-
Trinity transcript_abundance.pl error using RSEM
written by dimitrischat
-
RSEM Error - rsem-calculate-expression
written by anoo •
-
need help,please -->RNA-seq, STAR+RSEM wrong ,
written by dq18 •
-
rsem-calculate-expression error message
written by walker •
-
rsem-calculate-expression fails to open iso_res file
written by emyers •
