Question: RPKM and differential expression analysis
gravatar for harrietjamieson86
2.9 years ago by
harrietjamieson8620 wrote:

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

ADD COMMENTlink written 2.9 years ago by harrietjamieson8620

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

ADD REPLYlink written 2.9 years ago by Devon Ryan89k
gravatar for Sean Davis
2.9 years ago by
Sean Davis25k
National Institutes of Health, Bethesda, MD
Sean Davis25k wrote:

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

ADD COMMENTlink written 2.9 years ago by Sean Davis25k
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: 1617 users visited in the last hour