Question: rsem matrix to DESeq2
gravatar for Folder40g
4 months ago by
Folder40g110 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 • 438 views
ADD COMMENTlink modified 4 months ago by Kevin Blighe30k • written 4 months ago by Folder40g110
gravatar for Kevin Blighe
4 months ago by
Kevin Blighe30k
Republic of Ireland
Kevin Blighe30k 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 4 months ago by Kevin Blighe30k
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: 1259 users visited in the last hour