Question: rsem matrix to DESeq2
gravatar for Folder40g
8 days 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 • 83 views
ADD COMMENTlink modified 7 days ago by Kevin Blighe21k • written 8 days ago by Folder40g110
gravatar for Kevin Blighe
7 days ago by
Kevin Blighe21k
University College London Cancer Institute
Kevin Blighe21k 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 7 days ago by Kevin Blighe21k
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: 1378 users visited in the last hour