RSEM for DESeq2 without abundance matrix
1
0
Entering edit mode
le2336 ▴ 20
@le2336-10789
Last seen 4.2 years ago

I am working with RSEM outputs of RNA-seq data that I downloaded from a repository as .tsv files. I would like to apply VST to this data by first using tximport to import the RSEM outputs. However, the RSEM outputs from this repository only include the lengths and estimated counts, and are missing the abundance information. Would it still be possible to use these data with tximport and DESeq2?

Thanks very much.

tximport deseq2 rsem • 1.1k views
ADD COMMENT
2
Entering edit mode
@mikelove
Last seen 1 day ago
United States

No, tximport will give an error without the abundances.

If you have the gene-level counts and lengths, you can import the data manually and use DESeqDataSetFromMatrix with the counts matrix (rounded to integers). Then you can assign the lengths matrix as:

assays(dds)[["avgTxLength"]] <- ...

And this will then behave as if you had imported with tximport (one caveat is that you should make sure there are no genes with length of 0 -- tximport ensures this does not happen but it can occur with RSEM lengths).

ADD COMMENT
0
Entering edit mode

Thank you, this works great!

ADD REPLY

Login before adding your answer.

Traffic: 1004 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