Question: inconsistent results of rsem quantification
gravatar for gnilihzeux
17 months ago by
gnilihzeux0 wrote:

I find rsem-calculate-expression put out two different results after my twice running with default parameters.
I use the bowtie2 aligner and foud parameter -k will search for one or more alignments and RSEM has a default set --bowtie2-k 200.
So, I guess that rsem will assign multiple alignment reads randomly.
Anyone suffers this ?

rna-seq isoform rsem • 398 views
ADD COMMENTlink modified 15 months ago by Biostar ♦♦ 20 • written 17 months ago by gnilihzeux0

How different are the results? Could you do a log2 exp vs log2 exp of the two runs?

ADD REPLYlink written 17 months ago by kristoffer.vittingseerup3.4k

The -k 200 parameter does _not_ mean that "rsem will assign multiple alignment reads randomly." The -k 200 parameter means that Bowtie2 will report up to 200 distinct alignments for a read, which RSEM will then determine how to probabilistically allocate in accordance with its model. When there are > 200 distinct "valid" alignments (alignments with a sufficiently hight score to be found and reported by Bowtie2), then, for such reads, Bowtie2 will report only up to 200 alignments. However, it is worth noting that this parameter for use with RSEM was certainly chosen with some care, and very few reads are likely to produce > 200 distinct valid alignments.

ADD REPLYlink written 15 months ago by Rob4.5k
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1002 users visited in the last hour