Question: How calculate RPKM values for DESeq?
1
gravatar for barrypraveen
2.2 years ago by
barrypraveen20
barrypraveen20 wrote:

Hi, I want to calculate RPKM values. I did DESeq and i got the results. My next step is i have to identify RPKM values for each gene. My pipeline is Bowtie>tophat>HTSeq>DESeq .plesase guide me through right path.

rna-seq R software error • 3.4k views
ADD COMMENTlink modified 2.2 years ago by ddiez1.7k • written 2.2 years ago by barrypraveen20
1

barrypraveen : I had asked you previously to provide titles that made sense for questions you are asking. Please do not include entire error message as the title of the post (I have removed it).

If you don't change the title of this post to something suitable we may close this question.

ADD REPLYlink modified 2.2 years ago • written 2.2 years ago by genomax67k

sorry for the inconvenience.

ADD REPLYlink written 2.2 years ago by barrypraveen20
1

That looks like code I wrote to determine the GC content and length of each gene. That's completely and totally unrelated to what you're trying to do. Also, you need alignments or reads or something along those lines to get RPKMs, you can't somehow generate that from an annotation. Please post what your actual goals are (the answer should have absolutely nothing to do with RPKMs).

ADD REPLYlink written 2.2 years ago by Devon Ryan90k

Okay. Thanks for your answer. yes, you are right i copied from Seqanswers. Sorry for the inconvenience.

ADD REPLYlink written 2.2 years ago by barrypraveen20
1

Also, as someone has pointed out before, you shouldn't use RPKMs with DEseq. You need raw counts. Maybe you don't mean to use the RPKMs for DEseq analysis?

ADD REPLYlink written 2.2 years ago by ddiez1.7k

The problem is i have to cross check. That is why i need RPKM Values. I have been following two pipelines. First one is Bowtie>tophat>cufflinks>cuffmerge>cuffqunat>cuffdiff (sometimes i use hisat instead of Bowtie and Tophat). In this pipeline i can get RPKM values from cufflinks. But recently i am using another pipeline Bowtie>tophat>HTSeq>DESeq. I am very new of R programming that is why i am asking. Please do share your comments on that. If i did anything wrong.

ADD REPLYlink modified 2.2 years ago • written 2.2 years ago by barrypraveen20
2

That is why i need RPKM Values.

No you don't.
If you want to compare two pipelines, use both correctly. Don't mess with the input of DESeq2. There is no good reason to use FPKM here.

ADD REPLYlink written 2.2 years ago by WouterDeCoster39k
3
gravatar for ddiez
2.2 years ago by
ddiez1.7k
Japan
ddiez1.7k wrote:

You can compute rpkm from a DESeqDataSet in the way indicated in @igor's answer. However, RPKMs should only be used for downstream analysis and not for testing differential expression. This is explicitly mentioned in the documentation of DEseq2.

In order to test for differential expression, we operate on raw counts and use discrete distributions as described in the previous Section 1.4. However for other downstream analyses – e.g. for visualization or clustering – it might be useful to work with transformed versions of the count data.

Same is said in the documentation for the older package DEseq:

The count values must be raw counts of sequencing reads. This is important for DESeq’s statistical model to hold, as only the actual counts allow assessing the measurement precision correctly. Hence, please do do not supply other quantities, such as (rounded) normalized counts, or counts of covered base pairs – this will only lead to nonsensical results.

Other relevant discussion:

Answer by Gordon Smyth about the use of RPKMs with voom, edgeR and DEseq2.

ADD COMMENTlink modified 2.2 years ago • written 2.2 years ago by ddiez1.7k
2
gravatar for igor
2.2 years ago by
igor7.6k
United States
igor7.6k wrote:

If you are using DESeq2, there is an fpkm() function to get RPKM/FPKM.

https://rdrr.io/bioc/DESeq2/man/fpkm.html

ADD COMMENTlink modified 2.2 years ago • written 2.2 years ago by igor7.6k
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1020 users visited in the last hour