[BioC] DESeq cuont table normalization

Fatemehsadat Seyednasrollah fatsey at utu.fi
Wed Sep 19 11:22:40 CEST 2012


Hello,

I am following the DESeq vignette and repeating what is mentioned there for my dataset. 
There is a step in the vignette which tries to normalize the count table to gain a common scale :
> head( counts( cds, normalized=TRUE ) )

When I run this line and then agin run the command : 
head(counts(cds))

I see the amount of reads before normalization. Should not it change to the new amounts I mean the normalized amounts?


> print (head(counts(mine, normalized = TRUE)))
                  X3        X5          X4       X4.1        X2       X3.1
A1BG     186.0036964 96.798204 231.6504223 90.4352020 92.027309  40.852657
A1BG-AS1  22.3204436 29.143545  15.0666941 28.3718281 18.191445   8.881012
A1CF       0.9300185  1.040841   0.9416684  1.7732393  1.070085   0.000000
A2LD1     93.0018482 73.899704  92.2835016 86.0021039 63.135014 113.676959
A2M        4.6500924  5.204205  21.6583728  0.8866196  5.350425   5.328607
A2ML1      0.0000000  0.000000   0.0000000  0.0000000  0.000000   1.776202
                X1       X1.1
A1BG      66.40690  90.562533
A1BG-AS1  21.75399  14.916182
A1CF       0.00000   1.065442
A2LD1    100.75530 121.460339
A2M       11.44947   5.327208
A2ML1      0.00000   0.000000


> head(counts(mine))
          X3 X5  X4 X4.1 X2 X3.1 X1 X1.1
A1BG     200 93 246  102 86   46 58   85
A1BG-AS1  24 28  16   32 17   10 19   14
A1CF       1  1   1    2  1    0  0    1
A2LD1    100 71  98   97 59  128 88  114
A2M        5  5  23    1  5    6 10    5
A2ML1      0  0   0    0  0    2  0    0



More information about the Bioconductor mailing list