Question: non-integer counts for edgeR
5
gravatar for Gordon Smyth
7.0 years ago by
Gordon Smyth37k
Walter and Eliza Hall Institute of Medical Research, Melbourne, Australia
Gordon Smyth37k wrote:

Dear Mete,

No, you cannot use non-integer counts with edgeR.

If you must use non-integer counts, please use the voom() function in the limma package instead. This will do an analysis that is not too different from edgeR, just a little less powerful, and is not bothered by non-integer values.

Best wishes
Gordon

Edit: more recent versions of edgeR do now allow non-integer counts. However the values must be still be on the same scale as the counts, e.g., the values should add up approximately to the correct library sizes. So it is still not suitable for your weighted counts.

> Date: Mon, 21 May 2012 12:24:43 -0700
> From: Mete Civelek <mcivelek@mednet.ucla.edu>
> To: <bioconductor@r-project.org>
> Subject: [BioC] non-integer counts for edgeR
>
> Dear All,
>
> I am analyzing microRNAseq data with edgeR. Because some of the reads map
> to multiple locations, I weighed the counts based on the number of genomic
> locations that a read maps. For example, if a read maps to 3 locations,
> each location gets a count of 1/3. Of course, this means that the read counts
> for some miRNAs in some samples are not integers. Is it possible to use these
> counts to do differential expression analysis with a quantitative trait in
> edgeR? My understanding is "no" but I thought someone can suggest a
> solution.
>
> Thank you for your help.
>
> Best Regards,
> Mete
 

edger • 1.1k views
ADD COMMENTlink modified 5 months ago • written 7.0 years ago by Gordon Smyth37k
Answer: non-integer counts for edgeR
0
gravatar for Tim Triche
7.0 years ago by
Tim Triche4.2k
United States
Tim Triche4.2k wrote:

Hi Dr. Smyth,

Will you be writing / have you already written a paper summarizing what voom() does and how?

My understanding of voom() comes from

http://statsandgenomes.wordpress.com/2011/11/23/bioinformatics-seminar-professor-gordon-smyth-variance-models-for-rna-seq/

wherein it appears to be the current best general-purpose tool for treating RNAseq data as if it were microarray data. This would seem to bring up issues with rare transcripts and the estimation of their abundance, which is why I ask about a writeup.

Thanks as always for your patient explanations and high standards of work.

--t

ADD COMMENTlink modified 4.1 years ago by Gordon Smyth37k • written 7.0 years ago by Tim Triche4.2k

Hi Tim,

We have a few RNA-seq related papers in the pipeline, of which the voom paper is one.  We're investing a lot of time into trying to push edgeR and voom both along a similar pathway, to get the best use out of both.

Everyone in my lab would be happier if I didn't answer Bioconductor questions, because I'd get to their draft papers more quickly, and the voom one in particular ...

Regards
Gordon

Edit: voom paper now available at http://genomebiology.com/2014/15/2/R29

ADD REPLYlink modified 4.1 years ago • written 7.0 years ago by Gordon Smyth37k
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 16.09
Traffic: 357 users visited in the last hour