How to import RSEM results to DESeq2 for differential analysis
1
2
Entering edit mode
3.9 years ago
402374688 ▴ 30

The results of RSEM are expected counts of transcripts and genes and they are not integer. DESeq2, however, need a matrix of integer as input. The manual of RSEM said the results of RSEM can be used by DESeq. So my question is how to import non-integer RSEM counts to DESeq2 for downstream differential analysis? Thank you!

rsem deseq2 • 3.7k views
ADD COMMENT
0
Entering edit mode

Michael,

Thanks for your response. Could you please tell me why transcript-level counting is non-integer?

Thanks, Hadi

ADD REPLY
1
Entering edit mode

Please do not ask unrelated questions in old threads, and please use the ADD COMMENT button, not the answer field. It is not integer because it is abundance estimates, not counts. tximport will then produce gene level counts.

ADD REPLY
7
Entering edit mode
3.9 years ago
Michael Love ★ 2.3k

You can use tximport to import RSEM quants and then DESeqDataSetFromTximport() in DESeq2. The development version of tximport allows importing transcript-level abundances as well (previously just gene-level).


Edit 4.3.2020 by ATpoint since this thread gets a lot of views:

The manual contains example code on how to use RSEM output with tximport, which can then be processed with DESeqDataSetFromTximport()

ADD COMMENT

Login before adding your answer.

Traffic: 1849 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6