Question: rsem matrix to DESeq2
gravatar for Folder40g
15 months ago by
Folder40g120 wrote:

I know that I can pass the gene-level rsem matrix to limma and work with it.

But can I use it in DESeq?

I've a rsem matrix from firebrowse, gene-level ~20000 rows. DESeq requires raw counts, and as I only have the matrix with samples/gene I can't use tximport().


rsem • 1.8k views
ADD COMMENTlink modified 15 months ago by Kevin Blighe48k • written 15 months ago by Folder40g120
gravatar for Kevin Blighe
15 months ago by
Kevin Blighe48k
Kevin Blighe48k wrote:

You can use these with DESeq2 if you just round the numbers to whole integers and then input to DESeq2 with DESeqDataSetFromMatrix(). It is not ideal and using tximport would be preferred, as it does some adjustments for transcript length and transcript isoform abundances.

If you don't take my word, then take that of the DESeq2 developer: DESeq2 Following RSEM

Please be sure that you definitively have the un-normalised FireBrowse counts, though.


ADD COMMENTlink written 15 months ago by Kevin Blighe48k
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: 1664 users visited in the last hour