RPKM and differential expression analysis
1
2
Entering edit mode
6.5 years ago

Hello all,

I am trying to analyse differential expression in a dataset for which I only have RPKM values available to me. I usually use the R/BioC environment for RNA-seq analysis, and have read in various BioC documentation that using RPKM values in packages such as DESeq is really bad, from my understanding this is because the models used within the DESeq pipeline assume negative-binomial distribution of raw read counts, and RPKM doesn't fit at all into these assumptions (happy for anyone to correct my understanding there though).

So, I am trying to consider my options. Has anyone attempted something similar? I wondered if the best thing to do would be to just take the data out of the BioC environment entirely, and analyse with an appropriate GLM in R normal.

Any ideas would be appreicated!

Many thanks, Harriet

RNA-Seq rpkm differential expression R BioC • 3.0k views
ADD COMMENT
0
Entering edit mode

A GLM would seem to be the most appropriate approach. Obviously this isn't ideal, but it seems the only option.

ADD REPLY
2
Entering edit mode
6.5 years ago

Applying limma to log-transformed RPKM is an alternative approach.

ADD COMMENT
0
Entering edit mode

Hello Sean

Could you explain the rationale behind log-transforming RPKM? I have a similar problem to harrietjamieson86. I need to analyze public data and only RPKM is available.

ADD REPLY

Login before adding your answer.

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