[BioC] Normalization
Ryan C. Thompson
rct at thompsonclan.org
Thu Feb 28 06:20:05 CET 2013
To answer your first question, when you first create a DGEList object,
all the normalization factors are initially set to 1 by default. This
is equivalent to no normalization. Once you use calcNormFactors, the
normalization factors will be set appropriately.
I'm not sure about the second question. Could you provide an example of
how you are obtaining pseudocounts with edgeR?
On Wed 27 Feb 2013 05:12:27 PM PST, Vittoria Roncalli wrote:
> Hi, I am a edgeR user and I am a little bit confused on the normalization
> topic.
> I am using EdgeR to get different expressed genes within 3 conditions
> (RnaSeq) with 3 replicates each.
> I am following the user guide step:
>
> -update counts file (from mapping against reference transcriptome)
> - filter the low counts reads (1cpm)
> - reassess library size
> - estimate common dispersion
>
> Mi first question is related to the normalization. Why, after I import my
> file, next to the library size there is then column with norm.factors?
>
> $samples
>
> group lib.size norm.factors
>
> X48h_C_r1.sam CONTROL 10898526 1
>
> X48h_C_r2.sam CONTROL 7176817 1
>
> X48h_C_r3.sam CONTROL 9511875 1
>
> X48h_LD_r1.sam LD 11350347 1
>
> X48h_LD_r2.sam LD 14836541 1
>
> X48h_LD_r3.sam LD 12635344 1
>
> X48h_HD_r1.sam HD 11840963 1
>
> X48h_HD_r2.sam HD 17335549 1
>
> X48h_HD_r3.sam HD 10274526 1
>
>
>
> Is the normalization automated? What is the difference with the
> "calNormFactors?"
>
> Moreover, if I do not run the calNormFactors, what is into the
> pseudo.counts output?
>
>
> I am very confused about those points.
>
>
> Thanks in advance for your help.
>
>
> Looking forward to hearing from you.
>
>
> Vittoria
>
>
> _______________________________________________
> 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
More information about the Bioconductor
mailing list