Generating RPKM data from dds (DESeq2 pipeline)
1
0
Entering edit mode
@harryhorsnell17-15582
Last seen 5.9 years ago

Dear all,

I have just generated my  DESeqDataSet and I am curious about how I can generate RPKM data from this

dds <- DESeqDataSet(se, design = ~ Cell + Hour )

I am quite new to the DESeq2 workflow so any suggestions will be much appreciated 

Many thanks,

Harry Horsnell

deseq2 dds RPKM • 959 views
ADD COMMENT
1
Entering edit mode
@mikelove
Last seen 2 hours ago
United States

We have an fpkm() function, try this code below to pop up all the functions. But we can't generate this with just a count matrix. You need to provide the gene lengths. The best way to get abundance is to use a transcript quantifier like Salmon, then import with tximport. The abundances will be available as txi$abundance.

> help(package="DESeq2",help_type="html")
ADD COMMENT

Login before adding your answer.

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