non-integer counts for edgeR
1
5
Entering edit mode
@gordon-smyth
Last seen 8 hours ago
WEHI, Melbourne, Australia

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 • 2.0k views
ADD COMMENT
0
Entering edit mode
Tim Triche ★ 4.2k
@tim-triche-3561
Last seen 20 months ago
United States

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 COMMENT
0
Entering edit mode

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 REPLY

Login before adding your answer.

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