how to get normalized counts from edgeR
1
2
Entering edit mode
wang peter ★ 2.0k
@wang-peter-4647
Last seen 9.3 years ago
dear all: i have two questions 1. is the "size factor" from DESeq equally to the "norm factor" in the edgeR 2. i can get the normalized counts by counts(obj,normalized=T) how to get normalized counts from edgeR can i use the each col of count matrix divided by norm factor -- shan gao Room 231(Dr.Fei lab) Boyce Thompson Institute Cornell University Tower Road, Ithaca, NY 14853-1801 Office phone: 1-607-254-1267(day) Official email:sg839 at cornell.edu Facebook:http://www.facebook.com/profile.php?id=100001986532253
DESeq DESeq • 12k views
1
Entering edit mode
Mark Robinson ▴ 880
@mark-robinson-4908
Last seen 5.1 years ago
Hi Shan, On 09.07.2012, at 05:05, wang peter wrote: > dear all: > i have two questions > > 1. is the "size factor" from DESeq equally to the "norm factor" in the edgeR No. In edgeR, the library size and additional normalization scaling factors are separated. See the two different columns in the $samples element of a 'DGEList' object: > d$samples group lib.size norm.factors sample.1 1 12367 0.9746280 sample.2 1 12088 1.0243994 sample.3 2 12823 0.9897893 sample.4 2 12374 1.0119267 In all the downstream code, the lib.size and norm.factors are multiplied together to act as the effective library size; this (product) would be similar to DESeq's size factor. > 2. i can get the normalized counts by > > counts(obj,normalized=T) > > how to get normalized counts from edgeR > can i use the each col of count matrix divided by norm factor In edgeR, you can use cpm(). For example: d <- DGEList(counts=...) nc <- cpm(d, normalized.lib.sizes=FALSE) You might have a skim through the edgeR user's guide; many of these concepts are described in more detail there. Best, Mark > -- > shan gao > Room 231(Dr.Fei lab) > Boyce Thompson Institute > Cornell University > Tower Road, Ithaca, NY 14853-1801 > Office phone: 1-607-254-1267(day) > Official email:sg839 at cornell.edu > Facebook:http://www.facebook.com/profile.php?id=100001986532253 > > _______________________________________________ > Bioconductor mailing list > Bioconductor at r-project.org > https://stat.ethz.ch/mailman/listinfo/bioconductor > Search the archives: http://news.gmane.org/gmane.science.biology.informatics.conductor ---------- Prof. Dr. Mark Robinson Bioinformatics Institute of Molecular Life Sciences University of Zurich Winterthurerstrasse 190 8057 Zurich Switzerland v: +41 44 635 4848 f: +41 44 635 6898 e: mark.robinson at imls.uzh.ch o: Y11-J-16 w: http://tiny.cc/mrobin ---------- http://www.fgcz.ch/Bioconductor2012
0
Entering edit mode

By this method are the new 'nc' table values in counts per million?